The S-wave resonance contributions in the $B^0_s$ decays into $ \psi(2S,3S)$ plus pion pair

TL;DR

This paper analyzes the S-wave resonance contributions in $B^0_s$ decays to $\psi(2S,3S)$ plus pion pairs using perturbative QCD, accounting for resonances and nonresonant effects, and compares predictions with experimental data.

Contribution

It introduces a detailed perturbative QCD analysis of $B^0_s$ decays including next-to-leading order corrections and resonance effects, providing predictions consistent with experimental results.

Findings

Predicted $B^0_s ightarrow \psi(2S) \pi^+ \pi^-$ branching ratio matches LHCb data.

S-wave contributions estimated at $6.0 imes 10^{-5}$ for $\\psi(2S)$ decay.

Predictions for $B^0_s ightarrow \\psi(3S) \\pi^+ \\pi^-$ are of order $10^{-5}$.

Abstract

The three-body decays $B^0_s \rightarrow \psi(2S,3S) \pi^+ \pi^-$ are studied based on the perturbative QCD approach. With the help of the nonperturbative two-pion distribution amplitudes, the analysis is simplified into the quasi-two-body processes. Besides the traditional factorizable and nonfactorizable diagrams at the leading order, the next-to-leading order vertex corrections are also included to cancel the scale dependence. The $f_0(980)$, $f_0(1500)$ resonance contributions as well as the nonresonant contributions are taken into account using the presently known $\pi\pi$ time-like scalar form factor for the $s\bar{s}$ component. It is found that the predicted $B^0_s \rightarrow \psi(2S) \pi^+ \pi^-$ decay spectra in the pion pair invariant mass shows a similar behavior as the experiment. The calculated S-wave contributions to the branching ratio of $B^0_s \rightarrow \psi(2S) \pi^+ \pi^-$ is $6.0\times 10^{-5}$, which is in agreement with the LHCb data $\mathcal {B}(B^0_s \rightarrow \psi(2S) \pi^+ \pi^-)=(7.2\pm 1.2)\times 10^{-5} $ within errors. The estimate of $\mathcal {B}(B^0_s \rightarrow \psi(3S) \pi^+ \pi^-)$ can reach the order of $10^{-5}$, pending for the corresponding measurements.

Tables

Table 1. Table 1: The properties of ψ ​ ( 2 ​ S ) 𝜓 2 𝑆 \psi(2S) and ψ ​ ( 3 ​ S ) 𝜓 3 𝑆 \psi(3S) mesons.

Mesons	$I^G$	$J^{PC}$	Mass (MeV)	Width (keV)
$\psi \left(2S\right)$	$0^{-}$	$1^{--}$	$3686.097\pm 0.025$	$296\pm 8$
$\psi \left(3S\right)$	$0^{-}$	$1^{--}$	$4039\pm 1$	$800\pm 100$

Table 2. Table 2: Contributions to the branching ratios ( in unit of 10 − 5 superscript 10 5 10^{-5} ) of B s 0 → ψ ​ ( 2 ​ S , 3 ​ S ) ​ π + ​ π − → subscript superscript 𝐵 0 𝑠 𝜓 2 𝑆 3 𝑆 superscript 𝜋 superscript 𝜋 B^{0}_{s}\rightarrow\psi(2S,3S)\pi^{+}\pi^{-} decays from the factorizable amplitudes F 𝐹 F , from the nonfactorizable amplitudes M 𝑀 M , and from the interference between F 𝐹 F and M 𝑀 M .

Modes	FF	MM	FM	Total
$B_s^0\to \psi \left(2S\right)f_0\left(980\right)\to \psi \left(2S\right){\pi }^{+}{\pi }^{-}$	6.4	0.2	-1.5	5.1
$B_s^0\to \psi \left(2S\right)f_0\left(1500\right)\to \psi \left(2S\right){\pi }^{+}{\pi }^{-}$	0.34	0.04	-0.15	0.23
$B_s^0\to \psi \left(2S\right){\left({\pi }^{+}{\pi }^{-}\right)}_{\text{S-wave}}$	7.3	0.3	-1.6	6.0
$B_s^0\to \psi \left(3S\right)f_0\left(980\right)\to \psi \left(3S\right){\pi }^{+}{\pi }^{-}$	2.0	0.2	-0.5	1.7

Equations

$$\mathcal{A}(B_s^0\to \psi (2S,3S){\pi }^{+}{\pi }^{-})=H\otimes {\varphi }_B\otimes {\varphi }_{\psi }\otimes {\varphi }_{\pi \pi }^S,$$

$$p_B$$

$$k_B$$

$$p_1=\frac{M_B}{\sqrt{2}}((1-r^2)\zeta ,\eta (1-\zeta ),{\text{p}}_{1\mathrm{T}}),p_2=\frac{M_B}{\sqrt{2}}((1-r^2)(1-\zeta ),\eta \zeta ,{\text{p}}_{2\mathrm{T}}).$$

$${\Phi }_{B_s}(x,b)=\frac{i}{\sqrt{2N_c}}[({p/}_B+M_B){\gamma }_5{\varphi }_{B_s}(x,b)],$$

$${\varphi }_{B_s}(x,b)=N{x}^2(1-x{)}^2\exp [-\frac{x^2{M}_B^2}{2{\omega }_b^2}-\frac{{\omega }_b^2{b}^2}{2}],$$

$${\int }_0^1{\varphi }_{B_s}(x,b=0)dx=\frac{f_{B_s}}{2\sqrt{2N_c}}.$$

$${\Phi }_{\psi }^L=\frac{1}{2\sqrt{N_c}}[m{ϵ/}_L{\varphi }^L(x,b)+{ϵ/}_L{p/}_3{\varphi }^t(x,b)],$$

$${\varphi }_{\pi \pi }^S=\frac{1}{2\sqrt{N_c}}[p/{\Phi }_{v\nu =-}^{I=0}(z,\zeta ,{\omega }^2)+\omega {\Phi }_s^{I=0}(z,\zeta ,{\omega }^2)+\omega (n/v/-1){\Phi }_{t\nu =+}^{I=0}(z,\zeta ,{\omega }^2)],$$

$$\frac{d\mathcal{B}}{d\omega }=\frac{\tau \omega |\overrightarrow{p_1}||\overrightarrow{p_3}|}{32{\pi }^3{M}_B^3}|\mathcal{A}{|}^2,$$

$$|\overrightarrow{p_1}|=\frac{1}{2}\sqrt{{\omega }^2-4m_{\pi }^2},|\overrightarrow{p_3}|=\frac{1}{2\omega }\sqrt{[M_B^2-(\omega +m{)}^2][M_B^2-(\omega -m{)}^2]},$$

$$\mathcal{A}=V_{cb}^{\ast }{V}_{cs}(F^{LL}+M^{LL})-V_{tb}^{\ast }{V}_{ts}(F^{\prime LL}+F^{LR}+M^{\prime LL}+M^{SP}),$$

$$a_1$$

$$a_2$$

$$a_3$$

$$M_B$$

$$\lambda$$

$$f_{B_s}=259\pm 32,f_{\psi (2S)}=296_{-2}^{+3},f_{\psi (3S)}=187\pm 8.$$

$$\mathcal{B}(B_s^0\to \psi (2S)f_0(980)\to \psi (2S){\pi }^{+}{\pi }^{-})$$

$$\mathcal{B}(B_s^0\to \psi (2S)f_0(1500)\to \psi (2S){\pi }^{+}{\pi }^{-})$$

$$\mathcal{B}(B_s^0\to \psi (2S)({\pi }^{+}{\pi }^{-}{)}_{S-wave})=[6.0_{-1.4}^{+1.7}({\omega }_b{)}_{-1.4}^{+1.6}(f_{B_s}{)}_{-0.9}^{+0.9}(a_2^{I=0}{)}_{-0.5}^{+0.6}(t)]\times 10^{-5},$$

$$\frac{\mathcal{B}(B_s^0\to \psi (2S){\pi }^{+}{\pi }^{-})}{\mathcal{B}(B_s^0\to J/\psi {\pi }^{+}{\pi }^{-})}=0.34\pm 0.04(\text{stat})\pm 0.03(\text{syst})\pm 0.01(\mathcal{B}),$$

$$\mathcal{B}(B_s^0\to \psi (3S)f_0(980)\to \psi (3S){\pi }^{+}{\pi }^{-})=[1.7_{-0.3}^{+0.5}({\omega }_b{)}_{-0.4}^{+0.5}(f_{B_s}{)}_{-0.1}^{+0.3}(a_2^{I=0}{)}_{-0.0}^{+0.2}(t)]\times 10^{-5}$$

Full text

I Introduction

The decays of the neutral $B$ mesons to the charmonium state plus light meson pair have attracted lots of attentions recently. These decay processes could be used to study the spectroscopy in various meson pair systems, and they also have significant roles in understanding of the substructure of different resonant states. Up to now, several experimental collaborations, like LHCb and BaBar, have measured the decays $B_{s}^{0}\rightarrow J/\psi\pi^{+}\pi^{-}$ prd86052006 ; prd89092006 , $B_{s}^{0}\rightarrow J/\psi K^{+}K^{-}$ prd87072004 ; prl114041801 , $B^{0}\rightarrow J/\psi\pi^{+}\pi^{-}$prl90091801 ; prd87052001 ; prd90012003 ; plb74238 , $B^{0}\rightarrow J/\psi K^{+}K^{-}$ prd88072005 , $B^{0}_{s}\rightarrow J/\psi\phi\phi$ jhep03040 , $B^{0}_{s}\rightarrow\psi(2S)K^{+}\pi^{-}$ plb747484 .

On the theoretical side, the studies of $B$ mesons decay into three hadronic final states have been performed using different approaches plb56490 ; prd72094031 ; plb727136 ; plb726337 ; prd89074043 ; plb728579 ; prd91014029 ; HYCPRD1 ; HYCPRD2 ; HYCPRD3 ; HYCPRD4 ; prd87076007 ; epjc75536 ; prd89094007 ; pma58031001 ; npb899247 ; BEI1 ; BEI2 ; prd81094033 ; prd89094013 ; prd92054010 ; jhep02009 ; epjc75609 ; ps88035101 ; jhep09074 in recent years. The perturbative QCD approach (PQCD) pqcd1 ; pqcd2 has some unique features that are particularly suitable for dealing with the three-body $B$ meson decays. First, by introducing two-meson distribution amplitudes  fp42101 ; prl811782 ; prd62073014 ; epjc26261 ; npb555231 ; sjnp38289 ; tmp691109 , the analysis of three-body hadronic $B$ meson decays is simplified into the one for two-body decays. As pointed out in  plb561258 ; prd70054006 , it is not practical to make a direct calculation for the three-body processes due to the enormous number of diagrams which contain two gluon exchange at lowest order. Besides, its contribution is not important because the region with the two gluons being hard simultaneously is power suppressed. The dominant contributions come from the kinematic region where the two light mesons move almost parallelly for producing a resonance with an invariant mass below $O(\bar{\Lambda}M_{B})$ ($\bar{\Lambda}$ being the heavy-meson and heavy-quark mass difference), which can be catched by this new nonperturbative inputs. Second, the three-body decays of $B$ meson receive both resonant and nonresonant contributions, it is difficult to separate the two parts clearly  npb899247 . After absorbing the nonperturbative dynamics associated with the meson pair into the complex time-like form factors in the two-meson distribution amplitudes, both resonant and nonresonant contributions plb561258 ; prd70054006 are included in the PQCD approach. Third, in the above dominant contribution region, the end-point singularities are smeared by the two-meson invariant mass, which suggest that the PQCD approach has a good predictive power without any arbitrary cutoffs. Finally, different from the two-body cases, some possible final state interactions (FSIs) may be significant in the three-body decays prd89094013 ; fsi2 . According to fsi2 , there are two distinct FSIs mechanisms. One is the interactions between the meson pair in the resonant region associated with various intermediate states. The other is the rescattering between the third particle and the pair of mesons. In our opinion, the former can be factorized into two-meson distribution amplitudes while the latter is ignored in the quasi-two-body approximation. Over past few years, the braching ratios and direct $CP$ asymmetries of the three-body decays such as $B^{\pm}\rightarrow\pi^{\pm}(K^{\pm})\pi^{+}\pi^{-}$ prd89074031 , $B^{0}_{(s)}\rightarrow J/\psi\pi^{+}\pi^{-}$ prd91094024 , $B^{0}_{(s)}\rightarrow\eta_{c}\pi^{+}\pi^{-}$ 150906117 , $B^{+}_{c}\rightarrow D^{+}_{(s)}\pi^{+}\pi^{-}$ prd94034040 , $B_{(s)}\rightarrow(D_{(s)},\bar{D}_{(s)})\pi^{+}\pi^{-}$ 161108786 , $B_{(s)}\rightarrow P\pi^{+}\pi^{-}$ 161205934 and $B^{0}_{(s)}\rightarrow\eta_{c}(2S)\pi^{+}\pi^{-}$ 170101844 have been studied systematically in the PQCD approach. In a recent work, the P-wave resonance contributions have been calculated in the $B\rightarrow K\pi\pi$ decays 160904614 .

The three-body hadronic $B$ meson decays with the radially excited charmonium mesons in the final state have not received much attentions in the literature. As the LHC gathering more and more data, the processes of the $B$ mesons decays including excited charmonium states must have much possibilities to be found. Recently, the LHCb Collaboration have measured the three-body decay channels of $B^{0}_{(s)}\rightarrow\psi(2S)\pi^{+}\pi^{-}$ with $pp$ collision data collected at $\sqrt{s}=7$ TeV recently npb871403 . In the quark model, $\psi(2S)$ ($\psi(3S)$) is the first (second) radially excited vector charmonium with the radial quantum number $n=2(3)$ and the orbital angular momentum $l=0$. Both $\psi(2S)$ and $\psi(3S)$ were observed by the processes of the $e^{+}e^{-}$ annihilation into hadronic prl331453 ; dasp . The properties for the two high excited charmonium were updated in PDG 2016 pdg ; they are listed in Table 1.

We have previously studied the semileptonic and two-body nonleptonic decays of the $B_{c}$ meson to radially excited charmonium mesons in the PQCD approach by using the harmonicosillator wave functions for the charmonium states 160208918 ; epjc75293 . In the present work, we will extend our analysis to the three-body decays $B^{0}_{s}\rightarrow\psi(2S,3S)\pi^{+}\pi^{-}$ to provide an order of magnitude estimation. By introducing the two-pion distribution amplitudes, the S-wave contributions, which are the main contributions in the three-body decays $B^{0}_{s}\rightarrow\psi(2S,3S)\pi^{+}\pi^{-}$ npb871403 , could be described by the quasi-two-body processes $B^{0}_{s}\rightarrow\psi(2S,3S)f_{0}\rightarrow\psi(2S,3S)\pi^{+}\pi^{-}$ containing the S-wave resonant states $f_{0}$, $f_{0}(980)$ and $f_{0}(1500)$ as two examples. Following the steps in Ref. prd91094024 , the decay amplitude $\mathcal{A}(B^{0}_{s}\rightarrow\psi(2S,3S)\pi^{+}\pi^{-})$ can be written as the convolution

[TABLE]

where the hard function $H$ could be treated by PQCD including both factorizable and nonfactorizable contributions in the expansion at the leading order in $\alpha_{s}$ (single gluon exchange as depicted in Fig.1). The hadron wave functions $\phi_{B},\phi_{\psi}$ and $\phi^{S}_{\pi\pi}$ absorb the nonperturbative dynamics, and they can be extracted from experimental data or other nonperturbative methods.

Following the introduction, in Sect. II, we shall define kinematics and describe the wave functions of the initial and final states, then we will briefly review the related theoretical formulas. In Sect. III, we will calculate the $B^{0}_{s}\rightarrow\psi(2S,3S)\pi^{+}\pi^{-}$ decays in the PQCD approach with discussions. Finally we will close this paper with a conclusion.

II framework

It is convenient to work in the rest frame of the $B_{s}^{0}$ meson. Its momentum $p_{B}$, along with the charmonium meson momentum $p_{3}$, the pion pair momentum $p$ and other quark momenta $k_{i}$ in each meson, is chosen as prd91094024

[TABLE]

with the ratio $r=m/M_{B}$ and $m(M_{B})$ is the mass of the charmonium ($B^{0}_{s}$) meson. The factor $\eta=\omega^{2}/(M_{B}^{2}-m^{2})$ with the invariant mass squared $\omega^{2}=p^{2}$ for the pion pair. The $k_{iT}$, $x_{i}$ represent the transverse momentum and longitudinal momentum fraction of the quark inside the meson, respectively. If we choose the $\zeta=p_{1}^{+}/p^{+}$ as the $\pi^{+}$ meson momentum fraction, the two pion momenta $p_{1,2}$ can be written as

[TABLE]

Similar to the situation of the $B$ meson npb591313 , $B^{0}_{s}$ meson could also be treated as a heavy-light system, its wave function in the $b$ space can be expressed by prd63054008 ; prd65014007 ; epjc28515 ; ppnp5185

[TABLE]

where $b$ is the conjugate space coordinate of the transverse momentum $k_{\rm BT}$, and $N_{c}$ is the color factor. Here, we only consider one of the dominant Lorentz structures in our calculation. The distribution amplitude $\phi_{B_{s}}$ is adopt the same form as it in Refs. prd76074018 ; prd86-114025

[TABLE]

with shape parameter $\omega_{b}=0.50\pm 0.05$ GeV and the normalization constant $N$ being related to the decay constant $f_{B_{s}}$ through

[TABLE]

For the considered decays, the vector charmonium meson is longitudinally polarized. The longitudinal polarized component of the wave function is defined as 160208918 ; epjc75293

[TABLE]

with the longitudinal polarization vector $\epsilon_{L}=\frac{M}{\sqrt{2}m}(-r^{2},1-\eta,\textbf{0}_{T})$. For the twist-2 (twist-3) distribution amplitudes $\phi^{L}(\phi^{t})$ of 2S and 3S states, the same form and parameters are adopted as their in Refs. 160208918 ; epjc75293 .

According to prd91094024 , the S-wave two-pion distribution amplitudes are organized into

[TABLE]

where $n=(1,0,\textbf{0}_{T})$ and $v=(0,1,\textbf{0}_{T})$ are two dimensionless vectors. The asymptotic models for the twist-2 distribution amplitude $\Phi_{v\nu=-}$ and the twist-3 distribution amplitude $\Phi_{s},\Phi_{t\nu=+}$ and relevant time-like scalar form factor can be found in prd91094024 .

Now, we write down the differential branching ratio for $B^{0}_{s}\rightarrow\psi(2S,3S)\pi^{+}\pi^{-}$ decays,

[TABLE]

where $\tau=1.512\times 10^{-12}s$ is the lifetime of $B^{0}_{s}$ meson. The three-momenta of the pion in the pion pair center-of-mass system $\overrightarrow{p_{1}}$ and that of the charmonium $\overrightarrow{p_{3}}$ are given by

[TABLE]

with $m_{\pi}$ the pion mass.

The decay amplitude $\mathcal{A}$ is written as

[TABLE]

with $V_{ij}$ the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements. The detailed expressions of $F$ (the factorizable emission contributions) and $M$ (the nonfacorizable contributions) are the same as the $B^{0}_{s}\rightarrow J/\psi\pi^{+}\pi^{-}$ process in the appendix of Ref. prd91094024 , except for the replacement $J/\psi\rightarrow\psi(2S,3S)$. In this work, we also consider the vertex corrections to the factorizable amplitudes $F$ at the current known next-to-leading order (NLO) level. Their effects can be combined in the Wilson coefficients as usual npb591313 ; beneke1 ; beneke2 . In the NDR scheme we have

[TABLE]

The functions $f_{I}$ and $g_{I}$ arise from the vertex corrections; they can be found in prd63074011 .

III results

The mesons and quark masses (in units of GeV) and the Wolfenstein parameters are taken from the Particle Data Group pdg

[TABLE]

The decay constants needed (in units of MeV) are prl95212001 ; 160208918

[TABLE]

By using Eq. (9), we have the predictions of the branching ratios as

[TABLE]

where the errors are induced by the shape parameter $\omega_{b}=0.50\pm 0.05$ GeV, the decay constant $f_{B_{s}}=0.259\pm 0.032$, the Gegenbauer coefficient $a_{2}^{I=0}=0.2\pm 0.2$ for the two-pion system, and the hard scale $t$ which varies from $0.75t$ to $1.25t$, respectively. The errors from the uncertainty of the CKM matrix elements and the decay constants of charmonia are very small, and have been neglected. It is found that the main uncertainties in our approach come from the $B_{s}$ meson wave function (the first two errors), which can reach $30-50\%$ in size. The uncertainties caused by the Gegenbauer coefficient are less than $20\%$ which is similar to that of $B^{0}_{s}\rightarrow J/\psi\pi^{+}\pi^{-}$ prd91094024 . The scale-dependent uncertainty is largely reduced due to the inclusion of the next-to-leading order vertex corrections. Summing the two resonances $f_{0}(980)$ and $f_{0}(1500)$ in the strange scalar form factors, we have the total S-wave contribution for $B^{0}_{s}\rightarrow\psi(2S)\pi^{+}\pi^{-}$ decay as

[TABLE]

which is compatible with the LHCb data $(7.2\pm 1.2)\times 10^{-5}$ pdg when considering its uncertainty. The measured ratio between the branching ratios of $B^{0}_{s}\rightarrow\psi(2S)\pi^{+}\pi^{-}$ and $B^{0}_{s}\rightarrow J/\psi\pi^{+}\pi^{-}$ is npb871403

[TABLE]

which is consistent with the PQCD prediction $0.37^{+0.25}_{-0.18}$ (uncertainties added in quadrature), where the $\mathcal{B}(B^{0}_{s}\rightarrow J/\psi\pi^{+}\pi^{-})$ is read from prd91094024 . Because the mass of the $f_{0}(1500)$ resonance is beyond the $\pi^{+}\pi^{-}$ invariant mass spectra ($2m_{\pi}<\omega<M-m$) in $B^{0}_{s}\rightarrow\psi(3S)\pi^{+}\pi^{-}$ decay, there is only one resonant state $f_{0}(980)$ contributing to this decay. We have the branching ratio:

[TABLE]

In Fig. 2, we plot the differential branching ratio for $B^{0}_{s}\rightarrow\psi(2S,3S)\pi^{+}\pi^{-}$ as a function of the $\pi^{+}\pi^{-}$ invariant mass $\omega$. The red dashed and blue dotted curves represent the contributions from the resonances $f_{0}(980)$ and $f_{0}(1500)$, respectively. As expected, the $f_{0}(980)$ production is clearly dominant. From the branching ratios in Eqs. (III) and (16) one could find that the $f_{0}(980)$ resonance accounts for $85.0\%$ of the total branching ratio, the $f_{0}(1500)$ resonance $3.8\%$, and positive interference between the two terms is for $11.2\%$.

For a more direct comparison with the available experimental data npb871403 , we also present the S-wave $\pi^{+}\pi^{-}$ invariant mass distributions for the $B^{0}_{s}\rightarrow\psi(2S)\pi^{+}\pi^{-}$ decay in this work as well as the Fig. 4b in Ref. npb871403 . One observes an appreciable peak arising form the $f_{0}(980)$ resonance and a less strong, but clearly visible peak for the $f_{0}(1500)$ in Fig. 3 a. Comparing with the data, our distribution below $1.4$ GeV for the resonances agrees quite well, showing a similar behavior in this region.

Different from the fixed kinematics of the two-body decays, the decay amplitudes of the quasi-two-body decays are dependent on the $\pi^{+}\pi^{-}$ invariant mass. Therefore, it is more convenient to compare these contributions to the branching ratios, whose results are displayed separately in Table 2, where the labels $FF$, $MM$, $FM$ correspond to the contribution of the factorizable, nonfactorizable cases, and the interferences between them, respectively, while the label “Total” denotes the total contribution. It can be found that the dominant contributions to the branching ratios coming from the factorizable topology due to the vertex corrections, which are enhanced by the Wilson coefficient $C_{2}$ (see Eq. (II)). The interference contributions are of the same order as the factorizable ones with an opposite sign, which reflects the importance of nonfactorizable effects in the color-suppressed processes. This is similar to the case of the two-body $B$ meson decays into charmonia prd561615 ; prd71114008 ; prd68034004 .

IV conclusion

In this work, the quasi-two-body decays $B^{0}_{s}\rightarrow\psi(2S,3S)\pi^{+}\pi^{-}$ have been analyzed in the PQCD approach, assuming dominance of the S-wave resonance states $f_{0}(980)$ and $f_{0}(1500)$ in the invariant $\pi^{+}\pi^{-}$ mass distributions. Both the factorizable (including the vertex corrections) and nonfactorizable contributions are taken into account. We discussed theoretical uncertainties arising from the nonperturbative shape parameters, the decay constant, the Gegenbauer coefficient, and the scale dependence. It is found that the main uncertainties of the concerned processes come from the shape parameters and the decay constant of the $B_{s}$ meson. The predicted branching ratio and the invariant mass distributions for $B^{0}_{s}\rightarrow\psi(2S)\pi^{+}\pi^{-}$ decay are in agreement with the results from LHCb Collaboration. The decay mode $B^{0}_{s}\rightarrow\psi(3S)\pi^{+}\pi^{-}$ has not been measured yet, while the large value of the prediction $\mathcal{B}=1.7^{+0.8}_{-0.5}\times 10^{-5}$ for it in this work waits the the future measurements.

Acknowledgements.

The authors are grateful to Hsiang-nan Li for helpful discussions. This work was supported in part by the National Natural Science Foundation of China under Grants Nos. 11547020, 11605060, and 11547038, in part by the Natural Science Foundation of Hebei Province under Grant No. A2014209308, in part by the Program for the Top Young Innovative Talents of Higher Learning Institutions of Hebei Educational Committee under Grant No. BJ2016041, and in part by the Training Foundation of North China University of Science and Technology under Grant Nos. GP201520 and JP201512.