$B_{s}\to K \ell \nu_\ell$ and $B_{(s)} \to \pi (K) \ell^+\ell^-$ decays at large recoil and CKM matrix elements

We provide hadronic input for the $B$-meson semileptonic transitions to a light pseudoscalar meson at large recoil. The $B_s\to K$ form factor calculated from QCD light-cone sum rule is updated, to be used for a $|V_{ub}|$ determination from the $B_s\to K \ell \nu$ width. Furthermore, we calculate the hadronic input for the binned observables of $B \to \pi \ell^+ \ell^-$ and $B \to K \ell^+ \ell^-$. In addition to the form factors, the nonlocal hadronic matrix elements are obtained, combining QCD factorization and light-cone sum rules with hadronic dispersion relations. We emphasize that, due to nonlocal effects, the ratio of branching fractions of these decays is not sufficient for an accurate extraction of the $|V_{td}/V_{ts}|$ ratio. Instead, we suggest to determine the Wolfenstein parameters $A,\rho,\eta$ of the CKM matrix, combining the branching fractions of $B \to K \ell^+ \ell^-$ and $B \to \pi \ell^+ \ell^-$ with the direct $CP$-asymmetry in the latter decay. We also obtain the hadronic matrix elements for a yet unexplored channel $B_s \to K \ell^+ \ell^-$.


Introduction
Determination of CKM matrix elements from the semileptonic decays of B meson remains a topical problem. Most importantly, one has to clarify the origin of the tension between the |V ub | values extracted from the exclusive B → π ν and inclusive B → X u ν decays (see e.g., the review [1]). The B → π vector form factor f + Bπ (q 2 ) is the only theory input sufficient for the |V ub | determination from B → π ν . This hadronic matrix element is calculated in the lattice QCD at small recoil of the pion (at large q 2 ) or from QCD lightcone sum rules (LCSRs) at large recoil of the pion (at small and intermediate q 2 ).
Apart from increasing the accuracy of the form factor calculation, it is important to extend the set of "standard" exclusive processes used for |V ub | determination. The B s → K * (→ Kπ) ν decay, as one possibility, was discussed in [2]. A simpler process is the B s → K ν decay, where the data are anticipated from LHCb collaboration. Our first goal in this paper is to provide this decay mode with a hadronic input, updating the calculation of the B s → K form factors from LCSRs. This method [3][4][5] is based on the operator-product expansion (OPE) of a correlation function expressed in terms of lightmeson distribution amplitudes (DAs) with growing twist. The violation of the SU (3) f l symmetry in B s → K with respect to B → π transition emerges in LCSRs due to the s-quark mass effects in the correlation function, including the asymmetry between the sand {u, d}-partons in the kaon DAs. Earlier LCSR results on the B s → K form factors can be found in [6], where the NLO corrections to the correlation function computed in [7] were taken into account. In this paper we update the LCSRs for f + BsK (q 2 ) and also for the tensor form factor f T BsK (q 2 ). In particular, we correct certain terms in the subleading twist-3,4 contributions to LCSRs for both vector and tensor form factors. In parallel, we recalculate the B → K and B → π form factors using a common set of input parameters, e.g., the updated [8] 2-point QCD sum rule for the decay constants f B and f Bs . Importantly, the twist-5,6 corrections to the LCSRs estimated by one of us [9] are negligibly small, ensuring the reliability of the adopted twist ≤ 4 approximation.
The calculated form factors are then used to address the second goal of this paper: determination of CKM parameters from the flavour-changing neutral current (FCNC) decays B → K + − , B → π + − and B s → K + − . Recently, |V td |, |V ts | and their ratio were determined by LHCb collaboration [10] from the measured B → π + − and B → K + − partial widths. We suggest to make the extraction of CKM parameters from these decays more accurate and comprehensive. As well known, in addition to the semileptonic form factors, the hadronic input in FCNC decays includes also nonlocal hadronic matrix elements emerging due to the electromagnetic lepton-pair emission combined with the weak transitions. These hadronic matrix elements in the B → π + − decay amplitude are multiplied by the CKM parameters other than V td , making the determination of the latter not straightforward. We take into account the nonlocal hadronic effects in B → K + − and B → π + − , employing the methods used in [12], [13] and originally suggested in [11]. The nonlocal hadronic matrix elements are calculated at spacelike q 2 , using OPE, QCD factorization [14] and LCSRs, and are then matched to their values at timelike q 2 via hadronic dispersion relations. The results of this calculation are reliable at large hadronic recoil, below the charmonium region, that is, at q 2 < m 2 J/ψ . Here we also extend the calculation of nonlocal effects to the previously unexplored channel B s → K + − .
The binned widths and direct CP -asymmetries of FCNC semileptonic decays are then expressed in a form combining the CKM parameters with the quantities determined by the calculated hadronic input. Here we find it more convenient to switch to the Wolfenstein parametrization of the CKM matrix. In this form, three observables: the width of B → K + − , the ratio of B → π + − and B → K + − widths and the direct CP -asymmetry in B → π + − , are sufficient to extract the three Wolfenstein parameters A, η and ρ from experimental data, provided the parameter λ is known quite precisely. Two additional observables for the same determination are given by the yet unobserved B s → K + − decay. The current data on the B → K + − and B → π + − decays are not yet precise enough to yield the CKM parameters with an accuracy comparable to the other determinations. Hence, here we limit ourselves with the Wolfenstein parameters taken from the global CKM fit and predict the binned observables of all three FCNC decays in the optimal interval 1.0 GeV 2 < q 2 < 6.0 GeV 2 of the large recoil region.
In what follows, in Sect. 2 we specify and discuss the hadronic input and observables in the exclusive semileptonic B (s) decays. In Sect. 3 we present the numerical results and Sect. 4 is devoted to the final discussion. In the Appendices, we briefly recapitulate the calculation (A) of the form factors from LCSRs and (B) of the nonlocal hadronic matrix elements.

Observables in semileptonic B (s) decays and CKM parameters
The form factors of semileptonic transitions of B-meson to a light pseudoscalar meson P = π, K are defined in a standard way: where p µ and q µ are the four-momenta of the P -meson and lepton pair, respectively, and the vector and scalar form factors coincide at q 2 = 0, that is, f + BP (0) = f 0 BP (0). We start from the weak semileptonic decayB s → K + ν , where the hadronic input for = e, µ in the m = 0 approximation is given by the vector form factor f + BsK . We use the following quantity related to the differential width integrated over an interval 0 ≤ q 2 ≤ q 2 0 : is the 3-momentum of P meson in the rest frame of B meson. Our choice for the integration interval is q 2 0 = 12.0 GeV 2 , covering the region where the LCSRs used for the calculation of the form factors (see Appendix A) are valid. The same interval was adopted for the analogous quantity ∆ζ Bπ [0, q 2 0 ] for B → π ν calculated in [15,16]. The numerical estimate of ∆ζ BsK [0, q 2 0 ] presented in the next section can be directly used for |V ub | determination, provided the integrated branching fraction on the r.h.s. of Eq. (2.3) is measured.
Turning to semileptonic decays generated by the b → s(d) + − transitions ( = e, µ), we use a generic notationB → P + − for the three channels: B − → K − + − , B − → π − + − andB s → K 0 + − 1 , denoting the CP conjugated channels by B →P + − . The decay amplitude can be represented in the following form: where λ (q) p = V pb V * pq (p = u, c, t; q = d, s), m = 0, and we use unitarity of the CKM matrix, fixing hereafter λ u ). In Eq. (2.4) we introduce a compact notation: where m q is the mass of d or s-quark and In addition, we introduce the phase difference of the hadronic amplitudes defined above: BP (q 2 ) parametrize the nonlocal contributions to B → P + − , generated by the current-current, quark-penguin and chromomagnetic operators in the effective Hamiltonian, combined with an electromagnetically produced lepton pair. The definition of nonlocal amplitudes is given in Appendix B, where also the method of their calculation is briefly explained. In Refs. [11,12], this part of hadronic input was cast in the form of an effective (process-and q 2 -dependent) addition ∆C BP 9 (q 2 ) to the Wilson coefficient C 9 . Here, as in Ref. [13], we find it more convenient to separate the parts proportional to λ Squaring the amplitude (2.4) and integrating over the phase space, one obtains for the q 2 -binned branching fraction, defined as: the following expression: where the ratio of CKM matrix elements is parametrized in terms of its module and phase: and we use the following notation for the phase-space weighted and integrated parts of the decay amplitude squared: (2.14) The binned branching fraction for the CP -conjugated mode B →P + − is obtained from Eq. (2.9) by changing the sign at the term proportional to sin ξ q . Furthermore, we consider two binned observables: the CP -averaged branching fraction: (2.15) and the corresponding direct CP -asymmetry: . (2.16) In Eqs. (2.15) and (2.16) the CKM-dependent coefficients are conveniently separated from the quantities F BP , D BP , C BP , S BP , which contain the calculable hadronic matrix elements, Wilson coefficients and kinematical factors. In the next section we present numerical results for these quantities for a definite q 2 -bin in the large-recoil region.
Turning to the observables for the specific decay channels, we neglect λ (s) u , hence, put κ s = 0 and obtain for B → K + − : with vanishing CP asymmetry. For B − → π − + − and its CP -conjugated process both observables, and  (2.19), with Bπ replaced by B s K. Here we do not consider the decaysB 0 →K 0 + − and B 0 → π 0 + − , which are the isospin counterparts of, respectively, B − → K − + − and B − → π − + − and can be treated in a similar way (see [12,13]). We also postpone to a future study the time-dependent CP -asymmetry in theB s → K 0 + − decay.
Dividing Eq. (2.18) by Eq. (2.17), we notice that an accurate extraction of the ratio can only be achieved if the contributions of process-dependent nonlocal effects are taken into account for both decay modes. Moreover, this ratio depends also on the other CKM parameters, most importantly, on the V ub value 2 .
Here we suggest a different, more systematic way to extract the parameters of CKM matrix from the observables (2.17)- (2.19). First of all, we find it more convenient to switch to the four standard Wolfenstein parameters λ, A, ρ and η defined as in [17]. The relevant CKM factors can be represented as follows: where we neglect very small O(λ 4 ) corrections to these expressions 3 .
Hereafter, we suppose, that the parameter λ, precisely determined from the global CKM fit [17], is used as an input. Then, it is possible to extract all three remaining Wolfenstein parameters combining the three observables (2.17)-(2.19) for semileptonic FCNC decays. First, the parameter A is determined from the binned branching fraction of B → K + − , as follows after substituting Eq. (2.20) in Eq. (2.17): Then, combining the ratio of the B → π + − and B → K + − binned branching fractions with the CP -asymmetry of the pion mode, and employing Eqs.

Parameter
Ref.

Numerical Results
The most important input parameters used in our numerical analysis are listed in Table 1.
In particular, the electroweak parameters, the strong coupling and the meson masses are taken from [17]. For the quark masses in M S scheme, entering the correlation functions for QCD sum rules, we adopt, following, e.g., [8], the intervals covering the non-lattice determinations in [17]. We put m u,d = 0, except in the combination µ π(K) = m 2 π(K) /(m u + m d(s) ) entering the pion and kaon DAs. In LCSRs, parameters of the pion and kaon twist-2 DA's include the decay constants, and the Gegenbauer moments a π 2,4 and a K 1,2 . Normalization of the twist-3 DAs is determined by µ π,K , where the ChPT relations [20] between light-quark masses are used (see e.g., [18]). The remaining parameters of the twist-3 and twist-4 DAs, not shown in Table 1 for brevity, are taken from [21], they were also used in [11,15,18]. Furthermore, in LCSRs the renormalization scale µ and the Borel parameters M for the sum rules with B (B s ) interpolating current quoted in Table 1 are chosen, largely following [15]. The effective quark-hadron duality threshold is determined calculating the B (s) -meson mass from the differentiated LCSR. The decay constants f B and f Bs entering LCRSs are replaced by the two-point sum rules in NLO, their expressions and input parameters (in particular, the vacuum condensate densities) are the same as in [8]. The intervals obtained from these sum rules in NLO are MeV. Note that the above uncertainties are effectively smaller in LCSRs (Eq. (A.4) in Appendix A) due to the correlations of common parameters.
Using the input described above, we obtain the updated prediction for the B s → K vector and tensor form factors in the region 0 ≤ q 2 ≤ 12.0 GeV 2 where the OPE for LCSRs in the adopted approximation is reliable (see Appendix A). In parallel, we also recalculate the B → K and B → π form factors. For convenience, we fit the LCSR predictions for the B → P form factors in this region to the two-parameter BCL-version of z-expansion [22] in the form adopted in [18]:  Table 2. Note that, adopting a more complicated z-expansion with more slope parameters, only insignificantly changes the quality of the fit, and reveals strong correlations between these parameters. In any case the actual form of parametrization does not play a role as soon as we stay within the q 2 -region where the form factors are directly calculated from LCSRs. Our results for the form factors are also plotted in Fig. 1, where the error bands correspond to the uncertainties of the fitted parameters shown in Table 2. For comparison, we also show in the same figures the extrapolations of the recent lattice QCD results obtained at large q 2 (low hadronic recoil) and continued to the small q 2 region using the z-series parametrization. For the vector B s → K form factor this extrapolation was obtained by HPQCD Collaboration [23]. The same form factor was also calculated by ALPHA Collaboration [24] at a single large-q 2 value. For the vector and which is, as it should be, very close to the interval predicted in [16]. The latter interval is somewhat narrower than (3.5), reflecting the statistical (Bayesian) treatment applied in [16] which generally produces less conservative errors. In the future, when sufficiently accurate data on B s → K ν become available, a global statistical treatment of all B → P form factors is desirable. Comparing our results in Table 2 with the earlier LCSR calculation [6] of the B → K and B s → K form factors, we emphasize that, albeit the numerical results look close to ours, there are differences in the subleading twist-3,4 terms. We follow Ref. [18] where these terms have already been discussed and corrected. Also, as compared to [6], we use slightly different B (s) decay constants and twist-3 normalization parameter µ K .
Furthermore, the interval for our updated result for the B → K vector form factor in Table 2 lies somewhat above the previous LCSR prediction [11], f + BK (0) = 0.34 +0.05 −0.02 , mainly due to the smaller value of f B from the two-point sum rule used here and due to the slightly smaller value of the effective threshold in LCSR used in [11]. On the other hand, in the LCSR for f T BK (q 2 ) some minor corrections, implemented here in the subleading twist-4 terms, largely compensate the shift caused by the B-decay constant, so that our result in Table 2 is close to f T BK (0) = 0.39 +0.05 −0.03 obtained in [11]. Turning finally to the LCSR result for the vector B → π form factor, which was updated several times in past, let us mention that although we use the same analytical expressions as in Ref. [7], the input parameters such as µ π (determined by the light quark masses) and Gegenbauer moments a π 2 , a π 4 became more accurate, leading to a narrower interval of our prediction, compared to the interval f + Bπ (0) = 0.26 +0.04 −0.03 obtained in Ref. [7]. The central value of the latter is somewhat below the one we present in Table 2, since we use a smaller (larger) central input value of f B (of µ π ). In Ref. [7] one can also find a detailed comparison with the LCSR B → π form factor obtained earlier in Ref. [29].
We turn to the numerical analysis of B → P + − observables, where the B → P form factors obtained above are used. We recalculate the nonlocal amplitudes, following [12,13]. In Appendix B a brief outline of the calculational method is given. Here we need some additional input parameters. The most important are: the inverse moment λ B of the Bmeson DA (we assume λ Bs = λ B ) and the Borel and threshold parameters in π, K channel in the LCSRs for the soft-gluon emission contributions. They are displayed in Table 1. The same input parameters for the pion, kaon and B-meson DAs as the ones given in Table 1 serve as an input in the hard-gluon contributions for which we use the QCD factorization expressions [14] at spacelike q 2 .
The effective FCNC Hamiltonian (see Eq. (B.1) in Appendix B) is chosen as in [13] (see Table V there), with all Wilson coefficients C i taken at leading order in α s . This accuracy is sufficient for C 1−6 , C ef f 8 entering the nonlocal hadronic amplitudes, having in mind the overall accuracy of our method for these amplitudes. At the same time, the numerically large Wilson coefficients C 9 , C 10 and C ef f 7 of the FCNC operators multiplying the factorizable parts of the decay amplitudes, have a noticeable impact on the observables. Therefore, we adopt here the values of these coefficients at the next-to-leading order in α s (see Table 3).  Table 3. Wilson coefficients of the FCNC operators at next-to-leading (leading) order in α s used in our numerical analysis at various scales.
For completeness and future use, in Appendix B the numerical results for the separate nonlocal amplitudes H  Table 4, where the (uncorrelated) uncertainties are obtained by adding in quadrature the individual variations due to changes of input parameters.
Note that the binned quantities F BP are not much sensitive to the magnitude of the nonlocal amplitudes H (c) BP (q 2 ), which enter the numerically subleading contributions to the coefficients c BP (q 2 ). Hence, the differences between F BK , F Bπ and F BsK in Table 4    binned quantities D BP , C BP and S BP are essentially determined by the nonlocal effects in B → P + − . In particular, the large differences between D Bπ , C Bπ , S Bπ and D BsK , C BsK , S BsK emerge mainly due to the enhancement of the weak annihilation mechanism in the nonlocal amplitude H (u) [13]. The same mechanism does not play a role in the amplitude H (u) BsK (q 2 ) contributing toB s → K 0 + − , due to a different quark content of the initial B s meson, and due to a suppressed combination of Wilson coefficients.
As shown in the previous section, the binned quantities F BP , D BP , C BP , S BP can in principle be used for an independent determination of the Wolfenstein parameters A, η and ρ from the combination of observables measured in B → P + − decays. The important role in this determination is played by the direct CP -asymmetry in B → π + − which is not available yet in the large-recoil region bins. Hence, here we limit ourselves by an inverse procedure. Taking the values of all Wolfenstein parameters λ = 0.22506 ± 0.00050, A = 0.811 ± 0.026, from the global fit of CKM matrix [17] and using the calculated hadronic input from Table 4, we predict the values of the binned branching fractions presented in Table 5 and the binned direct CP -asymmetries: The numerical results for the B → K + − and B → π + − decays presented here update the previous ones obtained, respectively, in [12] 4 and [13].

Discussion
In this paper we updated the LCSR predictions for the B s → K form factors in the large recoil region of the kaon. We predicted the ratio of the integrated B s → K ν decay width and |V ub | 2 . Our result can be used to determine this CKM matrix element from the future data on B s → K ν in the kinematically dominant large recoil region. 4 Note that the branching fractions given in the literature are adjusted to our definition, which implies division by the width (q 2 2 − q 2 1 ) of the bin  Table 5. We also calculated the hadronic input for the branching fractions and direct CPasymmetries of B → P + − FCNC decays in the large recoil bin 1.0 ≤ q 2 ≤ 6.0 GeV 2 . Our results include the B → P form factors and nonlocal hadronic matrix elements, all obtained in the same framework and with a uniform input. The LCSRs used in this calculation take into account the soft-overlap nonfactorizable contributions to the form factors and nonlocal amplitudes. Extending the application of LCSRs to other nonlocal contributions represents an important task for the future. For example, as discussed in more detail in [13], the weak annihilation contribution which is important in the B → π + − decay can be obtained from LCSRs with B-meson DAs, alternative to QCD factorization and potentially including subleading effects.
Furthermore, we suggested a systematic way to extract the CKM matrix elements, cast in a form of the Wolfenstein parameters, from the combination of observables in B → P + − decays, independent of the other methods involving the nonleptonic B-decays and/or B −B mixing.
Note that an independent extraction of CKM parameters is also possible from other modes of FCNC exclusive B-decays, such as The corresponding combinations of observables demand, apart from B → V form factors, a dedicated calculation of all relevant nonlocal hadronic matrix elements. For this not yet accomplished task, a variety of methods combining QCD factorization with various versions of LCSRs may prove to be useful. In case of radiative decays the sum rules with photon and vector-meson DAs and heavy-meson interpolating currents can be also of use (for previous works in this direction see [35][36][37]).
In Table 5 we compare our results for the binned branching fractions 4 with the experimental measurements and lattice QCD predictions [28,38,39]. In the lattice QCD studies of B → P + − decays, as explained in detail in [39], the nonlocal contributions cannot be calculated in a fully model-independent way. Instead, the (continuum) QCD-factorization [14] in the timelike region of q 2 is employed. Let us also mention in this context the earlier estimates of B → K [40,41] and B → π [42] where the QCD-factorization approach was used combined with various inputs and extrapolations for the form factors.
As seen from Table 5, the theory predictions for the B → K + − branching fraction reveal some tension with the experimentally measured values, making this observable an important ingredient of the global fits of rare B decays (see e.g., [43]). Adding the characteristics of B → π + − and B s → K + − decays to the set of fitted observables will further extend the possibilities to test the Standard Model in the quark-flavour sector. The fact that these very rare B-decay modes are within the reach of LHCb experiment, makes this task realistic.

A LCSR calculation of the B → P form factors
The LCSRs for B → P (P = π, K) form factors at large recoil of P (parametrically, at q 2 m 2 b ) are derived from the correlation function of the weak flavour-changing current and B-interpolating quark current, sandwiched between the vacuum and on-shell P -state: where the quark-flavour combination q 1 = u, q 2 = s corresponds to theB s → K + weak transition; q 1 = s, q 2 = u and q 1 = d and q 2 = u (q 2 = s) correspond, respectively to the The invariant amplitudes F BP (q 2 , (p + q) 2 ) and F T BP (q 2 , (p + q) 2 ) in (A.1) are used to derive the LCSRs for the vector f + BP (q 2 ) and tensor f T BP (q 2 ) form factors, respectively. At q 2 m 2 b and (p + q) 2 m 2 b the OPE near the light-cone x 2 0 is applied for the correlation function (A.1) and the result is cast in a form of convolution, e.g.,:  ∼ O(m s ) accuracy is adopted. The factorizable parts of twist-5,6 contributions to LCSRs for B → P form factors were calculated by one of us [9] and their numerical impact on the total invariant amplitude was found negligible, < 0.1% of the total. This strengthens the argument for using a truncated twist expansion to the accuracy t = 4.
The NLO O(α s ) corrections to the twist-2 and (two-particle) twist-3 hard-scattering amplitudes T (2,3) 1 are taken into account. In the latter we neglect the s-quark mass, hence, the double suppressed O(α s m s /Λ) effects. We use the expressions for OPE derived in [7] extending them to the B → K and B s → K cases (see also [18]). We do not include the O(β 0 ) estimate of the twist-2 O(α 2 s ) contribution to the twist-2 hard-scattering amplitude calculated in [44], since the resulting effect in LSCR is very small and does not yet represent a complete NNLO computation of T  We equate the contribution of the excited and continuum B-states in the hadronic dispersion relation to the part of the above integral at s > s B 0 , where s B 0 is the effective, process-dependent threshold. The integral at s ≤ s B 0 is then equated to the contribution of the ground-state of B-meson. The subsequent Borel transformation with respect to the variable (p + q) 2 exponentiates denominators, so that, e.g., 1/[s − (p + q) 2 ] → e −s/M 2 . Here M 2 is the Borel parameter chosen so that M 2 ∼ Λm b ∼ µ 2 guarantees a power suppression of higher-twist contributions. One finally obtains the LCSRs for the B → P form factors: ds ImF The effective weak Hamiltonian of the b → q + − transitions (q = d, s) generating the B → P + − decays has the following form in the Standard Model (see e.g., the review [45]):   In the case of the B → K + − decay only the amplitude H (c) BK (q 2 ) contributes. The calculation of the nonlocal amplitudes following the method suggested in [11] proceeds in two stages. First, the amplitudes H (c,u) BP (q 2 ) are splitted in the contributions with different topologies, including c or u quark emission in LO, NLO factorizable corrections, nonfactorizable effects of soft gluon emission, hard-spectator and annihilation contributions. They are calculated one by one at spacelike q 2 < 0 where the light-cone OPE for the correlation function (B.2) is valid. For the hard-gluon NLO and spectator contributions we apply the QCD factorization and for the soft gluon emission the dedicated LCSRs. A detailed account of this calculation can be found in Refs. [12] and [13]. After that, the resulting functions H (c,u) BP (q 2 < 0) are fitted to the hadronic dispersion relations in the q 2 variable where the contributions from the lowest vector mesons V = ρ, ω, φ, J/ψ, ψ(2S) are isolated and the excited states and continuum contributions are modeled, employing the quark-hadron duality. Here we employ as an additional input the experimental data on branching fractions of the nonleptonic B → V P decays determining together with the vector meson decay constants the moduli of the residues in the pole terms of the dispersion relation. The phases of these contributions are included in the set of fit parameters. Since in this paper we are interested only in the large recoil (low q 2 ) region, the integral over hadronic spectral density at q 2 > 4m 2 D with no singularities in the large recoil region is modeled by a polynomial with complex parameters (see Refs. [12,13] for details). Indeed, for our purposes it is not necessary to use a more detailed hadronic representation, like the ansatz suggested in [46] and used in [47], where the broad charmonium resonances located above the open charm threshold are resolved with separate relative phases.
Having fitted the parameters of dispersion relations, we continue them to the positive values of q 2 in the large recoil region, where there is a minor influence of the modeldependent contributions. Finally, we note that in our approach the differences between the B s → K, B → K and B → π nonlocal amplitudes originate from the SU (3) f l -violating differences between the decay constants, parameters of light-meson DAs and nonleptonic amplitudes, as well as from the different spectator-quark flavours, determining the diagram content of these amplitudes.