Quasitwobody decays in the perturbative QCD approach
Abstract
In this work, we calculate the averaged branching ratios and direct violating asymmetries of the quasitwobody decays by employing the perturbative QCD (PQCD) factorization approach, where is a light pseudoscalar meson , and . The considered decay modes are studied in the quasitwobody framework by parametrizing the twopion distribution amplitude . The wave timelike form factor in the resonant regions associated with the and is estimated based on available experimental data. The PQCD predictions for the averaged branching ratios of the decays are in the order of . The branching ratios of the twobody decays are extracted from the corresponding quasitwobody decay modes. The whole pattern of the squared pion form factor measured by BABAR Collaboration could also be understood based on our studies. The PQCD predictions in this work will be tested by the precise data from the LHCb and the future Belle II experiments.
pacs:
13.25.Hw, 12.38.BxI Introduction
In recent years, prompted by a large number of experimental measurements bfbook (); prd78072006 (); prd79072004 (); lhcb0 (); prl111101801 (); jhep10143 (); prd90112004 (); prl112011801 (); prd95012006 (); hfag2016 (), threebody hadronic meson decays have been studied by using different theoretical frameworks plb56490 (); prd91014029 (); CY16 (); Wang2014a (); Wang2015a (); Wang2016 (); ly15 (); ly16 (); zhou17 (). For such threebody decays, both resonant and nonresonant contributions may appear, as well as the possible final state interactions prd89094013 (); 151209284 (); 89053015 (). The nonresonant contributions have been studied with the method of heavy meson chiral perturbation theory prd461148 (); prd452188 (); plb280287 () in Ref. CY16 (). Meanwhile, the resonant contributions are usually described with the isobar model prd113165 () in terms of the BreitWigner formalism BWmodel (). Based on the QCDimproved factorization prl83 (), such decays have been studied by many authors plb622207 (); prd74114009 (); B.E:2009th (); ST:15 (); CY01 (); prd76094006 (); CY16 (); prd89094007 (); prd87076007 (). By employing the perturbative QCD (PQCD) approach, the decays have also been investigated in Refs. Wang2014a (); Wang2015a (); Wang2016 (); ly15 (); ly16 (); zhou17 (); ma16 (); ma17 (); Chen:2002th (); chen:2004th ().
In the PQCD approach Chen:2002th (); chen:2004th () for the cases of a meson decaying into three final states, we restrict ourselves to the specific kinematical configurations, in which two energetic mesons are almost collimating to each other. The contribution from the region, where there is at least one pair of light mesons having an invariant mass below , being the meson and quark mass difference, as discussed in Refs. Wang2014a (); Wang2015a (); Chen:2002th (); chen:2004th (), is assumed dominant. The final state interactions are expected to be suppressed in such conditions. As a result, the dynamics associated with the meson pair could be factorized into a twomeson distribution amplitude MP (); MT01 (); MT02 (); MT03 (); MN (); Grozin01 (); Grozin02 (). The typical PQCD factorization formula for the decay amplitude can be written in the form of Chen:2002th ()
(1) 
The hard kernel describes the dynamics of the strong and electroweak interactions in the threebody hadronic decays in a similar way as the cases of the twobody decays, and and are the wave functions for the meson and the final state , which absorb the nonperturbative dynamics in the related processes.
In this work, we extend the previous studies Wang2016 (); ly16 () to the decays and in the PQCD approach with the help of the twopion distribution amplitudes , where the stands for the light pseudoscalar mesons, , or . For simplicity, in the following parts of this work, and will be adopted to take the place of and , respectively. The theoretical studies of the excited states will provide us with a deeper understanding of the internal structure of hadrons. For and , there are not many studies except Refs. prd60094020 (); prc79025201 (); prd77116009 (); 12056793 () in the frameworks of the quark model, the large limits, or the doublepole QCD sum rules. For the phenomenological study of the twobody decays and , we still lack the distribution amplitudes of the states and at present. Fortunately, we are allowed to single out the (and ) component according to the twopion distribution amplitudes as has been done in Ref. Wang2016 (). Following Ref. Wang2016 (), we here make an attempt to study the and decays in the quasitwobody framework based on the PQCD factorization approach. And the branching fractions for the twobody decays will be extracted from the quasitwobody processes .
This paper is organized as follows. In Sec. II, we give a brief introduction for the theoretical framework. The numerical values, some discussions, and the conclusions will be given in the last two sections.
Ii Framework
In the lightcone coordinates, the meson momentum , the total momentum of the pion pair, , the momentum of the final state meson , the momentum of the spectator quark in the meson, the momentum for the resonant state , and for the final state are chosen as
(2) 
where is the mass of meson, and the variable is defined as with the invariant mass squared . The parameter denotes the momentum fraction of the positive quark in each meson and runs from zero to unity. , and denote the transverse momentum of the positive quark, respectively. If we choose as one of the pion pair’s momentum fractions, the two pions momenta can be written as
(3) 
The twopion distribution amplitudes can be described in the same way as in Ref. Wang2016 () ,
(4) 
with
(5)  
(6)  
(7) 
where the Legendre polynomial .
After taking the  interference and excitedstate contributions into account, the timelike form factor in Eq. (5) can be written in the following form prd86032013 ():
(8) 
where is the twopion invariant mass square, and () is the decay width (mass) for the relevant resonance . The mass and width for these excited mesons, and the values of the complex parameters and in Eq. (8) can be found in Ref. prd86032013 (). The explicit expressions of the resonant state functions , and can be found, for example, in Ref. ly16 (). In this paper, we only consider the contributions from and . Following Ref. Wang2016 (), we also assume that
(9) 
for the form factors and that appeared in Eqs. (6) and (7). In the numerical calculations, we use the Gegenbauer moments
(10) 
for the twopion distribution amplitudes as used in Ref. ly16 ().
We here use the same wave functions for the and mesons as those in Refs. li2003 (); Xiao:2011tx (). The widely used distribution amplitude is of the form
(11) 
The normalization factor depends on the value of and , which is defined through the normalization relation . We set GeV and GeV li2003 (); Xiao:2011tx () in the numerical calculations. The wave function of the final state pseudoscalar meson (, or ) is of the form
(12) 
where is the chiral mass. The expressions of the relevant distribution amplitudes can also be found, for example, in Refs. ly16 (); ball99 (); ball9901 (); ball05 (); ball06 (); prd76074018 (); ly14 (); prd90114028 ().
The mesons and are considered as the mixtures from and through the relation
(13) 
with the and . We adopt the decay constants and mixing angle as prd58114006 (); plb499339 ()
(14) 
with .
Iii Numerical results
In numerical calculations, we use the following input parameters (in units of GeV except ) pdg2016 ():
(15) 
The values of the Wolfenstein parameters are the same as those given in Ref. pdg2016 (): , , .
For the decay , the differential branching ratio is written as pdg2016 ()
(16) 
with the mean lifetime of meson, and the invariant mass squared. The kinematic variables and denote one of the pion pair’s and ’s momentum in the centerofmass frame of the pion pair,
(17) 
Decay modes  Quasitwobody (in )  Twobody (in )  



Decay modes  Quasitwobody (in )  Twobody (in )  



By using the differential branching fraction, Eq. (16), and the decay amplitudes in the Appendix of Ref. ly16 (), we calculate the averaged branching ratios () and direct violating asymmetries () for the decays and list the results in Table 1. Meanwhile, and for the decays are shown in Table 2. The four errors of these PQCD predictions as listed in Tables 1 and 2 come from the uncertainties of , , , and , respectively.
For and the other three decay modes, the PQCD predictions for their branching ratios as listed in Table I are a little different from those as given previously in Table I of Ref. [16]. The reason is very simple: the Gegenbauer moments used here have been modified slightly from those in Ref. [16] as discussed in Ref. [18].
Taking the quasitwobody decay as an example, the PQCD prediction for its branching ratio and violating asymmetry are the following:
(18)  
(19) 
Here the individual errors as listed in Table 1 have been added in quadrature. Such a PQCD prediction for its branching ratio agrees well with the measured value from Collaboration within errors prd79072006 (). Furthermore, the PQCD prediction for this decay mode is also consistent with the measured value from prd79072006 ().
The width for the process was found to be MeV in Ref. ijmpa135443 (), which is consistent with the value MeV estimated from the annihilation experiments zpc62455 (). The branching fraction could be induced with the GeV zpc62455 (). The branching fraction, on the other hand, could be estimated from the relation EPJC2269 ()
(20) 
where the coupling is fetched from the component of the timelike form factor in Eq. (8) according to at . The decay constant GeV resulting from the data keV zpc62455 () is adopted in this work, which agrees with GeV from the doublepole QCD sum rules 12056793 (), GeV from the perturbative analysis in the large limit prd77116009 (), or GeV from the relativistic constituent quark model prd60094020 (). Utilizing Eq. (20), we find . From the definition of the decay rates between the quasitwobody and the corresponding twobody decay modes
(21) 
we then can find the PQCD predictions for , as listed in the third column of Table 1, where the individual errors have been added in quadrature.
For the cases of the considered quasitwobody and twobody decays involving instead of , in principle, one can obtain the PQCD predictions for the branching ratios and violating asymmetries in a similar way as the case for . But, there is not much reliable information about the properties of the meson except its mass and width ( GeV and GeV) pdg2016 (). What we can do here is to make some rough estimations of the branching ratios and violating asymmetries for the considered decays, and list the PQCD predictions in Table 2. For given keV zpc62455 (), we find the longitudinal decay constant GeV. And then can be obtained by using the same methods as for the decays involving . The errors of the PQCD predictions listed in the third column of Table 2 have been added in quadrature.
In Fig. 1(a), we show the dependence of the differential decay rate after the inclusion of the contributions from the resonant state , , and . One can see that there exists a clear dip near in Fig. 1(a). The position of this dip and the pattern of the whole curve do agree well with Fig. 45 of Ref. prd86032013 (), where the pion form factorsquared measured by are illustrated as a function of [i.e., ] in the region from to GeV. In fact, the differential decay rate does depend on the values of . In Fig. 1(a), we find the prominent peak, a shoulder around the and a clear dip followed by an enhancement (second a little lower and wide peak) in the region. The clear dip at is caused by the destructive interference between the resonant state and . We calculated numerically the interference terms between and amplitudes and found the large negative contribution to the branching ratios.
Taking