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

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.


I. INTRODUCTION
In the past two decades, many excited charmed states have been observed experimentally [1][2][3][4] 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 [5]. The conventional quark-potential models provide a successful description of most of those low-lying excitations, see Ref. [6] 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 π 0 invariant mass distribution and later confirmed by Belle and CLEO Collaborations [2][3][4]. It couples to the DK channel and decays mainly into the isospin breaking channel D s π 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 [7]. Were this interpretation true, one can learn much about the interaction between the D/D s mesons and π/K mesons from studying the D * s0 (2317) and related states. This path has been followed in Refs. [8][9][10][11][12][13][14][15] where the S-wave interaction between charmed D mesons and Goldstone bosons (denoted as φ hereafter) has been studied systematically up to the nextto-leading order (NLO) using chiral perturbation theory (ChPT) for heavy mesons [16][17][18] in combination with a unitarization procedure such as the one in Ref. [19].
In the meantime, significant progress has also been made in lattice QCD [20,21]. Using the Lüscher formalism and its extension to coupled channels (for early works on this topic, see e.g. Refs. [22,23]), scattering lengths and recently phase shifts for the Dφ interaction have been calculated at unphysical quark masses [24][25][26][27][28][29]. The first calculation only concerns the channels free of disconnected Wick contractions [24,25], i.e., Dπ with isospin I = 3/2, DK with I = 0, 1, D s K and D s π. The channels with disconnected Wick contractions such as Dπ with I = 1/2 and DK with I = 0 were calculated later in Refs. [26][27][28]. On the one hand, the lattice results can be used to determine the low-energy constants (LECs) in the chiral Lagrangian [25,[30][31][32][33]. 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. [25], 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. [34], and a remarkable agreement with the lattice results reported in Ref. [29] 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 [35] in fact corresponds to two states with poles located at 2105 +6 −8 − i 102 +10 −12 MeV and 2451 +36 −26 − i 134 +7 −8 MeV, respectively [34], similar to the well-known two-pole scenario of the Λ(1405) [19]. 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 [36]: all loop graphs containing internal matter field propagators start contributing at O (p 2 ). 1 Various approaches have been proposed to remedy the PCB issue, e.g. heavy baryon ChPT [37,38], infrared regularization [39], and the extended-on-mass-shell (EOMS) scheme [40]. 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 πN scattering up to O(p 3 ) [41] and O(p 4 ) [42], and up to leading one-loop order in the presence of ∆-resonance [43]. The first aim of this paper is to present a full calculation of Dφ 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φ scattering to one-loop was made in Ref. [44] in the framework of nonrelativistic heavy meson ChPT, which neglects sizeable recoil corrections [45]. 2 The first one-loop calculation within the covariant formalism was given in Ref. [45]. However, in that paper the NNLO contact terms of Dφ are not included and the kinetic term for the D * 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 [32], the Dφ 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 * 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φ 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 * -less case has been accomplished formally using the path integral formalism in Refs. [47,48]. 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 [19,32,49]. 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 1 The 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 [36]. 2 Such recoil corrections can be restored by using the extended heavy baryon propagator i/(v · k + k 2 /2m) UV-part and EOMS subtractions of the LECs are collected in Appendix A.

A. 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 O (p). However, the nonvanishing masses of the D and D * in the chiral limit introduce new scales M 0 and M * 0 , both counted as O (1). As a result, at low energies, the temporal components of the momenta of the D and D * are counted as O (1), while the spatial components are counted as O (p). Therefore, the virtuality q 2 − M ( * )2 0 in the propagators scales as O (p), and the propagators scale as O (p −1 ). The Goldstone boson propagators are counted as O (p −2 ) 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φ potentials up to leading one-loop order can be written as with the superscripts specifying the chiral dimension. The needed terms in the Goldstone sector read [50] where the trace in flavor space is denoted by · · · , F 0 is the pion decay constant in the chiral limit, and L 4,5 are LECs. Furthermore, χ = 2B 0 diag(m u , m d , m s ), with B 0 a constant related to the quark condensate, and The terms corresponding to interactions between the D = (D 0 , D + , D + s ) mesons and the Goldstone bosons are given by [12,[16][17][18]32] where h i and g j are LECs and the chiral building blocks are given by The covariant derivative is defined via and The socalled chiral connection in the covariant derivatives is defined as Γ µ = u † ∂ µ u + u∂ µ u † /2. Similarly, the relevant terms for the interaction between the D * and the Goldstone bosons are [16][17][18] L (1) withh 0,1 analogous to h 0,1 and F µν = (D µ D * ν − D ν D * µ ). Finally, the LO axial coupling has the form L (1) As pointed out in Refs. [43,51], the resonance-exchange contributions of O(p 2 ) and 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 O (p 2 ) and O (p 3 ) terms for the D * Dφ coupling.

B. 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 )φ 1 (p 2 ) → D 2 (p 3 )φ 2 (p 4 ) can be written as As usual, the Mandelstam variables are defined by s = (p 1 + p 2 ) 2 and t = ( The potentials at tree-level are given by where the involved coefficients corresponding to various scattering processes are shown in Table I. The functions in the D * -exchange potentials read The functions in the NLO potentials read while the one in the NNLO potentials is As for the one-loop potentials at NNLO, the parts without explicit D * mesons can be found in the appendix of Ref. [32] and the ones involving explicit D * states are too lengthy to be shown here. Note that V (Loop) 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 where 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β h i andβ g 0 are given in Appendix A 2.

C. 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 I by crossing and isospin symmetry, see Refs. [13,32] 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 ππ potentials are strictly restored within a range of energies [52], 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 [53,54]. 4 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 [53,54,56].
Before unitarization, the partial wave projection to a definite orbital angular momentum l should be performed where θ is the scattering angle between the incoming and outgoing particles in the centerof-mass frame, and the Mandelstam variable t is expressed as where λ(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 [19] T where the function G(s) encodes the two-body right-hand cut and is given by the two-point loop function The explicit expression for G(s) reads [56] where a(µ) is a subtraction constant with µ the renormalization scale and Note that the logarithmic scale dependence can be absorbed into the subtraction constance a(µ) and we do not distinguish the scale µ 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 [19,49] Eq. (23) is an algebraic approximation of the standard N/D method [56], and it should be understood in the matrix form for coupled-channels, for which G(s) = diag{G i (s)}, with i the channel index.

III. NUMERICAL ANALYSES
A. 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 [25][26][27][28][29] 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 a (S,I) 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: [25] whereM K ,M D andM Ds denote the masses in the two-flavor chiral limit (M 2 π (∝m) → 0 but with the fixed strange quark mass m s ). They have the form Using Eqs. (28) and (29), one gets Here,M D * andM D * s denote the corresponding two-flavor chiral limit masses, which can be estimated by the relations with M Phy D * and M Phy D * s denoting the corresponding physical masses, 2.008 GeV and 2.112 GeV, respectively. One hash 1 = h 1 andh 0 = h 0 in the heavy quark limit. 5 These relations as well as similar relations for other LECs will be employed in order to reduce the number of parameters. The DD * π axial coupling constant g 0 can be fixed by the decay width Γ D * + →D 0 π + . As discussed in Refs. [32], one gets g = (1.113±0.147) GeV for the renormalized coupling g, which contains the bare constants g 0 and one-loop chiral corrections. At the oneloop level, the pion decay constant has the form [50] In this paper, the scattering length data we use are taken from Ref. [25] and Refs. [26][27][28], respectively. The lattice simulation data for M K , M D , M Ds and F π , F K are taken from Refs. [25] and [58], which share the same ensembles (M007, M010, M020 and M030) with   [25]. The masses and decay constant in the chiral limit are in units of GeV. h 0 , h 1 , L r 4 and L r 5 are dimensionless. The asterisk indicates an input value.
Ref. [25]  Ref. [25]. The fit results are listed in the Table 3 in Ref. [32]. In addition, we also include the DK scattering length with (S, I) = (1, 0) at M π = 0.156 GeV [27], as discussed in Ref. [32]. 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 forM K ,M D andM Ds , as listed in Table II. Furthermore, in order to reduce the correlations between the LECs, we introduce the following redefinitions of the LECs [25,32] [25].
To compare with the result of Ref. [32] 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. [25], including 5 channels at pion masses 0.301 GeV, 0.364 GeV, 0.511 GeV and 0.617 GeV, as well as the isoscalar DK channel at the pion mass 0.156 GeV in Ref. [27], 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 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 π = 0.617 GeV. Results of both fits are listed in Table III.
In both fits, as in Ref. [32], the absolute value of the dimensionless LEC h ′ 5 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 χ 2 where χ 2 is the standard chi-squared, and χ 2 prior 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 ′ ) and UChPT-6(b ′ ) are obtained by minimizing χ 2 aug instead of χ 2 using the same data as in UChPT-6(a) and UChPT-6(b). Although the values of LECs are more natural, the χ 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. [32].
Compared to the NLO fits [13,33], the NNLO fits have larger χ 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 UChPT-6(a) UChPT-6(b) UChPT-6(a ′ ) UChPT-6(b ′ ) Among various fits, UChPT-6(b) has the smallest χ 2 , which is also true for the previous fits without D * [32]. 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 * on the quantities in question is small. 6 It is due to the nonvanishingM K in loops contributing to the D s π → D s π potential in the SU (2) chiral limit. Near the D s π 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 (1,1) Dsπ→Dsπ is negligible in the SU (2) chiral limit [25,33]. The results of the UChPT-6(b) fits to the lattice data of the scattering lengths. The filled circles are lattice results in Ref. [25], and the filled square (not included in the fits because it refers to N f = 2) and diamond are taken from Ref. [27].

B. Dynamically generated resonances
The unitary S-matrix could have poles in the complex energy ( √ 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 [53]  The results of the UChPT-6(b ′ ) fits to the lattice data of the scattering lengths. The filled circles are lattice results in Ref. [25], and the filled square (not included in the fits because it refers to N f = 2) and diamond are taken from Ref. [27].
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 (+, +, +, . . .), while (−, +, +, . . .), (−, −, +, . . .), (−, −, −, . . .) 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ǫ to s − iǫ between the thresholds thr n−1 and 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φ scattering processes is very small. We thus search for poles using the amplitudes without D * derived in Ref. [32]. In this way, the complexity of analytically continuing the threeand 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. IV. For the (S, I) = (1, 0) coupled-channel system, in addition to the pole at √ s = 2.317 GeV on the physical sheet, which corresponds to D * s0 (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, √ s = (2.439 ± i0.01) GeV. The only work which reported an analogous pole is Ref. [33], where a virtual state at √ s = 2.356 GeV below DK threshold on the second Riemann sheet was reported in an NLO calculation including DK, D s η and D s η ′ channels. We check whether such a pole exists using the parameters from the NLO fits in Ref. [25], 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 √ 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 nonnegligible 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φ scattering observables are needed to better pin down the LECs.
In addition, we also found a pair of poles √ s = (2.534 ± i0.097) GeV on the second Riemann sheet which are not included in Table IV. 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 IV since  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 KK → KK amplitude leads a violation of unitarity for the ππ scattering in the ππ-KK coupled-channel system [59,60]. 7 The same unitarity violation happens to the Dφ scattering with (S, I) = (0, 1/2), which has three coupled channels: Dπ, Dη and D sK . One of the left-hand cuts from the inelastic channel D sK → Dη amplitude, from (1.488 GeV) 2 to (2.318 GeV) 2 , overlaps with the right-hand cut starting from the Dπ threshold, which can be verified by the discontinuity across the real axis below the Dπ 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 ± i0.050) GeV are found on the first Riemann sheet for the coupledchannel (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π 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 V, 8 corresponding to the lower pole at (2.105 − i0.102) GeV of the two-pole structure of D * 0 (2400) advocated in Ref. [34]. 7 As pointed out by Refs [59,61], the unitarity violation is numerically small in the ππ-KK case, hence no serious problem was caused there. 8 Notice that the poles found in both the single-channel and coupled-channel unitarized NLO amplitudes are similar to each other in Ref. [13].   Table 4 of Ref. [32]. 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.  Table 4 of Ref. [32]. Physical masses are used to obtain the poles. The Riemann sheets on which the poles are located are indicated in the last column. 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. [62]. The M π 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 π . The trajectory of D * s0 (2317) is quite similar to the NLO fit result, as shown in Ref. [33]. On the contrary, the pion mass dependence trajectory of the pole around 2.1 GeV on the second Riemann sheet for the single-channel (S, I) = (0, 1/2) in Table V is quite complicated, as shown in Fig. 6. As the value of M π increases from M Phy π , both the real and imaginary parts of the pole tend to decrease on the second Riemann sheet. At some point around 3.2M Phy π , the real part of the pole becomes lower than the corresponding Dπ threshold. When M π increases to around 3.6M 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 Phy π . If we keep increasing M π , 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. [13,33] as well as the pion mass dependence of the f 0 (500) in Ref. [63]. 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φ 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. [32]. 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φ scattering in the threshold region is negligible.
Based on the small contribution of the D * to the Dφ 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 coupledchannel (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π 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. [33]. 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. [64].

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. (20) arē Here the involved scalar one-loop integrals stand for their finite parts only, which are obtained from the original ones, defined in Eq. (A1), by performing the MS − 1 subtraction.