The decays $B\to \Psi(2S)\pi(K),\eta_c(2S)\pi(K)$ in the pQCD approach beyond the leading-order

TL;DR

This paper calculates branching ratios for B meson decays involving excited charmonium states using pQCD, including NLO effects, and compares predictions with experimental data, highlighting areas for future tests and potential new physics signals.

Contribution

It provides next-to-leading-order pQCD predictions for B decays to excited charmonium states, incorporating vertex corrections and Wilson coefficients, and compares results with experimental data.

Findings

Predicted branching ratios align with current data when including NLO effects.

Decays involving K mesons have larger predicted branching ratios than those with pions.

No significant direct CP violation observed in these decay modes.

Abstract

Two body $B$ meson decays involving the radially excited meson $\psi(2S)/\eta_c(2S)$ in the final states are studied by using the perturbative QCD (pQCD) approach. We find that: (a) The branching ratios for the decays involving $K$ meson are predicted as $Br(B^+\to\psi(2S)K^+)=(5.37^{+1.85}_{-2.22})\times10^{-4}, Br(B^0\to\psi(2S)K^0)=(4.98^{+1.71}_{-2.06})\times10^{-4}, Br(B^+\to\eta_c(2S)K^+)=(3.54^{+3.18}_{-3.09})\times10^{-4}$, which are consistent well with the present data when including the next-to-leading-order (NLO) effects. Here the NLO effects are from the vertex corrections and the NLO Wilson coefficients. The large errors in the decay $B^+\to\eta_c(2S)K^+$ are mainly induced by using the decay constant $f_{\eta_c(2S)}=0.243^{+0.079}_{-0.111}$ GeV with large uncertainties. (b) While there seems to be some room left for other higher order corrections or the non-perturbative long distance contributions in the decays involving $\pi$ meson, $Br(B^+\to\psi(2S)\pi^+)=(1.17^{+0.42}_{-0.50})\times10^{-5}, Br(B^0\to\psi(2S)\pi^0)=0.54^{+0.20}_{-0.23}\times10^{-5}$, which are smaller than the present data. The results for other decays can be tested at the running LHCb and forthcoming Super-B experiments. (c) There is no obvious evidence of the direct CP violation being seen in the decays $B\to \psi(2S)\pi(K), \eta_c(2S)\pi(K)$ in the present experiments, which is supported by our calculations. If a few percent value is confirmed in the future , it would indicate new physics definitely.

Tables

Table 1. Table 1: Our LO and NLO predictions of the branching ratios for the decays B → ψ ​ ( η c ) ​ K , ψ ​ ( η c ) ​ π → 𝐵 𝜓 subscript 𝜂 𝑐 𝐾 𝜓 subscript 𝜂 𝑐 𝜋 B\to\psi(\eta_{c})K,\psi(\eta_{c})\pi . The fist uncertainty comes from the ω b = 0.4 ± 0.1 subscript 𝜔 𝑏 plus-or-minus 0.4 0.1 \omega_{b}=0.4\pm 0.1 for B 𝐵 B meson, the second and the third uncertainties are from the free parameter ω = 0.2 ± 0.1 𝜔 plus-or-minus 0.2 0.1 \omega=0.2\pm 0.1 and the c 𝑐 c quark mass 1.275 ± 0.025 plus-or-minus 1.275 0.025 1.275\pm 0.025 GeV for ψ ​ ( 2 ​ S ) / η c ​ ( 2 ​ S ) 𝜓 2 𝑆 subscript 𝜂 𝑐 2 𝑆 \psi(2S)/\eta_{c}(2S) meson. The last one comes from the decay constant f ψ ​ ( 2 ​ S ) = 0.296 − 0.002 + 0.003 subscript 𝑓 𝜓 2 𝑆 subscript superscript 0.296 0.003 0.002 f_{\psi(2S)}=0.296^{+0.003}_{-0.002} GeV/ f η c ​ ( 2 ​ S ) = 0.243 − 0.111 + 0.079 subscript 𝑓 subscript 𝜂 𝑐 2 𝑆 subscript superscript 0.243 0.079 0.111 f_{\eta_{c}(2S)}=0.243^{+0.079}_{-0.111} GeV. The data are listed in the second column pdg14 .

mode	data	LO	NLO
$B^{+}\to \psi \left(2S\right)K^{+}(\times {10}^{-4})$	$6.27\pm 0.24$	${2.39}_{-0.81-0.31-0.03-0.03}^{+1.11+0.18+0.04+0.05}$	${5.37}_{-1.17-1.89-0.03-0.07}^{+1.61+0.90+0.04+0.11}$
$B^0\to \psi \left(2S\right)K^0(\times {10}^{-4})$	$6.2\pm 0.5$	${2.22}_{-0.67-0.29-0.03-0.02}^{+1.00+0.16+0.03+0.04}$	${4.98}_{-1.09-1.75-0.03-0.06}^{+1.49+0.84+0.04+0.10}$
$B^{+}\to \psi \left(2S\right){\pi }^{+}(\times {10}^{-5})$	$2.44\pm 0.30$	${0.47}_{-0.15-0.04-0.00-0.00}^{+0.24+0.04+0.01+0.01}$	${1.17}_{-0.25-0.43-0.01-0.01}^{+0.36+0.22+0.02+0.03}$
$B^0\to \psi \left(2S\right){\pi }^0(\times {10}^{-5})$	$1.17\pm 0.17\pm 0.08$	${0.22}_{-0.07-0.02-0.00-0.00}^{+0.11+0.03+0.01+0.02}$	${0.54}_{-0.12-0.20-0.01-0.00}^{+0.17+0.10+0.01+0.02}$
$B^{+}\to {\eta }_c\left(2S\right)K^{+}(\times {10}^{-4})$	$3.4\pm 1.8$	${2.33}_{-0.69-0.17-0.02-1.64}^{+1.02+0.23+0.02+1.76}$	${3.54}_{-0.87-1.62-0.04-2.49}^{+1.23+1.18+0.05+2.68}$
$B^0\to {\eta }_c\left(2S\right)K^0(\times {10}^{-4})$	$--$	${2.16}_{-0.64-0.16-0.02-1.52}^{+0.95+0.21+0.02+1.63}$	${3.29}_{-0.81-1.51-0.05-2.32}^{+1.13+1.08+0.04+2.48}$
$B^{+}\to {\eta }_c\left(2S\right){\pi }^{+}(\times {10}^{-6})$	$--$	${5.82}_{-1.69-0.49-0.06-4.10}^{+2.51+0.70+0.06+4.40}$	${9.03}_{-2.20-2.56-0.11-6.37}^{+3.05+3.08+0.11+6.82}$
$B^0\to {\eta }_c\left(2S\right){\pi }^0(\times {10}^{-6})$	$--$	${2.72}_{-0.79-0.23-0.02-1.92}^{+1.18+0.33+0.03+2.06}$	${4.19}_{-1.02-1.97-0.05-2.96}^{+1.41+1.43+0.06+3.16}$

Equations

$$Br(B^{+}\to {\eta }_c(2S)K^{+})=(3.4\pm 1.8)\times 10^{-4},$$

$$Br(B^{+}\to \psi (2S)K^{+})=(6.27\pm 0.24)\times 10^{-4},$$

$$Br(B^{+}\to \psi (2S){\pi }^{+})=(2.44\pm 0.30)\times 10^{-5},$$

$$Br(B^0\to \psi (2S)K^0)=(6.2\pm 0.5)\times 10^{-4}.$$

$$A_{CP}(B^{+}\to \psi (2S)K^{+})$$

$$A_{CP}(B^{+}\to \psi (2S){\pi }^{+})$$

$$\ln \frac{m_B}{\Lambda }(1+{\rho }_H{e}^{i{\delta }_H}),\ln \frac{m_B}{\Lambda }(1+{\rho }_A{e}^{i{\delta }_A}),$$

$${\text{H}}_{eff}=\frac{G_F}{\sqrt{2}}\left[V_{cb}^{\ast }{V}_{c{q}^{\prime }}(C_1(\mu )O_1^c(\mu )+C_2(\mu )O_2^c(\mu ))-V_{tb}^{\ast }{V}_{t{q}^{\prime }}\sum _{i=3}^{10}{C}_i(\mu )O_i(\mu )\right],$$

$$O_1^c$$

$$O_3$$

$$O_5$$

$$O_7$$

$$O_9$$

$$F_{\psi P}^{V-A}$$

$$E_e(t)={\alpha }_s(t)\exp [-S_B(t)-S_{\pi }(t)],S_t(x)=\frac{2^{1+2c}\Gamma (3/2+c)}{\sqrt{\pi }\Gamma (1+c)}[x(1-x){]}^c,$$

$$h_e(x_1,x_3,b_1,b_3)$$

$$M_{\psi P}^{V-A}$$

$$E_d(t)$$

$$h_d(x_1,x_2,x_3,b_1,b_2)$$

$$M_{{\eta }_{c}P}^{V-A}$$

$$\text{M}(B\to \psi ({\eta }_c)\pi )$$

$$a_2$$

$$a_i$$

$$f_I$$

$$f_B$$

$$M_B$$

$$M_W$$

$$\Gamma (B\to \psi ({\eta }_c)P)=\frac{G_F^2}{32\pi m_B}(1-r_{\psi ({\eta }_c)}^2)|M(B\to \psi ({\eta }_c)P){|}^2.$$

$${\psi }^{L,v}(x,b)$$

$${\psi }^t(x,b)$$

$$\mathcal{T}(x)=1-4b^2{m}_c\omega x(1-x)+\frac{m_c(2x-1{)}^2}{\omega x(1-x)},$$

$$A_{CP}(B^{+}\to \psi (2S){\pi }^{+})$$

$$A_{CP}^{dir}=\frac{|{\lambda }_{CP}{|}^2-1}{1+|{\lambda }_{CP}{|}^2},A_{CP}^{mix}=\frac{2Im({\lambda }_{CP})}{1+|{\lambda }_{CP}{|}^2},$$

$${\lambda }_{CP}={\eta }_f\frac{V_{tb}^{\ast }{V}_{td}⟨f|H_{eff}|{\bar{B}}^0⟩}{V_{tb}{V}_{td}^{\ast }⟨f|H_{eff}|B^0⟩}={\eta }_f{e}^{-2i\alpha (\beta )}\frac{⟨f|H_{eff}|{\bar{B}}^0⟩}{⟨f|H_{eff}|B^0⟩},$$

Full text

I Introduction

It is well known that $\eta_{c}(2S)$ and $\Psi(2S)$ are the first radially excited states of the S-wave ground states $\eta_{c}(1S)$ and $J/\Psi$, respectively. These two excited charmonia states have been observed in the $B$ meson decays in the experiments pdg14 ,

[TABLE]

Furthermore, the direct CP-violating asymmetries of the two charged decay channels are given by PDG pdg14 , though with large uncertainties:

[TABLE]

In fact, $B$ meson exclusive decays into charmonia have been received a lot of attentions for many years. They are regarded as the golden channels in researching CP violation and exploring new physics. At the same time, they play the important roles in testing the unitarity of the Cabibbo-Kobayashi-Maskawa (CKM) triangle. Moreover, these decays are ideal modes to check the different factorization approaches, such as the naive factorization assumption (FA) ali , QCD-improved factorization (QCDF) beneke ; beneke2 , light-cone sum rules (LCSR) melic , and the perturbtive QCD (pQCD) approach chen . Most of these approaches can not work well when describing these decays, such as $B\to(J/\Psi,\eta_{c},\chi_{c0},\chi_{c1})K^{(*)}$. Because their predictions about the branching ratios are too small to explain the data. For example, the branching ratios of the decays $B^{+}\to J/\Psi K^{+}$ and $B^{0}\to J/\Psi K^{0}$ measured by Babar were about $(10.61\pm 0.50)\times 10^{-4}$ and $(8.69\pm 0.37)\times 10^{-4}$ babar , respectively, which imply a larger Wilson coefficient $a_{2}\approx 0.20\sim 0.30$ cheng . So the naive factorization assumption does not work with the Wilson coefficient $a_{2}\sim 0$. As for the QCDF approach, it has been found that the end point singularities exist in the amplitudes from the hard spectator scattering and annihilation diagrams, where the contributions are parametrized into

[TABLE]

respectively. These non-universal and uncontrollable parameters $\rho_{H,A}$ and $\delta_{H,A}$ will induce large theoretical uncertainties. What’s worse, when the emitted meson is heavy, such as $D$ and charmonia states $J/\Psi,\Psi(2S),\eta_{c}(2S)$, the QCDF will break down. For example, in the $\bar{B}^{0}\to D^{0}\pi^{0}$ decay, since the $D^{0}$ meson is not a compact object with small transverse extension, it will strongly interact with the $B\pi$ system, which makes the factorization fail. Fortunately, the transverse size of the charmonium is small in the heavy quark limit, so the authors of Ref.chay considered that the QCDF method could be used to the decay $B\to J/\Psi K$, while they found that the leading-twist (twist-2) contributions were too small to explain the data. Then the authors of Ref.cheng calculated the twist-3 contributions, where the divergent integral was parametrized as Eq.(7). But they still could not explain the experimental data. The end-point singularities were also found in the twist-3 amplitudes for the decays $B\to\eta_{c}(1S)K,\eta_{c}(2S)K$ song , where the QCDF approach was used. The LCSR approach was also insufficient to account for the data of these B meson decays into charmonia melic .

While under the pQCD approach, the spectator quark and other quarks are connected by one hard gluon (shown in Fig.1). Unlike the QCDF approach, the hard part of the pQCD approach consists of six quarks rather than four. Certainly, there also exist the soft and collinear gluon exchanges between quarks. So the double logarithms $\ln^{2}{Pb}$ will arise from the overlap of the soft and collinear divergences in radiative corrections to the meson wave functions, $P$ being the dominant light-cone component of a meson momentum, $b$ being the conjugate variable of parton transverse momentum $k_{T}$. One can use the $k_{T}$ resummation lihn to organize these leading double logarithms for all loop diagrams into a Sudakov factor, which suppresses the long distance contributions in the large $b$ region. When the end-point region with a momentum fraction $x\to 0,1$ is important for the hard amplitude, the corresponding large double logarithms $\alpha_{s}\ln^{2}x$ shall appear in the hard amplitude. One can use the threshold resummation lihn2 to organize this type of double logarithms for all loop diagrams into a jet function, which suppresses the end-point behavior of the hard amplitude. With the Sudakov factor and the jet function, one can evaluate all possible Feynman diagrams for the six-quark amplitude directly, including the nonfactorizable emission diagrams and annihilation type diagrams. But it is difficult to calculate these two kinds of contributions in QCDF approach. The pQCD approach has been used to B meson decays into charmonia such as $B\to J/\Psi K^{(*)}$ in Refs.chen ; li0 , where the consistent results with the experimental data were obtained.

So we would like to use the pQCD approach to study $B\to\psi(2S)P,\eta_{c}(2S)P$ (here $P$ refers to the pseudo-scalar meson $\pi$ or $K$) decays. Except the full leading-order (LO) contributions, the next-to-leading-order (NLO) contributions, which are mainly from the NLO Wilson coefficients and the vertex corrections to the hard kernel, are also included. Certainly, other NLO contributions, such as the quark loops and the magnetic penguin corrections, are available in the literatures li ; xiao , but they will not contribute to these considered decays.

We review the LO predictions for the $B\to\psi(2S)P,\eta_{c}(2S)P$ decays including those for the main NLO contributions, in Sec. II. We perform the numerical study in Section IV, where the theoretical uncertainties are also considered. Section V is the conclusion. It is noticed that we will use the abbreviation $\psi$ and $\eta_{c}$ to denote the mesons $\psi(2S)$ and $\eta_{c}(2S)$ unless specified in the following.

II the Leading-Order Predictions and the main next-to-leading order corrections

As previously stated, the pQCD factorization approach has been used to calculate many two-body charmed B meson decays, and has obtained consistent results with experiments. So we use this approach to consider the decays $B\to\psi(2S)\pi(K),\eta_{c}(2S)\pi(K)$, and the corresponding effective Hamiltonian can be written as:

[TABLE]

where $C_{i}(\mu)$ are Wilson coefficients at the renormalization scale $\mu$, $q^{\prime}=d$ ($q^{\prime}=s$) for $b\to d$ ($b\to s$) transition. $V$ represents for the Cabibbo-Kobayashi-Maskawa (CKM) matrix element, and the four fermion operators $O_{i}$ are given as:

[TABLE]

with $i,j$ being the color indices.

At the leading order, the relevant Feynman diagrams only include the factorizable and non-factorizable emission diagrams, as shown in Fig.1, where we take the decays $B^{0}\to\psi(\eta_{c})\pi^{0}$ as examples. If the emission particle is the vector meson $\psi(2S)$, then the amplitude for the factorizable emission diagrams Fig.1(a) and Fig.1(b) can be written as:

[TABLE]

where $r_{\psi}=m_{\psi}/m_{B},r_{P}=m^{P}_{0}/m_{B}(P=\pi,K$), the color factor $C_{F}=4/3$, and $\phi^{A(P,T)}_{P(P=\pi,K)}$ are the twist-2 (twist-3) distribution amplitudes for the meson $\pi$ or $K$. $x_{1}$ and $x_{3}$ are the light quark momentum fractions in the $B$ and $P$ mesons, respectively. The evolution factors evolving the Sudakov factors and the jet function $S_{t}(x)$ are listed as:

[TABLE]

with $c=0.3$, and $S_{B}(t),S_{\pi}(t)$ being the Sudakov factors, which can be found in Ref.zhang . The hard functions $h_{e}$ is given as:

[TABLE]

The amplitude for the non-factorizable spectator diagrams Fig.1(c) and Fig.1(d) is given as:

[TABLE]

where the twist-2(twist-3) distribution amplitudes $\psi^{L(t)}(x_{2},b_{2})$ can be found in the next section. $x_{2}$ is the $c$ quark momentum fraction in $\psi(2S)$ meson. The evolution factor $E_{d}(t)$ and the hard function $h_{d}$ are listed as:

[TABLE]

with the variable $A^{2}_{d}=r_{c}^{2}+(x_{1}-x_{2})(x_{2}r^{2}_{\psi}+x_{3}(1-r^{2}_{\psi}))$, and $K_{0},I_{0}$ and $H_{0}$ being the modified Bessel functions.

If the emission particle is the pseudo-scalar meson $\eta_{c}(2S)$, then the corresponding amplitude $F^{V-A}_{\eta_{c}P}$ can be obtained from $F^{V-A}_{\psi P}$ by the replacements of the parameters $f_{\psi}$ and $r_{\psi}$ with $f_{\eta_{c}}$ and $r_{\eta_{c}}$, respectively. While there are many differences between the nonfactorizable spectator amplitudes $M^{V-A}_{\eta_{c}P}$ and $M^{V-A}_{\psi P}$ because of the different Lorentz structures between the wave functions of $\eta_{c}(2S)$ and $\psi(2S)$. Here $M^{V-A}_{\eta_{c}P}$ is listed as following:

[TABLE]

where $r_{\eta_{c}}=m_{\eta_{c}}/m_{B}$, and the twist-3 distribution amplitudes of $\eta_{c}$ meson do not contribute to the amplitude. It is different from the case of $M^{V-A}_{\psi P}$. While the evolution factor $E_{d}(t)$ and the hard function $h_{d}$ are similar with those given in Eq.(18) and Eq.(19).

By combining the amplitudes from the different Feynman diagrams, one can get the total decay amplitude for the decays $B\to\psi(\eta_{c})\pi$:

[TABLE]

where $V_{ij}$ is the CKM matrix element and the combinations of Wilson coefficients $a_{2}=C_{1}+C_{2}/3,a_{i}=C_{i}+C_{i+1}/3$ with $i=3,5,7,9$. The amplitudes $F^{V-A}_{\psi(\eta_{c})\pi}$ and $M^{V-A}_{\psi(\eta_{c})\pi}$ are given in Eqs.(14) and (17), respectively. As for the decays $B\to\psi(\eta_{c})K$, the total amplitude can be obtained by replacing $F^{V-A}_{\psi(\eta_{c})\pi},M^{V-A}_{\psi(\eta_{c})\pi},V_{cd}$ and $V_{td}$ with $F^{V-A}_{\psi(\eta_{c})K},M^{V-A}_{\psi(\eta_{c})K},V_{cs}$ and $V_{ts}$, respectively in Eq.(21).

As stated before, the NLO corrections to the hard kernel for the decays $B\to\psi(2S)P$ and $B\to\eta_{c}(2S)P$ are simpler compared with other B meson decays such as $B\to\pi K,\rho K$. Here only the vertex corrections are need to be considered. Since the vertex corrections can reduce the dependence of the Wilson coefficients on the renormalization scale $\mu$, they play the important roles in the NLO analysis. It is well known that the nonfactorizable contributions are small li1 , we concentrate only on the vertex corrections to the factorizable amplitudes, as shown in Fig.2. Furthermore, the infrared divergences from the soft and collinear gluons in these Feynman diagrams can be canceled each other. That is to say, these corrections are free from the end-point singularity in collinear factorization theorem, so we can simply quote the QCDF expressions for the vertex corrections: their effects can be combined into the Wilson coefficients,

[TABLE]

with the function $f_{I}$ defined as liu :

[TABLE]

where we have neglected the terms proportional to $r^{4}_{\psi(\eta_{c})}$. Certainly, in the following numerical analysis, the NLO Wilson coefficients will be used in the NLO calculations.

III Numerical results and discussions

We use the following input parameters for the numerical calculations pdg14 ; zhou :

[TABLE]

For the CKM matrix elements, we adopt the Wolfenstein parametrization and the updated values $A=0.814,\lambda=0.22537,\bar{\rho}=0.117\pm 0.021$ and $\bar{\eta}=0.353\pm 0.013$ pdg14 . With the total amplitudes, one can write the decay width as:

[TABLE]

The wave functions of $B,\pi$ and $K$ have been well defined in many works, while those of the two excited charmonia states exist many uncertainties. Here we take the harmonic-oscillator model zhou :

[TABLE]

with

[TABLE]

where $f_{2S}$ refers to the decay constant $f_{\psi}$ or $f_{\eta_{c}}$, the free parameter $\omega=0.2\pm 0.1$ GeV and the $c$ quark mass $m_{c}=1.275\pm 0.025$ GeV. The main errors come from these parameters and $\omega_{b}=0.4\pm 0.1$ for $B$ meson wave funcion in our calculations.

Using the input parameters and the wave functions as specified in this section, we give the LO and NLO predictions for the considered decays in Table I. For these color-suppressed decays, the amplitudes associated with the Wilson coefficient $C_{1}+C_{2}/3$ usually play the dominant roles. It is instructive to check the contributions from these amplitudes at the leading-order and the next-to-leading-order. Here we take the decay $B^{0}\to\psi(2S)K^{0}$ as an example, the value of the factorizable color-suppressed amplitude $F^{V-A}_{\psi K}(a_{2})$ is about $-3.55\times 10^{-2}$, which will be partly canceled by the real part of nonfactorizable one $Re(M^{V-A}_{\psi K}(C_{2}))=1.16\times 10^{-2}$. And the imagine part $Im(M^{V-A}_{\psi K})$ is small and only about $-5.59\times 10^{-3}$. When the vertex corrections are included, the factorizable color-suppressed amplitude $F^{V-A}_{\psi K}(a_{2})$ becomes a complex number. It’s real part reduces to $-1.49\times 10^{-2}$, which will be largely canceled by $Re(M^{V-A}_{\psi K}(C_{2}))=1.09\times 10^{-2}$ (the difference from the leading order value is because of using the NLO Wilson coefficient). While the imagine part induced by the vertex corrections is large and about $-3.42\times 10^{-2}$. So the total amplitude from the tree operators increases after including the NLO contributions. From the numerical results, we find that the contribution (the square of the amplitudes) from the penguin operators is very small compared with that from the tree operators, about $1.2\%$ at the leading order and $0.05\%$ at the next-to-leading order. There is the similar situation in the decay $B^{0}\to J/\Psi K^{0}$ li0 ; Grossman . That is to say the penguin pollution is very small in these decays. From our results, we can find that the branching ratios of the channels $B\to\psi(2S)K,\eta_{c}(2S)K$ are larger than those of the decays $B\to\psi(2S)\pi,\eta_{c}(2S)\pi$ by one even two orders. This is mainly because of the CKM suppression factor $\lambda=0.22$ within the latter. Just like the argument given in Ref.beneke2 that the soft gluon contribution is suppressed by a factor $\Lambda_{QCD}/(m_{b}\alpha_{s})$ rather than $\Lambda_{QCD}/m_{b}$ in this type of decay, so the perturbative and power corrections can be sizeable. Our predictions shown in Table 1 support this argument: the NLO contributions can provide a $(52\sim 55)\%$ enhancement for the decays involving $\eta_{c}(2S)$ meson. And the enhancement for the decays involving $\psi(2S)$ meson is much more large. Furthermore, the charmed meson rescattering effects Colangelo in the decays $B\to\psi(2S)K,\eta_{c}(2S)K$ might not be very important. While the branching ratios for the decays $B^{+}\to\psi(2S)\pi^{+},B^{0}\to\psi(2S)\pi^{0}$ are still smaller than the data even with the vertex corrections and the NLO Wilson coefficients included. Maybe other possible higher order contributions or the contributions from the Glauber gluons glauber in the spectator diagrams play the important roles.

In the following we will discuss the CP-violating asymmetries of $B\to\psi(\eta_{c})\pi,\psi(\eta_{c})K$ decays. For the charged decays $B^{+}\to\psi(\eta_{c})K^{+}$, there is no weak phase in their decay amplitudes, so the direct CP asymmetries of these decays are zero, and it is in agreement with the experimental value $A_{CP}(B^{+}\to\psi(2S)K^{+})=(-2.4\pm 2.3)\%$ within $1\sigma$ error. The direct CP asymmetries of the charged decay $B^{+}\to\psi(2S)\pi^{+}$ predicted by pQCD approach are listed as:

[TABLE]

which are consistent with the experimental data $(2.2\pm 8.5\pm 1.6)\%$ determined by Belle belle and $(4.8\pm 9.0\pm 1.1)\%$ by LHCb lhcb with no evidence of direct CP violation being seen.

As for the neutral decay channels, the effects of $B^{0}-\bar{B}^{0}$ mixing should be considered. The direct and mixing induced CP violating asymmetries are defined as:

[TABLE]

where the CP-violating parameter $\lambda_{CP}$ is

[TABLE]

for $b\to d(b\to s)$ transition and $\eta_{f}$ is the CP-eigenvalue of the final states. From the left panel in Fig.3, one can find that the dependence of the mixing CP-asymmetry $A^{mix}_{CP}(B^{0}\to\psi(2s)\pi^{0})$ on the CKM weak phase $\alpha$. If taking the CKM weak phase $\alpha=(85.4^{+3.9}_{-3.8})^{\circ}$ pdg14 , we find that the value of $A^{mix}_{CP}(B^{0}\to\psi(2s)\pi^{0})$ is $(-31.3^{+26.4}_{-22.3})\%$. The mixing induced CP violating asymmetry for the channel $B^{0}\to\psi(2s)\pi^{0}$ is very sensitive to the angle $\alpha$. As to another CKM weak phase $\beta$, there are much more uncertainties: by including earlier $\sin(2\beta)$ measurements exper and recent results from LHCb lhcb2 and Belle belle2 , the Heavy Flavor Averaging Group (HFAVG) gives two possible solutions $2\beta=(43.8\pm 1.4)^{\circ}$ and $2\beta=(136.2\pm 1.4)^{\circ}$ hfavg , which is very consistent with our predictions shown in the right panel of Fig.3. If the assumption that $A^{dir}_{CP}=0$ is relaxed, then $A^{mix}_{CP}=-\eta_{f}\sqrt{1-A^{dir}_{CP}}\sin{2\beta}$. There are similar results between the decays $B^{0}\to\psi(2S)K^{0}_{S}(\pi^{0})$ and $B^{0}\to\eta_{c}(2S)K^{0}_{S}(\pi^{0})$ about the mixing induced CP violating asymmetries. As for the direct CP asymmetries of the decays $B^{0}\to\eta_{c}(2S)\pi^{0}$ and $B^{+}\to\eta_{c}(2S)\pi^{+}$, they are very small and at $10^{-5}$ order.

IV Summary

We study the B meson decays $B\to\psi(2S)K(\pi),\eta_{c}(2S)K(\pi)$ within the pQCD approach, where the radially excited charmonia states are involved. With the wave functions of these two mesons $\psi(2S)$ and $\eta_{c}(2S)$ derived from the harmonic-oscillator-model, we find that the branching ratios for the decays $B^{+}\to\psi(2S)K^{+},\eta_{c}(2S)K^{+}$ and $B^{0}\to\psi(2S)K^{0}$ can agree well with the data within errors after including the NLO corrections. While there is still some room left for other high order corrections or the non-perturbative long distance contributions for the decays $B^{+}\to\psi(2S)\pi^{+}$ and $B^{0}\to\psi(2S)\pi^{0}$. The pQCD predictions for the direct CP-violating asymmetries support the present experimental opinion, that is no evidence of direct CP violation being observed in these decays. If a few percent value is confirmed in the future, it would indicate new physics definitely.

Acknowledgment

This work is partly supported by the National Natural Science Foundation of China under Grant No. 11347030, by the Program of Science and Technology Innovation Talents in Universities of Henan Province 14HASTIT037.