Study of open-charm $0^+$ states in unitarized chiral effective theory with one-loop potentials

TL;DR

This paper derives and unitarizes chiral potentials for Goldstone boson and pseudoscalar charmed meson interactions, validating simplified models and exploring open-charm $0^+$ states as dynamically generated poles.

Contribution

It introduces a covariant chiral effective theory with explicit vector charmed mesons and demonstrates the validity of $D^*$-less potentials for scattering length calculations.

Findings

Negligible difference between models with and without explicit vector mesons for S-wave scattering.

Successfully unitarized one-loop potentials to fit lattice QCD data.

Identified poles corresponding to open-charm $0^+$ states and analyzed their trajectories.

Abstract

Chiral potentials are derived for the interactions between Goldstone bosons and pseudoscalar charmed mesons up to next-to-next-to-leading order in a covariant chiral effective field theory with explicit vector charmed-meson degrees of freedom. Using the extended-on-mass-shell scheme, we demonstrate that the ultraviolet divergences and the so-called power counting breaking terms can be properly absorbed by the low-energy constants of the chiral Lagrangians. We calculate the scattering lengths by unitarizing the one-loop potentials and fit them to the data extracted from lattice QCD. The obtained results are compared to the ones without an explicit contribution of vector charmed mesons given previously. It is found that the difference is negligible for $S$-wave scattering in the threshold region. This validates the use of $D^\ast$-less one-loop potentials in the study of the pertinent scattering lengths. We search for dynamically generated open-charm states with $J^P=0^+$ as poles of the $S$-matrix on various Riemann sheets. The trajectories of those poles for varying pion masses are presented as well.

Tables

Table 1. Table 1: The coefficients in the tree-level amplitudes for the ten relevant physical processes, with Δ K ​ π = M K 2 − M π 2 subscript Δ 𝐾 𝜋 superscript subscript 𝑀 𝐾 2 superscript subscript 𝑀 𝜋 2 \Delta_{K\pi}=M_{K}^{2}-M_{\pi}^{2} .

	Physical process	${𝒞}_{\text{LO}}$	${𝒞}_S$	${𝒞}_U$	${𝒞}_0^{\left(2\right)}$	${𝒞}_1^{\left(2\right)}$	${𝒞}_{24}^{\left(2\right)}$	${𝒞}_{35}^{\left(2\right)}$	${𝒞}_{1a}^{\left(3\right)}$	${𝒞}_{1b}^{\left(3\right)}$	${𝒞}_{23}^{\left(3\right)}$
1	$D^0{K}^{-}\to D^0{K}^{-}$	1	0	2	$M_K^2$	$-M_K^2$	1	1	$M_K^2$	0	1
2	$D^{+}{K}^{+}\to D^{+}{K}^{+}$	0	0	0	$M_K^2$	0	1	0	0	0	0
3	$D^{+}{\pi }^{+}\to D^{+}{\pi }^{+}$	1	0	2	$M_{\pi }^2$	$-M_{\pi }^2$	1	1	$M_{\pi }^2$	0	1
4	$D^{+}\eta \to D^{+}\eta$	0	$\frac{1}{3}$	$\frac{1}{3}$	$M_{\eta }^2$	$-\frac{1}{3}{M}_{\pi }^2$	1	$\frac{1}{3}$	0	0	0
5	$D_s^{+}{K}^{+}\to D_s^{+}{K}^{+}$	1	0	2	$M_K^2$	$-M_K^2$	1	1	$M_K^2$	0	1
6	$D_s^{+}\eta \to D_s^{+}\eta$	0	$\frac{4}{3}$	$\frac{4}{3}$	$M_{\eta }^2$	$\frac{4}{3}\left(M_{\pi }^2-2M_K^2\right)$	1	$\frac{4}{3}$	0	0	0
7	$D_s^{+}{\pi }^0\to D_s^{+}{\pi }^0$	0	0	0	$M_{\pi }^2$	0	1	0	0	0	0
8	$D^0\eta \to D^0{\pi }^0$	0	$\frac{1}{\sqrt{3}}$	$\frac{1}{\sqrt{3}}$	0	$-\frac{1}{\sqrt{3}}{M}_{\pi }^2$	0	$\frac{1}{\sqrt{3}}$	0	0	0
9	$D_s^{+}{K}^{-}\to D^0{\pi }^0$	$-\frac{1}{\sqrt{2}}$	$\sqrt{2}$	0	0	$-\frac{1}{2\sqrt{2}}\left(M_K^2+M_{\pi }^2\right)$	0	$\frac{1}{\sqrt{2}}$	$-\frac{1}{\sqrt{2}}{M}_K^2$	$\frac{1}{\sqrt{2}}{\Delta }_{K\pi }$	$-\frac{1}{\sqrt{2}}$
10	$D_s^{+}{K}^{-}\to D^0\eta$	$-\sqrt{\frac{3}{2}}$	$\sqrt{\frac{2}{3}}$	$-\sqrt{\frac{8}{3}}$	0	$\frac{1}{2\sqrt{6}}\left(5M_K^2-3M_{\pi }^2\right)$	0	$-\frac{1}{\sqrt{6}}$	$-\sqrt{\frac{3}{2}}{M}_K^2$	$\frac{-1}{\sqrt{6}}{\Delta }_{K\pi }$	$-\sqrt{\frac{3}{2}}$

Table 2. Table 2: Parameters for the chiral extrapolation for different configurations. h 0 subscript ℎ 0 h_{0} , L 4 r superscript subscript 𝐿 4 𝑟 L_{4}^{r} and L 5 r superscript subscript 𝐿 5 𝑟 L_{5}^{r} are fixed by the data in Ref.  Liu:2012zya . The masses and decay constant in the chiral limit are in units of GeV. h 0 subscript ℎ 0 h_{0} , h 1 subscript ℎ 1 h_{1} , L 4 r superscript subscript 𝐿 4 𝑟 L_{4}^{r} and L 5 r superscript subscript 𝐿 5 𝑟 L_{5}^{r} are dimensionless. The asterisk indicates an input value.

(37)

Table 3. Table 3: Values of the LECs from the 6-channel fits using the method of UChPT. The h i subscript ℎ 𝑖 h_{i} ’s are dimensionless, and the g 1 ′ superscript subscript 𝑔 1 ′ g_{1}^{\prime} , g 23 subscript 𝑔 23 g_{23} and g 3 ′ superscript subscript 𝑔 3 ′ g_{3}^{\prime} are in GeV -1 .

(50)

Table 4. Table 4: Poles in the coupled-channel amplitudes based on UChPT-6(b) in Table 4 of Ref.  Yao:2015qia . Physical masses and decay constants are used to obtain the poles. The Riemann sheets on which the poles are located are indicated in the last column.

(57)

Table 5. Table 5: Poles in the single-channel amplitudes based on UChPT-6(b) in Table 4 of Ref.  Yao:2015qia . Physical masses are used to obtain the poles. The Riemann sheets on which the poles are located are indicated in the last column.

(62)

Equations

$$\displaystyle\mathscr{L}_{\rm eff}=\sum_{i=1}^{2}\mathscr{L}_{\phi\phi}^{(2i)}+\sum_{j=1}^{3}\mathscr{L}_{D\phi}^{(j)}+\sum_{k=1}^{2}\mathscr{L}_{D^{\ast}\phi}^{(k)}+\sum_{l=1}^{3}\mathscr{L}_{D^{\ast}D\phi}^{(l)}\$$

$${\mathcal{L}}_{\varphi \varphi }^{(2)}$$

$${\mathcal{L}}_{\varphi \varphi }^{(4)}$$

$$\varphi =\left(\begin{array}{ccc}\frac{1}{\sqrt{2}}{\pi }^0+\frac{1}{\sqrt{6}}\eta & {\pi }^{+}& K^{+}\\ {\pi }^{-}& -\frac{1}{\sqrt{2}}{\pi }^0+\frac{1}{\sqrt{6}}\eta & K^0\\ K^{-}& {\bar{K}}^0& -\frac{2}{\sqrt{6}}\eta \end{array}\right).$$

$${\mathcal{L}}_{D\varphi }^{(1)}$$

$${\mathcal{L}}_{D\varphi }^{(2)}$$

$${\mathcal{L}}_{D\varphi }^{(3)}$$

$$\displaystyle u_{\mu}=i\bigg{[}u^{\dagger}\partial_{\mu}u-u\partial_{\mu}u^{\dagger}\bigg{]}\ ,\qquad U=u^{2}\ ,\qquad\chi^{\pm}=u^{\dagger}\chi u^{\dagger}\pm u\chi^{\dagger}u\ .$$

$${\mathcal{D}}_{\mu }H=H(\overleftarrow{{\partial }_{\mu }}+{\Gamma }_{\mu }^{†}),{\mathcal{D}}_{\mu }{H}^{†}=({\partial }_{\mu }+{\Gamma }_{\mu })H^{†},$$

$${\mathcal{L}}_{D^{\ast }\varphi }^{(1)}$$

$${\mathcal{L}}_{D^{\ast }\varphi }^{(2)}$$

$${\mathcal{L}}_{D^{\ast }D\varphi }^{(1)}=i{g}_0\left(D_{\mu }^{\ast }{u}^{\mu }{D}^{†}-D{u}^{\mu }{D}_{\mu }^{\ast †}\right).$$

$${\mathcal{V}}_{D_1{\varphi }_1\to D_2{\varphi }_2}(s,t)={\mathcal{V}}_{\mathrm{LO}}^{(\mathrm{WT})}+{\mathcal{V}}_{\mathrm{LO}}^{(\mathrm{EX})}+{\mathcal{V}}_{\mathrm{NLO}}^{(\mathrm{CT})}+{\mathcal{V}}_{\mathrm{NNLO}}^{(\mathrm{CT})}+{\mathcal{V}}_{\mathrm{NNLO}}^{(\mathrm{Loop})}.$$

$${\mathcal{V}}_{\mathrm{LO}}^{(\mathrm{WT})}(s,t)$$

$${\mathcal{V}}_{\mathrm{LO}}^{(\mathrm{EX})}(s,t)$$

$${\mathcal{V}}_{\mathrm{NLO}}^{(\mathrm{CT})}(s,t)$$

$${\mathcal{V}}_{\mathrm{NNLO}}^{(\mathrm{CT})}(s,t)$$

$${\mathcal{F}}_S(s,t)$$

$${\mathcal{F}}_U(s,t)$$

$${\mathcal{H}}_{24}(s,t,u)$$

$${\mathcal{H}}_{35}(s,t,u)$$

$${\mathcal{G}}_{23}(s,t,u)$$

$$M_0^2$$

$$M_0^{\ast 2}$$

$$h_i$$

$$g_j$$

$$h_i^r(\mu )$$

$$g_0^r(\mu )$$

$${\mathcal{V}}_l^{(S,I)}(s{)}_{D_1{\varphi }_1\to D_2{\varphi }_2}=\frac{1}{2}{\int }_{-1}^{+1}\mathrm{d}\cos \theta P_l(\cos \theta ){\mathcal{V}}_{D_1{\varphi }_1\to D_2{\varphi }_2}^{(S,I)}(s,t(s,\cos \theta )),$$

$$t(s,\cos \theta )$$

$$T(s)=[1-N(s)\cdot G(s){]}^{-1}\cdot N(s),$$

$$G(s)=i\int \frac{d^{4}q}{(2\pi {)}^4}\frac{1}{(q^2-M_D^2+iϵ)\left[(q+p{)}^2-M_{\varphi }^2+iϵ\right]},s\equiv p^2.$$

$$G(s)$$

$$N(s)={\mathcal{V}}_{\text{LO}}^{\text{(WT+EX)}}(s)+{\mathcal{V}}_{\text{NLO}}^{\text{(CT)}}(s)+{\mathcal{V}}_{\text{NNLO}}^{\text{(CT+Loop)}}(s)-{\mathcal{V}}_{\text{LO}}^{\text{(WT+EX)}}(s)\cdot G(s)\cdot {\mathcal{V}}_{\text{LO}}^{\text{(WT+EX)}}(s).$$

$$a_{D\varphi \to D\varphi }^{(S,I)}=-\frac{1}{8\pi (M_D+M_{\varphi })}{T}_{l=0}^{(S,I)}(s_{\text{thr}}{)}_{D\varphi \to D\varphi },s_{\text{thr}}=(M_D+M_{\varphi }{)}^2.$$

$$M_K=\sqrt{{\mathring{M}}_K^2+M_{\pi }^2/2},M_D={\mathring{M}}_D+(h_1+2h_0)\frac{M_{\pi }^2}{{\mathring{M}}_D},M_{D_s}={\mathring{M}}_{D_s}+2h_0\frac{M_{\pi }^2}{{\mathring{M}}_{D_s}},$$

$${\mathring{M}}_D^2=M_0^2+4h_0{\mathring{M}}_K^2,{\mathring{M}}_{D_s}^2=M_0^2+4(h_0+h_1){\mathring{M}}_K^2.$$

Full text

Chiral Lagrangians, Meson-meson interaction, Charmed mesons

I Introduction

In the past two decades, many excited charmed states have been observed experimentally Bevan:2014iga ; Aubert:2003fg ; Krokovny:2003zq ; Besson:2003cp and further experiments are intended either to investigate their properties more precisely or to search for more new states, e.g., by the LHCb Collaboration Aaij:2015vea . The conventional quark-potential models provide a successful description of most of those low-lying excitations, see Ref. Chen:2016spr for a recent review. However, quantum chromodynamics (QCD) at low energies has a much richer structure than quark models. There exist observed charmed mesons whose properties are in disagreement with the expectations from quark models, of which the most interesting one is the $D_{s0}^{*}(2317)$. It was first observed by the BABAR Collaboration in the inclusive $D_{s}^{+}\pi^{0}$ invariant mass distribution and later confirmed by Belle and CLEO Collaborations Aubert:2003fg ; Krokovny:2003zq ; Besson:2003cp . It couples to the $DK$ channel and decays mainly into the isospin breaking channel $D_{s}\pi$ due to its location below the $DK$ threshold. Many theoretical investigations were triggered consequently, attempting to reveal the nature of the $D_{s0}^{*}(2317)$ as well as other newly observed charmed states with $J^{P}=0^{+}$ and trying to reveal their internal structure. For instance, the $D_{s0}^{*}(2317)$ has been suggested to be a $DK$ bound state Barnes:2003dj . Were this interpretation true, one can learn much about the interaction between the $D/D_{s}$ mesons and $\pi/K$ mesons from studying the $D_{s0}^{*}(2317)$ and related states. This path has been followed in Refs. Kolomeitsev:2003ac ; Guo:2006fu ; Gamermann:2006nm ; Hofmann:2003je ; Guo:2008gp ; Guo:2009ct ; Cleven:2010aw ; Cleven:2014oka where the $S$-wave interaction between charmed $D$ mesons and Goldstone bosons (denoted as $\phi$ hereafter) has been studied systematically up to the next-to-leading order (NLO) using chiral perturbation theory (ChPT) for heavy mesons Burdman:1992gh ; Wise:1992hn ; Yan:1992gz in combination with a unitarization procedure such as the one in Ref. Oller:2000fj .

In the meantime, significant progress has also been made in lattice QCD Moir:2013ub ; Cichy:2015tma . Using the Lüscher formalism and its extension to coupled channels (for early works on this topic, see e.g. Refs. Liu:2005kr ; Lage:2009zv ), scattering lengths and recently phase shifts for the $D\phi$ interaction have been calculated at unphysical quark masses Liu:2008rza ; Liu:2012zya ; Mohler:2012na ; Mohler:2013rwa ; Lang:2014yfa ; Moir:2016srx . The first calculation only concerns the channels free of disconnected Wick contractions Liu:2008rza ; Liu:2012zya , i.e., $D\pi$ with isospin $I=3/2$, $D\bar{K}$ with $I=0,1$, $D_{s}K$ and $D_{s}\pi$. The channels with disconnected Wick contractions such as $D\pi$ with $I=1/2$ and $DK$ with $I=0$ were calculated later in Refs. Mohler:2012na ; Mohler:2013rwa ; Lang:2014yfa . On the one hand, the lattice results can be used to determine the low-energy constants (LECs) in the chiral Lagrangian Wang:2012bu ; Liu:2012zya ; Altenbuchinger:2013vwa ; Yao:2015qia ; Guo:2015dha . On the other hand, with these lattice calculations more insights into the nature of the $D_{s0}^{*}(2317)$ and other positive-parity charmed mesons are obtained. In particular, in Ref. Liu:2012zya , it is concluded that the lattice calculation of other channels performed there supports the interpretation that $D_{s0}^{*}(2317)$ is dominantly a $DK$ hadronic molecule. In addition, using the parameters fixed in that work, energy levels in the $I=1/2$ channel were computed in Ref. Albaladejo:2016lbb , and a remarkable agreement with the lattice results reported in Ref. Moir:2016srx was found. This agreement was taken to be as a strong evidence that the particle listed as $D_{0}^{*}(2400)$ in the Review of Particle Physics Olive:2016xmw in fact corresponds to two states with poles located at $\left(2105^{+6}_{-8}-i\,102^{+10}_{-12}\right)$ MeV and $\left(2451^{+36}_{-26}-i\,134^{+7}_{-8}\right)$ MeV, respectively Albaladejo:2016lbb , similar to the well-known two-pole scenario of the $\Lambda(1405)$ Oller:2000fj . In this scenario, the puzzle that the non-strange $D_{0}^{*}(2400)$ has a mass larger than the strange partner $D_{s0}^{*}(2317)$ can be easily understood. The poles were searched for in unitarized ChPT with the interaction kernel computed at NLO. In view of the phenomenological importance of the $D_{s0}^{*}(2317)$ and $D_{0}^{*}(2400)$, it is crucial to check the stability of the NLO predictions by extending to the next-to-next-to-leading order (NNLO), which is one of the purposes of this work.

When massive matter fields are included in ChPT, the nonvanishing matter-field mass in the chiral limit leads to the notable power counting breaking (PCB) issue Gasser:1987rb : all loop graphs containing internal matter field propagators start contributing at $\mathcal{O}\left(p^{2}\right)$.111The closed matter field loops are not taken into account since they are real below the two-matter-field threshold and can be absorbed by the redefinition of LECs Gasser:1987rb . Various approaches have been proposed to remedy the PCB issue, e.g. heavy baryon ChPT Jenkins:1990jv ; Bernard:1992qa , infrared regularization Becher:1999he , and the extended-on-mass-shell (EOMS) scheme Fuchs:2003qc . Recently, the EOMS scheme has been demonstrated to be a good solution to the PCB problem. It has been successfully applied to the study of $\pi N$ scattering up to $\mathcal{O}(p^{3})$ Alarcon:2012kn and $\mathcal{O}(p^{4})$ Chen:2012nx , and up to leading one-loop order in the presence of $\Delta$-resonance Yao:2016vbz .

The first aim of this paper is to present a full calculation of $D\phi$ scattering using the EOMS scheme within a manifestly Lorentz invariant chiral effective theory with explicit vector charmed mesons, to be denoted as $D^{*}$, up to NNLO, i.e. the leading one-loop order. The first study on $D\phi$ scattering to one-loop was made in Ref. Liu:2009uz in the framework of non-relativistic heavy meson ChPT, which neglects sizeable recoil corrections Geng:2010vw .222 Such recoil corrections can be restored by using the extended heavy baryon propagator $i/(v\cdot k+k^{2}/2m)$ instead of $i/v\cdot k$ Bernard:1993ry . The first one-loop calculation within the covariant formalism was given in Ref. Geng:2010vw . However, in that paper the NNLO contact terms of $D\phi$ are not included and the kinetic term for the $D^{\ast}$ is incomplete, both of which are necessary for a proper renormalization in the EOMS scheme. Furthermore, the scattering amplitudes are calculated perturbatively without considering resonant charmed mesons close-to or even below thresholds such as the $D_{s0}^{*}(2317)$ in the channel $(S,I)=(1,0)$, where $S$ and $I$ denote the strangeness and isospin, respectively. Thus, the results in those channels of such a calculation are incomplete and thus can not be considered significant. In our previous work Yao:2015qia , the $D\phi$ scattering amplitudes are presented up to NNLO explicitly in the absence of the $D^{*}$. Then, the scattering lengths are calculated based on the unitarized amplitudes, which are also used to fit to the lattice QCD data at unphysical pion masses. In order to judge the importance of the $D^{*}$, a selection of all the diagrams containing the $D^{*}$ is calculated and found to be negligible. However, a complete calculation including the $D^{\ast}$ mesons is still lacking and, furthermore, a systematic renormalization using the EOMS scheme is required. It is also important to check whether the full finite $c$-quark mass effects, corresponding to including the $D^{*}$ which are degenerate with the $D$ mesons in the heavy quark limit, are sizeable. These gaps will be closed in this paper.

Specifically, in Section II, we derive the covariant one-loop $D\phi$ scattering potentials with the $D^{*}$ resonances as dynamical degrees of freedom. We perform renormalization in the EOMS scheme and explicitly show that the ultraviolet (UV) divergences and PCB terms can be absorbed by redefining the LECs. The $D^{\ast}$-less case has been accomplished formally using the path integral formalism in Refs. Du:2016ntw ; Du:2016xbh . Then a unitarization procedure is taken to generate resonances not far from the corresponding thresholds. In such a scheme, one can also deal with larger pion masses as compared with a purely perturbative approach. However, in general crossing symmetry is no longer exactly fulfilled. When the unitarization is extended to the one-loop order, an additional subtraction in the potentials is needed to remove the right-hand cut in the $N(s)$ functions, see below Oller:2000fj ; Guo:2011pa ; Yao:2015qia . In Section III, by fitting the so-obtained scattering lengths to the lattice results, we determine the LECs in the effective Lagrangian. Then we search for poles in the unitarized amplitudes, and study their trajectories with varying pion mass. Section IV comprises a brief summary. The explicit UV-part and EOMS subtractions of the LECs are collected in Appendix A.

II Theoretical discussions of the $D\phi$ interactions

II.1 Effective Lagrangian

To set up the effective Lagrangian, we first specify the corresponding power counting rules. At low energies, the external momenta as well as the masses of the Goldstone bosons are counted as $\mathcal{O}\left(p\right)$. However, the nonvanishing masses of the $D$ and $D^{*}$ in the chiral limit introduce new scales $M_{0}$ and $M_{0}^{*}$, both counted as $\mathcal{O}\left(1\right)$. As a result, at low energies, the temporal components of the momenta of the $D$ and $D^{*}$ are counted as $\mathcal{O}\left(1\right)$, while the spatial components are counted as $\mathcal{O}\left(p\right)$. Therefore, the virtuality $q^{2}-M_{0}^{(*)2}$ in the propagators scales as $\mathcal{O}\left(p\right)$, and the propagators scale as $\mathcal{O}\left(p^{-1}\right)$. The Goldstone boson propagators are counted as $\mathcal{O}\left(p^{-2}\right)$ as usual. Based on the counting rules for the vertices and propagators, one can assign a chiral order for a given Feynman diagram, and thus for any physical quantity. However, for the specific Feynman graphs with loops, there exist terms with chiral order lower than the naive power counting order, which are called PCB terms. In the EOMS scheme, the PCB terms are absorbed into the redefinition of the LECs so that the resulting physical observables obey the power counting rules.

The effective Lagrangian relevant to our calculation of the $D\phi$ potentials up to leading one-loop order can be written as

[TABLE]

with the superscripts specifying the chiral dimension. The needed terms in the Goldstone sector read Gasser:1984gg

[TABLE]

where the trace in flavor space is denoted by $\langle\cdots\rangle$, $F_{0}$ is the pion decay constant in the chiral limit, and $L_{4,5}$ are LECs. Furthermore, $\chi=2B_{0}\,{\rm diag}(m_{u},m_{d},m_{s})$, with $B_{0}$ a constant related to the quark condensate, and $U=\exp\left({i\sqrt{2}\phi}/{F_{0}}\right)$, with

[TABLE]

The terms corresponding to interactions between the $D=(D^{0},D^{+},D_{s}^{+})$ mesons and the Goldstone bosons are given by Burdman:1992gh ; Wise:1992hn ; Yan:1992gz ; Guo:2008gp ; Yao:2015qia 333As discussed in Ref. Du:2016ntw , the dimension-three terms proportional to $g_{4,5}$ in that work do not contribute to the $D\phi$ scattering because their contribution from contact term to amplitudes will be canceled out by their contribution to the wave function renormalization.

[TABLE]

where $h_{i}$ and $g_{j}$ are LECs and the chiral building blocks are given by

[TABLE]

The covariant derivative is defined via

[TABLE]

and $\mathcal{D}^{\mu\nu\rho}=\{\mathcal{D}_{\mu},\{\mathcal{D}_{\nu},\mathcal{D}_{\rho}\}\}$, where $H\in\{D,D^{\ast}\}$ with $D^{\ast}=(D^{\ast 0},D^{\ast+},D_{s}^{\ast+})$. The so-called chiral connection in the covariant derivatives is defined as $\Gamma_{\mu}=\left(u^{\dagger}\partial_{\mu}u+u\partial_{\mu}u^{\dagger}\right)/2.$ Similarly, the relevant terms for the interaction between the $D^{\ast}$ and the Goldstone bosons are Burdman:1992gh ; Wise:1992hn ; Yan:1992gz

[TABLE]

with $\tilde{h}_{0,1}$ analogous to $h_{0,1}$ and $\mathcal{F}_{\mu\nu}=(\mathcal{D}_{\mu}D^{\ast}_{\nu}-\mathcal{D}_{\nu}D^{\ast}_{\mu})$. Finally, the LO axial coupling has the form

[TABLE]

As pointed out in Refs. Krebs:2009bf ; Yao:2016vbz , the resonance-exchange contributions of $\mathcal{O}(p^{2})$ and $\mathcal{O}(p^{3})$ can be taken into account by shifting the coupling in the LO resonance-exchange contribution and the LECs in the contact terms. This also holds true for our case. Thus, we do not need the $\mathcal{O}\left(p^{2}\right)$ and $\mathcal{O}\left(p^{3}\right)$ terms for the $D^{\ast}D\phi$ coupling.

II.2 Chiral potentials up to leading one-loop order

Up to NNLO, the Feynman diagrams needed for our calculation are displayed in Fig. 1. Accordingly, the chiral potential for the process $D_{1}(p_{1})\phi_{1}(p_{2})\to D_{2}(p_{3})\phi_{2}(p_{4})$ can be written as

[TABLE]

As usual, the Mandelstam variables are defined by $s=(p_{1}+p_{2})^{2}$ and $t=(p_{1}-p_{3})^{2}$, while $u$ can be obtained via $u=\sum_{i=1}^{2}(M_{D_{i}}^{2}+M_{\phi_{i}}^{2})-s-t$. The potentials at tree-level are given by

[TABLE]

where the involved coefficients corresponding to various scattering processes are shown in Table 1. The functions in the $D^{*}$-exchange potentials read

[TABLE]

The functions in the NLO potentials read

[TABLE]

while the one in the NNLO potentials is

[TABLE]

As for the one-loop potentials at NNLO, the parts without explicit $D^{\ast}$ mesons can be found in the appendix of Ref. Yao:2015qia and the ones involving explicit $D^{\ast}$ states are too lengthy to be shown here. Note that $\mathcal{V}^{\rm(Loop)}_{\rm NNLO}$ in Eq. (9) contains the contribution from wave function renormalization as well. We performed renormalization of the one-loop potentials using the so-called EOMS scheme. In this scheme, the UV divergence are absorbed by the counterterms when the bare LECs are expressed in terms of the renormalized ones via

[TABLE]

where $R=\frac{2}{d-4}+\gamma_{E}-1-\ln(4\pi)$, with $\gamma_{E}$ the Euler constant and $d$ the number of space-time dimensions. The $\beta$-functions are given in Appendix A.1. Here, $\mu$ is the scale introduced in dimensional regularization. Then additional subtractions are performed by splitting the UV-renormalized LECs via

[TABLE]

such that the PCB terms from the one-loop potentials are canceled. The remaining LECs $g_{1}$, $g_{2}$ and $g_{3}$ are untouched at the chiral order we are working. The coefficients $\bar{\beta}_{h_{i}}$ and $\bar{\beta}_{g_{0}}$ are given in Appendix A.2.

II.3 Partial waves and unitarization

In the present paper, we do not consider the effect of isospin violation. It is convenient to study the potentials in the isospin basis instead of the particle basis. All possible processes with definite strangeness $S$ and isospin $I$ can be obtained from the ten processes given in Table 1 by crossing and isospin symmetry, see Refs. Guo:2009ct ; Yao:2015qia for details.

Since the standard ChPT is organized in a double expansion in terms of small external momenta and light quark masses, it is expected to work well in the low-energy region. With increasing energy the convergence of the chiral series becomes worse. Especially, when the energy reaches the region where resonances appear, the perturbative chiral potentials start to violate unitarity largely and cannot be directly applied anymore. One way to restore unitarity is to unitarize the potentials, but usually at the price of violating the crossing symmetry. While the unitarity and analyticity of the single-channel $\pi\pi$ potentials are strictly restored within a range of energies Pelaez:2015qba , a rigorous solution for the coupled-channel case is still missing. A convenient approximation is to treat the right-hand cut nonperturbatively, while the cross-channel effects are incorporated in a perturbative manner Oller:1997ti ; Oller:1998hw .444A method of calculating the left-hand cut nonperturbatively was proposed very recently in Ref. Entem:2016ipb . The unitarization is equivalent to a resummation of the $s$-channel potentials, and can extend the applicable energy range of the perturbative amplitudes. For instance, the scattering data for the pion-kaon systems up to $1.2$ GeV can be well described Oller:1997ti ; Oller:1998hw ; Oller:1998zr .

Before unitarization, the partial wave projection to a definite orbital angular momentum $l$ should be performed

[TABLE]

where $\theta$ is the scattering angle between the incoming and outgoing particles in the center-of-mass frame, and the Mandelstam variable $t$ is expressed as

[TABLE]

where $\lambda(a,b,c)=a^{2}+b^{2}+c^{2}-2ab-2ac-2bc$ is the Källén function. We only deal with the $S$-wave scattering in this paper, and will drop the subscript $l=0$ for brevity.

The unitarized two-body scattering amplitude has the form Oller:2000fj

[TABLE]

where the function $G(s)$ encodes the two-body right-hand cut and is given by the two-point loop function

[TABLE]

The explicit expression for $G(s)$ reads Oller:1998zr

[TABLE]

where $a(\mu)$ is a subtraction constant with $\mu$ the renormalization scale and $\sigma(s)=\sqrt{[s-(M_{D}+M_{\phi})^{2}][s-(M_{D}-M_{\phi})^{2}]}$. Note that the logarithmic scale dependence can be absorbed into the subtraction constance $a(\mu)$ and we do not distinguish the scale $\mu$ with the one introduced by dimensional regularization in the perturbative one-loop potentials.

While the right-hand cut effect is collected in the $G(s)$ function, the $N(s)$ function is free of any two-body right-hand cut. However, it may include the left-hand cuts due to the crossed channels. Up to NNLO, the $N(s)$ function can be expressed as Oller:2000fj ; Guo:2011pa

[TABLE]

Eq. (23) is an algebraic approximation of the standard $N/D$ methodOller:1998zr , and it should be understood in the matrix form for coupled-channels, for which $G(s)=\text{diag}\{G_{i}(s)\}$, with $i$ the channel index.

III Numerical analyses

III.1 Fit to lattice data of the scattering lengths

Up to now, there is no experimental measurement on the light pseudoscalar mesons scattering off heavy bosons. We can only rely on lattice QCD results Liu:2012zya ; Mohler:2012na ; Mohler:2013rwa ; Lang:2014yfa ; Moir:2016srx to determine the relevant LECs. We will fit to the lattice results on the scattering lengths. For the channels with definite strangeness and isospin, the $S$-wave scattering lengths are obtained from the unitarized amplitudes ${T}(s)$ via

[TABLE]

Since the current lattice simulations are performed at unphysical pion masses with fixed charm and strange quark masses, in order to fit these lattice data, one needs to know the pion mass dependence of the scattering lengths, which is achieved by employing the following mass extrapolation formulae: Liu:2012zya

[TABLE]

where $\mathring{M}_{K}$, $\mathring{M}_{D}$ and $\mathring{M}_{D_{s}}$ denote the masses in the two-flavor chiral limit ($M_{\pi}^{2}(\propto\hat{m})\to 0$ but with the fixed strange quark mass $m_{s}$). They have the form

[TABLE]

Using Eqs. (28) and (29), one gets

[TABLE]

which is fixed as $h_{1}=0.427$ with the physical masses, i.e., $M_{\pi}^{\text{Phy}}=0.138~{}\text{GeV}$, $M_{K}^{\text{Phy}}=0.496~{}\text{GeV}$, $M_{D}^{\text{Phy}}=1.867~{}\text{GeV}$ and $M_{D_{s}}^{\text{Phy}}=1.968~{}\text{GeV}$. Similar to the case of the pseudoscalar charmed mesons, the pion mass dependence of the masses of the vector mesons read, consistent with the general expression derived in Ref. Bruns:2004tj ,

[TABLE]

Here, $\mathring{M}_{D^{*}}$ and $\mathring{M}_{D_{s}^{*}}$ denote the corresponding two-flavor chiral limit masses, which can be estimated by the relations

[TABLE]

with $M_{D^{*}}^{\text{Phy}}$ and $M_{D_{s}^{*}}^{\text{Phy}}$ denoting the corresponding physical masses, $2.008~{}\text{GeV}$ and $2.112~{}\text{GeV}$, respectively. One has $\tilde{h}_{1}=h_{1}$ and $\tilde{h}_{0}=h_{0}$ in the heavy quark limit.555Analogous to Eq. (30), it is easy to see that $\tilde{h}_{1}=\left(M_{D_{s}^{\ast}}^{2}-M_{D^{\ast}}^{2}\right)/\left[4\left(M_{K}^{2}-M_{\pi}^{2}\right)\right]=0.472$ , which is close to $h_{1}$ numerically. These relations as well as similar relations for other LECs will be employed in order to reduce the number of parameters. The $DD^{*}\pi$ axial coupling constant $g_{0}$ can be fixed by the decay width $\Gamma_{D^{*+}\to D^{0}\pi^{+}}$. As discussed in Refs. Yao:2015qia , one gets $g=(1.113\pm 0.147)~{}\text{GeV}$ for the renormalized coupling $g$, which contains the bare constants $g_{0}$ and one-loop chiral corrections. At the one-loop level, the pion decay constant has the form Gasser:1984gg

[TABLE]

with $\mu_{\phi}=\frac{M_{\phi}^{2}}{16\pi^{2}F_{0}^{2}}\ln\frac{M_{\phi}^{2}}{\mu^{2}}$.

In this paper, the scattering length data we use are taken from Ref. Liu:2012zya and Refs. Mohler:2012na ; Mohler:2013rwa ; Lang:2014yfa , respectively. The lattice simulation data for $M_{K}$, $M_{D}$, $M_{D_{s}}$ and $F_{\pi}$, $F_{K}$ are taken from Refs. Liu:2012zya and WalkerLoud:2008bp , which share the same ensembles (M007, M010, M020 and M030) with Ref. Liu:2012zya . The fit results are listed in the Table 3 in Ref. Yao:2015qia . In addition, we also include the $DK$ scattering length with $(S,I)=(1,0)$ at $M_{\pi}=0.156~{}\text{GeV}$ Mohler:2013rwa , as discussed in Ref. Yao:2015qia . It should be noticed that different lattice configurations usually take different values for both the strange and charm quark masses, which leads to different values for $\mathring{M}_{K}$, $\mathring{M}_{D}$ and $\mathring{M}_{D_{s}}$, as listed in Table 2.

Furthermore, in order to reduce the correlations between the LECs, we introduce the following redefinitions of the LECs Liu:2012zya ; Yao:2015qia

[TABLE]

where $\bar{M}_{D}$ stands for the average of the physical masses of the charmed mesons $D$ and $D_{s}$, $\bar{M}_{D}=(M_{D}^{\text{Phy}}+M_{D_{s}}^{\text{Phy}})/2$. The new parameters $h_{4}^{\prime}$, $h_{5}^{\prime}$, $h_{24}$ and $h_{35}$ are dimensionless, and $g_{1}^{\prime}$, $g_{3}^{\prime}$ and $g_{23}$ have the dimension of inverse mass. They are fixed by fitting to the lattice scattering lengths at varying pion masses. It is well known that the state $D_{s0}^{*}(2317)$ in the $(S,I)=(1,0)$ channel is produced as a bound state pole below the $DK$ threshold Liu:2012zya ; Guo:2015dha ; Guo:2009ct . Due to the large number of the parameters and the small number of data, we constrain further the subtraction constant $a(\mu)$ by requiring the existence of a bound state pole at $2.317~{}\text{GeV}$ in the $(S,I)=(1,0)$ amplitude when all the parameters take their physical values Liu:2012zya .

To compare with the result of Ref. Yao:2015qia where the $D^{*}$ mesons are not included, we utilize the same fit procedures. In the fit UChPT-6(a), we fit all of the data in Ref. Liu:2012zya , including 5 channels at pion masses $0.301~{}\text{GeV}$, $0.364~{}\text{GeV}$, $0.511~{}\text{GeV}$ and $0.617~{}\text{GeV}$, as well as the isoscalar $DK$ channel at the pion mass $0.156~{}\text{GeV}$ in Ref. Mohler:2013rwa , because all these refer to $N_{f}=3$. However, as we know, the standard ChPT only works well in the small pion mass and low energy region. It is naively expected that the unitarized approach has a larger convergence range, but the convergence for the pion mass larger than $0.6~{}\text{GeV}$ is still questionable. Therefore, to compare with the fit UChPT-6(a), another fit denoted as UChPT-6(b) is performed excluding the lattice data at $M_{\pi}=0.617~{}\text{GeV}$. Results of both fits are listed in Table 3.

In both fits, as in Ref. Yao:2015qia , the absolute value of the dimensionless LEC $h_{5}^{\prime}$ is much larger than 1, too large to be natural. We therefore use the same method as therein to perform further fits by minimizing the augmented $\chi^{2}$

[TABLE]

where $\chi^{2}$ is the standard chi-squared, and $\chi_{\rm prior}^{2}$ is a prior quantity constraining the LECs to take natural values. It is set to be the sum of squares of the free LECs . Two fits UChPT-6($a^{\prime}$) and UChPT-6($b^{\prime}$) are obtained by minimizing $\chi_{\rm aug}^{2}$ instead of $\chi^{2}$ using the same data as in UChPT-6(a) and UChPT-6(b). Although the values of LECs are more natural, the $\chi^{2}$ values, with the prior parts subtracted, become very large. As a result, the scattering lengths from the new fits have larger deviations from the lattice data. More details about the fit procedures can be found in Ref. Yao:2015qia .

Compared to the NLO fits Guo:2015dha ; Guo:2009ct , the NNLO fits have larger $\chi^{2}$ values, even though three more LECs $g_{i}\,(i=1,2,3)$ are included in the fits. This could be because the unitarization method we use works better for the tree-level potentials than one-loop ones. On the one hand, the left-hand cuts, stemming from the $t$- and $u$- channels, appear in the one-loop potentials, which would cause the problem of violation of right-hand unitarity in the region where the left- and right- cuts overlap, see more discussions in the next section. On the other hand, the off-shell effects are partially included in the unitarized amplitudes if the one-loop potentials are employed. Both of the above-mentioned effects have non-trivial analytical structures and could make the NNLO unitarization much more cumbersome than the NLO one. Related to this is the fact that the scattering length $a_{D_{s}\pi\to D_{s}\pi}^{(1,1)}$ remains sizeable in the $SU(2)$ chiral limit of $M_{\pi}\to 0$, as shown in Figs. 2 and 3. This was also the case in Ref. Yao:2015qia .666It is due to the nonvanishing $\mathring{M}_{K}$ in loops contributing to the $D_{s}\pi\to D_{s}\pi$ potential in the $SU(2)$ chiral limit. Near the $D_{s}\pi$ threshold, the $u$- and $t$-channel loops dominate in the $SU(2)$ chiral limit, which indicates that the left-hand cuts are non-negligible. In the NLO case, the scattering length $a_{D_{s}\pi\to D_{s}\pi}^{(1,1)}$ is negligible in the $SU(2)$ chiral limit Liu:2012zya ; Guo:2015dha .

Among various fits, UChPT-6(b) has the smallest $\chi^{2}$, which is also true for the previous fits without $D^{*}$ Yao:2015qia . In addition, the fits with and without dynamical $D^{*}$ have similar values of the chi-squared and the LECs, which indicates that the influence of the $D^{\ast}$ on the quantities in question is small.

III.2 Dynamically generated resonances

The unitary $S$-matrix could have poles in the complex energy ($\sqrt{s}$) plane in the region not far from the relevant thresholds. Bound states and resonances are poles located on the physical and unphysical Riemann sheets, respectively. Different Riemann sheets are characterized by the sign of the imaginary part of the loop function on the right branch cuts. Each loop function $G_{i}(s)$ has two sheets: the physical/first Riemann sheet and the unphysical/second Riemann sheet, denoted as $G^{i}_{I}(s)$ and $G^{i}_{II}(s)$, respectively. The expression in Eq. (25) defines the physical Riemann sheet, while the expression on the second sheet is given by analytic continuation via Oller:1997ti

[TABLE]

For the $n$-channel case, there exist $2^{n}$ Riemann sheets in total. Different sheets can be accessed by properly choosing the loop functions $G^{i}_{I/II}(s)$. We use the sign of the imaginary part of $G^{i}(s)$ above threshold to indicate the $G^{i}_{I/II}(s)$. In this convention, for the coupled-channel case, the first Riemann sheet is labelled as $(+,+,+,\ldots)$, while $(-,+,+,\ldots)$, $(-,-,+,\ldots)$, $(-,-,-,\ldots)$ and so on correspond to the second, third, fourth, …sheets, respectively. Normally, at a given energy $s$, only the sheet which can be reached from the physical one by crossing the branch cut from $s+i\epsilon$ to $s-i\epsilon$ between the thresholds $\text{thr}_{n-1}$ and $\text{thr}_{n}$, has a significant impact on physical observables.

As discussed earlier, we have found that the impact of the $D^{*}$ mesons on the $S$-wave $D\phi$ scattering processes is very small. We thus search for poles using the amplitudes without $D^{*}$ derived in Ref. Yao:2015qia . In this way, the complexity of analytically continuing the three- and four-point loops in Fig. 1 is avoided. The physical meson masses and decay constant are employed in the pole searching, and the obtained poles are listed in Table. 4.

For the $(S,I)=(1,0)$ coupled-channel system, in addition to the pole at $\sqrt{s}=2.317$ GeV on the physical sheet, which corresponds to $D_{s0}^{\ast}(2317)$ and was used as a condition to constrain the parameters, using the central values of the parameters we also found a pair of poles with a small but nonvanishing imaginary part on the second Riemann sheet, $\sqrt{s}=(2.439\pm i0.01)$ GeV. The only work which reported an analogous pole is Ref. Guo:2015dha , where a virtual state at $\sqrt{s}=2.356$ GeV below $DK$ threshold on the second Riemann sheet was reported in an NLO calculation including $DK,D_{s}\eta$ and $D_{s}\eta^{\prime}$ channels. We check whether such a pole exists using the parameters from the NLO fits in Ref. Liu:2012zya , and found that only part of the allowed parameter space allows for the pole on the unphysical Riemann sheet. Moreover, the effect of this virtual state pole located at $\sqrt{s}=2.356$ GeV on the physical amplitude is negligible in the NLO calculation, as can be seen from the left column of Fig. 4. However, the pole on the second Riemann sheet in the NNLO calculation can have a non-negligible effect on specific physical amplitudes, as shown in the right column of Fig. 4. These different behaviors are mainly due to different locations of the poles. Nevertheless, we see that the lattice data on the scattering lengths are insufficient to constrain the parameters, and as a result, calculations at different orders may even have a sizeable discrepancy in amplitudes not far from thresholds. More lattice data on $D\phi$ scattering observables are needed to better pin down the LECs.

In addition, we also found a pair of poles $\sqrt{s}=(2.534\pm i0.097)$ GeV on the second Riemann sheet which are not included in Table 4. They have a negligible effect on physical amplitudes and would disappear if the $u$- and $t$-channels are turned off. Likewise, we do not include the following poles in Table 4 since they are located far from the physical region and have little effect: poles at $\sqrt{s}=(2.448\pm i0.049)$ GeV and $(2.267\pm i0.099)$ GeV on the third Riemann sheet for $(S,I)=(1,0)$ and $(1,1)$, respectively; poles at $(2.257\pm i0.018)$ GeV on the second Riemann sheet in the $(-1,0)$ $D\bar{K}\to D\bar{K}$ channel.

It is well-known that the unitarization approach, relying on right-hand unitarity and the on-shell approximation, has the problem of violation of unitarity when the left-hand cut occurs in the on-shell potential. For instance, the left-hand cut in the $K\bar{K}\to K\bar{K}$ amplitude leads a violation of unitarity for the $\pi\pi$ scattering in the $\pi\pi$–$K\bar{K}$ coupled-channel system GomezNicola:2001as ; Dai:2011bs . 777As pointed out by Refs Guerrero:1998ei ; GomezNicola:2001as , the unitarity violation is numerically small in the $\pi\pi$–$K\bar{K}$ case, hence no serious problem was caused there. The same unitarity violation happens to the $D\phi$ scattering with $(S,I)=(0,1/2)$, which has three coupled channels: $D\pi$, $D\eta$ and $D_{s}\bar{K}$. One of the left-hand cuts from the inelastic channel $D_{s}\bar{K}\to D\eta$ amplitude, from ($1.488$ GeV)2 to ($2.318$ GeV)2, overlaps with the right-hand cut starting from the $D\pi$ threshold, which can be verified by the discontinuity across the real axis below the $D\pi$ threshold. Although this left-hand cut is not numerically important, its presence together with other left-hand cuts and right-hand cuts make the whole real axis nonanalytic. Since Eq. (23) was derived using the $N/D$ method neglecting the left-hand cuts, its continuation to the complex plane near the left-hand cut is untrustworthy. As a result, the coupled-channel amplitudes obtained from Eq. (23) do not have the correct analytic properties even in the relevant energy region. Consequently, a pair of pole at $(2.046\pm i0.050)~{}\text{GeV}$ are found on the first Riemann sheet for the coupled-channel $(S,I)=(0,1/2)$ amplitude. As we know, poles on the first Riemann sheet can only be located on the real axis below the lowest threshold, which are associated with bound states. A pole on the first sheet with a nonvanishing imaginary part or above the lowest threshold is inconsistent with causality. The appearance of the pole on the first sheet in the coupled-channel $(S,I)=(0,1/2)$ is due to the existence of the coupled-channel cut. The left-hand cuts stem from the one-loop potentials, and are absent in the NLO cases.

If we consider only the single-channel $D\pi$ for $(S,I)=(0,1/2)$, there is no such a problem as it comes from the left-hand cut of the inelastic channels. We searched for poles in the single-channel amplitude, and found a pair of poles in the second Riemann sheet given in Table 5,888Notice that the poles found in both the single-channel and coupled-channel unitarized NLO amplitudes are similar to each other in Ref. Guo:2009ct . corresponding to the lower pole at $(2.105-i0.102)$ GeV of the two-pole structure of $D_{0}^{*}(2400)$ advocated in Ref. Albaladejo:2016lbb .

In addition, we also investigated the pole movements with varying pion masses. The pion mass dependence trajectories of the poles can provide us with useful information about the properties of the different states, as discussed, e.g., in Ref. Hanhart:2014ssa . The $M_{\pi}$ trajectory for the pole corresponding to $D_{s0}^{*}(2317)$ is plotted in Fig. 5. The pole positions on the first Riemann sheet, which are identified as the pole mass, are shown as the solid line. The dotted line stands for the trajectory of the $DK$ threshold. From Fig. 5, one can see that the $D_{s0}^{*}(2317)$ always stays below the corresponding $DK$ threshold as a bound state for a wide range of $M_{\pi}$. The trajectory of $D_{s0}^{*}(2317)$ is quite similar to the NLO fit result, as shown in Ref. Guo:2015dha .

On the contrary, the pion mass dependence trajectory of the pole around $2.1~{}\text{GeV}$ on the second Riemann sheet for the single-channel $(S,I)=(0,1/2)$ in Table 5 is quite complicated, as shown in Fig. 6. As the value of $M_{\pi}$ increases from $M_{\pi}^{\text{Phy}}$, both the real and imaginary parts of the pole tend to decrease on the second Riemann sheet. At some point around $3.2M_{\pi}^{\text{Phy}}$, the real part of the pole becomes lower than the corresponding $D\pi$ threshold. When $M_{\pi}$ increases to around $3.6M_{\pi}^{\text{Phy}}$, the pair of poles hits the real axis below the threshold and these become two virtual states on the second Riemann sheet. If we further increase the pion mass, one of the virtual poles would move along the real axis away from the threshold, while the other one moves towards the threshold and becomes a bound state of the first Riemann sheet at $3.8M_{\pi}^{\text{Phy}}$. If we keep increasing $M_{\pi}$, both the virtual state and bound state move away from the threshold along the real axis. The behavior of the pole is similar to the corresponding ones in Refs. Guo:2009ct ; Guo:2015dha as well as the pion mass dependence of the $f_{0}(500)$ in Ref. Hanhart:2008mx . It is general for $S$-wave states.

IV Summary and conclusions

We have calculated the potential for the scattering of Goldstone bosons off charmed mesons, including the charmed vector mesons as explicit degrees of freedom, up to the NNLO in a framework of covariant ChPT. We explicitly show that the UV divergences and the so-called power counting breaking terms from the one-loop potentials can be absorbed by a redefinition of the LECs. In the EOMS scheme, we obtained the $D\phi$ scattering potentials possessing the good properties, i.e., they are free of UV divergences and power counting violating terms. In order to describe the $S$-wave scattering lengths at high pion masses and to study the possible dynamically generated resonances that are absent in the Lagrangian, e.g. the $D_{s0}^{*}(2317)$, at relatively high energies, which is the nonperturbative effect, a unitarization procedure is employed.

In order to determine the LECs, we performed fits to scattering lengths in a few channels computed in lattice QCD at various unphysical pion masses. Since the lattice simulations are performed with fixed charm and strange quark masses with varying up and down quark masses, we derived the corresponding pion mass dependence of the scattering lengths by extrapolation of the involved masses and the pion decay constant. For an easy comparison to the previous case without including $D^{*}$ explicitly, we used a similar fit procedure as in Ref. Yao:2015qia . It turns out that, for UChPT-6(b), the current fit result is quite close to the previous one which was done without explicit vector charmed mesons. It is thus a firm conclusion that the $D^{*}$ contribution to the $S$-wave $D\phi$ scattering in the threshold region is negligible.

Based on the small contribution of the $D^{*}$ to the $D\phi$ scattering potentials, we investigated possible dynamically generated resonances using the unitarized scattering amplitudes without explicit $D^{*}$ by analytic continuation. It is worth noticing that a pair of poles with nonvanishing imaginary parts are found on the physical Riemann sheet in the coupled-channel $(S,I)=(0,1/2)$ amplitude, which are at odds with causality. The issue is caused by the coupled-channel left-hand cut, which is not taken into account in the unitarization procedure we used. It may be avoided by using a single-channel potential if we only focus on the region near the $D\pi$ threshold. In the end, we studied the trajectories of the poles corresponding to the $D_{s0}^{*}(2317)$ and a resonance in $(S,I)=(0,1/2)$ channel with varying the pion mass. They exhibit similar behaviors as in the NLO case given by Ref. Guo:2015dha . To summarize, the LECs are badly determined due to the scarcity of available data. Thus to come to firmer conclusions, more lattice data are required, as also concluded from investigations based on the resonance-exchange model and $S$-matrix properties in Ref. Du:2016tgp .

Acknowledgements

We would like to thank Jiunn-Wei Chen, Jambul Gegelia and Zhi-Hui Guo for helpful discussions. MLD acknowledges the warm hospitality of the ITP of CAS where part of this work was done. This work is supported in part by DFG and NSFC through funds provided to the Sino-German CRC 110 “Symmetries and the Emergence of Structure in QCD” (NSFC Grant No. 11621131001, DFG Grant No. TRR110), by NSFC (Grant No. 11647601), by the Thousand Talents Plan for Young Professionals, by the CAS Key Research Program of Frontier Sciences (Grant No. QYZDB-SSW-SYS013), and by the CAS President’s International Fellowship Initiative (PIFI) (Grant No. 2017VMA0025).

Appendix A Renormalization of the LECs within EOMS scheme

In this appendix, we use the following notations for the chiral limit masses of $D$ and $D^{\ast}$: $m_{D}=M_{0}$ and $m_{D^{\ast}}=M_{0}^{\ast}$. Following Ref. Denner:1991kt , the $N$-point one-loop integrals are defined by

[TABLE]

The one-, two- and three-point one-loop scalar integrals are denoted by $A$, $B$ and $C$ as follows:

[TABLE]

A.1 $\beta$-functions

The $\beta$-functions in Eq. (II.2) read

[TABLE]

A.2 Coefficients of finite shifts

In this appendix, we express the EOMS subtractions in terms of the standard loop function. The explicit coefficients of the finite shifts in Eq. (II.2) are

[TABLE]

[TABLE]

Here the involved scalar one-loop integrals stand for their finite parts only, which are obtained from the original ones, defined in Eq. (A), by performing the $\overline{MS}-1$ subtraction.