Contributions from ${\Phi}_{B2}$ to the $B$ ${\to}$ $PP$ decays within the QCD factorization

TL;DR

This paper investigates the role of the ${\

Contribution

It introduces the ${\

Findings

${\Phi}_{B2}$ affects nonfactorizable annihilation amplitudes.

${\Phi}_{B2}$ contributions are small but may improve decay predictions.

Potential impact on pure annihilation $B$ decays.

Abstract

With the potential for the improvements of measurement precision,the refinement of theoretical calculation on hadronic $B$ weak decays is necessary. In this paper, we study the contributions of $B$ mesonic distribution amplitude ${\Phi}_{B2}$ within the QCD factorization approach, and find that ${\Phi}_{B2}$ contributes to only the nonfactorizable annihilation amplitudes for the $B$ ${\to}$ $PP$ decays ($P$ denotes the ground $SU(3)$ pseudoscalar mesons). Although small, the ${\Phi}_{B2}$ contributions might be helpful for improving the performance of the QCD factorization approach, especially for the pure annihilation $B_{d}$ ${\to}$ $K^{+}K^{-}$ and $B_{s}$ ${\to}$ ${\pi}^{+}{\pi}^{-}$ decays.

Tables

Table 1. Table 1: The input parameters pdg2018 .

$m_b$ $=$ $4.78\pm 0.06$ GeV,	${\overline{m}}_b\left({\overline{m}}_b\right)$ $=$ ${4.18}_{-0.03}^{+0.04}$ GeV,	${\overline{m}}_s\left(2\mathrm{GeV}\right)$ $=$ $95\pm 5$ MeV,
$\frac{{\overline{m}}_s\left(2\mathrm{GeV}\right)}{{\overline{m}}_{u,d}\left(2\mathrm{GeV}\right)}$ $=$ $27.3\pm 0.7$,	$f_{\pi }$ $=$ $130.2\pm 1.7$ MeV,	$f_K$ $=$ $155.6\pm 0.4$ MeV,
$m_{B_d}$ $=$ $5279.63\pm 0.15$ MeV,	$f_{B_d}$ $=$ $187.1\pm 4.2$ MeV,	${\tau }_{B_d}$ $=$ $1.520\pm 0.004$ ps,
$m_{B_s}$ $=$ $5366.89\pm 0.19$ MeV,	$f_{B_s}$ $=$ $227.2\pm 3.4$ MeV,	${\tau }_{B_s}$ $=$ $1.509\pm 0.004$ ps.

Table 2. Table 2: The C ​ P 𝐶 𝑃 CP -averaged branching ratios in the unit of 10 − 7 superscript 10 7 10^{-7} , where the theoretical uncertainties are from the input parameters listed in Table 1 . Different scenarios are explained in the text.

	our results				Ref. npb675.333
decay	S1		S2		S3
mode	$A_k^{i,B2}$ $=$ $0$	$A_k^{i,B2}$ $\ne$ $0$	$A_k^{i,B2}$ $=$ $0$	$A_k^{i,B2}$ $\ne$ $0$	$A_k^{i,B2}$ $=$ $0$	data
$B_s$ $\to$ ${\pi }^{+}{\pi }^{-}$	${3.13}_{-0.43}^{+0.56}$	${5.08}_{-0.86}^{+1.05}$	${3.44}_{-0.47}^{+0.62}$	${5.63}_{-0.97}^{+1.18}$	$1.49$	$6.7\pm 0.8$
$B_d$ $\to$ $K^{+}{K}^{-}$	${0.85}_{-0.14}^{+0.17}$	${1.01}_{-0.16}^{+0.20}$	${0.78}_{-0.13}^{+0.15}$	${0.91}_{-0.15}^{+0.18}$	$0.79$	$0.80\pm 0.15$

Table 3. Table 3: The values of A k i superscript subscript 𝐴 𝑘 𝑖 A_{k}^{i} in Eq.( 17 - 19 ). See text for explanations of different scenarios.

decay		$A_k^{i,B2}$ $=$ $0$			$A_k^{i,B2}$ $\ne$ $0$
mode	scenario	$A_1^i$	$A_2^i$	$A_3^i$	$A_1^i$	$A_2^i$	$A_3^i$
$B_s$ $\to$ ${\pi }^{+}{\pi }^{-}$	S1	$116$	$116$	$0$	$117$	$169$	$-20$
	S2	$103-i\mathrm{\hspace{0.17em}64}$	$103-i\mathrm{\hspace{0.17em}64}$	$0$	$104-i\mathrm{\hspace{0.17em}67}$	$139-i\mathrm{\hspace{0.17em}111}$	$-15+i\mathrm{\hspace{0.17em}17}$
$B_d$ $\to$ $K^{+}{K}^{-}$	S1	$115$	$115$	$0$	$116$	$170$	$-21$
	S2	$102-i\mathrm{\hspace{0.17em}64}$	$102-i\mathrm{\hspace{0.17em}64}$	$0$	$103-i\mathrm{\hspace{0.17em}67}$	$140-i\mathrm{\hspace{0.17em}113}$	$-15+i\mathrm{\hspace{0.17em}18}$

Equations

$${\langle}0{|}\bar{q}_{\alpha}(z)[...]b_{\beta}(0){|}\bar{B}(p){\rangle}\,=\,-\frac{i\,f_{B}}{4}\,\big{\{}\big{(}\!\!\not{p}+m_{b}\big{)}\,{\gamma}_{5}\big{\}}_{{\beta}{\gamma}}{\int}d{\xi}\,e^{-i{\xi}p_{+}z_{-}}\Big{[}{\Phi}_{B1}({\xi})+\not{n}_{-}{\Phi}_{B2}({\xi})\Big{]}_{{\gamma}{\alpha}},$$

$${\int }_0^{1}d\xi {\Phi }_{B1}(\xi )=1,$$

$${\int }_0^{1}d\xi {\Phi }_{B2}(\xi )=0.$$

$$⟨P_1{P}_2|{\hat{O}}_i|\bar{B}⟩$$

$$H_k(P_1,P_2)=H_k^{B1}(P_1,P_2)+H_k^{B2}(P_1,P_2),$$

$$A_k^i=A_k^{i,B1}+A_k^{i,B2},$$

$$A^{B2}=\pi {\alpha }_s{\int }_0^1\xi {\Phi }_{B2}(\xi )d\xi =\pi {\alpha }_s⟨\xi {⟩}_{B2},$$

$$\displaystyle A^{i,B2}_{1}\,=\,-A^{B2}\,{\int}_{0}^{1}\frac{dx}{\bar{x}}{\int}_{0}^{1}\frac{dy}{y}\Big{\{}2\,\frac{{\Phi}_{P2}(x)\,{\Phi}_{P1}(y)}{1-x\,\bar{y}}-r_{\chi}^{P1}\,r_{\chi}^{P2}\,\frac{\bar{x}}{y}\,\frac{{\Phi}_{P2}^{p}(x)\,{\Phi}_{P1}^{p}(y)}{1-x\,\bar{y}}\Big{\}},$$

$$\displaystyle A^{i,B2}_{2}\,=\,A^{B2}\,{\int}_{0}^{1}\frac{dx}{\bar{x}}{\int}_{0}^{1}\frac{dy}{y}\Big{\{}2\,\frac{{\Phi}_{P2}(x)\,{\Phi}_{P1}(y)}{\bar{x}\,y}-r_{\chi}^{P1}\,r_{\chi}^{P2}\,{\Phi}_{P2}^{p}(x)\,{\Phi}_{P1}^{p}(y)\,\Big{[}\frac{\bar{x}}{1-x\,\bar{y}}-\frac{x}{\bar{x}\,y}\Big{]}\Big{\}},$$

$$\displaystyle A^{i,B2}_{3}\,=\,-A^{B2}\,{\int}_{0}^{1}\frac{dx}{\bar{x}}{\int}_{0}^{1}\frac{dy}{y}\Big{\{}2\,r_{\chi}^{P2}\,\frac{x\,{\Phi}_{P2}^{p}(x)\,{\Phi}_{P1}(y)}{1-x\,\bar{y}}+r_{\chi}^{P1}\,{\Phi}_{P2}(x)\,{\Phi}_{P1}^{p}(y)\,\Big{[}\frac{y-\bar{y}}{1-x\,\bar{y}}+\frac{1}{\bar{x}\,y}\Big{]}\Big{\}},$$

$$\mathcal{B}(B_s\to {\pi }^{+}{\pi }^{-})=(6.7\pm 0.8)\times 10^{-7},$$

$$\mathcal{B}(B_d\to K^{+}{K}^{-})=(8.0\pm 1.5)\times 10^{-8}.$$

$${\int }_0^1\frac{du}{u}\to X_A,$$

$${\int }_0^1\frac{du}{u^2}\to X_L,$$

$${\int }_0^{1}du\frac{\ln u}{u}\to -\frac{1}{2}{X}_A^2,$$

$${\int }_0^{1}du\frac{\ln u}{u^2}\to X_L-X_A-X_L{X}_A.$$

$${A}_{1}^{i,B2}\,=\,-A^{B2}\,\Big{\{}12\,({\pi}^{2}-6)-r_{\chi}^{P1}\,r_{\chi}^{P2}\,\Big{[}\frac{{\pi}^{2}}{6}+\frac{1}{2}\,X_{A}^{2}+X_{L}X_{A}+X_{A}-X_{L}\Big{]}\Big{\}},$$

$${A}_{2}^{i,B2}\,=\,A^{B2}\,\Big{\{}72\,(X_{A}-1)^{2}-r_{\chi}^{P1}\,r_{\chi}^{P2}\,\Big{[}\frac{{\pi}^{2}}{6}+\frac{1}{2}\,X_{A}^{2}+X_{L}X_{A}-X_{L}^{2}\Big{]}\Big{\}},$$

$${A}_{3}^{i,B2}\,=\,-A^{B2}\,6\,\Big{\{}\Big{[}\frac{{\pi}^{2}}{6}+\frac{1}{2}\,X_{A}^{2}-X_{A}\Big{]}(2\,r_{\chi}^{P2}-r_{\chi}^{P1})+r_{\chi}^{P1}\,(X_{A}X_{L}-X_{L}+1)\Big{\}}.$$

$$X_A=(1+{\rho }_A{e}^{i{\varphi }_A})\ln \frac{m_B}{{\Lambda }_h},$$

$$X_L=(1+{\rho }_A{e}^{i{\varphi }_A})\frac{m_B}{{\Lambda }_h},$$

$${\langle}{\xi}{\rangle}_{B2}\,=\,\frac{1}{2}\,\Big{(}{\langle}{\xi}{\rangle}_{+}-{\langle}{\xi}{\rangle}_{-}\Big{)}\,=\,\frac{\bar{\Lambda}}{3\,m_{b}}\,,$$

$${\phi}_{B_{q}}^{+}({\xi})\,=\,N^{+}\,{\xi}\,{\exp}\Big{(}-\frac{{\xi}\,m_{B_{q}}}{{\omega}_{B_{q}}}\Big{)},$$

$${\phi}_{B_{q}}^{-}({\xi})\,=\,N^{-}\,{\exp}\Big{(}-\frac{{\xi}\,m_{B_{q}}}{{\omega}_{B_{q}}}\Big{)},$$

$${\cal A}(B_{s}{\to}{\pi}^{+}{\pi}^{-})\,=\,i\,\frac{G_{F}}{\sqrt{2}}\,f_{B_{s}}\,f_{\pi}^{2}\Big{\{}V_{ub}^{\ast}V_{us}\,\Big{(}b_{1}+2\,b_{4}+\frac{1}{2}\,b_{4,{\rm EW}}\Big{)}+V_{cb}^{\ast}V_{cs}\,\Big{(}2\,b_{4}+\frac{1}{2}\,b_{4,{\rm EW}}\Big{)}\Big{\}},$$

$${\cal A}(B_{d}{\to}K^{+}K^{-})\,=\,i\,\frac{G_{F}}{\sqrt{2}}\,f_{B_{d}}\,f_{K}^{2}\Big{\{}V_{ub}^{\ast}V_{ud}\,\Big{(}b_{1}+2\,b_{4}+\frac{1}{2}\,b_{4,{\rm EW}}\Big{)}+V_{cb}^{\ast}V_{cd}\,\Big{(}2\,b_{4}+\frac{1}{2}\,b_{4,{\rm EW}}\Big{)}\Big{\}},$$

$$b_1=\frac{C_F}{N_c^2}{C}_1{A}_1^i,$$

$$b_{4}\,=\,\frac{C_{F}}{N_{c}^{2}}\,\Big{[}C_{4}\,A_{1}^{i}+C_{6}\,A_{2}^{i}\Big{]},$$

$$b_{4,{\rm EW}}\,=\,\frac{C_{F}}{N_{c}^{2}}\,\Big{[}C_{10}\,A_{1}^{i}+C_{8}\,A_{2}^{i}\Big{]},$$

Full text

Because of successive impetus from both experiments and theoretical improvements, the study of nonleptonic $B$ meson weak decays has been one of the hot topics of particle physics. Most of the two-body hadronic $B$ decays with branching ratio larger than $10^{-6}$ have been investigated thoroughly and carefully at the BaBar and Belle experiments pdg2018 ; epjc74.3026 in the past years. A huge amount of $B$ meson experimental data will be accumulated at the high luminosity colliders in the near future, about $50$ $ab^{-1}$ by the Belle-II detector at the $e^{+}e^{-}$ SuperKEKB collider 1808.10567 and about $300$ $fb^{-1}$ by the LHCb Upgrade II detector at the hadron HL-LHC collider 1808.08865 ; 1812.07638 . With the advent of a new age of $B$ physics at the intensity frontier, besides some new phenomena, the unprecedented precision will offer a much more rigorous test on the standard model of elementary particles. The prospective experimental sensitivities for $B$ mesons require more and more accuracy of theoretical calculation.

As is well known, the participation of the strong interactions make it very complicated to calculate the $B$ meson weak decays, especially for the nonleptonic cases. Based on power-counting rules in the heavy quark limits and perturbative QCD theory, some phenomenological models, such as QCD factorization (QCDF) prl83.1914 ; npb591.313 ; npb606.245 ; plb488.46 ; plb509.263 ; prd64.014036 , perturbative QCD (pQCD) approach prd63.074006 ; plb504.6 ; prd63.054008 ; prd63.074009 and so on, have been developed and employed to compute the hadronic matrix elements (HMEs) describing the transformations between the initial $B$ meson and final hadrons through local quark interactions. However, the nonperturbative contributions to HMEs bring theoretical results on branching ratios with many and large uncertainties, particularly for the internal $W$-boson emission and the neutral current processes. To reduce theoretical uncertainties and satisfy the precision requirements of experimental analysis, a careful and comprehensive examination of all possible nonperturbative factors within a phenomenological model is necessary. In this paper, the contributions from the $B$ meson wave functions will be reassessed in detail within the theoretical framework of QCDF.

Wave functions (WFs) or distribution amplitudes (DAs) of the $B$ meson are the essential ingredients of the master formulas in QCDF npb591.313 and pQCD plb504.6 approaches to evaluate the nonfactorizable contributions to HMEs, such as the spectator scattering amplitudes. However, the knowledge of the $B$ mesonic WFs and DAs is still limited so far. It is intuitive that the component quarks of a hadron should move with the same velocity to form a color singlet, and thus the valence quarks would share momentum fractions according to their masses. It is expected that the $B$ mesonic DAs should be very asymmetric with ${\xi}$ at the scales of order $m_{b}$ or smaller, if the light spectator quark carries a longitudinal momentum fraction ${\xi}$ ${\sim}$ ${\cal O}({\Lambda}_{\rm QCD}/m_{b})$, where ${\Lambda}_{\rm QCD}$ and $m_{b}$ are respectively the characteristic QCD scale and the mass of $b$ quark. Generally, the $B$ meson is described by two scalar functions up to the leading power in $1/m_{b}$ prd55.272 ; npb592.3 ; npb625.239 ; npb642.263 , which is written as npb591.313

[TABLE]

where the dots denote the path-ordered exponential gauge factor; the light spectator quark moves along the light-like $z_{-}$ line; $n_{-}$ $=$ $(1,0,0,-1)$ is a null vector; and the normalization conditions of DAs are npb591.313

[TABLE]

[TABLE]

According to the conventions of Refs.prd55.272 ; npb592.3 , ${\Phi}_{B1}$ $=$ ${\phi}_{B}^{+}$ and ${\Phi}_{B2}$ $=$ $({\phi}_{B}^{+}$ $-$ ${\phi}_{B}^{-})/2$. Generally, the two functions ${\phi}_{B}^{\pm}$ are not identical, ${\phi}_{B}^{+}$ ${\neq}$ ${\phi}_{B}^{-}$, and satisfy the relation ${\phi}_{B}^{+}({\xi})$ $+$ ${\xi}\,{\phi}_{B}^{-{\prime}}({\xi})$ $=$ [math] npb592.3 . So, ${\Phi}_{B2}$ ${\neq}$ [math]. The contributions of ${\Phi}_{B2}$ part are suppressed by the power factor of ${\Lambda}_{\rm QCD}/m_{b}$, compared with those of ${\Phi}_{B1}$. In the actual calculations for the $B$ ${\to}$ $PP$ decays with the QCDF approach ($P$ denotes the light $SU(3)$ ground pseudoscalar meson), for example in Ref.npb606.245 , only the contributions from ${\Phi}_{B1}$ part are considered appropriately, while those from ${\Phi}_{B2}$ part are not included explicitly. It should be pointed out that the value of ${\Lambda}_{\rm QCD}/m_{b}$ is not a negligible number, because the mass of the $b$ quark is finite rather than infinite. It has been shown in Refs.npb625.239 ; npb642.263 ; prd71.034018 ; epjc28.515 ; prd74.014027 that there is a large contribution of ${\Phi}_{B2}$ to the hadronic $B$ ${\to}$ ${\pi}$ transition formfactors within the pQCD approach, and its share could reach up to ${\sim}$ $30\%$ with some specific inputs epjc28.515 ; prd74.014027 . This means that the contributions of ${\Phi}_{B2}$ to branching ratios for the $W$ emission processes can reach up to ${\sim}$ $70\%$ for some cases. The ${\Phi}_{B2}$ contribution that were neglected in most cases should be given due attention with the QCDF approach, which is the focus of this paper.

Here, it should be pointed out that a possibly large contribution of ${\Phi}_{B2}$ to formfactors is present only with the pQCD approach rather than the QCDF approach, due to different understandings on the nature of the hadronic transition formfactors. With the pQCD approach prd63.074006 ; plb504.6 ; prd63.054008 ; prd63.074009 , it is assumed that the light quark with a soft momentum of ${\cal O}({\Lambda}_{\rm QCD})$ in the initial $B$ meson should interact with a hard gluon, so it could receive a large boost in order to form a colorless final state with a light energetic quark originating from the $b$ quark decaying interaction point. It is therefore arguable that the hadronic transition formfactors are computable perturbatively with the help of the Sudakov factor regulation on soft contributions. The hadronic transition formfactors are written as the convolution of wave functions of both the $B$ meson and final hadron. Contrarily, it is argued npb591.313 ; npb592.3 with the QCDF approach that the hard and soft contributions to the heavy-to-light formfactors have the same scaling behavior, and the hard contributions are suppressed by one power of ${\alpha}_{s}$ compared with the soft contributions. Because of the dominance of soft contributions, the formfactors for the transition between $B$ meson and light hadron are not fully calculable with the perturbative QCD theory. So, the formfactors are regarded as nonperturbative inputs with the QCDF approach, and therefore have nothing to do with the $B$ mesonic wave functions.

We will concentrate on the $B$ ${\to}$ $PP$ decays for the moment. Up to power corrections of $1/m_{b}$, the general QCDF formula of HMEs for an effective operator $\hat{O}_{i}$ is written as npb591.313 ,

[TABLE]

where $F_{0}^{B{\to}P_{i}}$ denotes the formfactor; $T^{I}$, $H^{I}$ and $T^{II}$ are hard scattering kernels; the mesonic DAs, ${\phi}_{P_{2}}(x)$ and ${\phi}_{P_{1}}(y)$, are the functions of longitudinal momentum fractions $x$ and $y$ of light quarks.

For the two terms in the first line of Eq.(4), soft contributions are assumed to be embodied in the formfactors $F_{0}^{B{\to}P_{i}}$ and DAs. Contributions of $T^{I}$ and $H^{I}$ are dominated by hard gluon exchange. So these contributions, which are irrelevant to $B$ mesonic wave functions, are considered as perturbative corrections to the naive factorization formula, which involve only decay constants and formfactors, but no DAs.

The term in the second line of Eq.(4) corresponds to nonfactorizable contributions. The spectator scattering interactions (see Fig.1) entangle the initial $B$ meson with the final hadrons, which make separating one hadron from others impossible. Therefore, the spectator scattering amplitudes are usually written as the convolution integral of the hard kernels $T^{II}$ and all participating DAs. The hard spectator scattering amplitudes contain the contributions from both ${\Phi}_{B1}$ and ${\Phi}_{B2}$, and can be written as

[TABLE]

where $P_{1}$ is the emitted meson; $P_{2}$ is the recoiled meson that incorporates the spectator quark from $B$ meson into itself; $H_{k}^{B1}$ ($H_{k}^{B2}$) is the contribution from ${\Phi}_{B1}$ (${\Phi}_{B2}$); the subscript $k$ on $H_{k}$ refers to the possible Dirac current structure ${\Gamma}{\otimes}{\Gamma}$ of an operator $\hat{O}$, namely, $k$ $=$ $1$, $2$ and $3$ correspond to ${\Gamma}{\otimes}{\Gamma}$ $=$ $(V-A){\otimes}(V-A)$, $(V-A){\otimes}(V+A)$ and $-2(S-P){\otimes}(S+P)$ respectively. After the straightforward calculation, we find that considering the $SU(3)$ flavor symmetry, the expressions of $H_{1}^{B1}$ and $H_{2}^{B1}$ are entirely consistent with Eq.(47) and Eq.(48) of Ref.npb675.333 , and $H_{3}^{B1}$ $=$ [math]. Our calculations also show that $H_{k}^{B2}$ corresponding to Fig.1(a) and Fig.1(b) are nonzero. Moreover, the terms of both ${\int}_{0}^{1}\frac{{\Phi}_{P1}(y)}{\bar{y}^{2}}dy$ and ${\int}_{0}^{1}\frac{{\Phi}_{P1}^{p}(y)}{\bar{y}^{2}}dy$ appear in $H_{k}^{B2}$, where ${\Phi}_{P1}(y)$ and ${\Phi}_{P1}^{p}(y)$ are the leading twist (twist-2) and twist-3 DAs of the emitted meson $P_{1}$ and $\bar{y}$ $=$ $1$ $-$ $y$. It is clearly seen that with the asymptotic forms of ${\Phi}_{P1}(y)$ $=$ $6\,y\,\bar{y}$ and ${\Phi}_{P1}^{p}(y)$ $=$ $1$, the integrals of ${\int}_{0}^{1}\frac{{\Phi}_{P1}(y)}{\bar{y}^{2}}dy$ and ${\int}_{0}^{1}\frac{{\Phi}_{P1}^{p}(y)}{\bar{y}^{2}}dy$ exhibit logarithmic and linear infrared divergences. Fortunately, because of the opposite sign between the emitted quark and antiquark propagators plus the condition of Eq.(3), the contributions of $H_{k}^{B2}$ exactly cancel each other out. The total contributions from ${\Phi}_{B2}$ to spectator scattering amplitudes are zero.

Compared with the leading contributions, the weak annihilation (WA) contributions are thought to be suppressed by one power of ${\Lambda}_{\rm QCD}/m_{b}$ npb591.313 . However, the WA contributions are significant and can not be ignored in practical application of the QCDF approach to the hadronic $B$ decays npb606.245 ; npb675.333 ; prd65.074001 ; prd65.094025 ; npb774.64 . Therefore, the QCDF master formula of Eq.(4) is generalized to estimate the WA contributions. The WA interactions have two types of topologies within the QCDF approach. The nonfactorizable and factorizable topologies respectively correspond to gluon emission from the initial $B$ meson and final quarks, see Fig.2. The factorizable WA amplitudes can be written as the product of the time-like [math] ${\to}$ $P_{1}P_{2}$ formfactors and the integral of $B$ mesonic WFs, see Fig.2(a) and (b). With the normalization condition of Eq.(3), it is clearly seen that ${\Phi}_{B2}$ contributes nothing to the factorizable WA amplitudes $A^{f}_{k}$, where the superscript $f$ means factorizable, i.e., gluon emission from the final quarks; the subscript $k$ has the same meaning as that of $H_{k}$ in Eq.(5). The nonfactorizable WA amplitudes, corresponding to Fig.2(c) and (d), can be written as the convolution integral of all participating hadronic DAs, and contain the contributions from both ${\Phi}_{B1}$ and ${\Phi}_{B2}$.

[TABLE]

where the superscript $i$ means gluon emission from the initial $B$ meson; $A^{i,B1}_{k}$ ($A^{i,B2}_{k}$) is the contribution from ${\Phi}_{B1}$ (${\Phi}_{B2}$). The expressions of $A^{i,B1}_{k}$ have been explicitly given by Eq.(62) of Ref.npb606.245 and Eq.(54) of Ref.npb675.333 . Here, we will give the new components $A^{i,B2}_{k}$.

[TABLE]

[TABLE]

[TABLE]

[TABLE]

where the factor $r_{\chi}^{P}$ $=$ $\frac{2\,m_{P}^{2}}{\bar{m}_{b}\,(\bar{m}_{q_{1}}+\bar{m}_{q_{2}})}$.

It is easy to find that contributions from ${\Phi}_{B2}$ to the WA amplitudes are nonzero, because the moment parameter ${\langle}{\xi}{\rangle}_{B2}$ is nonzero. Hence, ${\Phi}_{B2}$ may present nontrivial effects on the observables of hadronic $B$ decays, especially for the WA dominant ones.

In order to better investigate the ${\Phi}_{B2}$ contributions and eliminate other pollution, the pure WA decays $B_{d}$ ${\to}$ $K^{+}K^{-}$ and $B_{s}$ ${\to}$ ${\pi}^{+}{\pi}^{-}$ will be restudied in this paper. Although their branching ratios are tiny, they have been measured accurately by now hfag .

[TABLE]

[TABLE]

With the asymptotic twist-2 and -3 DAs, ${\Phi}_{P}(u)$ $=$ $6\,u\,\bar{u}$ and ${\Phi}_{P}^{p}(u)$ $=$ $1$, the integrals in Eq.(8-10) exhibit logarithmic and linear infrared divergences. For an estimation of the WA contributions from ${\Phi}_{B2}$, these divergent endpoint integrals will be parameterized by the commonly used notations within the QCDF approach npb606.245 ; npb675.333 ; npb774.64 .

[TABLE]

[TABLE]

[TABLE]

[TABLE]

The phenomenological parameters $X_{A}$ and $X_{L}$, are usually treated as universal for hadronic $B$ decays in previous literatures npb606.245 ; npb675.333 ; prd65.074001 ; prd65.094025 ; npb774.64 .

With the above parameterization scheme, the WA amplitudes can be rewritten as

[TABLE]

[TABLE]

[TABLE]

The parameters of $X_{A}$ and $X_{L}$ including part of strong phases are complex, and are usually parameterized as npb606.245 ; npb675.333 ; prd65.074001 ; prd65.094025 ; npb774.64

[TABLE]

[TABLE]

where ${\Lambda}_{h}$ $=$ $0.5$ GeV npb606.245 ; npb675.333 , and ${\phi}_{A}$ is an undetermined strong phase. In addition, according to the relations given by Refs.prd55.272 ; npb592.3 , the moment parameter in Eq.(7) is

[TABLE]

with ${\langle}{\xi}{\rangle}_{+}$ $=$ $2\,{\langle}{\xi}{\rangle}_{-}$ $=$ $\frac{4}{3}\,\frac{\bar{\Lambda}}{m_{b}}$ and $\bar{\Lambda}$ $=$ $m_{B}$ $-$ $m_{b}$ ${\approx}$ $0.55$ GeV prd55.272 . Using the exponential type model for $B$ meson DAs

[TABLE]

[TABLE]

where $N^{\pm}$ is the normalization constant determined via ${\int}_{0}^{1}{\phi}_{B_{q}}^{\pm}({\xi})d{\xi}$ $=$ $1$, one can obtain ${\langle}{\xi}{\rangle}_{B2}$ $=$ $0.042{\pm}0.01$ with the shape parameter ${\omega}_{B_{s}}$ $=$ $0.45{\pm}0.10$ GeV for $B_{s}$ meson epjc73.2437 , and ${\langle}{\xi}{\rangle}_{B2}$ $=$ $0.039{\pm}0.01$ with ${\omega}_{B_{d}}$ $=$ $0.42{\pm}0.10$ GeV for $B_{d}$ meson prd74.014027 , which are basically in agreement with the estimation of Eq.(22).

Using the commonly used notations in the QCDF approach npb606.245 ; npb675.333 ; prd65.074001 ; prd65.094025 ; npb774.64 , the amplitudes for the pure WA decays $B_{d}$ ${\to}$ $K^{+}K^{-}$ and $B_{s}$ ${\to}$ ${\pi}^{+}{\pi}^{-}$ are written as

[TABLE]

[TABLE]

where the Fermi weak coupling constant $G_{F}$ ${\simeq}$ $1.166{\times}10^{-5}$ ${\rm GeV}^{-2}$ pdg2018 ; $f_{B_{q}}$, $f_{\pi}$ and $f_{K}$ are decay constants; $V_{ij}$ ($i$ $=$ $u$, $c$ and $j$ $=$ $d$, $s$, $b$) is the Cabibbo-Kobayashi-Maskawa (CKM) matrix element. The definition of parameter $b_{i}$ is

[TABLE]

[TABLE]

[TABLE]

where $C_{F}$ $=$ $4/3$ is the color factor; $N_{c}$ $=$ $3$ is the number of colors; $C_{i}$ is the Wilson coefficient; $A^{i}_{k}$ is the amplitude building block of Eq.(6).

To provide a quantitative estimate of the ${\Phi}_{B2}$ contributions, the inputs listed in Table.1 are used in our numerical calculation. Their central values will be regarded as the default inputs unless otherwise specified.

The constraints on annihilation parameters from data are illustrated in Fig.3. It is clearly seen from Fig.3(a) that it is impossible to accommodate simultaneously $B_{d}$ ${\to}$ $K^{+}K^{-}$ and $B_{s}$ ${\to}$ ${\pi}^{+}{\pi}^{-}$ decays within $2\sigma$ errors with the same values of ${\rho}_{A}$ and ${\phi}_{A}$ when the ${\Phi}_{B2}$ contributions are overlooked. Other studies of $B$ decays, such as Refs.npb675.333 ; prd88.014043 , have uncovered similar results. It seems not easy to clarify discrepancies between data and the QCDF results with the same set of parameters ${\rho}_{A}$ and ${\phi}_{A}$. To clam down this situation, the factorizable and nonfactorizable annihilation parameters corresponding to different topologies are introduced in Refs.prd90.054019 ; plb740.56 . However, more annihilation parameters make the method uneconomical and unsatisfactory. Interestingly, by including the ${\Phi}_{B2}$ contributions, Fig.3(b) shows overlapping areas of annihilation parameters, which implies that the ${\Phi}_{B2}$ contributions are nontrivial for accommodating the tension between data and QCDF predictions for ${\cal B}(B_{d}{\to}K^{+}K^{-})$ and ${\cal B}(B_{s}{\to}{\pi}^{+}{\pi}^{-})$. In addition, if theoretical uncertainties from inputs are taken into account, the overlapping bands will be inevitably enlarged. The same annihilation parameters suitable for pure WA hadronic $B$ decays might be obtained with the QCDF approach.

As is shown by Fig.3(b), strict limits on annihilation parameters ${\rho}_{A}$ and ${\phi}_{A}$ can not be obtained only from experimental data on ${\cal B}(B_{d}{\to}K^{+}K^{-})$ and ${\cal B}(B_{s}{\to}{\pi}^{+}{\pi}^{-})$. In principle, considering more $B$ decays, such as a global fit on nonleptonic $B$ decays in Refs.prd90.054019 ; plb740.56 , is helpful for extracting the informations of annihilation parameters. However, for many hadronic $B$ decays, other contributions, such as spectator scattering interactions, will complicate the determination of annihilation parameters. How to get annihilation parameter spaces as compact as possible from data is beyond the scope of this paper.

It is seen from Fig.3(b) that, in general, the value of ${\rho}_{A}$ increase with the increasing value of ${|}{\phi}_{A}{|}$. A large value of parameter ${\rho}_{A}$ will spoil the self-consistency and confidence level of the QCDF approach, and ${\rho}_{A}$ ${\leq}$ $1$ is proposed in Refs.npb606.245 ; npb675.333 . The strong phase ${\phi}_{A}$ describes the rescattering among hadrons and relates closely to $CP$ violation of nonleptonic $B$ decays. Focusing on the pure WA decays of $B_{d}$ ${\to}$ $K^{+}K^{-}$ and $B_{s}$ ${\to}$ ${\pi}^{+}{\pi}^{-}$, to roughly estimate branching ratios, two scenarios based on Fig.3(b) are considered in our numerical calculation. Scenario S1 is with parameters ${\rho}_{A}$ $=$ $1$ and ${\phi}_{A}$ $=$ $0^{\circ}$, and scenario S2 is with ${\rho}_{A}$ $=$ $1.2$ and ${\phi}_{A}$ $=$ $-40^{\circ}$. Practically, for the scenario S1, it is intuitive that zero strong phase ${\phi}_{A}$ seems a little unnatural. Trying to combine the value of ${\rho}_{A}$ as close to one as possible with a nonzero ${\phi}_{A}$, the scenario S2 is considered. In addition, the scenario S2 is comparable with the scenario S3 of Ref.npb675.333 , where the “universal annihilation” parameters ${\rho}_{A}$ $=$ $1$ and ${\phi}_{A}$ $=$ $-45^{\circ}$ are used.

Using such inputs, we list the QCDF results for ${\cal B}(B_{d}{\to}K^{+}K^{-})$ and ${\cal B}(B_{s}{\to}{\pi}^{+}{\pi}^{-})$ with and without considering the ${\Phi}_{B2}$ contributions in Table 2, in which the theoretical predictions of scenario S3 of Ref.npb675.333 and experimental data are also listed for convenience of comparison. In order to show the effects of ${\Phi}_{B2}$ much more clearly, we collect the numerical results of $A_{i}^{k}$ in Table 3.

From Table 2, it can be found that: (i) The experimental data for both $B_{d}$ and $B_{s}$ decays can not be well explained simultaneously by QCDF approach without considering the ${\Phi}_{B2}$ contributions; (ii) The numerical difference between the case for $A_{k}^{i,B2}$ $=$ [math] of scenario S2 and scenario S3 of Ref.npb675.333 arises from different inputs, such as decay constants, the CKM parameters and so on, besides parameters ${\rho}_{A}$ and ${\phi}_{A}$. (iii) With the scenario S2, the ${\Phi}_{B2}$ contributions present about $60\%$ and $20\%$ corrections to ${\cal B}(B_{d}{\to}K^{+}K^{-})$ and ${\cal B}(B_{s}{\to}{\pi}^{+}{\pi}^{-})$, respectively, which significantly improve the QCDF predictions and can explain the data within uncertainty.

The results in Table 2 show that ${\Phi}_{B2}$ contributions to nonfactorizable WA amplitude building blocks $A^{i}_{k}$ are small, due to the small moment ${\langle}{\xi}{\rangle}_{B2}$. In addition, according to the conventions of Refs.npb606.245 ; npb675.333 , building block $A^{i}_{3}$ is always accompanied by the small value of Wilson coefficient $C_{5}$. Hence, on one hand, the dominant contributions to WA amplitudes come from ${\Phi}_{B1}$ part; on the other hand, to some certain extent, the ${\Phi}_{B2}$ contributions present un-negligible correction to the amplitude especially for the pure annihilation decay modes and can improve the performances of the QCDF approach.

In summary, the improvements of measurement precision with the running Belle-II and LHCb experiments call for the refinements of theoretical calculation on hadronic $B$ weak decays. For the $B$ mesons, there are two scalar DAs ${\Phi}_{B1}$ and ${\Phi}_{B2}$. The ${\Phi}_{B2}$ contributions to formfactors and branching ratios can be significant for some cases with the pQCD approach. In this paper, we study the ${\Phi}_{B2}$ contributions with the QCDF approach, and find that for the $B$ ${\to}$ $PP$ decays, they can be safely neglected in the spectator scattering amplitudes, and contribute to only the nonfactorizable WA amplitudes. The ${\Phi}_{B2}$ contributions to WA amplitudes are small compared with the dominant ${\Phi}_{B1}$ contributions, due to the small moment ${\langle}{\xi}{\rangle}_{B2}$. However, the participation of ${\Phi}_{B2}$ plays a positive role in accommodating the pure WA decays $B_{d}$ ${\to}$ $K^{+}K^{-}$ and $B_{s}$ ${\to}$ ${\pi}^{+}{\pi}^{-}$ to data with the universal annihilation parameters ${\rho}_{A}$ and ${\phi}_{A}$. The values of annihilation parameters ${\rho}_{A}$ and ${\phi}_{A}$ with the QCDF approach have been under discussion for a long period. More information about WA parameters ${\rho}_{A}$ and ${\phi}_{A}$ could be obtained by a comprehensive study on nonleptonic $B$ decays.

Acknowledgments

This work is supported by the National Natural Science Foundation of China (Grant Nos. 11875122, 11705047, U1632109 and 1191101296) and the Program for Innovative Research Team in University of Henan Province (Grant No.19IRTSTHN018).