CP violation induced by the double resonance for pure annihilation decay process in perturbative QCD

In a perturbative QCD approach we study the direct CP violation in the pure annihilation decay process of B¯s0→π+π-π+π-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\bar{B}^0_{s}\rightarrow \pi ^+\pi ^-\pi ^+\pi ^-$$\end{document} induced by the ρ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho $$\end{document} and ω\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\omega $$\end{document} double resonance effect. Generally, the CP violation is small in the pure annihilation type decay process. However, we find that the CP violation can be enhanced by double ρ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho $$\end{document}–ω\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\omega $$\end{document} interference when the invariant masses of the π+π-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pi ^+\pi ^-$$\end{document} pairs are in the vicinity of the ω\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\omega $$\end{document} resonance. For the decay process of B¯s0→π+π-π+π-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\bar{B}^0_{s}\rightarrow \pi ^+\pi ^-\pi ^+\pi ^-$$\end{document}, the CP violation can reach ACP(B¯s0→π+π-π+π-)=27.20-0.15-0.31-6.11+0.05+0.28+7.13%\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A_{CP}(\bar{B}^0_{s}\rightarrow \pi ^+\pi ^-\pi ^+\pi ^-)=27.20^{+0.05+0.28+7.13}_{-0.15-0.31-6.11}\%$$\end{document}.


Introduction
CP violation is an important area in searching new physics signals beyond the standard model (SM). It is generally believed that the B meson system provides rich information as regards CP violation. Theoretical work has been done in this direction in the past few years. CP violation arises from the weak phase in the Cabibbo-Kobayasgi-Maskawa (CKM) matrix [1,2] in SM. Meanwhile, it is remarkable that CP violation can still be produced by the interference effects between the tree and penguin amplitudes. Since the kinematic suppression, the strong phase associated with long distance rescattering is generally neglected during the past decades. The appearance of the ρ and ω resonance is associated with complex strong phase which is responsible for the CP violation except for the weak phase. The CP violation can be enhanced in the decay process of B ± → π ± π + π − via a ρ-ω mixing mechanism [3]. Recently, the LHCb Collaboration found the large CP violation in the three-body decay channels of B ± → π ± π + π − and B ± → K ± π + π − [4][5][6]. Hence, the nonleptonic B meson decay from the three-body a e-mail: ganglv66@sina.com b e-mail: luye189@163.com and four-body decay channels has been become an important area in searching for CP violation via the ρ and ω resonance.
A mixing between the u and d flavor leads to the breaking of isospin symmetry for the ρ-ω system. The chiral dynamics has been shown to restore the isospin symmetry [7]. The ρω mixing matrix element ρω gives rise to isospin violation which refers to the contribution of the mixing of ω → ρ → 2π . The magnitude has been extracted by the pion form factor through the cross section of e + e − → π + π − from relatively broad ρ resonance region [8][9][10][11][12][13]. The appearance of the ρ and ω resonance is associated with a complex strong phase. Especially, there is perhaps larger strong phase from double ρ and ω interference. The CP violation origins from the weak phase difference and the strong phase difference. Hence, the decay process ofB 0 s → π + π − π + π − is a great candidate for studying the origin of the CP violation.
Meanwhile, it is well known that the CP violation is extremely tiny from the pure annihilation decay process in experiment. There is relatively large error in dealing with the decay amplitudes from the QCD factorization approach [14,15]. The perturbative QCD (PQCD) factorization approach [16][17][18][19][20][21][22][23][24] is based on k T factorization. The amplitude can be divided into the convolution of the Wilson coefficients, the light-cone wave function, and hard kernels by the low energy effective Hamiltonian. The endpoint singularity can be eliminated by introducing the transverse momentum. However, the transverse momentum integration leads to the double logarithm term, which is resummed into the Sudakov form factor. The nonperturbative dynamics is included in the meson wave function which can be extracted from experiment. The perturbative contribution can be calculated by perturbation theory.
The remainder of this paper is organized as follows. In Sect. 2 we present the form of the effective Hamiltonian. In Sect. 3 we give the calculation formalism and details of CP violation from ρ-ω mixing in theB 0 s → ρ 0 (ω)ρ 0 (ω) → π + π − π + π − decay. In Sect. 4 we show the input parame-ters. We present the numerical results in Sect. 5. A summary and discussion are included in Sect. 6. The related functions defined in the text are given in the appendix.

The effective hamiltonian
With the operator product expansion, the effective weak Hamiltonian can be written as [25] where q = (d, s), G F represents Fermi constant, C i (i = 1, . . . , 10) are the Wilson coefficients, V q 1 q 2 (q 1 and q 2 represent quarks) is the CKM matrix element, and O i is the four quark operator. The operators O i have the following forms: where α and β are color indices, and q = u, d, s, c or b quarks. In Eq.

Formalism
The amplitudes A σ of the processB s ( p) where σ is the helicity of the vector meson. 1 ( p 1 ) and 2 ( p 2 ) are the polarization vectors (momenta) of V 1 and V 2 , respectively. m 1 and m 2 refer to the masses of the vector mesons V 1 and V 2 . The invariant amplitudes a, b, c are associated with the amplitude A i (i refers to the three kind of polarizations, longitudinal (L), normal (N) and transverse (T)). Then we have The longitudinal H 0 , transverse H ± of helicity amplitudes can be expressed The interaction of the photon and the hadronic matter can be described by the vector meson dominance model (VMD) [27]. The photon can couple to the hadronic field through a ρ meson. The mixing matrix element ρω (m 2 ρ ) is precisely extracted by fitting to recent data from the cross section for e + e − → π + π − , which does not include the non-resonant contribution of direct ω → π + π − [12,13]. We can express ρω (m 2 ρ ) = Re ρω (m 2 ρ ) + Im ρω (m 2 ρ ) for the real part of Re ρω (m 2 ρ ) and the imaginary part of Im ρω (m 2 ρ ) at the ρ mass.
One can find the mixing parameters is the momentum dependence absorbing the non-resonant contribution of direct decay ω → π + π − . We introduce the momentum dependence of the ρ-ω mixing parameters ρω (s) where √ s is the invariant mass of π + π − . It is sensible to devote one's energies to search the mixing process at the ω mass where the two pions can be produced. We write ρω (s) = Re ρω (m 2 ω ) + Im ρω (m 2 ω ) as follows: Re ρω (m 2 ω ) = −4900 ± 400 MeV 2 , Im ρω (m 2 ω ) = −6500 ± 3100 MeV 2 .
The formalism of the CP violation is presented for thē B 0 s meson decay process in the following. The amplitude A (Ā) for the decay processB 0

can be written as
where H T and H P refer to the tree and penguin operators in the Hamiltonian, respectively. We define the relative magnitudes and phases between the tree and penguin operator contributions as follows: where δ and φ are strong and weak phases, respectively. The weak phase difference φ can be expressed as a combination of the CKM matrix elements: The parameter r is the absolute value of the ratio of tree and penguin amplitudes: The parameter of CP violating asymmetry, A cp , can be written as where and the T i (i = 0, +, −) are the tree-level helicity amplitudes. The r j ( j = 0, +, −) refer to the absolute value of the ratio of tree and penguin amplitude for the three kinds of polarizations, respectively. The sin δ k (k = 0, +, −) represent the relative strong phases between the tree and penguin operator contributions from three kinds of helicity amplitudes. We can see explicitly from Eq. (15) that both weak and strong phase differences are responsible for CP violation. ρ-ω mixing introduces the strong phase difference and is well known for the three-body decay processes of the bottom hadron [28][29][30][31][32][33][34]. Due to ρ-ω interference from u and d quark mixing, we can write the formalism in an approximate form in terms of the first order isospin violation: where t ρρ ( p ρρ ) and t ρω ( p ρω ) are the tree (penguin) amplitudes forB s → ρ 0 ρ 0 andB s → ρ 0 ω, respectively. g ρ is the coupling for ρ 0 → π + π − . ρω refers to the effective ρ-ω mixing amplitude which also effectively includes the direct coupling ω → π + π − . s V , m V and V (V =ρ or ω) is the inverse propagator, mass and decay rate of the vector meson V , respectively.

Calculation details
We can decompose the decay amplitude for the decay pro-cessB 0 s → ρ 0 (ω)ρ 0 (ω) in terms of tree-level and penguinlevel contributions depending on the CKM matrix elements of V ub V * us and V tb V * ts . Due to Eqs. (15), (20) and (21), we can calculate the amplitudes t ρρ , t ρω , p ρρ and p ρω in a perturbative QCD approach. The F and M function associated with the decay amplitudes can be found in the appendix from the perturbative QCD approach.
There are four types of Feynman diagrams contributing toB s → M 2 M 3 (M 2 , M 3 = ρ or ω) annihilation decay mode at leading order. The pure annihilation type process can be classified into factorizable diagrams and non-factorizable diagrams [39,40]. Through calculating these diagrams, we can get the amplitudes A (i) , where i = L , N , T stand for the longitudinal and two transverse polarizations. Because these diagrams are the same as those of B → K * φ and B → K * ρ decays [39,40], the formulas ofB s → ρρ orB s → ρω are similar to those of B → K * φ and B → K * ρ. We just need to replace some corresponding wave functions, Wilson coefficients and corresponding parameters.
With the Hamiltonian (1), depending on CKM matrix elements of V ub V * us and V tb V * ts , the decay amplitudes The tree-level amplitude t ρρ can written as where f B s refers to the decay constant ofB s meson. The penguin-level amplitude are expressed as follows: The decay amplitude forB 0 s → ρ 0 ω can be written as 2 We can give the tree-level contribution as follows: and the penguin-level contributions are given as follows: Based on the definition of (21), we can get where

Input parameters
The CKM matrix, which elements are determined from experiments, can be expressed in terms of the Wolfenstein parameters A, ρ, λ and η [37,38,41]: where O(λ 4 ) corrections are neglected. The latest values for the parameters in the CKM matrix are [42]: wherē The other parameters and the corresponding references are listed in Table 1.

The numerical results of CP violation in
In the numerical results, we find that the CP violation can be enhanced via double ρ-ω mixing for the pure annihilation type decay channelB 0 s → ρ 0 (ω)ρ 0 (ω) → π + π − π + π − when the invariant mass of π + π − is in the vicinity of the ω resonance within perturbative QCD scheme. The CP violation depends on the weak phase difference from CKM matrix elements and the strong phase difference. The CKM matrix elements, which relate to ρ, η, λ and A, are given in Eq. (35). The uncertainties due to the CKM matrix elements come from ρ, η, λ and A. In our numerical calculations, we let ρ, η, λ and A vary among the limiting values. The numerical results are shown from Figs. 1, 2 and 3 with the different parameter values of the CKM matrix elements. The dash line, dot line and solid line corresponds to the maximum, middle, and minimum CKM matrix element for the decay channel of B 0 s → ρ 0 (ω)ρ 0 (ω) → π + π − π + π − , respectively. We find the results are not sensitive to the values of ρ, η, λ and A. In Fig. 1, we give the plot of the CP violating asymmetry as a function of √ s. From Fig. 1, one can see that the CP violation parameter is dependent on √ s and changes rapidly due to ρω mixing when the invariant mass of π + π − is in the vicinity of the ω resonance [see Eq. (15)]. From the numerical results, it is found that the maximum CP violating parameter reaches 27.25% in the case of (ρ mini , η mini ).
From Eqs. (15) and (22), one can see that the CP violating parameter depends on both sin δ and r . The plots of sin δ and r as a function of √ s are shown in Figs. 2, and 3, respectively. It can be seen that sin δ 0 (sin δ − and sin δ + ) vary sharply in the  Fig. 1 The CP violating asymmetry, A cp , as a function of √ s for different CKM matrix elements. The dash line, dot line and solid line corresponds to the maximum, middle, and minimum CKM matrix element for the decay channel ofB 0 s → ρ 0 (ω)ρ 0 (ω) → π + π − π + π − , respectively range of the resonance in Fig. 2. One can see that r changes largely in the vicinity of the ω resonance.

Summary and conclusion
In this paper, we study the CP violation for the pure annihilation type decay process ofB 0 s → π + π − π + π − in perturbative QCD. It has been found that the CP violation can be enhanced greatly in the neighborhood of the ρ-ω resonance. There is the double resonance effect via ρ-ω mixing which can produce large strong phase in this decay process. The maximum CP violation value can reach 27.25% due to double ρ and ω resonance in Fig. 1. Hence, this decay process may be a good channel to test the ρ-ω mixing mechanism. The LHC experiment may search the large CP violation by reconstructing the π + π − pairs at the region of the ρ and ω mass. It is impossible to produce the decay process of s GeV r Fig. 3 Plot of r as a function of √ s corresponding to central parameter values of the CKM matrix elements forB 0 s → ρ 0 (ω)ρ 0 (ω) → π + π − π + π − . The dash line, dot line and solid line corresponds to r 0 , r + and r − , respectively B s → ρ(ω) → π + π − . So we could not obtain large CP violation by this decay process via ρ-ω mixing.
The theoretical errors are large, which follows from the uncertainties of the results. Generally, power corrections beyond the heavy quark limit give the major theoretical uncertainties. This implies the necessity of introducing 1/m b power corrections. Unfortunately, there are many possible 1/m b power suppressed effects and they are generally nonperturbative in nature and hence not calculable by the perturbative method. There are more uncertainties in this scheme. The first error refers to the variation of the CKM parameters, which are given in Eq. (35). The second error comes from the hadronic parameters: the shape parameters, form factors, decay constants, and the wave function of the B s meson. The third error corresponds to the choice of the hard scales, which vary from 0.75 to 1.25t, characterizing the size of next-to-leading order QCD contributions. Therefore, the results for CP violating asymmetry of the decay process B 0 s → π + π − π + π − is given as follows:

Appendix: Related functions defined in the text
In this appendix we present explicit expressions of the factorizable and non-factorizable amplitudes with perturbative QCD in Eqs. (24) and (27) [22][23][24]45,46]. The factorizable amplitudes F L L,i ann (a i ), and F S P,i ann (a i ) (i = L , N , T ) are written as with the color factor C F = 3/4 and f B s refer to the decay constant of theB s meson and a i represents the corresponding Wilson coefficients for annihilation decay channels. In the above functions, The hard scale t is chosen as the maximum of the virtuality of the internal momentum transition in the hard amplitudes, including 1/b i : The hard functions h are written as [47] h e (x 1 , where J 0 and Y 0 are Bessel functions with H (1) The threshold re-sums factor S t follows the parameterized [48] where the parameter c = 0.4. In the non-factorizable contributions, S t (x) has a very small numerical effect on the amplitude [49]. Therefore, we drop S t (x) in h n and h na . The evolution factors E ( ) e and E ( ) a entering in the expressions for the matrix elements are given by in which the Sudakov exponents are defined as where γ q = −α s /π is the anomalous dimension of the quark. The explicit form for the function s(Q, b) is where the variables are defined bŷ and the coefficients A (i) and β i are with n f is the number of the quark flavors and γ E the Euler constant. We will use the one-loop expression of the running coupling constant.
In this study, we use the model function where the share parameter ω b = 0.5 ± 0. uu + dd . The distribution amplitudes of a vector meson (v=ρ or ω), φ ρ , φ ω , φ t V , φ s V , φ v V , and φ a V , are calculated using the light-cone QCD sum rule [50,51]: where t = 2x − 1. Here f V is the decay constant of the vector meson with longitudinal polarization, whose values are shown in Table 1.