$\bar{B}_s \to f_0(980)$ form factors and the width effect from light-cone sum rules

TL;DR

This paper calculates the $ar{B}_s o f_0(980)$ form factors using light-cone sum rules, highlighting the significance of higher-twist contributions and exploring the $S$-wave $ar{B}_s o KK$ form factors with a resonance model.

Contribution

It provides a detailed calculation of $ar{B}_s o f_0(980)$ form factors including higher-twist effects and introduces a resonance model for $ar{B}_s o KK$ form factors.

Findings

Higher-twist contributions account for about 50% correction to form factors.

The background effect from $f'_0+f''_0$ resonance is approximately 5%.

The strong coupling $|g_{f_0 KK}|$ is extracted as 1.08 GeV with uncertainties.

Abstract

In this paper we calculate the $\bar{B}_s \to f_0(980)$ form factors from light-cone sum rules with $B$ meson DAs. With adopting the quark-antiquark configuration of light scalar mesons,the high twist two-particle and the three-particle contributions are found to be $\sim 25\%$ individual, and totally they give about $50\%$ correction to certain form factors in the considered energy regions. We further explore the light-cone sum rules approach to study the $S-$wave $\bar{B}_s \to KK$ form factors,the $f_0+f'_0+f''_0$ resonance model is proposed and the result shows that the background effect from $f'_0+f''_0$ accounts $\sim 5\%$. As a by-product, we extract the strong coupling $|g_{f_0 KK }| = 1.08^{+0.05}_{-0.14}$ GeV with taking the $\bar{B}_s \to f_0(980)$ form factors calculated previous under the narrow width approximation.

Tables

Table 1. Table 1: B ¯ s → f 0 ​ ( 980 ) → subscript ¯ 𝐵 𝑠 subscript 𝑓 0 980 \bar{B}_{s}\to f_{0}(980) form factors at q 2 = 0 superscript 𝑞 2 0 q^{2}=0 ( q 2 = 4 ​ m π 2 superscript 𝑞 2 4 superscript subscript 𝑚 𝜋 2 q^{2}=4m_{\pi}^{2} for ElBennichXY ) predicted in different methods.

$\mathrm{Methods}$	${ℱ}_{+}$	${ℱ}_{-}$	${ℱ}_T$
$\mathrm{PQCD}$ LiTK	$0.70$	$\equiv 0$	$0.40$
$\mathrm{CLFD}$ ElBennichXY	$0.80$	$\equiv 0$	$-$
$\mathrm{CQM}$ IssadykovIBA	$0.254$	$-$	$0.285$
$\mathrm{QCDSRs}$ GhahramanyZZ	$0.12$	$-0.07$	$-0.08$
$\mathrm{LCSRs}$ ColangeloBG	$0.37$	$\equiv 0$	$0.228$
$\mathrm{LCSRs}-\mathrm{chiral}$ SunNV	$0.44$	$-0.44$	$0.58$
$\mathrm{LCSRs}$ WangVRA	$0.90$	$0.14$	$0.60$
$\mathrm{this}\mathrm{work}$	$0.52$	$0.04$	$0.21$

Table 2. Table 2: The fitting result for ℱ i subscript ℱ 𝑖 \mathcal{F}_{i} in f 0 + f 0 ′ + f 0 ′′ subscript 𝑓 0 subscript superscript 𝑓 ′ 0 subscript superscript 𝑓 ′′ 0 f_{0}+f^{\prime}_{0}+f^{\prime\prime}_{0} model.

$\left|g_{S\pi \pi }\right|{ℱ}_j$	$\left|g_{S\pi \pi }\right|{ℱ}_{+}$	$\left|g_{S\pi \pi }\right|{ℱ}_0$
${\kappa }_{{ℱ}_j}^{f_0}\left(\mathrm{GeV}\right)$	${0.56}_{-0.08}^{+0.02}$	${0.40}_{-0.05}^{+0.01}$
${\eta }_{{ℱ}_j}^{f_0}\left(\mathrm{GeV}\right)$	$-{0.97}_{+0.13}^{-0.01}$	${0.19}_{-0.19}^{+0.15}$
${\rho }_{{ℱ}_j}^{f_0}\left(\mathrm{GeV}\right)$	$-{12.5}_{-1.42}^{+2.60}$	${0.61}_{+0.57}^{-0.60}$
${\gamma }_{{ℱ}_j}^{f_0^{\prime }}$	${0.49}_{-0.17}^{+0.45}$	${0.48}_{-0.16}^{+0.26}$
${\gamma }_{{ℱ}_j}^{f_0^{\prime \prime }}$	${0.66}_{-0.20}^{+0.80}$	${0.62}_{-0.26}^{+0.58}$
	${ℱ}_{+}^{{\overline{B}}_s\to f_0}\left(0\right)$	${ℱ}_0^{{\overline{B}}_s\to f_0}\left(0\right)$
	$0.52\pm 0.10$	$0.37\pm 0.06$

Equations

$${\Pi }_{\text{2pSRs}}(q)=i\int d^{4}x{e}^{iqx}⟨0|T\{J^s(x),J^s(0)\}|0⟩.$$

$$|f_0(980)⟩=|s\bar{s}⟩\cos \theta +|n\bar{n}⟩\sin \theta ,$$

$$|\sigma (500)⟩=-|s\bar{s}⟩\sin \theta +|n\bar{n}⟩\cos \theta .$$

$$⟨f_0|\bar{u}u|0⟩=\frac{1}{\sqrt{2}}{m}_{f_0}{\bar{f}}_{f_0}^n,⟨f_0|\bar{s}s|0⟩=m_{f_0}{\bar{f}}_{f_0}^s,$$

$$⟨\sigma |\bar{u}u|0⟩=\frac{1}{\sqrt{2}}{m}_{\sigma }{\bar{f}}_{\sigma }^n,⟨\sigma |\bar{s}s|0⟩=m_{\sigma }{\bar{f}}_{\sigma }^s.$$

$$2\text{Im}{\Pi }_{2pSRs}(q)$$

$$\sum _{S=\sigma ,f_0}\frac{m_S^2({\bar{f}}_S^s{)}^2}{(m_S^2-q^2)}=\frac{1}{\pi }{\int }_0^{s_0}ds\frac{\text{Im}{\Pi }_{\text{2pSRs}}^{\text{OPE}}(s)}{s-q^2},$$

$$\displaystyle\sum_{S=\sigma,f_{0}}m_{S}^{2}(\bar{f}_{S}^{s}(\mu_{0}))^{2}e^{-m_{S}^{2}/M^{2}}\Big{(}\frac{\alpha_{s}(\mu_{0})}{\alpha_{s}(M)}\Big{)}^{2/\beta_{1}}$$

$$\displaystyle f(n)=1-e^{-\frac{s_{0}}{M^{2}}}\Big{[}1+\frac{s_{0}}{M^{2}}+\frac{1}{2}\Big{(}\frac{s_{0}}{M^{2}}\Big{)}^{2}+\frac{1}{n!}\Big{(}\frac{s_{0}}{M^{2}}\Big{)}^{n}\big{]}\,,$$

$$I(n)={\int }_{e^{-\frac{s_0}{M^2}}}^{1}dt{\ln }^{n}t\ln (-\ln t),$$

$$\displaystyle h(n)=\frac{\alpha_{s}(M)}{\pi}\Big{(}\frac{17}{3}+2\frac{I(1)}{f(1)}-2\ln\frac{M^{2}}{\mu^{2}}\Big{)}\,.$$

$$\displaystyle\alpha_{s}(\mu)=\frac{\pi}{2\beta_{1}\log(\mu/\Lambda)}\Big{(}1-\frac{\beta_{2}}{\beta_{1}^{2}}\frac{\log(2\log(\mu/\Lambda))}{2\log(\mu/\Lambda)}\Big{)}\,,$$

$${\Lambda }^{(n_f)}=\mathrm{Which}[n_f=3,0.332,n_f=4,0.292,n_f=5,0.210].$$

$$⟨\bar{u}u⟩=(-0.25\mathrm{GeV}{)}^3,⟨\bar{s}s⟩=0.8⟨\bar{u}u⟩,$$

$$⟨g_s\bar{q}\sigma \cdot Gq⟩=-0.8⟨\bar{q}q⟩,$$

$$⟨{\alpha }_s/\pi G_{\mu \nu }^a{G}^{a\mu \nu }⟩=0.012{\mathrm{GeV}}^4,$$

$$⟨\bar{q}{\gamma }_{\mu }{\lambda }^{a}q\bar{q}{\gamma }^{\mu }{\lambda }^{a}q⟩=-\frac{16}{9}⟨\bar{q}q{⟩}^2,$$

$$⟨\bar{q}{\sigma }_{\mu \nu }{\lambda }^{a}q\bar{q}{\sigma }^{\mu \nu }{\lambda }^{a}q⟩=-\frac{3}{4}⟨\bar{q}q{⟩}^2.$$

$$m_s=0.125\mathrm{GeV}.$$

$$\displaystyle m_{q}(\mu)=m_{q}(\mu_{0})\Big{(}\frac{\alpha_{s}(\mu_{0})}{\alpha_{s}(\mu)}\Big{)}^{-1/\beta_{1}}\,,$$

$$\displaystyle\langle\bar{q}q\rangle_{\mu}=\langle\bar{q}q\rangle_{\mu_{0}}\Big{(}\frac{\alpha_{s}(\mu_{0})}{\alpha_{s}(\mu)}\Big{)}^{1/\beta_{1}}\,,$$

$$⟨{\alpha }_s{G}^2{⟩}_{\mu }=⟨{\alpha }_s{G}^2{⟩}_{{\mu }_0},$$

$$\displaystyle\langle g_{s}\bar{q}\sigma\cdot Gq\rangle_{\mu}=\langle g_{s}\bar{q}\sigma\cdot Gq\rangle_{\mu_{0}}\Big{(}\frac{\alpha_{s}(\mu_{0})}{\alpha_{s}(\mu)}\Big{)}^{-1/(6\beta_{1})}\,.$$

$$\sum _{S=\sigma ,f_0}\frac{m_S^4({\bar{f}}_S^s({\mu }_0){)}^2{e}^{-m_S^2/M^2}}{m_S^2({\bar{f}}_S^s({\mu }_0){)}^2{e}^{-m_S^2/M^2}}$$

$$m_{f_0}=(985\pm 122)\mathrm{MeV},$$

$$m_{\sigma }=(439\pm 304)\mathrm{MeV},$$

$${\bar{f}}_{f_0}^s(1\mathrm{GeV})=(358\pm 4)\mathrm{MeV},$$

$${\bar{f}}_{\sigma }^s(1\mathrm{GeV})\sim 0.$$

$${\Pi }_{\nu }(p,q)=i\int d^{4}x{e}^{ip\cdot x}⟨0|T\{J^s(x),J_{\nu }^{\mathcal{I}}(0)\}|{\bar{B}}_s(p+q)⟩,$$

$${\tilde{\Pi }}_{\nu }(p,\tilde{q})=i\int d^{4}x{e}^{ip\cdot x}⟨0|T\{J^s(x),{\tilde{J}}_{\nu }^{\mathcal{I}}(0)\}|{\bar{B}}_{s,v}(p+\tilde{q})⟩.$$

$$s_s(x,0,m_s)=\int \frac{d^{4}p}{(2\pi {)}^4}{e}^{-ipx}{\int }_0^{1}du{G}_{\mu \nu }(ux)$$

$$\displaystyle\cdot\Big{[}\frac{ux^{\mu}\gamma^{\nu}}{p^{2}-m_{s}^{2}}-\frac{(p/+m_{s})\sigma_{\mu\nu}}{2(p^{2}-m_{s}^{2})^{2}}\Big{]}\,,$$

$$⟨0|{\bar{s}}_{\alpha }(x)h_{v\beta }(0)|{\bar{B}}_{s,v}⟩$$

$$\displaystyle+\left(\frac{x_{\rho}\gamma_{\delta}-x_{\delta}\gamma_{\rho}}{v\cdot x}\right)\,\frac{x/}{v\cdot x}\,Z(\omega,\zeta)\,\Big{\}}\gamma_{5}\Big{]}_{\beta\alpha}\,,$$

$$⟨S(p)|J_{\nu }^A(0)|{\bar{B}}_s(p+q)⟩$$

$$=$$

$${\mathcal{F}}_1(q^2)=\frac{{\mathcal{F}}_{+}}{2}(q^2),$$

$$=$$

$$m_{f_0}{\bar{f}}_{f_0}^s{\mathcal{F}}_{+}^{{\bar{B}}_s\to f_0}(q^2)e^{-m_{f_0}^2/M^2}$$

Full text

1 Introduction

Form factor is a fundamental physical quantum in effective field theory (EFT), with including both the long distance (LD) and short distance (SD) physics. In order to extract the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements, reliably, in the semileptonic $B$ decay processes, the precise calculations of the relevant form factors are inevitable. The heavy-to-light form factors are calculated by different approaches, in which the lattice QCD (LQCD) gives the reliable simulations in the low recoiled regions HorganPVA ; LatticeTIA , while in the large and full recoiled regions, the QCD-based analytical approaches like the light-cone sum rules (LCSRs) BraunIJ ; BraunKW ; ColangeloDP ; BallYE ; BallRG ; StraubICA ; KhodjamirianST ; DuplancicIX ; WangVGV ; ChengSMJ ; Descotes-GenonBUD ; RusovCHR and the perturbative QCD (PQCD) LiNU ; LiIU ; KurimotoZJ ; LuNY ; ChenGV ; LiNK ; ChengFWA are applied.

There are many successful calculations for the heavy-to-light transition form factors with the final state being a pseudoscalar ($P$) or a vector ($V$) meson. While to our knowledge, the transition form factors with light scalar ($S$) meson final state are not understood well so far due to the unclear underlying structure and the large width effect. It has been suggested that the scalar mesons with masses below or near $1$ GeV (the isoscalar $\sigma/f_{0}(500)$ and $f_{0}(980)$, the isodoublet $\kappa$, and the isovector $a_{0}$) form a $SU(3)$ flavor nonet, and the scalar mesons with masses around $1.5$ GeV ($f_{0}(1370)$, $a_{0}(1450)$, $K_{0}^{\ast}(1430)$ and $f_{0}(1500)$) form the other one. The combined analysis based on orbital angular momentum JaffeIG ; ChengNB and data CloseZU ; AchasovHM ; AchasovFH implies that the heavier nonet states favour the quark-antiquark configuration replenished with some possible gluon content. From the spectral analysis, there is not a general agreement on the underlying assignments of the scalar mesons in the lighter nonet, like $f_{0}(980)$. Pictures like tetra-quark JaffeIG ; Agaev:2017cfz ; Agaev:2018sco , gluonball WeinsteinGC , hybrid state WeinsteinGD and molecule state WeinsteinGU are all discussed, in which the tetra-quark assignment is more favorite nowadays. The case is different in the $B\to f_{0}(980)$ decays when $f_{0}(980)$ is energetic and the process happens with large recoiling, where the conventional quark-antiquark assignment is the favorite one since the possibility to form a tetra-quark state is power suppressed with comparing to the state of quark pair ChengNB . So in this work, with the main purpose to calculate the heavy-quark transition form factors, we would take the usual quark-antiquark nature of $f_{0}(980)$, and postpone the tetra-quark study somewhere else.

Some attempts are carried out to calculate the $B\to S$ transition form factors from LCSRs with scalar mesons distribution amplitudes (DAs) GhahramanyZZ ; SunNV ; ColangeloBG ; WangVRA . We comment that these work considered only the scalar mesons in the heavier nonet, and their accuracy is debatable since some important informations of the input DAs, such as the standard conformal partial expansion and the width effect, are still missing. In this paper, with taking the $\bar{s}s$ configuration of $f_{0}(980)$, we suggest to study the $\bar{B}_{s}\to f_{0}(980)$ form factor from the alternative LCSRs with $B$ meson DAs. Although the width of $f_{0}(980)$ is suppressed by the phase space111In fact the width of $f_{0}$ is smaller than it of $\rho$ meson. It should be stressed that the width effect in $B\to\rho$ transition is usually neglected because in the experimental analysis the $\rho$ meson is identified by the P-wave $\pi\pi$ signal when the dipion invariant mass locates in the $\rho-$pole region BehrensVV ; AdamPV . This makes the narrow-width treatment for $B\to\rho$ form factor from LCSRs being consistent with the experimental measurement., we would like to access the width effect by applying the approach proposed to calculate $B\to\pi\pi,\,K\pi$ form factors ChengSMJ ; Descotes-GenonBUD , with substituting the isovector $\pi\pi$ state by the scalar isoscalar $KK$ state.

The rest of this paper is organised as follows. In section 2 we revisit and update the mass and the decay constant of $f_{0}(980)$ in the two-point QCD sum rules. Section 3 and section 4 are the main parts of this paper, where we present the LCSRs calculation for $\bar{B}_{s}\to f_{0}(980)$ form factors and generalizes it to study the $S-$wave $\bar{B}_{s}\to KK$ form factors, respectively. We summary in section 5. The coefficients in three-particle corrections from $B$ meson LCDAs are complemented in Appendix B.

2 $\sigma$ and $f_{0}(980)$ in the QCD sum rules revisited

QCD sum rules ShifmanBX is a powerful tool to study hadron spectrum. For the meson with $I^{G}(J^{PC})=0^{+}(0^{++})$, the scalar isoscalar currents include both $J^{n}=\bar{n}n=\frac{1}{\sqrt{2}}(\bar{u}u+\bar{d}d)$ and $J^{s}=\bar{s}s$. We start with the two-point correlation function

[TABLE]

Although the data of $D_{s}^{+}\to f_{0}\pi^{+}$ indicates $f_{0}(980)$ may be dominated by the $\bar{s}s$ component, much more measurements GokalpNY ; AnisovichZP ; ChengAI ; KaminskiQG ; MennessierXG ; ZhangKM ; LiSW ; ZhangQVQ ; LiuYMI ; FeldmannSH , especially the comparable branching ratios between $B\to f_{0}(980)\to\pi\pi$ and $B\to f_{0}(980)\to KK$ TanabashiOCA , support a mixing between $f_{0}$ and $\sigma$:

[TABLE]

The mixing implies that $f_{0}$ and $\sigma$ should be treated separately, and in the basis of flavour, two decay constants are needed to describe each of them,

[TABLE]

The neutral scalar meson can not be produced via the vector current because $f_{S}$ is vanished in the $SU(3)$/isospin limit with the charge conjugation invariance and the conservation of vector current.

The basic idea of QCD sum rules is to calculate independently for the correlation function in twofold ways: the QCD calculation at quark-gluon level in the Euclidean momenta space, and the summing of intermediate states from the view of hadron. The QCD calculation in the negative half plane of $q^{2}$ is guaranteed by the operator-product-expansion (OPE) technology, and the correlation function is then written in terms of various vacuum condensates. On the other hand, the average distance between two coordinate points ([math] and $x$ in Eq.1) grows when $q^{2}$ shifting from large negative to positive values, and then the LD quark-gluon interaction forms the hadrons ColangeloDP . In this way, the correlation function can be expressed as the sum of contributions from all possible intermediate states in the positive half-plane, with possible subtractions. The accuracy of LCSRs approach is mainly depended on how to match the QCD calculation to the hadron spectral analysis, in more word is how to take the quark-hadron duality in the dispersion relation to eliminate the contributions from excited and continuum states.

By inserting a complete set of intermediate states $|n\rangle$, the unitarity relation ($q^{2}>0$) of the correlation function reads as,

[TABLE]

where $d\tau_{n}$ denotes the phase space of each state $|n\rangle$, like $\sigma\,,f_{0}(980)$ and their excited states. We are now interesting in $f_{0}(980)$ which enters in the two-point sum rules as the first excited state with the ground state $\sigma$, so we include both them in the hadron inserting and take the threshold $s_{0}$ to truncate the higher excited states. After applying the dispersion relation and employing the quark-hadron duality, the matching between quark amplitude and the hadron spectral analysis is taken as

[TABLE]

where $\Pi^{\text{OPE}}_{\text{2pSRs}}(s)$ is the OPE result for the correlation function. In order to improve the convergence of OPE calculation and suppress the contributions from high excited states and continuum spectrums, we apply the Borel transformation to both sides of Eq.5. The result is quoted GovaertsUA ; BecchiVZ as follow, with $\alpha_{s}$ to one-loop order and the vacuum condensate terms up to dimension six,

[TABLE]

The functions defined in the perturbative terms are

[TABLE]

We use the two-loop expression for strong coupling TanabashiOCA ,

[TABLE]

with the evolution kernels $\beta_{1}=(33-2n_{f})/12$ and $\beta_{2}=(153-19n_{f})/24$. The hadronic scale $\Lambda$ is set to reproduce $\alpha_{s}(m_{Z})=0.118$ and $\alpha_{s}(1\,\mathrm{GeV})=0.474$, and written in terms of step function

[TABLE]

The input values for the nonperturbative vacuum condensates at default scale $1\,\mathrm{GeV}$ are taken as LeutwylerQG ; IoffeEE :

[TABLE]

The light quark masses are estimated as the ”current-quark” masses in the $\overline{\mathrm{MS}}$ ($\mu=1\,\mathrm{GeV}$) scheme TanabashiOCA ,

[TABLE]

We also consider the running of parameters with one-loop accuracy ChernyakAS ; LepageZB ; VermaserenFQ ,

[TABLE]

Differentiating both sides of Eq.6 by the Borel mass we obtain an auxiliary sum rules,

[TABLE]

which is further used, together with the renormalization-improved sum rules in Eq.6, to fit the masses and decay constants of $\sigma$ and $f_{0}$. The terms proportional to quark mass on the right hand side are neglected in Eq.13.

The Borel mass is fixed by the rule of thumb that the contribution from high dimension condensate terms is no larger than twenty percents in the truncated OPE, and simultanously the contribution from excited and continuum states is smaller than thirty percents when summing up the hadrons. The threshold value $s_{0}$ is usually close to the outset of the first higher excited state with the same quantum number, then a certain vicinity can be expected, we determine it with considering the maximal stability of physical quantities once the Borel mass has been set down. Within the interval $M^{2}=1.0\pm 0.1$ GeV2 at the fixed threshold value $s_{0}=2.0\pm 0.2$ GeV2 which is slightly larger than the one used in Agaev:2017cfz ; Agaev:2018sco because we are discussing the $s\bar{s}$ current, we do the combined quadratic fit to both sides of Eqs.(6,13), and obtain

[TABLE]

The $s-$flavor decay constant of $f_{0}$ agrees with the prediction $\bar{f}^{s}_{f_{0}}=(370\pm 20)\,\mathrm{MeV}$ obtained under the assumption that $f_{0}(980)$ and $f_{0}(1500)$ are the lowest scalar states with $\bar{s}s$ assignment ChengNB . The nearly zero $s-$flavor decay constant of $\sigma$ indicates that $f_{0}(980)$ is the lowest state in the channel with scalar isoscalar current $J^{s}$, standing by which we reevaluate the sum rules in Eqs.(6, 13) by considering only the $f_{0}$ state, and obtain the same result for $m_{f_{0}}$ and $\bar{f}^{s}_{f_{0}}$ as listed in Eq.14.

3 $\bar{B}_{s}\to f_{0}$ form factors from the LCSRs

The approach of LCSRs with B meson DAs was proposed to calculate the $B\to P,V$ form factors KhodjamirianST , in this section we implement it to calculate the $\bar{B}_{s}\to f_{0}$ form factors.

As we demonstrated in the last section that the scalar isoscalar current $\bar{s}s$ coupling to $\sigma$ is nearly zero, so it would be reasonable to consider $f_{0}(980)$ as the lowest scalar state in $B_{s}$ decays. This consideration is also supported by the fact that no evidence of $\sigma$ is observed in the $B$ decays so far AaijEMV . Let’s consider another correlation function

[TABLE]

with the weak current $J_{\nu}^{\mathcal{I}}=\bar{s}\Gamma_{\nu}^{\mathcal{I}}b$. The indicator $\mathcal{I}=A,T$ correspond to the gamma matrices $\Gamma_{\nu}^{A}=\gamma_{\nu}\gamma_{5}$ and $\Gamma_{\nu}^{T}=\sigma_{\nu\mu}\gamma_{5}q^{\mu}$, respectively222We use the convention $\sigma_{\mu\nu}=\frac{i}{2}(\gamma_{\mu}\gamma_{\nu}-\gamma_{\nu}\gamma_{\mu})$.. The heavy-to-light current is reduced to the light quark current after transiting to the heavy quark effective theory (HQET), and the correlation function is modified to

[TABLE]

We use the notations $p$ to denote the momentum carried by the scalar isoscalar current $J^{s}$, and $\tilde{q}=q-m_{b}v$ for the effective current $\tilde{J}^{\mathcal{I}}_{\nu}=\bar{s}\Gamma^{\mathcal{I}}_{\nu}h_{v}$. In the rest frame the effective $b$-quark field is defined by $h_{v}(x)=e^{im_{b}vx}b(x)$ with the unit vector $v=(1,0,0,0)$, and $|\bar{B}_{s}(p+q)\rangle=|\bar{B}_{s,v}(p+\tilde{q})\rangle$ holds up to the $\mathcal{O}(1/m_{b})$ accuracy. We would discuss in the intervals $|p^{2}|,|\tilde{q}^{2}|\gg\Lambda_{QCD}^{2},(m_{B_{s}}-m_{b})^{2}$ where the correlation function does not fluctuate violently and the OPE calculation works well.

In the correlation function, the same flavour quark fields with small displacement can be contracted in the form of quark propagator,

[TABLE]

in which the first term is the freedom part in the QCD limit, and the second term respects the soft one-gluon correction. Two-particle and three-particle DAs of $B_{s}$ meson are defined as GeyerFB ,

[TABLE]

respectively. Two variables $\omega$ and $\zeta$ are introduced to represent the plus components of light quark and the gluon momentum, respectively. The path-ordered gauge factor is always underlied in the matrix element sandwiched between meson state and vacuum. Recently, the renormalization group equations (RGE) are resolved in the $N_{C}$ limit for three-particle DAs BraunLIQ , and the models are suggested for the higher-twist $B$ meson DAs BenekeWJP , following which the power suppressed correction are supplemented to $B\to P,V,\gamma$ form factors Descotes-GenonBUD ; BenekeWJP ; GubernariWYI ; LuCFC ; ShenABS . The Lorentz definition in Eqs.(3,19) should not be confused with the definite twist definition of LCDAs, we collect their relations, as well as the general model for the later one in appendix A.

The definition of $\bar{B}_{s}\to S$ transition form factors is quoted as WirbelJ ; BauerBM

[TABLE]

The following relations are suggested for the form factors associated with axial-vector current in Eq.3,

[TABLE]

We calculate the correlation function in Eq.16 under the narrow width approximation, obtain the LCSRs result for $\bar{B}_{s}\to f_{0}(980)$ form factors,

[TABLE]

The contributions proportional to light quark mass $m_{s}$ are not shown explicitly above. The dimensionless variable $\sigma\equiv\omega/m_{B_{s}}$ is the longitudinal momentum fraction of the light quark inside $\bar{B}_{s}$ meson, and the virtuality of internal quark is $s_{q}=m_{B_{s}}^{2}\sigma-(q^{2}\sigma-m_{s}^{2})/\bar{\sigma}$. To process the calculation, we have defined an auxiliary distribution

[TABLE]

with the boundary conditions $\overline{\phi}_{\pm}(0)=\overline{\phi}_{\pm}(\infty)=0$. The derivation is $g^{\prime}_{+}(\sigma m_{B_{s}})=\left(d/d\sigma\right)g_{+}(\sigma m_{B_{s}})$. We give several comments in orders:

• (i)

With taking into account the quark mass effect, our results shown in Eqs.(3,3) consist with the calculations of $B\to D_{0}^{\ast}$ form factors ShenMM .

• (ii)

$\mathcal{F}_{T,p}(q^{2})$ and $\mathcal{F}_{T,q}(q^{2})$ are obtained by matching the coefficients associated with different Lorentz structures (say, $p_{\nu}$ and $q_{\nu}$, respectively) with the matrix element sandwiched by the tenser current. In principle, they should be equal to each other when considering only the two-particle DAs333Because the coefficients of three-particle correlation in $\Delta\mathcal{F}_{T,p}$ may have the $1/q^{2}$ factor, we would take $\Delta\mathcal{F}_{T,q}$ for this part contribution in numerical analysis., so we use the unified notation $\mathcal{F}_{T}(q^{2})$.

• (iii)

In the heavy quark limit $\sigma\to 0$, we reproduce the relations $\mathcal{F}_{+}(0)=2m_{B_{s}}/(m_{B_{s}}+m_{S})\,\mathcal{F}_{T}(0)$ and $\mathcal{F}_{-}(0)=0$ at the full recoiled point.

Multiplying both sides of Eq.3 by $q^{\nu}$, we derive the matrix elements deduced by the pseudo-scalar current $J^{P}=im_{b}\bar{s}\gamma_{5}b$,

[TABLE]

which suggests another sum rules,

[TABLE]

We remark that Eq.28 is established on the heavy quark limit, so the new LCSRs in Eq.29 can be used to estimate how well does the heavy quark expansion work by comparing it with the original LCSRs in Eqs.(3,3).

The contributions from three-particle DAs of $B$ meson are arranged in an universal form

[TABLE]

where the dimensionless variable $\eta=(\bar{\sigma}^{2}m_{B_{s}}^{2})/(\bar{\sigma}^{2}m_{B_{s}}^{2}-q^{2}+m^{2})$ can be understood as the ratio between the minimal virtuality of the $b$ quark field and the maximal virtuality carried by the internal light quark. The integral over the three-particle DAs is written as

[TABLE]

here another two auxiliary distributions are introduced for $\overline{X}_{A}\,,\overline{Y}_{A}\,,\overline{W}\,,\overline{\widetilde{X}}_{A}\,,\overline{\widetilde{Y}}_{A}$ and $\overset{=}{W}\,,\overset{=}{Z}$,

[TABLE]

The lower indicator $N=1,2,3$ stands for the power of Borel mass $M^{-2(N-1)}$ premultiplied with the integrals, the coefficients $C_{N,i}$ associated to each three-particle DA are presented in appendix B.

We take $m_{b}(m_{b})=4.2$ GeV for $b$ quark mass in the typical $\overline{\mathrm{MS}}$ scheme, and use $f_{B_{s}}=0.242$ GeV obtained from lattice QCD BazavovAA and two-point QCD sum rules GelhausenWIA . The inverse moment of $\bar{B}_{s}$ DAs is chosen in the interval $\lambda_{\bar{B}_{s}}=450\pm 50$ MeV, a little bit smaller then the conventional value $500\pm 50$ MeV, by considering the possible next-to-leading-order radiative correction effects444Where $\lambda_{B}=358^{+38}_{-30}(343^{+22}_{-20})$ MeV is obtained by comparing the $B\to\pi(\rho)$ form factor from LCSRs with pion KhodjamirianUB (rho StraubICA ) and $B$ meson DAs. WangVGV ; GaoLTA . The same ballpark of the LCSRs parameters, say, $M^{2}=1.0\pm 0.1$ GeV2 and $s_{0}=2.0\pm 0.2$ GeV2, are used here as in the two-point sum rules in the last section, with which the OPE convergence is automatically manifested by the relative small three-partical DAs contribution $\mathcal{F}^{3p}(Q^{2})/\mathcal{F}^{2p}(Q^{2})\lesssim 30\%$ for $Q^{2}\in[0,5]$ GeV2. We plot in figure 1 for the $\bar{B}_{s}\to f_{0}$ form factors obtained under the narrow width approximation, where the lightgray shadows reveal the total uncertainty came from the inverse moment $\lambda_{\bar{B}_{s}}$ and the LCSRs parameters. The high twist $B$ meson DAs with two-particle configuration ($g_{+}$) give about $25\%$ decrease to $\mathcal{F}_{+}$ and $\mathcal{F}_{T}$, about $20\%$ increase to $\mathcal{F}_{-}$, the tiny change of $\mathcal{F}_{0}$ from $g_{+}$ can be understood by the interplay between the negative correction for $\mathcal{F}_{+}$ and the positive correction for $\mathcal{F}_{-}$. The three-particle $\bar{B}_{s}$ DAs bring another $25\%$ decrease and increase to $\mathcal{F}_{T}$ and $\mathcal{F}_{-}$, respectively, while its corrections to $\mathcal{F}_{+}$ and $\mathcal{F}_{0}$ are tiny, indicating that the heavy quark limit is broken with considering the three-particle DAs correction. We compare our result for $\bar{B}_{s}\to f_{0}$ form factors with other methods in table 1.

4 The width effect and the $\bar{B}_{s}\to KK$ form factors

To investigate the width effect of intermediate states, let’s look back to the dispersion relation of correlation function in Eq.15,

[TABLE]

The imaginary part in the numerator, corresponding to the physical regions ($q^{2}>0$), can be obtained by interpolating a complete set of intermediate states between two local currents in the correlation function,

[TABLE]

where $p$ and $d\tau_{n}$ denote the momentum and phase space volume of each state $|S_{n}\rangle$, respectively. In the narrow width approximation, $|S_{n}\rangle$ are single mesons and the dispersion relation reduces to

[TABLE]

In Eq.35, only the contribution from ground state is singled out while the rest parts are retained in the integral, this is exactly what we did in the last section.

A straightforward way to consider the width effect is to substitute the interpolation of single mesons by stable multi-meson states, such as the $KK,\,\pi\pi,\,\eta\eta$ and their continuum states ChengSMJ ,

[TABLE]

Here we only write out explicitly the term contributed from $KK$ state, the ellipsis denotes the contributions from $\pi\pi$, $4\pi$, $\eta\eta$, $\eta\eta^{\prime}$ and their excited states. Because it is much more harder to generate a $\bar{s}s$ pair than a $\bar{n}n$ pair from vacuum, the contribution from $\eta\eta$, $\eta\eta^{\prime}$ channels can be expected to be small. The contributions from $\pi\pi$ and $4\pi$ channels are also small since they are produced via $KK\to\pi\pi,4\pi$ rescattering, which exceeds the scope we are discussing here. In the follow calculation we consider $KK$ state as the ground state which gives dominant contribution to the correlation function, while the contributions from high excited multi-meson states are suppressed by the Borel exponent $e^{-s/M^{2}}$ (with the Borel mass $M^{2}\lesssim 1.1$ GeV2).

4.1 Formalism

The scalar isoscalar kaon form factor is defined as DonoghueXH

[TABLE]

There is no experiment measurement for $\Gamma_{K}^{s}$ as yet555For the scalar isoscalar pion form factor $\Gamma_{\pi}$ defined with $\bar{n}n$ source current, the $N$-subtracted omnés representation gives a good description for the data up to $k^{2}=1.5^{2}$ GeV2 with considering the generalized Watson theorem of the $\pi\pi\to\pi\pi$ phase YndurainCM ; PelaezVS ; KaminskiQE ; GarciaMartinCN ; KangJAA ; YuanFOA . , so our calculation relies on the theoretical input. In the chiral perturbative theory (CHPT), $\varGamma_{K}^{s}$ has been derived to next-to-leading-order MeissnerBC with supplementing by the unitarity constraint LahdeWR , however, the unitarized CHPT works only at low energies, say, $s<1.1$ GeV2 DoringWKA ; MeissnerHYA . In order to include the high energy behaviour, one is forced to employ a model Hanhart:2012wi ; RopertzSTK and/or to adopt the perturbative QCD approximation ChengKHI .

From the view of hadron, $\varGamma_{K}^{s}$ is relevant to the $T$ matrix elements of $\pi\pi\to KK$ and $KK\to KK$ scatterings Hanhart:2012wi ; RopertzSTK via

[TABLE]

where $M_{i}$, in principle, is an analytic term describing the transition from the scalar isoscalar source to the channel $i$ ($i=1$ represents for $\pi\pi$ and $i=2$ for $KK$), and $G_{ii}$ is the free propagation of the particles in channel $i$. The $\pi\pi$ channel does not contribute to the $\bar{s}\Gamma_{\nu}^{\mathcal{I}}b$ transition at Born level, which is one of the reasons we can obtain the approximate dispersion relation in Eq.36, with retaining only the $KK$ channel. To maintain the self-consistency of the intrepolation, we take only the $KK$ channel in Eq.38 too, and the transition from source current $J^{s}$ to $KK$ channel happens as $G_{KK}M_{K}=1$ with subsequently the free propagating . In this way, a simple relation is obtained as

[TABLE]

with the normalisation $M_{K}(0)=m_{K}^{2}-m_{\pi}^{2}/2$.

We can read from Eqs.(38,39) that the phase of $s$ flavoured kaon form factor is only determined by the $KK\to KK$ scattering amplitude, whose expression with definite isospin ($I$) and partial wave ($J$) is parameterised as

[TABLE]

In the above equation, $\beta_{K}(k^{2})=\sqrt{1-4m_{K}^{2}/k^{2}}$ is the phase space of $KK$ system, $\eta^{I}_{J}$ and $\delta^{I}_{J}$ are the inelasticity and phase shift, respectively, and $\phi^{I}_{J}$ is the phase. We would use the result of $T^{0}_{0}(k^{2})$ obtained from the amplitude analysis DaiZTA ; DaiUAO as input for the scalar form factor666We thanks Ling-yun Dai for sharing us with the original result of their global fit analysis., which, as they claimed, can be extrapolated to a high energy $\sim 5$ GeV2. We show in figure 2 for the result in the energy regions $k^{2}\in[4m_{\pi}^{2},2.0]$ GeV2, these curves consist with the result obtained from CHPT in the low energy regions $k^{2}<1.1$ GeV2 DoringWKA , and also consist with the fully $K$-matrix description in the high energies RopertzSTK . The amplitude analysis result we adopted here considered all the measured data, the $\pi\pi-K\bar{K}$ final state interaction, the mass difference between the charged and neutral kaon, and the low energy Roy equation BuettikerPP . In the amplitude analysis demonstrated by coupled-channel treatment and combined fitting, the source current with $s\bar{s}$ configuration is overwhelming coupled to $KK$ channel through $f_{0}$, while the coupling to $\pi\pi$ is tiny.

The remaining matrix element in Eq.36 are defined in terms of $\bar{B}_{s}\to KK$ transition form factors FallerDWA , for the axial-vector current $j_{\nu}^{A}$ we have

[TABLE]

The kinematics of $\bar{B}_{s}$ decay to $KK$ state are described by three independent variables: $k^{2}$, $q^{2}$ and $\theta_{K}$, denoting the invariant mass of $KK$ system, the squared momentum transfer in the weak decay and the angle between the 3 momentum of $K^{-}(k_{2})$ and $\bar{B}_{s}$ meson in the $K^{+}K^{-}$ rest frame, respectively. The dot products are

[TABLE]

with the kinematic $\mathrm{K\ddot{a}ll\acute{e}n}$ function $\lambda_{B}\equiv\lambda(m_{B_{s}}^{2},q^{2},k^{2})=m_{B_{s}}^{4}+k^{4}+q^{4}-2(m_{B_{s}}^{2}k^{2}+m_{B_{s}}^{2}q^{2}+k^{2}q^{2})$. By the way, the matrix element in Eq.41 can also be defined by the helicity amplitudes,

[TABLE]

The helicity definition provides a possibility to study the contributions from different partial waves, because $H_{\lambda}$, with $\lambda=t,0,+,-$, can be expanded in terms of the associated Legendre polynomials. To study the partial waves contributions within the convenient definition in orthogonal Lorentz structures, we translate the partial wave expansion from the helicity amplitudes $H_{\lambda}$ to the form factors $F_{i}$,

[TABLE]

Considering the decomposition of isoscalar $KK$ state,

[TABLE]

substituting Eq.41 and Eq.44 into Eq.36, we obtain the $S-$wave contribution to the imaginary part with interpolating scalar isoscalar $KK$ state,

[TABLE]

In fact, the phase space $d\tau_{2K}$ plays as a $S-$wave projector for the timelike-helicity form factors $F_{t,0}$, which means that only the $S-$wave component $F_{t,0}^{(l=0)}$ survives after integrating over the angle $\theta_{K}$. For the form factor $F_{\parallel}$, the role of $S-$wave projector vanishes and the contribution starts from $D$-wave component ($l=2n,n=1,2,3\cdots$), which part is expected tiny in the $\bar{B}_{s}\to KK$ transtion and would not be discussed in this paper.

We take the global duality to eliminate the contributions beyond $KK$ state with the threshold $s_{0}^{2K}$,

[TABLE]

and arrive at the LCSRs result for $S-$wave $\bar{B}_{s}\to KK$ transition,

[TABLE]

[TABLE]

Multiplying both sides of Eq.46 by $q^{\nu}$, we obtain another independent LCSRs for the timelike-helicity form factor $F_{t}^{(l=0)}(q^{2},s)$,

[TABLE]

4.2 Models

Eqs.(48,49,50) are the main results in this section. Due to the convoluted integral, we can not solve out the form factors $F_{t/0}^{(l=0)}$ in terms of the $B$ meson DAs. To propel the calculation, one way we can try is to introduce the parameterisation of $S-$wave $\bar{B}_{s}\to KK$ form factors, and the first candidate coming into our mind is the single resonance ($f_{0}$) model in the generalized Breit-Wigner formula777We drop the $\sigma$ with the same reasons as described in section 3, which, phenomenologically, is further supported by the fact that no any signal is found for $KK$ coupling to $\sigma$ TanabashiOCA ..

[TABLE]

The strong coupling $g_{f_{0}\pi\pi}$ is normalized as

[TABLE]

An underlying condition implied in Eqs.(48,49,50) is the reality of expressions on the left hand side, and we take the more strict local reality at each point of invariant mass,

[TABLE]

To fulfil this requirement, a strong phase $\phi_{f_{0}}$ is introduced to compensate the phase difference between the kaon form factor and the modeled $\bar{B}_{s}\to KK$ form factors. Generally speaking, $\phi_{f_{0}}$ should depend on both the two variables $s$ and $q^{2}$, while in the single $f_{0}$ model the $q^{2}-$dependence disappears,

[TABLE]

The simple model in Eqs.(51-53) is inspired by the physics that the sum rules obtained for the $\bar{B}_{s}\to KK$ form factors, in the narrow width approximation, should reproduce the sum rules for the form factors of $\bar{B}_{s}\to f_{0}$ transition. To check this, let’s consider the energy-dependent width of $f_{0}$ with including the loop effects of two kaons coupling,

[TABLE]

The width of $f_{0}$ is usually parameterized under the Flatt$\acute{e}$ model FlatteRZ ; FlatteXV with considering the location of $f_{0}$ in the invariant mass, say, below or above the threshold. While in the $\bar{B}_{s}$ decays, the case is different because the invariant mass is alway above the threshold888 Another reason for us not using the Flatt$\acute{e}$ model is that it gives a smaller value of $I^{f_{0}}=4.44^{+0.81}_{-0.98}$, which is close to $I^{f^{\prime}_{0}}$ and damages the contribution hierarchy from different resonances, as we would see in Eqs.(68,69)., then we can take the conventional form of the width as described in Eq.57, In the single $f_{0}$ model, the $s$ flavoured kaon form factor is written as

[TABLE]

Substituting Eqs.(51,58) into Eq.48, the left hand side of Eq.48 becomes

[TABLE]

Similarly, we can reproduce the LCSRs for form factors $\mathcal{F}^{\bar{B}_{s}\to f_{0}}_{-}$ and $\mathcal{F}^{\bar{B}_{s}\to f_{0}}_{0}$ with taking into account the resonance models in Eq.52 and Eq.53, respectively. What’s more, the relation defined in Eq.23 also holds in the resonance models.

4.3 Numerics

We employ the $z-$series expansion for heavy-to-light transition form factors BourrelyZA , with $j=+,-,0$,

[TABLE]

$\mathcal{F}^{\bar{B}_{s}\to f_{0}}_{j}(0)$ is the value at the full recoiled energy, the parameters $b_{\mathcal{F}_{j}},c_{\mathcal{F}_{i}}$ indicate the coefficients associated with the $\zeta-$functions,

[TABLE]

with the definitions $t_{\pm}\equiv(m_{B_{s}}\pm m_{f_{0}})^{2}$ and $t_{0}\equiv t_{+}(1-\sqrt{1-t_{-}/t_{0}})$.

The $S-$wave $\bar{B}_{s}\to KK$ form factors, under the single $f_{0}$ model, is rearranged in a general formula as

[TABLE]

where for the sake of brevity we introduce the following notations:

[TABLE]

The integral coefficient on the left hand side reads as

[TABLE]

$I^{\mathrm{OPE}}_{j}$ represents the OPE calculations on the right hand side of Eqs.(48,49,50). There is no physical requirement that the threshold value $s_{0}^{2K}$ should be equal to $s_{0}$, we fixed it in an independent way by considering the correlation function in Eq.1 with $KK$ interpolating. The 2pSRs is then written in terms of the scalar isoscalar kaon form factor,

[TABLE]

We then determine the value $s_{0}^{2K}=2.0$ GeV2, closing to it taken in the case of single meson interpolating.

Besides the bound state $f_{0}$, it is nature to question what’s the roles of the excited states $f^{\prime}_{0},f^{\prime\prime}_{0}$ in the $\bar{B}_{s}\to KK$ transition999 Hereafter we take $f_{0}$, $f^{\prime}_{0}$ and $f^{\prime\prime}_{0}$ to denote the meson states $f_{0}(980)$, $f_{0}(1500)$ and $f_{0}(1710)$, respectively. We do not consider $f_{0}(1370)$ as a separated resonance due to the weak coupling of $f_{0}(1370)$ to $KK$ state.. To include these effects, we suggest the $f_{0}+f^{\prime}_{0}+f^{\prime\prime}_{0}$ model by appending $f^{\prime}_{0}$ and $f^{\prime\prime}_{0}$ states to Eq.(53),

[TABLE]

We tacitly assume that the strong phase $\phi_{S}$ associated to each intermediate state is only dependent on the invariant mass of dikaon state, which means we do not consider the interaction effect between different resonances when introducing the strong phases to satisfy Eq.55. By this way, the $\bar{B}_{s}\to S$ form factors for different intermediate resonances are linear to each other: $\mathcal{F}_{j}^{f^{\prime}_{0}}(q^{2})=\mathcal{\gamma}^{f^{\prime}_{0}}_{\mathcal{F}_{j}}\mathcal{F}_{j}^{f_{0}}(q^{2})$ and $\mathcal{F}_{j}^{f^{\prime\prime}_{0}}(q^{2})=\mathcal{\gamma}^{f^{\prime\prime}_{0}}_{\mathcal{F}_{j}}\mathcal{F}_{j}^{f_{0}}(q^{2})$. The dimensionless parameters $\mathcal{\gamma}^{f^{\prime}_{0}}_{\mathcal{F}_{j}}$ and $\mathcal{\gamma}^{f^{\prime\prime}_{0}}_{\mathcal{F}_{j}}$ indicate the relative size of $\bar{B}_{s}\to f^{\prime}_{0}$ and $f^{\prime\prime}_{0}$ form factors comparing to the $\bar{B}_{s}\to f_{0}$ form factors, respectively. The LCSRs in Eq.63, in the case of three resonance model, is modified to

[TABLE]

We quote the values of the integral coefficients (in unit of $10^{-2}$)

[TABLE]

In these integrals with taking the bound limit at $4m_{K}^{2}$, only the right half of the peaking region in figure 2 is taken into account, so the relative sizes of integral coefficients $I^{f^{\prime}_{0}},I^{f^{\prime\prime}_{0}}$ to $I^{f_{0}}$ ($I^{f^{\prime}_{0}}/I^{f_{0}},I^{f^{\prime\prime}_{0}}/I^{f_{0}}$) can be expected to be double of that in the $B\to\pi\pi$ case ChengSMJ .

The fitting result in the $f_{0}+f^{\prime}_{0}+f^{\prime\prime}_{0}$ model are presented in table 2, where the errors come from the sum rules parameters. For the total widths of the intermediate resonances we choose $\Gamma_{f_{0}}^{\mathrm{tot}}=0.055$ GeV, $\Gamma^{\mathrm{tot}}_{f^{\prime}_{0}}=0.112$ GeV and $\Gamma^{\mathrm{tot}}_{f^{\prime\prime}_{0}}=0.123$ GeV TanabashiOCA . We note that varying the width of $f_{0}$ in $[0.01,0.1]$ GeV brings another uncertainty to $I^{f_{0}}$ by $|^{+2.55}_{-1.40}$, while the width effects of $f^{\prime}$ ($f^{\prime\prime}$) to $I^{f^{\prime}_{0}}$ ($I^{f^{\prime\prime}_{0}}$) is negligible since their widths are much smaller. In the fit we also use the condition at $q^{2}=0$

[TABLE]

and also the hierarchy anstz of different resonances in the left hand side of Eq.68, say, $0<\mathcal{\gamma}^{f^{\prime}}_{\mathcal{F}_{j}},\mathcal{\gamma}^{f^{\prime\prime}}_{\mathcal{F}_{j}}<1$. At the bottom of Tab.2, for the comparison we supplement the $\bar{B}_{s}\to f_{0}$ form factors calculated in section 3 under the narrow width approximation. In principle, the strong coupling $|g_{f_{0}KK}|$ can be extracted out in case we have the reliable prediction for $\bar{B}_{s}\to f_{0}$ form factors, $\mathcal{F}_{+}$ and $\mathcal{F}_{0}$. With the result obtained in the narrow width approximation as we demonstrated in the last section, we can estimate $|g_{f_{0}KK}|=1.08^{+0.05}_{-0.14}$ GeV, but keep in mind that the width/non-resonant effect is sizeable and we reserve another $\sim 50\%$ uncertainty.

The contribution from each resonance to the OPE result (Eq.68) is listplotted in Fig.3, from which the expected leading role of $f_{0}$ is confirmed. We plot in Fig.4 for the $S-$wave $\bar{B}_{s}\to KK$ form factors in the $f_{0}+f^{\prime}_{0}+f^{\prime\prime}_{0}$ model, for convenience we also show the part of contributions from $f_{0}$ in blue dashed curves. It is easy to see the overwhelming role of $f_{0}$, while the contributions from $f^{\prime}_{0}$ and $f^{\prime\prime}_{0}$ account only $\sim 5\%$. The result at the full recoiled energy $\sqrt{q^{2}}F_{t}^{(l=0)}(1,0)/m_{B_{s}}=54.0^{+4.0}_{-7.0}$ is much larger than the result for $S-$wave $B\to\pi\pi$ form factors obtained in the LCSRs with $2\pi$DAs ChengHPQ ; HambrockAOR ; ChengSFK , with the asymptotic prediction $\sqrt{q^{2}}F_{t,\mathrm{asy}}^{(l=0)}(4m_{\pi}^{2},0)/m_{B}=5.40\pm 1.00$, this discrepancy is explained by the strong threshold effect of $f_{0}$ in the $\bar{B}_{s}\to KK$ decay.

5 Conclusion

In this paper we calculate the $\bar{B}_{s}\to f_{0}(980)$ form factor from the light-cone sum rules with $B-$meson DAs, and investigate the $S-$wave $\bar{B}_{s}\to KK$ form factors to study the width effect, basing on the assumption that $f_{0}$ is dominated by the $\bar{s}s$ configuration. With taking the conventional quark-antiquark assignment, we revisit the 2pSRs for the mass and decay constant of $f_{0}(980)$. For the $\bar{B}_{s}\to f_{0}(980)$ form factors, we find that the high twist two-particle and the three-particle $B-$meson DAs give $25\%$ correction separately, and their total correction to certain form factors can be about $50\%$. In order to investigate the width effect, we suggest the three resonance states model to parameterize the $S-$wave $\bar{B}_{s}\to KK$ form factors, the fitting result shows the dominant role of $f_{0}$, and as a by-product, suggest a new way to determine the strong coupling $|g_{f_{0}KK}|$. The residual uncertainty of our prediction mainly comes from the freedom to choose the widths of $f_{0}$.

Acknowledgments

We are grateful to Hai-yang Cheng and Ling-yun Dai for helpful discussions. This work is supported by the National Science Foundation of China under Grant No.118050 60 and No.11905056, ”the Fundamental Research Funds for the Central Universities” under Grant No.531118010176 and No.531118010258. S. C. is grateful to the High energy theory group at Institute of Physics, Academia Sinica for hospitality and for financial support where this work was finalized.

Appendix A $B$ meson LCDAs

Several models have been suggested for the LCDAs with definite twists, which incorporate the correct low-momentum behaviour and satisfy the (tree-level) equation of motion (EOM) constraints GrozinPQ ; BraunLIQ ; LuCFC . In BenekeWJP , a more general ansatz is proposed with comprising all of these models as particular cases.

For the two-particle $B_{s}$ meson DAs demonstrated in Eq.3, $\phi_{+}(\omega)$ and $\phi_{-}(\omega)$ are the leading and subleading twist DAs, and $g_{+}(\omega)$ and $g_{-}(\omega)$ are DAs at twist-4 and twist-5, respectively. The general model are quoted as BenekeWJP

[TABLE]

The general function $f(\omega)$ is normalized as $\int_{0}^{\infty}d\omega\omega f(\omega)=1$ and decreases sufficiently fast at $\omega\to\infty$. In Eqs.(72,74), three auxiliary functions are introduced to simplify the expression,

[TABLE]

The normalization constant $\varkappa$ is determined from the leading twist function $f(\omega)$ via the EOM relations

[TABLE]

Concerning the three-particle $B_{s}$ meson DAs, the definitions in Eq.19 by Lorentz structures should not be confused with the definitions by means of definite twists, and they have the following relations BraunLIQ

[TABLE]

The general model for the DAs with definite twists are BenekeWJP

[TABLE]

$\lambda_{E}$ and $\lambda_{H}$ are the parameters entered in the normalization conditions

[TABLE]

It is known that the EOM, as shown in Eqs.(78-80), imply the connections between the two-particle and three-particle LCDAs,

[TABLE]

with taking the general functions

[TABLE]

The normalization constants in these two particular models are

[TABLE]

We use the exponential models in our numerical evaluation.

Appendix B Coefficients in the three-particle correction

B.1 Correction coefficients to $\mathcal{F}_{+}(q^{2})$

[TABLE]

B.2 Correction coefficients to $\mathcal{F}_{-}(q^{2})$

[TABLE]

B.3 Correction coefficients to $\mathcal{F}_{T,p}(q^{2})$

[TABLE]

B.4 Correction coefficients to $\mathcal{F}_{T,q}(q^{2})$

[TABLE]

B.5 Correction coefficients to $\mathcal{F}_{0}(q^{2})$

[TABLE]