Study of the quasi-two-body decays Bs0→ψ(3770)(ψ(3686))π+π-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B^{0}_{s} \rightarrow \psi (3770) (\psi (3686))\pi ^+\pi ^-$$\end{document} with perturbative QCD approach

In this note, we study the contributions from the S-wave resonances, f0(980)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_{0}(980)$$\end{document} and f0(1500)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_{0}(1500)$$\end{document}, to the Bs0→ψ(3770)π+π-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B^{0}_{s}\rightarrow \psi (3770)\pi ^ {+}\pi ^{-}$$\end{document} decay by introducing the S-wave ππ\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} distribution amplitudes within the framework of perturbative QCD approach. Both the resonant and the non-resonant contributions are contained in the scalar form factor in the S-wave distribution amplitude ΦππS\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Phi ^S_{\pi \pi }$$\end{document}. Since the vector charmonium meson ψ(3770)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\psi (3770)$$\end{document} is a S–D wave mixed state, we calculate the branching ratios of S-wave and D-wave, respectively, and the results indicate that f0(980)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_{0}(980)$$\end{document} is the main contribution of the considered decay, and the branching ratio of the ψ(2S)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\psi (2S)$$\end{document} mode is in good agreement with the experimental data. We also take the S–D mixed effect into the Bs0→ψ(3686)π+π-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B^{0}_{s}\rightarrow \psi (3686)\pi ^ {+}\pi ^{-}$$\end{document} decay. Our calculations show that the branching ratio of Bs0→ψ(3770)(ψ(3686))π+π-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B^{0}_{s}\rightarrow \psi (3770)(\psi (3686))\pi ^ {+}\pi ^{-}$$\end{document} can be at the order of 10-5\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$10^{-5}$$\end{document}, which can be tested by the running LHCb experiments.


Introduction
In the present decade, the B mesons' three-body hadronic decays have drawn a lot of attention on both the experimental and the theoretical sides, since they can test the standard model (SM) and help us to have a better understanding of the scheme of QCD dynamics. The three-body decays of B mesons are more complicated than the two-body cases because they include both resonant and non-resonant contributions. We have two-gluon exchange in the decay amplitude within our framework of theoretical calculation, which will lead to a hard kernel that involves three bodies, and we may introduce possible final-state interactions [1][2][3][4]. So just as stated in Ref. [5], because of the interference between the resonance and non-resonance, it is difficult to make a direct a e-mail: Liangzr@email.swu.edu.cn b e-mail: dfbdfbok@163.com c e-mail: yuxq@swu.edu.cn calculation of the resonance and non-resonance contributions separately.
In the pQCD approach which is based on the k T factorization theorem, the three-body decay can be simplified into a two-body case by bringing in two-hadron distribution amplitudes [43,44], which contain the information of both resonance and non-resonance. The dominant contributions come from the region where the pair of two light mesons move parallel with an invariant mass below O(¯ M B ), wherē = M B − m b is the mass difference of the B meson and b-quark. Therefore, one can express the typical pQCD factorization formula of the B meson's three-body decay amplitude as [35,36] where the hard decay kernel H describes the contribution from the region with only a one gluon exchange diagram at leading order, which can be calculated in perturbation theory. The terms φ B , φ h 1 h 2 and φ h 3 , which can be extracted from experiment or calculated by several nonperturbative approaches, and can be regarded as nonperturbative input, are the distribution amplitude of the B meson, the h 1 h 2 pair and h 3 , respectively. The B 0 s → ψ(2S)π + π − decay was first observed by the LHCb Collaboration [7]; the data was displayed based on the integrated luminosity of 1 f b −1 in the pp collisions at a center-of-mass energy of √ s = 7 TeV. It is found that f 0 (980) is the main source of the decay rate by a method called the sPlot technique. The B 0 s → ψ(3770)π + π − decay has not been observed yet, so it is desirable to make a theoretical prediction for the branching ratios of this decay mode, for testing the three-body decay's mechanism and the mixing scheme of ψ(3770) as well. In this work, we will calculate the branching ratio of the quasi-two-body decay mode B 0 s → ψ(3770)π + π − . Since the vector charmonium meson ψ(3770), the lowest-lying charmonium state just above the DD threshold, is mainly regarded as a S-D mixture, we will take the S-wave and D-wave contribution into account, respectively. Here, ψ(2S) is the first radially excited charmonium meson, and the pure 1D state indicates the principal quantum number n = 1 and the orbital quantum number l = 2. The S-D mixing angle θ can be obtained from the ratio of the leptonic decay widths of ψ(3686) and ψ(3770) [45]. In Ref. [46], the authors make tentative calculations for different mixing solutions of the B meson exclusive decay B → ψ(3770)K with the QCD factorization, and they drew the conclusion that when taking account of higher-twist effects and adopting the S-D mixing angle θ = −(12 ± 2) • (which the widely accepted value), the branching ratio of the decay B → ψ(3770)K can fit the experimental data well. Also, in Refs. [47][48][49], the authors provided two sets of mixing schemes within the nonrelativistic potential model: θ = −(12 ± 2) • or θ = (27 ± 2) • . Here, the charmonium mesons ψ(3686) and ψ(3770) may in a good approximation be described as [47][48][49][50] ψ(3686) = cos θ |cc(1D) + sin θ |cc(2S) , The contents of this paper is as follows. After the introduction, we describe the theoretical framework and the wave function of the excited charmonium mesons ψ(2S) and ψ(1D) in Sect. 2. In Sect. 3, we list the decay amplitude of the considered decay modes. The numerical results and analysis as regards the results we have got are presented in Sect. 4. Finally we will finish this paper with a brief summary.

The theoretical framework and the wave function
For the quasi-two-body B 0 s → ψ f 0 (→ π + π − ) decays, the relevant weak effective Hamiltonian can be written as [51] (a) (b) (c) (d) Fig. 1 The lowest order Feynman diagrams for the B 0 where V * cb V cs and V * tb V ts are Cabibbo-Kobayashi-Maskawa (CKM) factors, O i (μ) is a local four-quark operator and C i (μ) is the corresponding Wilson coefficient.
It is convenient for us to choose the light-cone coordinates for simplicity. In these coordinates, we choose the B 0 s meson at rest, and let the ππ meson pair and ψ(2S, 1D) meson move along with the direction of n = (1, 0, 0 ) and v = (0, 1, 0 ), respectively, and the Feynman diagrams are described in Fig. 1. So the momenta of the B 0 s ( p B ), ππ ( p), and ψ(2S, 1D) ( p 3 ) are written as Meanwhile, the corresponding light quark's momentum in each meson reads as follows: where M B 0 s is the mass of B 0 s , and r = the pion-pair invariant mass ω 2 and its momentum p satisfy the relation ω 2 = p 2 and p = p 1 + p 2 . The quantities x 1 , z, and x 3 indicate the momentum fractions of the spectator quark inside the meson; they are in the range of 0 ∼ 1. By introducing the kinematic variables, ζ , of the pion pair, we define ζ = p + 1 / p + as the π + meson momentum fraction; the other component's kinematic variables of the pion pair can be expressed as In our calculations, the hadron B 0 s is usually treated as a heavy-light system, and the wave function of which can be found in Refs. [52][53][54]. We have where the distribution amplitude the normalization factor N B can be calculated by the normalization relation being the number of colors. Here, we choose the shape parameter ω B s = 0.50 ± 0.05 GeV [55].
The vector charmonium meson ψ(3770), as mentioned above, is commonly regarded as a S-wave and D-wave mixing state. We adopt the wave function form of this vector charmonium meson on the basis of a harmonic-oscillator potential, which has been applied to the charmonium state successfully, such as J/ψ, ψ(2S), ψ(3S) and so on [50,[56][57][58]. The theoretical results agree well with the measured experimental data, which indicates the reasonableness of adopting this form of the function. For the wave function of the pure 2S state, ψ(2S), and the pure 1D state, ψ(1D), the longitudinal polarized component is defined as [57,58] where p 3 is the momentum of the charmonium mesons ψ(2S) and ψ(1D), with the longitudinal polarization vec- Here the ψ L and ψ t correspond to twist-2 and twist-3 distribution amplitudes (DAs). The explicit forms are [50,57] with I( for ψ(2S) and for ψ(1D). For the shape parameter ω 1D in the DAs of the ψ(1D), we choose ω 1D = 0.5±0.05 GeV, for the reason we have discussed in Ref. [50], and ω 2S = 0.2 ± 0.1 GeV [57].
is the normalization constant, which satisfies the normalization conditions: and the decay constants of the radially excited state ψ(2S) and the angular excitation state ψ(1D) are given in Table  1. Both the wave functions of Eq. (10) and of Eq. (11) are symmetric under x ↔ x.
The differential branching ratios for the B 0 s → ψ(2S, 1D) π + π − decay in the B 0 s meson rest frame can be written as [63] with in the pion-pair center-of-mass system and with the B 0 s meson lifetime τ B 0 s .

The decay amplitudes
In the pQCD factorization approach, the B 0 s → ψ(2S)π + π − decay amplitude A can be expressed in form of where the explicit forms of and  (16)  In our work, we also take vertex corrections into account in the factorization diagrams, and the Wilson coefficients are combined in the NDR scheme [64][65][66] as follows: a 4 = C 4 + C 10 , The hard scattering functions f I and g I are given in Ref. [67], the renormalization scale μ is chosen at the order of m b . For the B 0 s → ψ(1D)π + π − decay, the amplitude is similar to the decay amplitude of B 0 s → ψ(2S)π + π − , just replacing the DAs of ψ(2S) with the corresponding DAs of ψ(1D) in Eq. (15).
As for the decay amplitude of the B 0 s → ψ(3770) (ψ(3686))π + π − decay, we give the expression based on the idea of the S-D mixing scheme:

Numerical results and discussions
In our numerical calculation, the input parameters are listed in Table 1, where the masses of the involved mesons, the lifetime of meson and the Wolfenstein parameters are obtained from the 2018 PDG [63]. The decay constant of ψ(2S) is calculated by the leptonic decay process ψ(2S) → e + e − [58] and the decay constant of ψ(1D) was calculated in Ref.
[68]. The masses of the b quark and c quark are running masses which are calculated under the modified minimal subtraction scheme at the renormalization scale μ, which is equal to the quark mass. By using the differential branching ratio formula Eq. (14), first we make predictions of the branching ratios of the decay mode B 0 s → ψ(2S)π + π − for different intermediate states, including the f 0 (980) and f 0 (1500) resonances, and the numerical results are listed as follows: where the three main errors come from the shape parameter ω B s of the wave function of B 0 s meson, the hard scale t, which varies from 0.9t ∼ 1.1t (not changing 1/b i , i = 1, 2, 3), and the Gegenbauer moment a 2 = 0.2 ± 0.2 [11] in the ππ distribution amplitude, respectively. The other errors from the uncertainty of the input parameters, for example, the decay constants of the B 0 s and charmonium mesons and the Wolfenstein parameters, are tiny and can be neglected safely. We see that the input parameter ω B s of the B 0 s meson is the primary source of the uncertainties, which range approximately in 9.7-17.1%, and then the Gegenbauer moment and the hard scale t, which characterizes the size of the next-leading-order contribution. When we consider the total S-wave contributions of the f 0 (980) and f 0 (1500), we can get which is in agreement with the new experiment data (7.1 ± 1.3) × 10 −5 within the allowed errors [63]. Comparing with previous work [39], we find that our calculation of the branching ratio of the B 0 s → ψ(2S)π + π − is closer to the latest experimental results, for which the main reason is that we adopt the new input parameters in 2018 PDG, and the latest parameters lead to our uncertainties being smaller.
From the numerical results, we can see that f 0 (980) is the principal contribution, which reaches a percentage of 94.7%, just as the experiment observed, and the f 0 (1500) is 1.2%, while the constructive interference between these two resonances can contribute nearly 4.1% to the total branching ratio.
In experiment, the calculated ratio of the branching fraction has been given in Ref. [7], which is By using the previous prediction as regards the branching ratio of the decay mode B 0 s → J/ψ(π + π − ) S [34], we obtain the ratio B(B 0 s → ψ(2S)π + π − )/B(B 0 s → J/ψπ + π − ) = 0.46 +0. 16 −0.18 , which is consistent with the experiment's measurement, and which indicates that the harmonic-oscillator wave function for excited charmonium is applicable and reasonable. Besides the decay mode B 0 s → ψ(2S)π + π − , we make a calculation for the part of 1D, and we consider similar contributions from the containing S-wave resonance state, f 0 (980) and f 0 (1500). The reason is that these two resonances masses are also within the scope of the ππ invariant mass spectra, which is 2m π < ω < M B 0 s − M ψ . After taking the integral over ω, the results are Also, the total S-wave contribution of the B 0 s → ψ (1D)π + π − decay is In Figs. 2a and 3, we plot the differential branching ratio of the B 0 s → ψ(2S, 1D)π + π − decay as a function of the ππ invariant mass ω, in which we can clearly see that the peak arises from f 0 (980), while f 0 (1500) is unsharp, which also makes a contribution to the decay. For comparison, at the same time, we present the experiment data from LHCb [7] in Fig. 2b, which shows a basic agreement with our predictive results. Comparing the results of ψ(2S) and ψ(1D), it is easy to see that the results of ψ(1D) are more sensitive to the Gegenbauer moment, a 2 = 0.2 ± 0.2. This means that, although the value is in good agreement with many decay modes, there is still a necessity to explore more accurate data to facilitate a better understanding of the nonperturbative hadron dynamics. In the ψ(2S) and ψ(1D) modes, since the f 0 (1500) mass is near the maximum of the ππ invariant mass, the corresponding contributions are very small compared to the total contributions of the S-wave. We note that the branching ratio of the ψ(1D) is smaller than that of the ψ(2S), which should be attributed to the dependence of the corresponding wave function and the decay constant. Fig. 3 The S-wave differential branching ratio of the B 0 s → ψ(1D)π + π − decay Furthermore, we calculate the branching fraction of the mode B 0 s → ψ(3770)(ψ(3686))π + π − based on the S-D mixing scheme, whose two sets of mixing angles have been introduced in Sect. 1, and we list the computational results in Table 2. Table 2 Branching ratios of the quasi-two-decay B 0 s → ψ(3686, 3770) f 0 (→ π + π − ) in the pQCD approach based on two sets of S-D mixing angles, where the uncertainties are similar to the previous ones except the last one, which shows dependence of the mixing angle Comparing with the pure D-wave state, we can notice that the branching ratio of the S-D mixing state for B 0 s → ψ(3770)π + π − will be increased approximately by a factor 2 when the mixing angle is −12 • ; the reason is mainly the small decay constant of ψ(1D), which is compatible with the summaries in Refs. [46,49,[68][69][70]. Moreover, we can observe that the results of the B 0 s → ψ(3686)π + π − change a little comparing with the pure 2S mode when taking the mixing effect into account, so ψ(3686) may be regarded as ψ(2S) state. Considering the size of the data collected in LHCb, we may expect the measurement of this decay mode to come in the near future; this will help us to understand the structure of ψ(3770) and the three-body decay mechanism.

Summary
In this work, we have calculated the contributions from the S-wave resonances, f 0 (980) and f 0 (1500), to the B 0 s → ψ(3770)(ψ(3686))π + π − decay by introducing the S-wave ππ distribution amplitudes within the framework of the perturbative QCD approach. Due to the 2S-1D mixing scheme character of ψ(3770), we calculate the branching ratios of Swave and D-wave, respectively, and the results indicate that f 0 (980) is the main contribution of the considered decay, and the differential result of the ψ(2S) mode is in good agreement with the experimental data. We also analyzed the theoretical uncertainties in this paper, to find that the result of ψ(1D) is sensitive to the Gegenbauer coefficient, because of which we need more accurate data to understand the nonperturbative hadron dynamics. Finally, by introducing the mixing angle θ = −12 • and θ = 27 • , we make a further calculation of B 0 s → ψ(3770)(ψ(3686))π + π − , and our calculations show that the branching ratio may be of the order of 10 −5 based on the small mixing angle θ = −12 • , which will be tested by the running LHCb experiments.
where J 0 is the Bessel function and K 0 , I 0 are modified Bessel function with H (1) 0 (x) = J 0 (x) + iY 0 (x). The threshold resummation factor S t (x) has been parameterized in [72],