Research Article  Open Access
Study on the Weak Decays
Abstract
Motivated by the prospects of the potential particle at highluminosity heavyflavor experiments, we studied the weak decays, where = , , . The nonfactorizable contributions to hadronic matrix elements are taken into consideration with the QCDF approach. It is found that the CKMfavored decay has branching ratio of , which might be measured promisingly by the future experiments.
1. Introduction
The first evidence for upsilons, the bound states of , was observed in collisions of protons with a stationary nuclear target at Fermilab in 1977 [1, 2]. From that moment on, bottomonia have been a subject of intensive experimental and theoretical research. Some of the salient features of upsilons are as follows [3]. In the upsilons rest frame, the relative motion of the quark is sufficiently slow. Nonrelativistic Schrödinger equation can be used to describe the spectrum of bottomonia states and thus one can learn about the interquark binding forces. The particles, with the radial quantum number , and , decay primarily via the annihilation of the quark pairs into three gluons. Thus the properties of the invisible gluons and of the gluonquark coupling can be gleaned through the study of the upsilons decay. Compared with the , , and light quarks, the relatively large mass of the heavy quark implies a nonnegligible coupling to the Higgs bosons, making upsilons one of the best hunting grounds for light Higgs particles. By now, our knowledge of the properties of bottomonia comes primarily from annihilation.
The particle is the ground state of the vector bottomonia with quantum number of [4]. The mass of the particle, [4], is about three times heavier than the mass of the particle (the ground state of charmonia with the same quantum number of ). On one hand, compared with the decay, much richer decay channels could be accessed by the particle. On the other hand, the coupling constant for the decay is smaller than that for the decay due to the QCD nature of asymptotic freedom, which results in hadronic partial width , although the possible phase space in the decay is larger than that in the decay. In addition, the squared value of the quark charge, , is less than that of the quark charge, , which results in electromagnetic partial width . So the full decay width of the particle, , is less than that of the particle, [4]. Furthermore, one of the outstanding properties of all upsilons below threshold is their narrow decay width of tens of keV [4].
The and particles share the similar decay mechanism. As is the case for the particle, strong decays of the particle are suppressed by the phenomenological OZI (OkuboZweigIizuka) rules [5–7], so electromagnetic interactions and radiative transitions become competitive. It is expected that, at the lowest order approximation, the decay modes of the particle could be subdivided into four types. The lion’s share of the decay width is the hadronic decay via the annihilation of the quark pairs into three gluons, that is, some via [4]. The partial width of the electromagnetic decay via the annihilation of the quark pairs into a virtual photon could be written approximately as , where the value of is the ratio of inclusive production of hadrons to the pair production rate at the energy scale of and is the partial width of decay into dileptons. Branching ratio of the radiative decay is about [4]. The uptodate research for light Higgs bosons in the radiative decay has been performed by CLEO [8], Belle [9], and BaBar [10] collaborations. The magnetic dipole transition decay, , is very challenging to experimental physicist due to the very soft photon and pollution from other processes, such as [11]. The experimental signal for has not been discovered until now. Besides, the particle could also decay via the weak interactions, although the branching ratio for a single or quark decay is tiny, about [4]. In this paper, we will estimate the branching ratios for the flavorchanging nonleptonic weak decays with the QCD factorization (QCDF) approach [12, 13], where , and . The motivation is listed as follows.
From the experimental point of view, there is plenty of upsilons at the highluminosity dedicated bottomonia factories, for example, over at Belle (see Table 1). Upsilons are also observed by the onduty ALICE [17], ATLAS [18], CMS [19], and LHCb [20] experiments at LHC. It is hopefully expected that more upsilons could be accumulated with great precision at the running upgraded LHC and forthcoming SuperKEKB. The huge data samples will provide good opportunities to search for the weak decays which in some cases might be detectable. Theoretical studies on the weak decays are just necessary to offer a ready reference. For the twobody decay, final states with opposite charges have definite energies and momenta in the centerofmass frame of the particle. Particularly, identification of a single charged meson in the final state, which is free from inefficiently double tagging of the bottomed hadron pairs occurring above the threshold, would provide an unambiguous signature of the weak decay. Of course, small branching ratios make the observation of the weak decays extremely difficult, and evidences of an abnormally large production rate of single mesons in the decay might be a hint of new physics beyond the standard model.
From the theoretical point of view, the bottomchanging upsilon weak decay could permit overconstraining parameters obtained from meson decay, but few studies are devoted to the nonleptonic upsilon weak decay in the past. For the decay, the amplitude is usually treated as the factorizable product of two independent factors: one describing the transition between heavy quarkonium and and the other depicting the production of the state from the vacuum. Previous works, such as [21] based on the spin symmetry and nonrecoil approximation, [22] based on the heavy quark effective theory, and [23] based on the BauerStechWirbel (BSW) model [24], concentrated mainly upon the transition form factors which are related to the space integrals of the meson wave functions. As is well known, there exist hierarchical scales with nonrelativistic quantum chromodynamics (NRQCD) [25–27] which is an approach to deal with the heavy quarkonium; that is, , where is the mass of heavy quark with typical velocities . The physical contributions at scales of are absorbed into wave functions of the and particles; thus the transition should be dominated by the nonperturbative dynamics. The nonfactorizable contributions above scales of have not been taken seriously in previous works. About 2000, Beneke et al. proposed the QCDF approach [12, 13], where nonfactorizable contributions could be estimated systematically with the perturbation theory based on collinear factorization approximation and power countering rules in the heavy quark limit [13], and the QCDF approach has been widely applied to nonleptonic meson decays. So it should be very interesting to study the weak decays by considering nonfactorizable contributions with the attractive QCDF approach.
This paper is organized as follows. In Section 2, we will present the theoretical framework and the amplitudes for the nonleptonic twobody weak decays with the QCDF approach. Section 3 is devoted to numerical results and discussion. The last section is our summary.
2. Theoretical Framework
2.1. The Effective Hamiltonian
The low energy effective Hamiltonian responsible for the decays is [28]where the Fermi coupling constant [4]. Using the Wolfenstein parameterization, the CabibboKobayashiMaskawa (CKM) factors can be expanded as a power series in the small parameter [4],The Wilson coefficients summarize the physical contributions above scales of . The Wilson coefficients are calculable with the perturbation theory and have properly been evaluated to the nexttoleading order (NLO) with the renormalization group (RG) equation. The numerical values of the Wilson coefficients at scales of in naive dimensional regularization scheme are listed in Table 2. The local tree fourquark operators are defined as follows:where and are color indices and the sum over repeated indices is understood. To obtain the decay amplitudes, the remaining and the most intricate works are to calculate accurately hadronic matrix elements of local operators.

2.2. Hadronic Matrix Elements
Phenomenologically, the simplest treatment with hadronic matrix elements of fourquark operators is approximated by the product of two current matrix elements with color transparency ansatz [29] and naive factorization (NF) scheme [30, 31], and current matrix elements are further parameterized by decay constants and transition form factors. For example, previous studies on the decay [22, 23] were based on NF approach.
As is well known, NF’s defect is the disappearance of the renormalization scale dependence, the strong phases, and the nonfactorizable contributions from hadronic matrix elements, resulting in nonphysical amplitudes and no violating asymmetries. To remedy NF’s deficiencies, Beneke et al. proposed that hadronic matrix elements could be written as the convolution integrals of hard scattering kernels and light cone distribution amplitudes with the QCDF approach [12, 13].
For the decay, the spectator quark is the heavy bottom (anti)quark. According to the QCDF’s power counting rules [13], contributions from the spectator scattering are power suppressed. With the QCDF master formula, hadronic matrix elements could be written aswhere transition form factor and light cone distribution amplitudes of the emitted meson are nonperturbative input parameters and hard scattering kernels are computable order by order with the perturbation theory in principle.
The leading twist twovalenceparticle distribution amplitudes of pseudoscalar and longitudinally polarized vector meson are defined in terms of Gegenbauer polynomials [32, 33]:where , is the Gegenbauer moment, and .
After calculation, the decay amplitudes could be written as
The coefficient in (6), including nonfactorizable contributions from QCD radiative vertex corrections, is written as [34]
For the transversely polarized vector meson, the factor is zero beyond leading twist (twist2) contributions. For the pseudoscalar and longitudinally polarized vector meson, with the modified minimal subtraction () scheme, the factor is written as [34]whereand the relations are as follows:
The numerical values of coefficient at scales of are listed in Table 2.
2.3. Decay Constants and Form Factors
The matrix elements of current operators are defined as follows:where and are the decay constants of pseudoscalar and vector mesons, respectively, and and denote the mass and polarization of vector meson, respectively.
The transition form factors are defined as follows [22–24]:where and is required compulsorily to cancel singularities at the pole . There is a relation among these form factors:
It is clearly seen that there are only three independent form factors, and , at the pole for the decays. The form factors at the pole could be written as the overlap integrals of wave functions of mesons [24]; that is,where is a Pauli matrix acting on the spin indices of the decaying bottom quark and and denote the fraction of the longitudinal momentum and the transverse momentum carried by the nonspectator quark, respectively.
Using the separation of momentum and spin variables, the wave functions of mesons can be written aswith the normalization condition,where denote the spin of valence quark in meson; ; and for the and particles, respectively.
For the ground states of heavy quarkonia and particles, we will take the solution of the Schödinger equation with nonrelativistic threedimensional scalar harmonic oscillator potential,where the parameter determines the average transverse quark momentum; that is, . According to the power counting rules of NRQCD [25], the characteristic magnitude of the moment is order of and . So we will take in our calculation. Employing the substitution ansatz [35],where and is the mass of valence quark. Setting , we can obtainwhere is a normalization factor.
Using the aforementioned convention, we getwhere the uncertainties come from variation of valence quark mass . In addition, according to the NRQCD argument, the relativistic corrections and highertwist effects might give uncertainties of , about 10%~30%. Values at could, in principle, be extrapolated by assuming the form factors dominated by a proper pole which is unknown or calculated with other methods, such as the approach using the BetheSalpeter wave functions with the help of the nonrelativistic instantaneous approximation and the potential model based on the Mandelstam formalism [36]. Here, we will follow the common practice for nonleptonic decays with the QCDF approach. Values of form factors at are taken to offer an order of magnitude estimation, because both and are heavy quarkonium and the recoil effects might be not so significant.
2.4. Decay Amplitudes
With the aforementioned definition of hadronic matrix elements, the decay amplitudes of decays can be written as
For the decays, the hadronic matrix elements in (6) can also be expressed as [37]
The definition of helicity amplitudes iswhere invariant amplitudes , , and and variable areThe scalar amplitudes , , and describe the , , and wave contributions, respectively. Clearly, compared with the wave amplitude, the and wave amplitudes are suppressed by a factor .
3. Numerical Results and Discussion
In the rest frame of particle, branching ratio for nonleptonic weak decays can be written aswhere the decay width [4]; the momentum of final states is
The input parameters, including the CKM Wolfenstein parameters, masses of and quarks, decay constants, and Gegenbauer moments of distribution amplitudes in (5), are collected in Table 3. If not specified explicitly, we will take their central values as the default inputs. Our numerical results on the averaged branching ratios for the decays are displayed in Table 4, where theoretical uncertainties of the QCDF results come from the CKM parameters, the renormalization scale , masses of and quarks, and hadronic parameters (decay constants and Gegenbauer moments), respectively. For the sake of comparison, previous results of [22, 23] are reevaluated with coefficient , where the scenario of the flavor dependent parameter in [23] is taken. The following are some comments.(1)The QCDF’s results fall in between those of [22, 23], because the form factors in our calculation fall in between those of [22, 23].(2)There is a clear hierarchical relationship, . These are two dynamical reasons. One is that the CKM factor responsible for the decay is suppressed by a factor of relative to the CKM factor responsible for the , decays. The other is that the orbital angular momentum .(3)The CKMfavored dominated decay has the largest branching ratio, ~10^{−10}, which should be sought for with high priority and firstly observed at the running LHC and forthcoming SuperKEKB. For example, the production cross section in Pb collision can reach up to a few with the ALICE detector at LHC [38]. Therefore, per 100 fb^{−1} data collected at ALICE, over particles are in principle available, corresponding to tens of events if with about 10% reconstruction efficiency.(4)There are many uncertainties on the QCDF’s results. The first uncertainty, about 7~8%, from the CKM factors could be lessened with the improvement on the precision of the Wolfenstein parameter . The second uncertainty from the renormalization scale should, in principle, be reduced by inclusion of higher order corrections to hadronic matrix elements. The third uncertainty is due to the fact that masses of and quark affect the shape lines of wave functions and hence the magnitude of form factors and branching ratios. The fourth uncertainty from hadronic parameters is expected to be reduced with the relative ratio of branching ratios.(5)Other factors, such as the contributions of higher order corrections to hadronic matrix elements, relativistic effects, and dependence of form factors, which are not considered in this paper, deserve the dedicated study. Our results just provide an order of magnitude estimation.
4. Summary
With the sharp increase of the data sample at highluminosity dedicated heavyflavor factories, the bottomchanging weak decays are interesting in exploring the underlying mechanism responsible for transition between heavy quarknoia, investigating perturbative and nonperturbative effects and overconstraining parameters from decays. The weak decays are allowable within the standard model, though their branching ratios are expected to be tiny in comparison to the conventional strong and electromagnetic decays. In this paper, we studied the nonleptonic weak decays, which are dominated based on the low energy effective theory and hence should have large branching ratios among weak decay modes. Considering the nonfactorizable contributions to hadronic matrix elements with the QCDF approach, we estimated the branching ratios of the weak decays, where transition form factors are obtained by the integrals of wave functions with the nonrelativistic isotropic harmonic oscillator potential. The prediction on branching ratios for the decays is the same order as previous works [22, 23]. The CKMfavored decay has relatively large branching ratio, ~10^{−10}, and might be detectable in future experiments.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The authors thank the referees for their helpful comments. The work is supported by the National Natural Science Foundation of China (Grant nos. 11475055, 11275057, U1232101, and U1332103). Dr. Wang is thankful for the support from CCNUQLPL Innovation Fund (QLPL201411).
References
 S. Herb, D. C. Hom, L. M. Lederman et al., “Observation of a dimuon resonance at 9.5 GeV in 400GeV protonnucleus collisions,” Physical Review Letters, vol. 39, p. 252, 1977. View at: Publisher Site  Google Scholar
 W. Innes, J. A. Appel, B. C. Brown et al., “Observation of structure in the Υ region,” Physical Review Letters, vol. 39, no. 20, pp. 1240–1242, 1977. View at: Publisher Site  Google Scholar
 P. Franzini and J. LeeFranzini, “Upsilon resonances,” Annual Review of Nuclear and Particle Science, vol. 33, pp. 1–28, 1983. View at: Publisher Site  Google Scholar
 Particle Data Group, “Review of particle physics,” Chinese Physics C, vol. 38, no. 9, Article ID 090001, 2014. View at: Publisher Site  Google Scholar
 S. Okubo, “$\phi $meson and unitary symmetry model,” Physics Letters, vol. 5, no. 2, pp. 165–168, 1963. View at: Publisher Site  Google Scholar
 G. Zweig, CERNTH401, 402, 412, 1964.
 J. Iizuka, “A systematics and phenomenology of meson family,” Progress of Theoretical Physics Supplement, vol. 37, pp. 21–34, 1966. View at: Publisher Site  Google Scholar
 W. Love, V. Savinov, H. Mendez et al., “Search for very light CPodd Higgs boson in radiative decays of ${\rm Y}(1S)$,” Physical Review Letters, vol. 101, Article ID 151802, 2008. View at: Publisher Site  Google Scholar
 M.C. Chang, “Search for light Higgs and dark photons at BaBar and Belle,” in Proceedings of the European Physical Society Conference on High Energy Physics (EPSHEP '13), p. 270, Stockholm, Sweden, July 2013. View at: Google Scholar
 J. Lees, V. Poireau, V. Tisserand et al., “Search for a light Higgs resonance in radiative decays of the $\mathrm{{\rm Y}}$(1S) with a charm tag,” Physical Review D, vol. 91, no. 7, Article ID 071102, 9 pages, 2015. View at: Publisher Site  Google Scholar
 C. Patrignani, T. Pedlar, and J. Rosner, “Recent results in bottomonium,” Annual Review of Nuclear and Particle Science, vol. 63, no. 1, pp. 21–44, 2013. View at: Publisher Site  Google Scholar
 M. Beneke, G. Buchalla, M. Neubert, and C. T. Sachrajda, “QCD factorization for $B\to \pi \pi $ decays: strong phases and $CP$ violation in the heavy quark limit,” Physical Review Letters, vol. 83, no. 10, pp. 1914–1917, 1999. View at: Publisher Site  Google Scholar
 M. Beneke, G. Buchalla, M. Neubert, and C. T. Sachrajda, “QCD factorization for exclusive nonleptonic Bmeson decays: general arguments and the case of heavylight final states,” Nuclear Physics B, vol. 591, no. 12, pp. 313–418, 2000. View at: Publisher Site  Google Scholar
 J. Brodzicka, T. Browder, P. Chang et al., “Physics achievements from the Belle experiment,” Progress of Theoretical and Experimental Physics, Article ID 04D001, 2012. View at: Publisher Site  Google Scholar
 J. P. Lees, V. Poireau, V. Tisserand et al., “Study of $\mathrm{{\rm Y}}(3S,2S)\to \eta \mathrm{{\rm Y}}(1S)$ and $\mathrm{{\rm Y}}(3S,2S)\to {\pi}^{+}{\pi}^{}\mathrm{{\rm Y}}(1S)$ hadronic transitions,” Physical Review D, vol. 84, no. 9, Article ID 092003, 2011. View at: Publisher Site  Google Scholar
 R. Briere, T. Ferguson, G. Tatishvili et al., “Comparison of particle production in quark and gluon fragmentation at $\sqrt{s}~10$ GeV,” Physical Review D, vol. 76, no. 1, Article ID 012005, 12 pages, 2007. View at: Publisher Site  Google Scholar
 B. Abelev, J. Adam, D. Adamová et al., “Measurement of quarkonium production at forward rapidity in pp collisions at $\sqrt{s}=7$ TeV,” The European Physical Journal C, vol. 74, article 2974, 2014. View at: Publisher Site  Google Scholar
 G. Ada, T. Abajyan, B. Abbott et al., “Measurement of upsilon production in 7 TeV $pp$ collisions at ATLAS,” Physical Review D, vol. 87, no. 5, Article ID 052004, 31 pages, 2013. View at: Publisher Site  Google Scholar
 V. Khachatryan, A. M. Sirunyan, S. Chatrchyan et al., “Measurements of the ${\rm Y}(1S)$, ${\rm Y}(2S)$, and ${\rm Y}(3S)$ differential cross sections in pp collisions at $\sqrt{s}=7$ TeV,” Physics Letter B, vol. 749, pp. 14–34, 2015. View at: Google Scholar
 R. Aaij, C. Abellan Beteta, B. Adeva et al., “Production of J/ψ and Υ mesons in pp collisions at $\sqrt{s}=8$ TeV,” Journal of High Energy Physics, vol. 2013, no. 6, article 064, 2013. View at: Publisher Site  Google Scholar
 M. A. SanchisLozano, “On the search for weak decays of heavy quarkonium in dedicated heavyquark factories,” Zeitschrift für Physik C Particles and Fields, vol. 62, no. 2, pp. 271–279, 1994. View at: Publisher Site  Google Scholar
 K. K. Sharma and R. C. Verma, “Rare decays of $\mathrm{\Psi}$ and $\mathrm{\gamma}$,” International Journal of Modern Physics A, vol. 14, no. 6, p. 937, 1999. View at: Publisher Site  Google Scholar
 R. Dhir, R. C. Verma, and A. Sharma, “Effects of flavor dependence on weak decays of $\mathrm{J}/\mathrm{\psi}$ and $\mathrm{{\rm Y}}$,” Advances in High Energy Physics, vol. 2013, Article ID 706543, 12 pages, 2013. View at: Publisher Site  Google Scholar
 M. Wirbel, B. Stech, and M. Bauer, “Exclusive semileptonic decays of heavy mesons,” Zeitschrift für Physik C: Particles and Fields, vol. 29, no. 4, pp. 637–642, 1985. View at: Publisher Site  Google Scholar
 G. Legage, L. Magnea, C. Nakhleh, U. Magnea, and K. Hornbostel, “Improved nonrelativistic QCD for heavyquark physics,” Physical Review D, vol. 46, no. 9, pp. 4052–4067, 1992. View at: Publisher Site  Google Scholar
 G. T. Bodwin, E. Braaten, and G. P. Legage, “Rigorous QCD analysis of inclusive annihilation and production of heavy quarkonium,” Physical Review D, vol. 51, no. 3, pp. 1125–1171, 1995. View at: Publisher Site  Google Scholar
 N. Brambilla, A. Pineda, J. Soto, and A. Vairo, “Effectivefield theories for heavy quarkonium,” Reviews of Modern Physics, vol. 77, no. 4, pp. 1423–1496, 2005. View at: Publisher Site  Google Scholar
 G. Buchalla, A. J. Buras, and M. E. Lautenbacher, “Weak decays beyond leading logarithms,” Reviews of Modern Physics, vol. 68, no. 4, article 1125, 1996. View at: Publisher Site  Google Scholar
 J. Bjorken, “Topics in Bphysics,” Nuclear Physics B—Proceedings Supplements, vol. 11, pp. 325–341, 1989. View at: Publisher Site  Google Scholar
 N. Cabibbo and L. Maiani, “Twobody decays of charmed mesons,” Physics Letters B, vol. 73, no. 45, pp. 418–422, 1978. View at: Publisher Site  Google Scholar
 D. Fakirov and B. Stech, “F and Ddecays,” Nuclear Physics B, vol. 133, no. 2, pp. 315–326, 1978. View at: Publisher Site  Google Scholar
 P. Ball and G. Jones, “Twist3 distribution amplitudes of ${K}^{*}$ and $\varphi $ mesons,” Journal of High Energy Physics, vol. 2007, no. 3, article 69, 2007. View at: Publisher Site  Google Scholar
 P. Ball, V. Braun, and A. Lenz, “Highertwist distribution amplitudes of the K meson in QCD,” Journal of High Energy Physics, vol. 2006, no. 5, article 004, 2006. View at: Google Scholar
 J. Sun, G. Xue, Y. Yang, G. Lu, and D. Du, “Study of ${B}_{\stackrel{}{c}}\to J/\psi {\pi}^{}$, ${\eta}_{c}{\pi}^{}$ decays with QCD factorization,” Physical Review D, vol. 77, no. 7, Article ID 074013, 2008. View at: Publisher Site  Google Scholar
 B. Xiao, X. Qin, and B. Ma, “The kaon form factor in the lightcone quark model,” The European Physical Journal A, vol. 15, no. 4, pp. 523–527, 2002. View at: Publisher Site  Google Scholar
 C. Chang and Y. Chen, “Decays of the ${B}_{c}$ meson,” Physical Review D, vol. 49, pp. 3399–3411, 1994. View at: Publisher Site  Google Scholar
 G. Valencia, “Angular correlations in the decay,” Physical Review D, vol. 39, no. 11, pp. 3339–3345, 1989. View at: Publisher Site  Google Scholar
 B. Abelev, J. Adam, D. Adamova et al., “Production of inclusive ${\rm Y}(1S)$ and ${\rm Y}(2S)$ in p–Pb collisions at $\sqrt{{s}_{\text{N}\text{N}}}=5.02$ TeV,” Physics Letters B, vol. 740, pp. 105–117, 2015. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2015 Junfeng Sun et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The publication of this article was funded by SCOAP^{3}.