Branching ratios, $CP$ asymmetries and polarizations of $B\rightarrow \psi(2S) V$ decays

We analyzed the nonleptonic decays $B/B_s\to \psi(2S) V $ with $V=(\rho, \omega, K^{*}, \phi)$ by employing the perturbative QCD (PQCD) factorization approach. Here the branching ratios, the $CP$ asymmetries and the complete set of polarization observables are investigated systematically. Besides the traditional contributions from the factorizable and nonfactorizable diagrams at the leading order, the next-to-leading order (NLO) vertex corrections could also provide considerable contributions. The PQCD predictions for the branching ratios of the $B_{(s)}\to \psi(2S)K^{*}, \psi(2S) \phi$ decays are consistent with the measured values within errors. As for $B\to \psi(2S) \rho, \psi(2S) \omega$ decays, the branching ratios can reach the order of $10^{-5}$ and could be measured in the LHCb and Belle-II experiments. The numerical results show that the direct $CP$ asymmetries of the considered decays are very small. Thus the observation of any large direct $CP$ asymmetry for these decays will be a signal for new physics. The mixing induced $CP$ asymmetries in the neutral modes are very close to $\sin 2\beta_{(s)}$, which suggests that these channels can give a cross-check on the measurement of the Cabbibo-Kobayashi-Maskawa (CKM) angle $\beta$ and $\beta_s$. We found that the longitudinal polarization fractions $f_0$ are suppressed to $\sim 50\%$ due to the large nonfactorizable contributions. The magnitudes and phases of the two transverse amplitudes $\mathcal {A}_{\parallel}$ and $\mathcal {A}_{\perp}$ are roughly equal, which is an indication for the approximate light quark helicity conservation in these decays. The overall polarization observables of $B\to \psi(2S) K^{*0}$ and $B_s\to \psi(2S) \phi$ channels are also in good agreement with the experimental measurements as reported by LHCb and BaBar. Other results can also be tested by the LHCb and Belle-II experiments.


I. INTRODUCTION
Studies of the decays of B mesons into a charmonium meson plus a light vector meson contribute a lot to our knowledge about the CP violation and mixing in the B meson system [1], and also provide a particularly important place to look for the physics beyond the standard model (SM). For example, the mode B 0 s → J/ψφ is the so-called "golden mode" for measuring β s , which is extracted from the angular analysis of the time-dependent differential decay rate [2][3][4]. The counterpart phase β in the B meson system can also be extracted in B 0 → J/ψK * decay [5]. The decay B 0 → J/ψφ, which would proceed mainly via a Cabibbo-suppressed and color-suppressed transition (bd →ccdd) with rescattering of dd into ss, provide useful information for understanding rescattering mechanisms [6][7][8]. In addition, combining the decays B 0 → J/ψφ and B 0 → J/ψω can be helpful to study the ω − φ mixing [9]. These decays are dominated by tree diagrams and the contributions from penguin diagrams are expected to be small. With continuously increasing high-precision measurements, the penguin effects, which play an important role in the extraction of the above phases, can be measured by means of an analysis of the angular distribution of B 0 → J/ψρ 0 [10] and B 0 s → J/ψK * 0 [11]. In the framework of SM, these decay modes are induced by transitions b → qcc with q = d, s. In principle, any mode involving various excitations of the cc mesons such as B → ψ(2S)V decays could be an alternative to that for J/ψ analogues, and give additional and complementary information. Experimentally, the ψ(2S) meson can be reconstructed in the decay channels ψ(2S) → µ + µ − and ψ(2S) → J/ψπ + π − , with the J/ψ meson decaying into two muons [12]. Nowadays, several experimental Collaborations have measured the decays B 0 s → ψ(2S)φ [13], B → ψ(2S)K * (892) [14][15][16][17][18], B 0 → ψ(2S)π 0 [19], B + → ψ(2S)π + [20], B s → ψ(2S)η ( ′ ) [21,22]. Some relative ratios of the branching ratios for B meson decays into ψ(2S) and J/ψ mesons are also measured by several experiments [23][24][25][26].
On the theory side, these B → ψ(2S)V modes do have some special properties. Since there are three possible values of the total angular momentum with different CP eigenvalues (L = 1 is odd, while L = 0, 2 are even), the angular analysis is needed to separate the contributions from the CP -even and CP -odd partial waves. Therefore the final state can be decomposed into three helicity amplitudes (one longitudinal and two transverse components ). The information about the phases of the transverse decay amplitudes provides a test of the factorization hypothesis [27]. Furthermore, B → ψ(2S)V are the color-suppressed modes and therefore a significant impact of nonfactorizable contributions is expected. Both improvements in the accuracy of the experimental measurements and the observation of new modes, especially involving ψ(2S) in the final states, can be helpful in understanding the role of any nonfactorizable corrections [28][29][30] and differentiating various theory approaches. Nowadays there exist several theoretical approaches as described in Refs. [31][32][33][34][35][36][37][38][39][40][41][42][43][44][45][46] which shed more light on the S-wave ground state charmonium decays of B mesons. The nonleptonic B decays with radially excited charmonium mesons in final state, however, have received less attention in the literature.
Based on the k T factorization theorem, The perturbative QCD (PQCD) approach [47,48] is suitable for describing different types of heavy hadron decays. After including the parton transverse momentum k T , which is not negligible in the end-point region, both factorizable and nonfactorizable contributions are calculable without endpoint singularity. The Sudakov resummation has also been introduced to suppress the long-distance contributions effectively. Therefore, the PQCD approach is a self-consistent framework and has a good predictive power. In our previous works [49][50][51], the semi-leptonic, two-body and three-body non-leptonic decays of the B c (B) mesons to ψ(2S) are studied in the PQCD framework. Here, furthermore, we will extend our previous analysis to the B → ψ(2S)V decays. In a recent work [52], The authors applied the PQCD approach to study B → J/ψV decays and also obtained the theoretical predictions in good agreement with currently available data. Therefore we have good reasons to believe that it is appropriate to analyze B → ψ(2S)V in this framework.
This work is organized as follows. In Sec. II we present some basic formulas such as the effective Hamiltonian and kinetic conventions, then briefly review the pQCD approach. The technical formulas of the calculation and the nonperturbative meson wave functions are summarised in Appendix A and B, respectively. Section III devoted to numerical calculation and discussion. Our conclusions are left for Sec. IV.

II. ANALYTIC FORMULAS AND PERTURBATIVE CALCULATIONS
For nonleptonic charmonium B decays, both the tree operators and the penguin operators of the standard effective weak Hamiltonian contribute, which is given by [53] with the CKM matrix element ξ c(t) = V * c(t)b V c(t)q . O i (µ) and C i (µ) are the effective four quark operators and their QCD corrected Wilson coefficients at the renormalization scale µ, respectively. Their explicit form can be found in Ref. [53].
At quark level, when theb →qcc decay occurs through the four quark operators, a cc state ψ(2S) is created while the other light anti-quarkq is flying away. Since the heavy b quark in B meson carry most of the energy of B meson, the spectator quark of the B meson is soft. A hard gluon is exchanged so that the spectator quark gets energy from the four quark operator and then form a fast moving vector meson with its partner anti-quark. This makes the perturbative calculations into a six-quark interaction, which involves the four quark operator and the spectator quark connected by a hard gluon. The relevant Feynman diagrams are shown in Fig. 1.
In the PQCD approach, the decay amplitudes are expressed as the convolution of the hard kernels H with the relevant meson wave functions Φ i where k i are the momentum of the quark in each meson, and "Tr" denotes the trace over all Dirac structure and color indices. C(t) is the short distance Wilson coefficients at the hard-scale t. The meson wave functions Φ, including all nonperturbative components in the k T factorization, can be extracted from experimental data or other nonperturbative methods. The hard kernel H(k i , t) describes the four quark operator and the spectator quark connected by a hard gluon, which can be perturbatively calculated including all possible Feynman diagrams without endpoint singularity. In the following, we start to compute the decay amplitudes of B → ψ(2S)V decay.
which satisfy the normalization (ǫ L 2,3 ) 2 = (ǫ T 2,3 ) 2 = −1 and the orthogonality ǫ L 2 · P 2 = ǫ L 3 · P 3 = 0. The decay amplitude can be decomposed into three parts of the polarizations amplitudes as follows: with the null vectors n = (1, 0, 0 T ) and v = (0, 1, 0 T ). The subscript L, N, T correspond to the longitudinal, normal and transverse polarization states, respectively. According to Eq. (1), the three different polarization amplitudes have the following expressions, where F (M) describes the contributions from the factorizable (nonfactorizable ) diagrams. The superscript LL, LR, and SP refers to the contributions from and (S − P ) ⊗ (S + P ) operators, respectively. These explicit factorization formulas are all listed in Appendix A. 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 [54]. In the NDR scheme, the vertex corrections are included by the modifications to the combinations a i 1 of the Wilson coefficients C i associated with the factorizable The functions f h I arise from the vertex corrections, which are given in Ref. [31].

III. NUMERICAL RESULTS
To be used in our numerical calculations, those parameters such as meson mass, the Wolfenstein parameters, decay constants, and the lifetime of B (s) mesons [55] are given in Table I, while the input wave functions and various parameters of the light vectors are shown in Appendix B. We now use the method previously illustrated to estimate the physical observables (such as the CP averaged branching ratios, direct and mixing CP violations, polarization fractions, and relative phases ) of the considered decays.
A. The CP averaged branching ratios For B → ψ(2S)V decays, the branching ratios can be written as where the terms A 0 , A , A ⊥ denote the longitudinal, parallel, and perpendicular polarization amplitude in the trasversity basis, respectively, which are related to A L,N,T of Eq. (6) via Here A 0 and A are the CP even amplitudes whereas A ⊥ correspond to CP odd ones. Note that an additional minus sign in A 0 (see Ref. [56] ) make our definitions of the relative phase between A (⊥) and A 0 takes the value of π in the heavy-quark limit. The CP averaged branching ratios for the B → ψ(2S)V decays are shown in Table II together with some of the experimental measurements. Some dominant uncertainties are considered in our calculations. The first error in these entries is caused by the hadronic parameters in the B (s) meson wave function: (1) the shape parameters: ω b = 0.40 ± 0.04 for the B meson, and ω b = 0.50 ± 0.05 for the B s meson; (2) the decay constants, which are given in Table I. The second error is from the uncertainty of the heavy quark masses. In the evaluation, we vary the values of m c(b) within a 10% range. The last one is caused by the variation of the hard scale from 0.8t to 1.2t, which characterizes the size of the NLO QCD contributions. It is found that the main uncertainties in our approach come from the B meson wave function, which can reach 20 − 30% in magnitude. The scale-dependent uncertainty is less than 20% due to the inclusion of the NLO vertex corrections. We have checked the sensitivity of our results to the choice of the shape parameter ω c (see Eq. (B7)) in charmonia meson wave function. The variation of ω c in the range 0.18 ∼ 0.22 will result in a small change of the branching ratio, say less than 10%. In addition, the uncertainties related to the light vector mesons, such as the vector meson decay constants and the Gegenbauer moments shown in Table VI, are only several percent. Therefore they have been neglected in our calculations. For the color-suppressed decays, it is expected that the factorizable diagram contribution is suppressed due to the cancellation of Wilson coefficients C 1 + C 2 /3. After the inclusion of the vertex corrections, the factorizable diagram contributions become comparable with the nonfactorizable ones. Some important features of the numerical results collected in Table II In the quark model, the difference between B 0 → ψ(2S)ω and B 0 → ψ(2S)ρ 0 decays comes from the sign of dd component, which only appears in penguin operators, so their difference should be relatively small. The branching ratio . This is a consequence of the fact that the ω vector and tensor decay constants are smaller than those of the ρ 0 according to Table VI; (III) The value of B(B s → ψ(2S)K * 0 ) have a tendency to be smaller than 2B(B 0 → ψ(2S)ρ 0 ). Although the K * and B s meson decay constants are larger than those of the ρ 0 and B 0 meson, the SU(3) breaking effects in the twist-2 distribution amplitudes, parametrized by the first Gegenbauer moment a 1K * (see Eq. (B10)), gives a negative contribution to the B s → ψ(2S)K * 0 decay, which induces the smaller branching ratio.
(IV) For the first four B (s) → ψ(2S)V decays as listed in Table II, one can see that the PQCD predictions for their branching ratios agree well with the world averaged values as given in HFAG 2016 and PDG 2016 [55,59] within one standard deviation. For B s → ψ(2S)K * 0 decay, the central value of our theoretical prediction for its branching ratio is slightly smaller than that of the PDG number [55]. But we know that the PDG result is obtained by multiplying the best value B(B 0 → ψ(2S)K * 0 ) with the measured ratio [60]. We hope the future experiment will provide a direct measurement to this mode.
(V) As for the channels with ρ and ω as the final state V meson, they have not been measured yet. The pQCD predictions for the decay rates of these three channels are at the order of 10 −5 , measurable in the future LHCb and Belle-II experiments.
For a more direct comparison with the available experimental measurements of the relative rates of B (s) meson decays into ψ(2S) and J/ψ mesons, we recalculated the corresponding B (s) decays to J/ψV by using the same input parameters as in this paper but with the replacement ψ(2S) → J/ψ, and we found numerically that where the errors have the same meaning as those for B (s) → ψ(2S)V decays. The above results are well consistent with the previous PQCD calculations [52] and also the present data [55]. Finally, as a cross-check, using the PQCD predictions as given in Table II and Eq. (10) we can estimate the relative ratios where all uncertainties are added in quadrature. Since the parameter dependences of the PQCD predictions for the branching ratios are largely canceled in their relative ratios, the total theoretical error of R V are only a few percent, much smaller than those for the branching ratios. Fortunately, two of these five ratios have been measured by LHCb [24] D0 [25], and CDF [14,26] experiments: where the third uncertainty is from the ratio of the ψ(2S) and J/ψ branching fractions to µ + µ − . It is easy to see that our PQCD predictions for both R φ and R K * 0 agree very well with the measured values.

B. CP ASYMMETRIES
Studying CP asymmetries is an important task in B physics. For the charged B decays, the CP asymmetries arise from the interference between the penguin diagrams and tree diagrams. The direct CP violation asymmetry including three polarization are defined by whereĀ 0, ,⊥ is the CP -conjugate amplitude of A 0, ,⊥ . For the neutral B 0 (s) decays, because of the B 0 (s) −B 0 (s) mixing, it is required to include time-dependent measurements in CP violation asymmetries. If the final states are CP eigen states, the time-dependent CP asymmetry is defined as where ∆m is the mass difference of the two mass eigenstates of the neutral B meson and f is a two-body final state. The direct CP asymmetry C 0, ,⊥ f and mixing-induced CP asymmetry S 0, ,⊥ f are referred to as  in the B 0 → ψ(2S)(ρ 0 , ω, φ) decays. The error arises from the hard scale t. Modes The parameter λ 0, ,⊥ f = η f e −2iβ (s)Ā 0, ,⊥ A 0, ,⊥ describes CP violation in the interference between mixing and decay. η f is the CP eigenvalue (±1) of the polarization state. β (s) is the CKM angle defined as usual [55]. Note that the final states of ψ(2S)K * 0 and its CP conjugate are flavor-specific, for example, the kaon and pion charges of K * 0 → K + π − andK * 0 → K − π + depend on whether we had a B andB meson in the initial state, and the time-dependent angular distributions do not show CP violation due to interference between mixing and decay. Therefore, we only calculate the direct CP asymmetry for B 0 s → ψ(2S)K * 0 and B 0 → ψ(2S)K * 0 decays. The PQCD predictions for the CP asymmetry parameters A dir 0, ,⊥ are listed in Table III and IV. Unlike the branching ratios, the direct CP asymmetry is not sensitive to the wave function parameters and heavy quark masses, but suffer from large uncertainties due to the hard scale t. In order to reduce the large scale dependence effectively, one has to know the complete NLO corrections, which are in fact not yet available now, thus beyond the scope of this paper.
Since the direct CP asymmetry is proportional to the interference between the tree and penguin contributions, while the Wilson coefficients of the penguin diagram are loop suppressed when compared with those tree contributions. Therefore, the direct CP asymmetry parameters of these processes are rather small ( only 10 −3 ∼ 10 −4 ), and the mixing-induced CP asymmetry for neutral B decays is almost proportional to the sin 2β (s) from Eq. (15). The mixinginduced CP asymmetry parameters S ψ(2S)ρ 0 (ω) and S ψ(2S)φ in Table IV are very close to the current world average values − sin 2β = −0.691 ± 0.017 and −2β s = −0.0376 +0.0008 −0.0007 [59], respectively. That is to say, these modes can serve as an alternative places to extract CKM angle β (s) . Furthermore, the large mixing-induced CP asymmetry S f for b → d transition can confront with future experimental results. It can also be seen that the CP asymmetry parameters for three polarization states are slightly different because the strong phases coming from the non-factorizable diagrams and vertex corrections are polarization-dependent [31]. On experimental side, so far only the charge asymmetries of B + → ψ(2S)K * + process was measured by BaBar Collaboration [55]: Of course, the statistical uncertainty is too large to make any statement. Any observation of large direct CP asymmetry for the considered decays B (s) → ψ(2S)V decays will be a signal for new physics. Besides, the precise measurements of these mixing-induced CP asymmetries serve to determine the CP phases related to the B 0 −B 0 and B 0 s −B 0 s mixing amplitudes.

C. Polarization fractions and relative phases
In experimental analyses, we usually define five observables corresponding to three polarization fractions f 0 , f , f ⊥ , and two relative phases φ , φ ⊥ , where  [15] 0.40 ± 0.14 ± 0.07 ----Belle [58] 0.455 +0.031+0.014 with normalisation such that f 0 + f + f ⊥ = 1. The polarization fractions as well as relative phases are shown in Table V, where the sources of the errors in the numerical estimates have the same origin as in the discussion of the branching ratios in Table II. It is easy to see that the most important theoretical uncertainties are caused by the heavy quark masses. From Eqs. (7) and (A10), we can see the mass terms m b and m c associated vertex corrections and nonfactorizable amplitudes, respectively. It can numerically change the real and imaginary parts of these contributions and have a significant effect on the polarization fractions, especially for the relative phases. The uncertainties from the wave function parameters are very small because they mainly give an overall change of all polarization amplitudes and the parameter dependence can be canceled out in Eq. (17). From Table V, both the B + → ψ(2S)(K * + , ρ + ) and B 0 → ψ(2S)(K * 0 , ρ 0 ) modes have the same polarization fractions and relative phases, since they differ only in the lifetimes or isospin factor in our formalism. Comparing the three polarization fractions, the perpendicular polarization fractions f ⊥ are less than 25% shows that the CP even component dominates in these decays. According to the power counting rules in the factorization assumption, the longitudinal polarization dominates the decay ratios and the transverse polarizations are suppressed [61] due to the helicity flips of the quark in the final state hadrons. However, the situation is very different for the color-suppressed decays, where the contributions from the nonfactorizable tree diagrams in Figs. 1(c) and 1(d) are comparable with those of the color-suppressed tree diagrams although the latter are enhanced by the involving vertex corrections. With an additional gluon, the transverse polarization in the nonfactorizable diagrams does not encounter helicity flip suppression, therefore numerically we get a longitudinal polarization fraction (f 0 ) smaller than 50%, which are compatible with those currently available data. The fact that the nonfactorizable diagrams can give a large transverse polarization contribution is also observed in the B c → J/ψD * + (s) decays [62]. There are another equivalent set of helicity amplitudes (A 0 , A + , A − ), which are related to the spin amplitudes (A 0 , A , A ⊥ ) introduced in Eq. (9) by while A 0 is common to both bases. It is expected that |A 0 | 2 > |A + | 2 > |A − | 2 if the two final-states are both light vector mesons. The larger the mass of the vector-meson daughters, the weaker the inequality. In B → ψ(2S)V decays with light V being a recoiled meson and heavy ψ(2S) an ejected one. The positive-helicity amplitude is suppressed by m ψ(2S) /M (almost of order unity) due to one of the quark helicities in ψ(2S) has to be flipped, while the negative-helicity one is subject to a further chirality suppression of order m V /M [61]. Therefore, A + and A 0 can be comparable and larger than A − . Using values of Table V and Eq. (18), the pQCD predictions do favor the hierarchy pattern |A 0 | 2 ∼ |A + | 2 > |A − | 2 .
The angular analysis of B 0 → ψ(2S)K * 0 and B 0 s → ψ(2S)φ has been carried out by BaBar [18] and LHCb [13], respectively. The obtained polarization observables are also summarized in Table V. As expected under SU (3)-flavor symmetry, both decay modes have similar magnitudes and phases of the amplitudes. Our results of polarization fractions can accommodate the data well within uncertainties, while the predicted relative phases are a bit smaller than the data. One can find a shift from π at the 6−7σ level in φ and φ ⊥ shows the existence of final-state interaction. However, the f − f ⊥ is about 5% and the difference between φ and φ ⊥ dose not exceed 0.3 radians, which suggest that our solutions are consistent with approximate s-quark helicity conservation despite substantial strong phases.
For the B 0 s → ψ(2S)K * 0 channel, the LHCb Collaboration [60] has reported the longitudinal polarisation fraction f 0 as 0.524 ± 0.056 ± 0.029, but a thorough angular analysis is still missing. As for other modes, we obtain reasonably accurate results, which could be tested by future experimental measurements.

IV. CONCLUSION
In this paper we have investigated the seven B → ψ(2S)V decay modes carefully by employing the PQCD factorization approach. Besides the color-suppressed factorizable diagrams, the nonfactorizable diagrams and the vertex correction diagrams can also be evaluated in this approach.
The predicted branching ratios and the relative rates of B meson decays into ψ(2S) and J/ψ mesons are compared with experiments wherever available. Our results indicate that the direct CP asymmetries in these channels are very small due to the suppressed penguin contributions as we mentioned above. The mixing-induced CP asymmetries are not far away from sin 2β (s) , these channels can therefore play an important role in the extraction of the CKM angle β (s) .
Finally, we made a comprehensive polarization analysis of the considered decays. The predicted polarization fractions and relative phases of B 0 → ψ(2S)K * 0 and B 0 s → ψ(2S)φ decays are consistent with data. Due to the large mass of ψ(2S) and the dominant contributions from the nonfactorizable diagrams, we obtain an equal amount of transverse and longitudinal polarization. The pattern of f ≈ f ⊥ , φ ≈ φ ⊥ favor the conservation of light quark helicity. The deviations from π at several standard deviations in φ and φ ⊥ indicate the existence of the still unknown final-state interaction.
We also discussed theoretical uncertainties arising from the hadronic parameters in B meson wave function, heavy quark masses and hard scale t. The total uncertainties are acceptable, around 30% in magnitude. The uncertainties from the hadronic parameters can give sizable effects on the PQCD predictions for branching ratios, while the CP asymmetries suffer a large error from the hard scale t. The further studies at the completely NLO level are certainly required to improve the accuracy of the theoretical predictions. Furthermore, the polarization observables f 0, ,⊥ and φ ,⊥ are more sensitive to the heavy quark masses, which suggest that the color-suppressed type decays may be more sensitive to the vertex corrections and nonfactorizable contributions. Our results and findings will be further tested by the LHCb and Belle-II experiments in the near future.
with r c = m c /M and m c is the charm quark mass; C f = 4/3 is a color factor; f ψ is the decay constant of the ψ(2S) meson. The coefficient (−) 1 √ 2 appears for B → ψ(2S)(ρ 0 )ω decay, because only the d quark component of the (ρ 0 )ω meson is involved. We neglect terms higher than r 2 v orders, since the vector light cone wave functions derived from sum rules are expanded to this order [63]. The functions h come from the Fourier transform of virtual quark and gluon propagators. They are defined by where J 0 is the Bessel function and K 0 , I 0 are modified Bessel function with K 0 (ix) = π 2 (−N 0 (x) + iJ 0 (x)). α e and β a,b,c,d are the virtuality of the internal gluon and quark, respectively. Their expressions are The hard scale t is chosen as the maximum of the virtuality of the internal momentum transition in the hard amplitudes, including 1/b i (i = 1, 2, 3): The Sudakov factors can be written as where the function s(Q, b) is given in [64]. γ q = −α s /π is the anomalous dimension of the quark. The threshold resummation factor S t (x) is adopted from [65], with a running parameter c(Q 2 ) = 0.04Q 2 − 0.51Q + 1.87 [66] and Q 2 = M 2 (1 − r 2 ) [67].

Appendix B: THE WAVE FUNCTIONS
In the PQCD approach, the necessary inputs contain the light-cone distribution amplitudes (LCDAs) which are constructed by the nonlocal matrix elements. The B u,d,s meson light-cone matrix element are decomposed into the following two Lorentz structures [54]: with the color factor N c . As usual the former Lorentz structure in above equation is the dominant contribution in the numerical calculations, while the latter Lorentz structure is negligible [68]. In impact coordinate space the B meson wave function can be expressed by [65,69] where b is the conjugate variable of the transverse momentum of the valence quark of the meson. The distribution amplitude φ B (x, b) as being used in Refs. [56,65] are adopted here with the shape parameter ω b and the normalization constant N being related to the decay constant f B by normalization: The shape parameter ω b = 0.40 ± 0.04 GeV for the B u,d mesons and ω b = 0.50 ± 0.05 GeV for the B s meson. For the ψ(2S) meson, the longitudinally and transversely polarized LCDAs up to twist-3 are defined by [49,50] ψ(2S)(P 2 , ǫ L 2 )|c(z) α c(0) β |0 = The asymptotic models for the twist-2 distribution amplitudes ψ L,T and the twist-3 distribution amplitudes ψ V,t are extracted from the correspond Schrödinger states for the harmonic-oscillator potential. Their expressions have been derived as [49]. with where the parameter ω c = 0.20 ± 0.02 GeV. N i (i = L, T, t, V ) are the normalization constants and the normalization conditions: For a light vector meson, the light-cone wave function for longitudinal (L) and transverse (T) polarization are written as [63] Φ respectively, where ǫ 0123 = 1 in our convention. Note that v is the moving direction of vector particle. The twist-2 distribution amplitudes are given by and those of twist-3 ones are with t = 2x − 1. The vector (tensor) decay constants f V (f T V ) together with the Gegenbauer moments [70] are shown numerically in Table VI. Note that positive a ,⊥ 1 refer to aK * 0 containing an s quark, while for a K * + (K * 0 ) with an s quark, a ,⊥ 1 changes sign [71].