Anatomy of $B_s \to VV $ decays and effects of nexttoleading order contributions in the perturbative QCD factorization approach
Abstract
By employing the perturbative QCD (PQCD) factorization approach, we calculated the branching ratios, CPviolating asymmetries, the longitudinal and transverse polarization fractions and other physical observables of the thirteen charmless hadronic decays with the inclusion of all currently known nexttoleading order (NLO) contributions. We focused on the examination of the effects of all those currently known NLO contributions and found that: (a) for the measured decays and , the NLO contributions can provide to enhancements to the leading order (LO) PQCD predictions of their CPaveraged branching ratios, and consequently the agreement between the PQCD predictions and the measured values are improved effectively after the inclusion of the NLO contributions; (b) for the measured decays, the NLO corrections to the LO PQCD predictions for and are generally small in size, but the weak penguin annihilation contributions play an important role in understanding the data about their decay rates, and ; (c) the NLO PQCD predictions for above mentioned physical observables do agree with the measured ones and the theoretical predictions from the QCDF, SCET and FAT approaches; (d) for other considered decays, the NLO PQCD predictions for their decay rates and other physical observables are also basically consistent with the theoretical predictions from other popular approaches, future precision measurements could help us to test or examine these predictions.
pacs:
13.25.Hw, 12.38.Bx, 14.40.NdKey Words: meson decays; The PQCD factorization approach; Branching ratios; Polarization fractions; Relative phases
I Introduction
During the past three decades, the twobody charmless hadronic decays, with being the light vector mesons and , have been studied by many authors based on rather different factorization approaches qcdf2 (); qcdf07 (); qcdf09 (); qcdfwa (); qcdfll (); gf99 (); ali07 (); jpg06 (); pqcd2 (); scet (); scetv (). Several such decay modes, such as decay, have been observed by CDF and LHCb experiments cdf (); lhcbks (); lhcbks1 (); lhcbphi (); lhcbphi2 (); lhcbphiks (); lhcbrhophi (); lhcb0 (); pdg2016 (); hfag2016 (). When compared with the similar (here , and ) decays, are indeed much more complicated due to the fact that more helicity amplitudes should be taken into account. The decays can offer, consequently, rich opportunities for us to test the Stand Model (SM) and to search for the exotic new physics beyond the SM.
Experimentally, a large transverse polarization fraction of was firstly observed in 2003 by BABAR and Belle Collaborations babe (). The new world averages of as given by HFAG2016 hfag2016 () for decays, for example, are the following:
(1)  
(2) 
These measured values were in strong confliction with the general expectation in the naive factorization ansatz naive1 (), which is the socalled “polarization puzzle” 2004pft (); 2004pa (); lipa (); lipa1 (). The similar deviations also be observed later for and decays Aubert:2008bc (); hfag2016 ().
For the charmless decays studied in this paper, the similar puzzles have also been observed by CDF and LHCb Collaboration for , and decay modes cdf (); lhcbks (); lhcbks1 (); lhcbphi (); lhcbphi2 (); lhcbrhophi (). The new world averages of and as given by HFAG2016 hfag2016 () for these three decay modes are the following:
(3)  
(4)  
(5) 
More measurements are expected in the near future.
Theoretically, a number of strategies were proposed to resolve the above mentioned ”polarization puzzle” within and/or beyond the SM. For example, the weak penguin annihilation contributions in QCD factorization (QCDF) approach was proposed by Kagan Kann (), the final state interactions were considered in Refs. final (); fsi (); 2004pft (), the formfactor tuning in the perturbative QCD (PQCD) approach was suggested by Li lipa (), and even the exotic new physics effects have been studied by authors in Refs. np (); np1 (). Obviously, it is hard to get a good answer to this seemingly longstanding puzzle at present. However, according to the statement in Ref. lipa1 (), the complicated QCD dynamics involved in such decays should be fully explored before resorting to the possible new physics beyond the SM. Therefore, the QCDF approach qcdf2 (); qcdf07 (); qcdf09 (); qcdfwa (); qcdfll (); qcdfv (), the softcollinear effective theory (SCET) scet (); scetv () and the PQCD approach ali07 (); Li2002 (); jpg06 (); pqcd2 (); pqcdv (); pqcdv1 (), have been adopted to investigate these kinds of decays systematically.
The twobody charmless hadronic decays have been systematically studied in the PQCD approach at leading order (LO) in 2007 ali07 (). Recently, the authoers of Ref. pqcd2 () made improved estimations for the modes by keeping the terms with the higher power of the ratios in the PQCD approach, with and being the masses of the initial and final states. However, there still existed some issues to be clarified, e.g. the measured large decay rates for and decays, the latest measurement of a smaller for decay, etc.
Therefore, we would like to revisit those twobody charmless decays by taking into account all currently known nexttoleading order (NLO) contributions in the PQCD factorization approach. We will focus on the effects of the NLO contributions arising from various possible sources, such as the QCD vertex corrections (VC), the quark loops (QL), and the chromomagnetic penguins Li05 (); nlo05 () in the SM. As can be seen from Refs. fan2013 (); xiao08b (); nlo05 (); Li05 (); xiao2014 (), the NLO contributions do play an important role in understanding the known anomalies of physics such as the amazingly large decay rates fan2013 (); xiao08b (), the longitudinalpolarization dominated Li05 () and the evidently nonzero , i.e., the famous “puzzle” nlo05 (); xiao2014 (), and so forth. Very recently, we extend these calculations to the cases such as decays xiao14a (), decays xiao14b () and decays xiao17 (). We found that the currently known NLO contributions can interfere with the LO part constructively or destructively for those considered meson decay modes. Consequently, the agreement between the PQCD predictions and the experimental measurements of the CPaveraged branching ratios, the polarization fractions and CPviolating asymmetries was indeed improved effectively due to the inclusion of the NLO contributions.
This paper is organized as follows. In Sec. II, we shall present various decay amplitudes for the considered decay modes in the PQCD approach at the LO and NLO level. We show the PQCD predictions and several phenomenological analyses for the branching ratios, CPviolating asymmetries and polarization observables of thirteen decays in Sec III. A short summary is given in Sec. IV.
Ii Decay amplitudes at LO and NLO level
We treat the meson as a heavylight system and consider it at rest for simplicity. By employing the lightcone coordinates, we define the meson with momentum , the emitted meson with the momentum along the direction of , and the recoiled meson with the momentum in the direction of (Here, and are the lightlike dimensionless vectors), respectively, as the following,
(6) 
The polarization vectors of the final states can then be parametrized as:
(7) 
with being the longitudinal(transverse) polarization vector.
The momenta carried by the light antiquark in the initial and final mesons are chosen as follows:
(8) 
The integration over and will lead conceptually to the decay amplitudes in the PQCD approach,
(9)  
in which, is the conjugate space coordinate of transverse momentum , stands for the Wilson coefficients evaluated at the scale , and denotes the hadron wave functions, which are nonperturbative but universal inputs, of the initial and final states. The kernel describes the hard dynamics associated with the effective ”sixquark interaction” exchanged by a hard gluon. The Sudakov factors and together suppress the soft dynamics in the endpoint region effectively li2003 ().
ii.1 Wave functions and decay amplitudes
Without the endpoint singularities in the evaluations, the hadron wave functions are the only input in the PQCD approach. These nonperturbative quantities are process independent and could be obtained with the techniques of QCD sum rule and/or Lattice QCD, or be fitted to the measurements with good precision.
For meson, its wave function could be adopted with the Lorentz structure ali07 (); pqcd2 ()
(10) 
in which the distribution amplitude is modeled as
(11) 
with being the shape parameter. We take GeV for the meson based on the studies of lattice QCD and lightcone sum rule wbs1 (); wbs2 (); wbs3 (). The normalization factor will be determined through the normalization condition: with the decay constant GeV.
For the vector meson, the longitudinally and transversely polarized wave functions up to twist3 are given by wbs3 (); pball98 ()
(12) 
where and are the momentum and the mass of the light vector mesons, and is the corresponding longitudinal(transverse) polarization vector Li2002 (). Here is LeviCivita tensor with the convention .
The twist2 distribution amplitudes and can be written in the following form wbs3 (); pball98 ()
(13)  
(14) 
where , is the decay constants of the vector meson with longitudinal(transverse) polarization. The Gegenbauer moments here are the same as those in Refs. wbs3 (); pball07 (); pball98 ():
(15) 
ii.2 Example of the LO decay amplitudes
In the SM, for the considered decays induced by the transition with , the weak effective Hamiltonian can be written asburas96 (),
(17) 
where the Fermi constant GeV, and is the CabbiboKobayashiMaskawa(CKM) matrix element, are the Wilson coefficients and are the local fourquark operators. For convenience, the combinations of the Wilson coefficients are defined as usual ali07 (); pqcd2 ():
(18) 
where the upper(lower) sign applies, when is odd(even).
At leading order, as illustrated in Fig. 1, there are eight types of Feynman diagrams contributing to the decays, which can be classified into three types: the factorizable emission diagrams ( Fig. 1(a) and 1(b)); the nonfactorizable emission diagrams (Fig. 1(c) and 1(d)); and the annihilation diagrams (Fig. 1(e)1(h)). As mentioned in the Introduction, the considered thirteen modes have been studied at LO in the PQCD approach ali07 (); pqcd2 (). The factorization formulas of decay amplitudes with various topologies have been presented explicitly in Ref. ali07 (). Therefore, after the confirmation by our independent recalculations, we shall not collect those analytic expressions here for simplicity. In this work, we aim to examine the effects of all currently known NLO contributions to the considered decay modes in the PQCD approach to see whether one can improve the consistency between the theory and the experiment in the SM or not, which would be help for us to judge the necessity of the exotic new physics beyond the SM.
For decays, both of the longitudinal and transverse polarizations will contribute. Then, the decay amplitudes can be decomposed into three partspqcd2 ():
(19) 
where and correspond to the longitudinally, normally and transversely polarized amplitudes, respectively, whose detailed expressions can be inferred from Refs. ali07 (); pqcd2 ().
ii.3 NLO contributions
In the framework of the PQCD approach, many twobody charmless decays have been investigated by including currently known NLO contributions,for example, in Refs. fan2013 (); nlo05 (); xiao14a (); xiao14b (); xiao17 (); xiao08b (); xiao2014 (); zhang09 (); zhou12 (). Of course, some NLO contributions are still not known at present, as discussed in Ref. fan2013 (). The currently known NLO corrections to the LO PQCD predictions of decays are the following:

(a) The NLO Wilson coefficients (NLOWC), the renormalization group running matrix at NLO level and the strong coupling constant at twoloop level as presented in Ref. buras96 ();

(d) The NLO contributions from the chromomagnetic penguin (MP) operator o8g2003 (); Li05 (); nlo05 () as illustrated in Figs. 2(g)2(h).
In this paper, we adopt directly the formulas for all currently known NLO contributions from Refs. nlo05 (); Li05 (); o8g2003 (); xiao08b (); fan2013 (); xiao14a (); xiao14b (); xiao2014 () without further discussions about the details. Moreover, some essential comments should be given for those still unknown NLO corrections to the nonfactorizable emission amplitudes and the annihilation amplitudes as follows:

(a) For the nonfactorizable emission diagrams as shown in Fig. 1, since the hard gluons are emitted from the upper quark line of Fig. 1(c) and the upper antiquark line of Fig. 1(d) respectively, the contribution from these two figures will be strongly cancelled each other, the remaining contribution is therefore becoming rather small in magnitude. In NLO level, another suppression factor will appear, the resultant NLO contribution from the hardspectators should become much smaller than the dominant contribution from the ¡±tree¡± emission diagrams (Fig. 1(a) and 1(b)).
Therefore, it is reasonable for us to expect that those still unknown NLO contributions in the PQCD approach are in fact the higher order corrections to the already small LO pieces, and should be much smaller than the dominant contribution for the considered decays, say less than 5% of the dominant ones.
According to Refs. Li05 (); nlo05 (), the vertex corrections can be absorbed into the redefinition of the Wilson coefficients by adding a vertexfunction to them.
(20) 
where denotes the vector meson emitted from the weak vertex ( i.e. the in Fig. 2(a)2(d)). The expressions of the vertexfunctions with both longitudinal and transverse components can be found easily in Refs. qcdfpppv (); nlovc ().
The NLO “QuarkLoop” and “MagneticPenguin” contributions are in fact a kind of penguin corrections with insertion of the fourquark operators and the chromomagnetic operator , respectively, as shown in Figs. 2(e,f) and 2(g,h). For the transition, for example, the corresponding effective Hamiltonian and can be written in the following form:
(21)  
(22) 
where is the invariant mass of the gluon which attaches the quark loops in Figs. 2(e,f), and the functions can be inferred from Refs. Li05 (); nlo05 (); fan2013 (); xiao08b (); xiao2014 (). The in Eq. (22) is an effective Wilson coefficient with the definition of buras96 ().
With explicit evaluations, we find the following three points:

For the pure annihilation decays of , and , they do not receive the NLO contributions from the vertex corrections, the quarkloop and the magneticpenguin diagrams. The NLO correction to these decay modes comes only from the NLOWCs and the the strong coupling constant at the twoloop level.

For the and channels with only transition and no annihilation diagrams, the ”quarkloop” and ”magneticpenguin” diagrams cannot contribute to these two decay modes. The related NLO contributions are mainly induced by the vertex corrections to the emitted or mesons.

For the remaining seven decay modes, besides the LO decay amplitudes, all of the currently known NLO contributions should be taken into account as follows:
(23) where and the terms stand for the LO amplitudes, while and are the NLO ones, which describe the NLO contributions arising from the uploop, charmloop, QCDpenguinloop, and magneticpenguin diagrams, respectively.
Now, we can calculate the decay amplitudes and in the PQCD approach. As mentioned in Eq. (19), for decays, there are three individual polarization amplitudes . For the longitudinal components, the NLO decay amplitudes and can be written as:
(24)  
(25)  
with .
The transverse components and of the corresponding decay amplitudes can be written in the form of
(26)  