This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

\UseRawInputEncoding

Chiral excitations of open-beauty systems

Xiao-Yu Guo1,2, Matthias F.M. Lutz2,3 1 Institute of Theoretical Physics, Faculty of Science,
Beijing University of Technology, Beijing 100124, China
2 GSI Helmholtzzentrum für Schwerionenforschung GmbH,
Planck Str. 1, 64291 Darmstadt, Germany
3 Technische Universität Darmstadt, D-64289 Darmstadt, Germany
Abstract

We study the scattering of open-beauty mesons and Goldstone bosons as predicted by the chiral SU(3) Lagrangian. The impact of subleading-order chiral interactions to systems with JP=0+J^{P}=0^{+} and JP=1+J^{P}=1^{+} quantum numbers is worked out. We estimate the relevant low-energy coefficients from the open-charm sector, for which their values have been determined previously from sets of QCD lattice data. The leading-order heavy-quark symmetry-breaking effects are estimated by matching the BB-meson ground-state chiral mass formula to the mass formula from the heavy-quark effective theory. We make refined predictions for the flavor antitriplet and sextet resonances that are generated dynamically by coupled-channel interactions.

chiral Lagrangian, heavy quark effective theory, open-beauty mesons, coupled-channel scattering
pacs:
12.38.−t,12.38.Bx,12.39.Fe,12.39.Hg,14.65.Fy,21.30.Fe

I Introduction

How to understand the spectrum of QCD is the long-lasting challenge in contemporary theoretical physics. While the strong interactions lead to a myriad of hadronic resonances, QCD, the fundamental theory, fails to describe these phenomena if tackled with a perturbative expansion in the gauge coupling constant. In order to unravel nonperturbative aspects of QCD, effective field theories are commonly employed. They are constructed to respect the asymptotic symmetries of QCD in its low-energy limit. Those symmetries include the chiral SU(3) symmetry (which holds exactly when the quarks uu, dd, and ss are massless) and the heavy-quark symmetry (which asymptotically holds as the quarks cc or bb tend to be infinitely massive). An open-charm or open-beauty system, which is composed of one heavy quark cc or bb and light quarks, possesses both of these two symmetries and therefore plays a crucial role in studies of QCD.

Chiral SU(3) symmetry drives the interactions between an open heavy-flavor meson and Goldstone bosons. Such interactions can be described by an effective chiral Lagrangian. The leading order s-wave interaction may be attractive or repulsive and resonant states can thus be dynamically generated (see e.g. Lutz:2001yb ). According to such an interaction, the nature of the open-charm meson Ds0(2317)D_{s0}^{*}(2317) can be successfully explained. While within a quark model approach DiPierro:2001dwf ; Albaladejo:2018mhb such state poses a puzzle, coupled-channel scattering studies between Goldstone bosons and JP=0J^{P}=0^{-} DD-meson ground states suggest that it is a bound state below the DKDK threshold Kolomeitsev:2003ac ; Guo:2006fu ; Lang:2014yfa . It comes as a part of a flavor antitriplet, which is completed by a broad state with isospin-strangeness (I,S)=(1/2,0)(I,S)=(1/2,0) quantum numbers. The heavy-quark spin symmetry implies that a heavy-flavor hadron always comes with almost degenerate spin partners. This reflects the marginal importance of the spin orientation of its heavy-quark content. For instance, the Ds1(2460)D_{s1}(2460) can be understood as the spin 1 partner of the Ds0(2317)D_{s0}^{*}(2317). In addition, the heavy-quark flavor symmetry suggests that open-beauty partners of the open-charm states should exist. We note, however, that so far no experimental evidence for the open-beauty partner of the Ds0(2317)D_{s0}^{*}(2317) has been found. How stable are these theoretical predictions? Answering this question requires systematic studies that take into account chiral and heavy-quark symmetry-breaking effects as based on effective Lagrangians suitably linked to QCD. In our work we use chiral perturbation theory (ChPT) and heavy-quark effective theory (HQET).

Chiral symmetry-breaking effects are encoded into low-energy constants (LEC) of the chiral Lagrangian Lutz:2007sk ; Guo:2008gp ; Liu:2012zya ; Altenbuchinger:2013vwa ; Guo:2018tjx . In a recent study Guo:2018kno , the LECs relevant in the open-charm chiral Lagrangian have been well estimated from QCD lattice data on ensembles with a variety of unphysical quark masses. It was shown that the open-charm mesons Ds0(2317)D_{s0}^{*}(2317) and Ds1(2460)D_{s1}(2460) are stable against the next-to-leading (NLO) chiral corrections Guo:2018kno ; Guo:2018gyd . Moreover, the same work Guo:2018kno predicted the behavior of the (1/2,0)(1/2,0) antitriplet partner of the Ds0(2317)D_{s0}^{*}(2317) at unphysical quark masses, which has been recently confirmed by lattice simulations Gayer:2021xzv .

Motivated by our success in the open-charm sector, we would like to derive more quantitative predictions for the open-beauty partner systems. In this article, we will study the chiral symmetry-breaking effects in the open-beauty systems, working out systematically heavy-quark spin and flavor symmetry-breaking effects. Such a projection requests a matching of ChPT to HQET.

We will first recall the part of the chiral Lagrangian responsible for the leading-order symmetry-breaking LECs. The heavy-quark mass dependence of those LECs is determined by matching the chiral formula of the BB-meson masses to the HQET result. Based on such LECs, we will derive the pole positions of resonance states in the complex plane as generated by coupled-channel scattering from the chiral Lagrangian. Special attention will be paid to the antitriplet states. Predictions for the open-beauty partners of Ds0(2317)D_{s0}^{*}(2317) and Ds1(2460)D_{s1}(2460) will be made. A brief discussion on the flavor-exotic sextet states will follow Kolomeitsev:2003ac ; Hofmann:2003je . In particular, we will comment on the disputed isospin and strangeness one X(5568)X(5568) state D0:2016mwd ; Aaij:2016iev ; Aaltonen:2017voc ; Sirunyan:2017ofq ; Aaboud:2018hgx ; Agaev:2016mjb ; Lang:2016jpk ; Burns:2016gvy ; Lu:2016kxm ; Chen:2016mqt ; Liu:2016xly ; Yang:2016sws ; Ke:2018stp .

II The chiral Lagrangian with heavy-quark flavor symmetry

The pseudoscalar BB-meson ground states form a flavor antitriplet B=(B0,B+,Bs+)B=(B^{0},-B^{+},B_{s}^{+}). They are the heavy-quark flavor partners of the open-charm antitriplets D=(D0,D+,Ds+)D=(D^{0},-D^{+},D_{s}^{+}) states. In both flavor sectors the heavy-quark spin symmetry predicts almost degenerate JP=1J^{P}=1^{-} partner states. According to the heavy-quark flavor symmetry, the chiral Lagrangian for these two types of mesons take the same form. Therefore it is useful to introduce generic heavy-meson fields HH and HμνH_{\mu\nu}, that refer either to a DD or a BB-meson antitriplet field. The kinetic terms of the heavy-flavored mesons constitute the leading-order chiral Lagrangian

kin=(^μH)(^μH¯)(M¯34Δ)2HH¯\displaystyle\mathcal{L}_{\mathrm{kin}}=(\hat{\partial}_{\mu}H)(\hat{\partial}^{\mu}\bar{H})-\Big{(}\bar{M}-\frac{3}{4}\Delta\Big{)}^{2}\,H\,\bar{H}
(^μHμα)(^νH¯να)+12(M¯+14Δ)2HμαH¯μα,\displaystyle\ \ -(\hat{\partial}_{\mu}H^{\mu\alpha})(\hat{\partial}^{\nu}\bar{H}_{\nu\alpha})+\,\frac{1}{2}\Big{(}\bar{M}+\frac{1}{4}\Delta\Big{)}^{2}H^{\mu\alpha}\,\bar{H}_{\mu\alpha}, (1)

where following Guo:2018kno ; Guo:2018gyd the tensor field representation is adopted for the vector mesons. The LECs M¯34Δ\bar{M}-\frac{3}{4}\Delta and M¯+14Δ\bar{M}+\frac{1}{4}\Delta are the chiral-limit masses of the heavy pseudoscalar and vector mesons respectively. They depend on the heavy-quark mass MQM_{Q}. The M¯=M¯(MQ)\bar{M}=\bar{M}(M_{Q}) scales as M¯MQ\bar{M}\sim M_{Q} in the heavy-quark mass limit. The hyperfine splitting Δ=Δ(MQ)\Delta=\Delta(M_{Q}) is caused by O(1/MQ)O(1/M_{Q}) heavy-quark spin symmetry-breaking effect Wise:1992hn ; Goity:1992tp .

The chiral covariant derivative ^μ\hat{\partial}_{\mu} in (1)

^μH¯=μH¯+ΓμH¯,^μH=μHHΓμ.\displaystyle\hat{\partial}_{\mu}\bar{H}=\partial_{\mu}\,\bar{H}+\Gamma_{\mu}\,\bar{H}\,,\qquad\;\;\hat{\partial}_{\mu}H=\partial_{\mu}\,H-H\,\Gamma_{\mu}\,.
Γμ=12eiΦ2fμe+iΦ2f+12e+iΦ2fμeiΦ2f,\displaystyle\Gamma_{\mu}={\textstyle\frac{1}{2}}\,e^{-i\,\frac{\Phi}{2\,f}}\,\partial_{\mu}\,e^{+i\,\frac{\Phi}{2\,f}}+{\textstyle\frac{1}{2}}\,e^{+i\,\frac{\Phi}{2\,f}}\,\partial_{\mu}\,e^{-i\,\frac{\Phi}{2\,f}}, (2)

involves the flavor octet Goldstone-boson fields as encoded in the 3×33\times 3 matrix Φ\Phi. The parameter ff is the chiral limit of the pion-decay constant for which we choose f=92.4f=92.4 MeV. The leading-order, so called Tomozawa-Weinberg, interaction between the BB mesons and the Goldstone bosons is implied by the kinetic terms (1) via its covariant derivative and the chiral connection Γμ\Gamma_{\mu} .

At the NLO the chiral Lagrangian has terms that are proportional to the masses of the uu-, dd-, and ss quarks. Such terms break the chiral symmetry explicitly and define corrections to the heavy-flavor ground-state masses. We recall such terms in the chiral Lagrangian

χ=(4c02c1)HH¯trχ+2c1Hχ+H¯\displaystyle\mathcal{L}_{\chi}=-\big{(}4\,c_{0}-2\,c_{1}\big{)}\,H\,\bar{H}\,{\rm tr\,}\chi_{+}-2\,c_{1}\,H\,\chi_{+}\,\bar{H}
+(2c~0c~1)HμνH¯μνtrχ++c~1Hμνχ+H¯μν,\displaystyle\quad\;\;+\,\big{(}2\,\tilde{c}_{0}-\tilde{c}_{1}\big{)}\,H^{\mu\nu}\,\bar{H}_{\mu\nu}\,{\rm tr\,}\chi_{+}+\tilde{c}_{1}\,H^{\mu\nu}\,\chi_{+}\,\bar{H}_{\mu\nu}\,,
χ±=12(e+iΦ2fχ0e+iΦ2f±eiΦ2fχ0eiΦ2f),\displaystyle\chi_{\pm}={\textstyle\frac{1}{2}}\left(e^{+i\,\frac{\Phi}{2\,f}}\,\chi_{0}\,e^{+i\,\frac{\Phi}{2\,f}}\pm e^{-i\,\frac{\Phi}{2\,f}}\,\chi_{0}\,e^{-i\,\frac{\Phi}{2\,f}}\right), (3)

where the quark-mass dependence is embodied in the diagonal matrix χ0=2B0diag(mu,md,ms)\chi_{0}=2\,B_{0}\,{\rm diag}(m_{u},m_{d},m_{s}), with the low-energy constant B0B_{0}. Throughout this work we assume perfect isospin symmetry with mu=md=mm_{u}=m_{d}=m. Those parameters contribute to the JP=0J^{P}=0^{-} and JP=1J^{P}=1^{-} meson masses at NLO

MH2={(M¯34Δ)2+(4c02c1)ΠH(2),0+2c1ΠH(2),1+ΠHHOifH[0](M¯+14Δ)2+(4c~02c~1)ΠH(2),0+2c~1ΠH(2),1+ΠHHOifH[1],\displaystyle M_{H}^{2}=\left\{\begin{aligned} &\Big{(}\bar{M}-\frac{3}{4}\Delta\Big{)}^{2}+(4\,c_{0}-2\,c_{1})\,\Pi_{H}^{(2),0}+2\,c_{1}\Pi_{H}^{(2),1}+\Pi_{H}^{\rm HO}\qquad{\rm if}\qquad H\in[0^{-}]\\ &\Big{(}\bar{M}+\frac{1}{4}\Delta\Big{)}^{2}+(4\,\tilde{c}_{0}-2\,\tilde{c}_{1})\,\Pi_{H}^{(2),0}+2\,\tilde{c}_{1}\Pi_{H}^{(2),1}+\Pi_{H}^{\rm HO}\qquad{\rm if}\qquad H\in[1^{-}]\end{aligned}\right.\,, (4)

where the LECs ci=ci(MQ)c_{i}=c_{i}(M_{Q}) and c~i=c~i(MQ)\tilde{c}_{i}=\tilde{c}_{i}(M_{Q}) depend on the heavy-quark mass MQM_{Q}. Throughout this work we will use an superscript (c)(c) and (b)(b) with e.g.

ci(c)ci(MQ=Mc),c~i(b)c~i(MQ=Mb),\displaystyle c_{i}^{(c)}\equiv c_{i}(M_{Q}=M_{c})\,,\qquad\tilde{c}_{i}^{(b)}\equiv\tilde{c}_{i}(M_{Q}=M_{b})\,, (5)

for the open-charm and open-beauty systems, respectively. The chiral terms ΠH(2),i\Pi_{H}^{(2),i} are linear combinations of the light-quark masses

ΠH(2),0=2B0(2m+ms),\displaystyle\Pi_{H}^{(2),0}=2\,B_{0}(2\,m+m_{s})\,, (6)
ΠH(2),1={2B0mifH{D,D,B,B}2B0msifH{Ds,Ds,Bs,Bs}.\displaystyle\Pi_{H}^{(2),1}=\left\{\begin{array}[]{ll}2\,B_{0}\,m\quad\!&{\rm if}\!\!\quad H\in\{D,\,\,D^{*},B,\,\,B^{*}\}\\ 2\,B_{0}\,m_{s}\quad\!&{\rm if}\!\!\quad H\in\{D_{s},D_{s}^{*},B_{s},B_{s}^{*}\}\end{array}\right.\,. (9)

that contribute to chiral order Qχ2Q_{\chi}^{2}, with

QχB0mq,q=u,d,s.\displaystyle Q_{\chi}\sim\sqrt{B_{0}m_{q}},\qquad q=u,d,s. (10)

The higher-order chiral corrections ΠHHO\Pi_{H}^{\rm HO} start from O(Qχ3)O(Q_{\chi}^{3}) at the one-loop level.

The parameters cic_{i} and c~i\tilde{c}_{i} scale with the heavy-quark mass MQM¯M_{Q}\sim\bar{M} in the heavy-quark mass limit Wise:1992hn ; Goity:1992tp . It is useful to make this more explicit. We factor out their heavy-quark mass-independent part, denoted by CiΛχ1C_{i}\sim\Lambda_{\chi}^{-1}. They are of dimension 1-1 with Λχ\Lambda_{\chi} the chiral symmetry-breaking scale. Higher-order corrections account for the heavy-quark symmetry breaking effects. At order 1/M¯1/\bar{M}, those effects enter. They consist of an overall shift to CiC_{i} and a hyperfine splitting between the cic_{i} and c~i\tilde{c}_{i} Boyd:1994pa ; Brambilla:2017hcq . We introduce dimensionless parameters ζi\zeta_{i} and ηi\eta_{i} responsible for these two kinds of effects, and arrive at the following representation

ci(MQ)=M¯(MQ)(Ci+ζiM¯(MQ)34ηi(MQ)M¯(MQ)),\displaystyle c_{i}(M_{Q})=\bar{M}(M_{Q})\Big{(}C_{i}+\frac{\zeta_{i}}{\bar{M}(M_{Q})}-\frac{3}{4}\frac{\eta_{i}(M_{Q})}{\bar{M}(M_{Q})}\Big{)},
c~i(MQ)=M¯(MQ)(Ci+ζiM¯(MQ)+14ηi(MQ)M¯(MQ)).\displaystyle\tilde{c}_{i}(M_{Q})=\bar{M}(M_{Q})\Big{(}C_{i}+\frac{\zeta_{i}}{\bar{M}(M_{Q})}+\frac{1}{4}\frac{\eta_{i}(M_{Q})}{\bar{M}(M_{Q})}\Big{)}\,. (11)

The parameters ηi=ηi(MQ)\eta_{i}=\eta_{i}(M_{Q}) depend on MQM_{Q}. The remaining four parameters CiC_{i}, ζi\zeta_{i} turn out to be independent on MQM_{Q}. This will be seen in the next section, in which an explicit matching with HQET is worked out.

In application of the LECs cic_{i} and c~i\tilde{c}_{i} in the charm sector, the CiC_{i} can be determined modulo an unknown ζi\zeta_{i}-dependence

Ci=14M¯(c)(ci(c)+3c~i(c)4ζi),\displaystyle C_{i}=\frac{1}{4\bar{M}^{(c)}}\Big{(}c_{i}^{(c)}+3\,\tilde{c}_{i}^{(c)}-4\,\zeta_{i}\Big{)}\,, (12)

where we apply our notation ci(c)ci(Mc)c^{(c)}_{i}\equiv c_{i}(M_{c}). Similarly, the value of ηi\eta_{i} at MQ=McM_{Q}=M_{c} can be easily derived as

ηi(c)=(c~i(c)ci(c)),\displaystyle\eta_{i}^{(c)}=\big{(}\tilde{c}_{i}^{(c)}-c_{i}^{(c)}\big{)}\,, (13)

with again ηi(c)ηi(Mc)\eta^{(c)}_{i}\equiv\eta_{i}(M_{c}).

In this work, the one-loop result is employed for the higher-order chiral correction terms ΠHHO\Pi_{H}^{\rm HO}. At the one-loop level, they receive contributions from bubble and tadpole diagrams and their corresponding counterterms

ΠHHO=ΠHbubble+ΠHtadpole+ΠHCT.\displaystyle\Pi_{H}^{\rm HO}=\Pi_{H}^{\rm bubble}+\Pi_{H}^{\rm tadpole}+\Pi_{H}^{\rm CT}\,. (14)

For explicit expressions, Refs. Guo:2018kno ; Guo:2018gyd are referred to. While more LECs are involved, they have been determined in the charm sector in Ref Guo:2018kno . Using the leading-order scaling behavior, ΠHHOM¯(MQ)\Pi_{H}^{\rm HO}\sim\bar{M}(M_{Q}), the LECs in the bottom sector can be well estimated.

III Matching the chiral Lagrangian and the HQET

So far we have not yet fully specified the LEC c0,1(b)c^{(b)}_{0,1} and c~0,1(b)\tilde{c}^{(b)}_{0,1} in (11). While at leading order in the 1/MQ1/M_{Q} expansion their values may be inferred from c0,1c_{0,1} and c~0,1\tilde{c}_{0,1} at MQ=McM_{Q}=M_{c} this is no longer true at subleading order. Since such counterterms contribute at NLO in the chiral counting scheme it appears reasonable to consider the effect of order 1/MQ1/M_{Q}. Note that our approach considers chiral N2LO effects in the heavy pseudoscalar and vector-meson masses. We now set up a more detailed matching with HQET, in which the heavy-meson masses take the form Falk:1992wt

MH(MQ)={MQ+Λ¯(H)+μπ(H)22MQμG(H)22MQifH[0]MQ+Λ¯(H)+μπ(H)22MQ+μG(H)26MQifH[1].\displaystyle M_{H}(M_{Q})=\left\{\begin{aligned} &M_{Q}+\bar{\Lambda}_{(H)}+\frac{\mu^{2}_{\pi(H)}}{2\,M_{Q}}-\frac{\mu_{G(H)}^{2}}{2\,M_{Q}}\qquad{\rm if}\qquad H\in[0^{-}]\\ &M_{Q}+\bar{\Lambda}_{(H)}+\frac{\mu^{2}_{\pi(H)}}{2\,M_{Q}}+\frac{\mu_{G(H)}^{2}}{6\,M_{Q}}\qquad{\rm if}\qquad H\in[1^{-}]\end{aligned}\right.\,. (15)

The quantities Λ¯\bar{\Lambda} and μπ2\mu^{2}_{\pi}, μG2\mu^{2}_{G} come with an explicit index (H)(H), that resolves the specifics of the light-quark content. The Λ¯\bar{\Lambda} is the contribution from light degrees of freedom, and therefore MQM_{Q} independentManohar:2000dt . The μπ2\mu_{\pi}^{2} term accounts for the kinetic energy of the heavy quark in the meson’s rest frame. Due to reparametrization invariance, it is MQM_{Q} independent as well Luke:1992cs . Finally, the μG2\mu_{G}^{2} is a chromomagnetic moment which leads to a hyperfine splitting between the 00^{-} and 11^{-} B-mesons. It depends on the heavy-quark mass. We assume that μG2\mu_{G}^{2} can be factorized as a product of the high-energy and low-energy contributions,

μG(H)2=C^cm(MQ)μ^G(H)2,\displaystyle\mu_{G(H)}^{2}=\hat{C}_{\rm cm}(M_{Q})\,\hat{\mu}_{G(H)}^{2}\,, (16)

where the factor μ^G\hat{\mu}_{G} accounts for low-energy contributions. The high-energy contributions are incorporated in the renormalization-group (RG) invariant Wilson coefficient C^cm(MQ)\hat{C}_{\rm cm}(M_{Q})Neubert:1993mb . The RG evolution starts at a scale close to the heavy-quark mass μMQ\mu\sim M_{Q}, where the value of the Wilson coefficient is determined by matching the chromomagnetic moment from HQET to the multiloop calculations from the QCD Lagrangian. For our purposes it suffices to know the ratio C^cm(Mb)/C^cm(Mc)\hat{C}_{\rm cm}(M_{b})/\hat{C}_{\rm cm}(M_{c}). This ratio has been derived at the one-loop and two-loop level in Falk:1990pz and Amoros:1997rx ; Czarnecki:1997dz respectively. The averaged result is

RC^cm(Mb)C^cm(Mc)0.80(4),\displaystyle R\equiv\frac{\hat{C}_{\rm cm}(M_{b})}{\hat{C}_{\rm cm}(M_{c})}\simeq 0.80(4)\,, (17)

where the uncertainty is estimated by the difference of the one-loop and the two-loop results. In the latest calculation, a poor convergence pattern has been claimed at the three-loop level Grozin:2007fh . Therefore we refrain from using the three-loop result here.

The expansion moments Λ¯\bar{\Lambda}, μπ2\mu_{\pi}^{2}, and μ^G\hat{\mu}_{G} depend on physical scales significantly lighter than MQM_{Q}: the light-quark masses mqm_{q} and an intrinsic nonperturbative QCD scale Λ\Lambda. While Λ\Lambda was once commonly regarded as the Landau scale ΛQCD\Lambda_{\rm QCD} in literature Isgur:1989vq ; Falk:1990pz ; Kitazawa:1993bk , we follow the way of Refs. Grinstein:1993ys ; Boyd:1994pa and identify ΛΛχ\Lambda\sim\Lambda_{\chi}. It has been demonstrated in Randall:1992ww that only via such an assignment, ChPT and HQET can be matched convincingly at the loop level. Moreover, the size of Λ\Lambda should be comparable to the mass difference M¯MQ\bar{M}-M_{Q} Falk:1990pz ; Boyd:1994pa . Using numerical values, this mass difference is 0.8\sim 0.8 GeV for both charm and bottom systems, indeed consistent with estimates of Λχ\Lambda_{\chi}.

In the chiral limit, the Λ¯\bar{\Lambda} scales as Λ\sim\Lambda whereas the μπ2\mu_{\pi}^{2} and μ^G2\hat{\mu}_{G}^{2} scale as Λ2\sim\Lambda^{2}. We can expand the components Λ¯\bar{\Lambda}, μπ2\mu_{\pi}^{2} and μ^G2\hat{\mu}_{G}^{2} in powers of the light-quark masses around their chiral limit. The corrections are suppressed by powers of B0mq/Λχ\sqrt{B_{0}m_{q}}/\Lambda_{\chi}. The small-scale expansion scheme entails the scaling behavior Δμ^G2/M¯Qχ\Delta\sim\hat{\mu}_{G}^{2}/\bar{M}\sim Q_{\chi}, and therefore the expansion parameters in the chiral and heavy-quark expansions are comparable

B0mqΛχΛMQ.\displaystyle\frac{\sqrt{B_{0}m_{q}}}{\Lambda_{\chi}}\sim\frac{\Lambda}{M_{Q}}\,. (18)

By matching the mass formula (15) with the chiral result (4), we can recover the relations between the M¯,Δ\bar{M},\Delta and the heavy-quark moments Wise:1992hn ; Brambilla:2017hcq . It is emphasized, that the chiral structure of (4) restricts the structure of the O(Qχ2)O(Q_{\chi}^{2}) corrections to the heavy-quark expansion moments. From the matching we obtain that, CiC_{i} and ζi\zeta_{i} are involved in the O(Qχ2)O(Q_{\chi}^{2}) corrections of Λ¯\bar{\Lambda} and μπ2\mu_{\pi}^{2} respectively. And they are indeed heavy-quark mass independent. In addition ηi\eta_{i} contributes to the O(Qχ2)O(Q_{\chi}^{2}) corrections of μG2\mu_{G}^{2}, and its scaling behavior is proportional to C^cm\hat{C}_{\rm cm}. We summarize,

M¯(b)Δ(b)M¯(c)Δ(c)=ηi(b)ηi(c)=R.\displaystyle\frac{\bar{M}^{(b)}\Delta^{(b)}}{\bar{M}^{(c)}\Delta^{(c)}}=\frac{\eta_{i}^{(b)}}{\eta_{i}^{(c)}}=R\,. (19)

Using (19) with R0.80R\simeq 0.80 together with (13) we are left with three unknown parameters M¯(b)\bar{M}^{(b)} and ζ0,1\zeta_{0,1}, given the charm-sector LECs M¯(c)\bar{M}^{(c)}, Δ(c)\Delta^{(c)}, ci(c),c~i(c)c_{i}^{(c)},\tilde{c}_{i}^{(c)} as inputs. The MQM_{Q} independent parameters CiC_{i} follow from (12).

The three unknown parameters M¯(b)\bar{M}^{(b)} and ζ0,1\zeta_{0,1} are determined by a fit with our chiral mass formula to the empirical values of the four BB-meson ground-state masses. Here we admit a residual systematic uncertainty of 55 MeV in the heavy-meson masses. Such a value was used in our previous open-charm system studies Guo:2018kno . It reflects the accuracy level at which we expect our one-loop chiral formula to hold. The results of M¯(b)\bar{M}^{(b)}, ζ0,1\zeta_{0,1} are shown in Tab. 1. In this table, we also show the parameters involved in the expansion (11) together with the associated LECs. In the fits, we used the results of Guo:2018kno , for the inputs of the LECs in the charm sector. In the charm sector, four sets of fitted results are determined according to the lattice data on DD-meson ground-state masses and πD\pi D s-wave scattering process. They are named as Fits 1-4. We will recall some of the fitting details in the following discussion.

Fit 1 Fit 2 Fit 3 Fit 4
M¯(b)\bar{M}^{(b)}[GeV] 5.37435.3743 4.85404.8540 5.33035.3303 5.36665.3666
ζ0\zeta_{0} 0.09210.0921 1.5072-1.5072 0.0839-0.0839 0.05230.0523
ζ1\zeta_{1} 0.16890.1689 0.12330.1233 0.15850.1585 0.16780.1678
C0[GeV1]C_{0}[{\rm GeV}^{-1}] 0.06020.0602 0.87770.8777 0.17740.1774 0.11450.1145
C1[GeV1]C_{1}[{\rm GeV}^{-1}] 0.23760.2376 0.39160.3916 0.33820.3382 0.34450.3445
η0(b)\eta_{0}^{(b)} 0.0145-0.0145 0.0302-0.0302 0.0176-0.0176 0.0170-0.0170
η1(b)\eta_{1}^{(b)} 0.0238-0.0238 0.03180.0318 0.0276-0.0276 0.0238-0.0238
Δ(b)\Delta^{(b)}[GeV] 0.05620.0562 0.06430.0643 0.05630.0563 0.05680.0568
c0(b)c_{0}^{(b)} 0.42620.4262 2.77572.7757 0.87500.8750 0.67970.6797
c~0(b)\tilde{c}_{0}^{(b)} 0.41170.4117 2.74552.7455 0.85740.8574 0.66280.6628
c1(b)c_{1}^{(b)} 1.46371.4637 2.00022.0002 1.98201.9820 2.03452.0345
c~1(b)\tilde{c}_{1}^{(b)} 1.43991.4399 2.03202.0320 1.95441.9544 2.01072.0107
χ2/N\chi^{2}/N 0.940.94 0.100.10 0.890.89 0.920.92
Table 1: The low-energy parameters in (11), corresponding to Fits 1-4 in the charm sector Guo:2018kno . The χ2/N\chi^{2}/N is the chi-square per data point, with the number of data points N=4N=4 and an ad hoc systematic error estimate of 55 MeV.

Consider first the scenarios of Fit 1, 3, 4. The masses of the BB-meson ground states can be reproduced within the systematic error of 5 MeV. All of the 3 fits give modest heavy-quark corrections to the leading-order expectations of c0,1c_{0,1} and c~0,1\tilde{c}_{0,1}. At leading order, c0,1c_{0,1} and c~0,1\tilde{c}_{0,1} are about 2.5 times larger at MQ=MbM_{Q}=M_{b} as compared to their values at MQ=McM_{Q}=M_{c}. For convenience we recall the ranges c0c~0(0.20.3)c_{0}\sim\tilde{c}_{0}\sim(0.2-0.3) and c1c~1(0.60.9)c_{1}\sim\tilde{c}_{1}\sim(0.6-0.9) at MQ=McM_{Q}=M_{c} from Guo:2018kno .

Scenarios 3 and 4 show quite similar values for the LEC. This is not the case for scenario 1. Here we recall a decisive distinction. Both Fit 1 and Fit 2 did not consider QCD lattice data on the πD\pi D s-wave scattering process Moir:2016srx . While Fit 1, nevertheless, appears reasonably consistent with the πD\pi D phase shift and inelasticity parameters as given in Moir:2016srx , this is not the case for Fit 2. The key feature of Fit 3 and also Fit 4 is their compatibility with the lattice data on the ηD\eta D phase shift. Such data play a crucial role in the determination of the LEC. Based on this observation we would disfavor scenarios 1 and 2. In this context it is amusing to observe that Fit 2 should be rejected also based on an unnaturally large value of the ζ0\zeta_{0} parameter as shown in Tab. 1. This is so despite the fact that it comes with the best chi-square value for the reproduction of the B-meson masses. The corresponding c0c~02.7c_{0}\sim\tilde{c}_{0}\sim 2.7 at MQ=MbM_{Q}=M_{b}, are nearly ten times larger than their charmed counterparts c0c~00.3c_{0}\sim\tilde{c}_{0}\sim 0.3. This implied a serious violation of the leading-order scaling behavior. Such large LECs lead to unnaturally large higher order corrections. Therefore, altogether we exclude Fit 2 from our further analysis.

IV Scatterings with coupled-channel dynamics

The chiral Lagrangian predicts the formation of JP=0+J^{P}=0^{+} and JP=1+J^{P}=1^{+} resonance state as a consequence of coupled-channel final-state interactions of the Goldstone bosons with the ground-state heavy mesons with JP=0J^{P}=0^{-} and JP=1J^{P}=1^{-} quantum numbers. Such states are an unavoidable consequence of the flavor SU(3) chiral Lagrangian. We focus on the resonances generated by the ss-wave scatterings with open-beauty quantum numbers. Here the leading order coupled-channel interaction is predicted by the Tomozawa-Weinberg theorem in terms of the ”known” parameter ff.

A resonance is dynamically generated from a scattering process when the reaction amplitude contains a pole in the complex ss plane. The scattering amplitude with manifest ss-channel unitarity is obtained from a self-consistent summation

Tab(s)=Vab(s)+c,dVac(s)Jcd(s)Tdb(s),\displaystyle T_{ab}(s)=V_{ab}(s)+\sum_{c,\,d}V_{ac}(s)\,J_{cd}(s)\,T_{db}(s)\,, (20)

with given out- and in-going two-body states aa and bb. The TT-matrix exhibits poles in the complex ss plane that we can determine by extending the definition of T(s)T(s) into the higher Riemann sheets of the ss plane. The potential Vab(s)V_{ab}(s) is obtained with an on-shell condition Lutz:2001yb , and set equivalent to the scattering amplitude as derived from the chiral Lagrangian at the matching point s=μM\sqrt{s}=\mu_{M}. It receives tree-level chiral symmetry-breaking contributions from the Lagrangian (3). Additional terms are implied by the LEC that imply the one-loop ΠHHO\Pi_{H}^{\rm HO} structures decomposed in (14). The latter will also contribute to Vab(s)V_{ab}(s) as tree-level chiral symmetry-preserving corrections, see Ref Guo:2018kno . Following Kolomeitsev:2003ac , we set μM=MB()\mu_{M}=M_{B^{(*)}} for S=0,2S=0,2 and μM=MBs()\mu_{M}=M_{B_{s}^{(*)}} for S=±1S=\pm 1 scatterings with the total quantum numbers JP=0+(1+)J^{P}=0^{+}(1^{+}). The diagonal analytic matrix J(s)J(s) function is universal as it leads to a scattering amplitude T(s)T(s) that is consistent with the coupled-channel unitarity condition and the microcausality condition. Such an approach can be justified if short-range forces largely dominate the system. Along the real axis from each threshold, there is a branch cut defining the doorway for the higher Riemann sheets. For an nn-dimensional coupled-channel system there are 2n2^{n} Riemann sheets, and we use the signature (±,,±)(\pm,\cdots,\pm) as introduced in Guo:2018gyd to label a specific one.

(I,S)=(1/2,0)(I,S)=(1/2,0) (I,S)=(0,1)(I,S)=(0,1)
JP=0+J^{P}=0^{+}
Fit 1 5.5202209+2820.0923254+463i5.5202^{+282}_{-209}-0.0923^{+463}_{-254}\,i 5.6296395+4315.6296^{+431}_{-395}
Fit 3 5.5137225+1790.1073387+623i5.5137^{+179}_{-225}-0.1073^{+623}_{-387}\,i 5.5689544+6145.5689^{+614}_{-544}
Fit 4 5.5126196+1540.1120384+595i5.5126^{+154}_{-196}-0.1120^{+595}_{-384}\,i 5.5755535+6105.5755^{+610}_{-535}
TW 5.5207190+2600.0905243+456i5.5207_{-190}^{+260}-0.0905^{+456}_{-243}\,i 5.6495310+3535.6495^{+353}_{-310}
JP=1+J^{P}=1^{+}
Fit 1 5.5652208+2830.0915252+460i5.5652^{+283}_{-208}-0.0915^{+460}_{-252}\,i 5.6763391+4285.6763^{+428}_{-391}
Fit 3 5.5595224+1790.1071384+624i5.5595^{+179}_{-224}-0.1071^{+624}_{-384}\,i 5.6179542+6105.6179^{+610}_{-542}
Fit 4 5.5586198+1560.1115382+597i5.5586^{+156}_{-198}-0.1115^{+597}_{-382}\,i 5.6242533+6055.6242^{+605}_{-533}
TW 5.5658190+2600.0903242+455i5.5658_{-190}^{+260}-0.0903^{+455}_{-242}\,i 5.6959311+3545.6959^{+354}_{-311}
Table 2: Complex pole masses (in GeV) of the flavor antitriplet states with JP=0+J^{P}=0^{+} and 1+1^{+}. The relevant Riemann sheets are (,+,+)(-,+,+) and (+,+)(+,+) for (I,S)=(1/2,0)(I,S)=(1/2,0) and (0,1)(0,1).

We start with a discussion of the flavor antitriplet channels. Poles in the complex ss-plane are found from the parameter sets 1, 3 and 4. The complex pole masses are compared with the results from the Tomozawa-Weinberg interaction in Tab. 2. To estimate the theoretical error, we allow a deviation of the matching points from their natural values for |ΔμM|=0.1|\Delta\mu_{M}|=0.1 GeV. We first look at the states with isospin and strangeness (I,S)=(0,1)(I,S)=(0,1). A pole below the BKBK threshold is found on the physical Riemann sheet in the JP=0+J^{P}=0^{+} scattering amplitude always. It is the open-beauty partner of Ds0(2317)D_{s0}^{*}(2317) . The pole mass is 5.59(8)5.59(8) GeV. Comparing to the leading-order result at 5.65(3)5.65(3) GeV, the higher-order chiral corrections slightly reduce the pole mass. This result is somewhat lower than previous predictions with values above 5.7 GeV Colangelo:2012xi ; Altenbuchinger:2013vwa ; Lang:2015hza ; Du:2017zvv . In the axial-vector sector, the open-beauty partner of Ds1(2460)D_{s1}(2460) is found as a bound state at 5.64(8)5.64(8) GeV, which should be compared with previous predictions at above 5.75 GeV Colangelo:2012xi ; Altenbuchinger:2013vwa ; Lang:2015hza ; Du:2017zvv . Besides the bound states in the (I,S)=(0,1)(I,S)=(0,1) channels, broad resonances are found in the (I,S)=(1/2,0)(I,S)=(1/2,0) channels with poles in the unphysical Riemann sheet denoted by (,+,+)(-,+,+). Their broad charmed partners were extensively discussed in previous theoretical studies Kolomeitsev:2003ac ; Hofmann:2003je ; Lutz:2007sk ; Altenbuchinger:2013vwa ; Yao:2015qia ; Du:2017zvv . Our prediction of their pole masses are (5.52(3)0.12(5)i)\big{(}5.52(3)-0.12(5)\,i\big{)} GeV for the 0+0^{+} state and (5.57(3)0.12(5)i)\big{(}5.57(3)-0.12(5)\,i\big{)} GeV for the 1+1^{+} state. Both of them are in agreement with previous theoretical predictions Kolomeitsev:2003ac ; Guo:2006rp ; Altenbuchinger:2013vwa ; Du:2017zvv .

(I,S)=(1/2,0)(I,S)=(1/2,0) (I,S)=(1,1)(I,S)=(1,1)
JP=0+J^{P}=0^{+}
Fit 1 5.8104099+0950.0148+010023i5.8104^{+095}_{-099}-0.0148^{-023}_{+010}\,i 5.8045205+2610.0730064+127i5.8045^{+261}_{-205}-0.0730^{+127}_{-064}\,i
Fit 3 5.7274144+2580.0595+158230i5.7274^{+258}_{-144}-0.0595^{-230}_{+158}\,i 5.7945154+1470.1871219+266i5.7945^{+147}_{-154}-0.1871^{+266}_{-219}\,i
Fit 4 5.7422171+2850.0517+186233i5.7422^{+285}_{-171}-0.0517^{-233}_{+186}\,i 5.7868134+1330.1938190+234i5.7868^{+133}_{-134}-0.1938^{+234}_{-190}\,i
TW 5.7730149+1480.0201+722i5.7730^{+148}_{-149}-0.0201^{-22}_{+7}\,i 5.7900245+3190.071952+107i5.7900^{+319}_{-245}-0.0719^{+107}_{-52}i
JP=1+J^{P}=1^{+}
Fit 1 5.8592096+0910.0153+011025i5.8592^{+091}_{-096}-0.0153^{-025}_{+011}\,i 5.8558208+2640.0730063+127i5.8558^{+264}_{-208}-0.0730^{+127}_{-063}\,i
Fit 3 5.7747141+2540.0638+160237i5.7747^{+254}_{-141}-0.0638^{-237}_{+160}\,i 5.8406150+1460.1913211+256i5.8406^{+146}_{-150}-0.1913^{+256}_{-211}\,i
Fit 4 5.7894166+2790.0564+186241i5.7894^{+279}_{-166}-0.0564^{-241}_{+186}\,i 5.8333132+1330.1972184+227i5.8333^{+133}_{-132}-0.1972^{+227}_{-184}\,i
TW 5.8192150+1480.0198+722i5.8192^{+148}_{-150}-0.0198^{-22}_{+7}\,i 5.8362246+3200.070851+105i5.8362^{+320}_{-246}-0.0708^{+105}_{-51}i
Table 3: Complex pole masses (in GeV) of two flavor sextet states with JP=0+J^{P}=0^{+} and 1+1^{+}. The relevant Riemann sheets are (,+,+)(-,+,+) and (,+)(-,+) for (I,S)=(1/2,0)(I,S)=(1/2,0) and (1,1)(1,1).

In the charmed and beauty sector, there are further poles belonging to a flavor sextet. In Tab. 3, we listed the complex pole masses found in the sextet channels with (I,S)=(1/2,0)(I,S)=(1/2,0) and (1,1)(1,1) from Fits 1,3, and 4. The other sextet channel, with (I,S)=(0,1)(I,S)=(0,-1), shows a pole in the vicinity of the scattering threshold within the range of our theoretical uncertainty always. The resonances with (I,S)=(1/2,0)(I,S)=(1/2,0) were also obtained in previous works Guo:2006rp ; Du:2017zvv with masses (5.845.85)(5.84-5.85) GeV and (5.885.91)(5.88-5.91) GeV, respectively, for JP=0+J^{P}=0^{+} and 1+1^{+}. Our predicted values are significantly smaller than those. The JP=0+J^{P}=0^{+} member with (I,S)=(1,1)(I,S)=(1,1) is of particular interest. It has been speculated that the controversial X(5568)X(5568) state has such quantum numbers (see e.g. Chen:2016spr ; Brambilla:2019esw ). In Kolomeitsev:2003ac , it was shown that the LO Tomozawa-Weinberg interaction implies the existence of an exotic state with JP=0+J^{P}=0^{+} and (I,S)=(1,1)(I,S)=(1,1) quantum numbers at s(5.79(3)0.07(1)i)\sqrt{s}\simeq\big{(}5.79(3)-0.07(1)\,i\big{)} GeV. We demonstrate that this prediction is quite stable against higher-order corrections from our Fits 1,3, and 4. The pole mass comes at (5.80(3)0.14(8)i)\big{(}5.80(3)-0.14(8)\,i\big{)} GeV as shown in Tab. 3. We conclude that the X(5568)X(5568) cannot be a chiral excitation, i.e. it is not explained convincingly by chiral coupled-channel dynamics. This supports previous such claims Albaladejo:2016eps ; Sun:2016tmz .

V Summary

We studied open-beauty mesons with JP=0+J^{P}=0^{+} and 1+1^{+} quantum numbers. Such states were predicted as a consequence of coupled-channel interactions based on the chiral SU(3) Lagrangian. Already the leading order Tomozawa-Weinberg interaction implies attractive forces in the flavor antitriplet and sextet channels between the Goldstone bosons and the heavy-meson ground states with JP=0J^{P}=0^{-} and 11^{-} quantum numbers.

In this article the role of the next-to-leading order chiral interactions in the s-wave open-beauty meson scattering processes was scrutinized. The LECs are derived mainly from corresponding LECs as obtained previously from global fits to the QCD lattice dataset in the charm sector Guo:2018kno . Where possible additional direct data form the beauty sector were taken into account. The heavy-quark scaling behavior is constrained by the RG-invariant Wilson coefficient for the chromomagnetic moment. We employ the results of the Wilson coefficient calculated at the 2-loop level.

We find that in the antitriplet but also in the exotic flavor sextet, the chiral correction terms lead to minor effects in the JP=0+J^{P}=0^{+} and 1+1^{+} pole masses only. This confirms the semiquantitative predictions made almost two decades ago by one of the authors. Our refined values should be used in ongoing experimental searches and QCD lattice simulations. It is noted, however, that like in the open-charm sector, we expect the light-quark mass dependence in the flavor antitriplet and sextet states to be significant. This should be investigated further, in particular for the πB\pi B s-wave phase shifts.

VI Acknowledgements

Xiang-Dong Gao and Daniel Mohler are acknowledged for stimulating discussions.

References