Can the two-pole structure of the D ∗ 0 ( 2300 ) be understood from recent lattice data?

It was demonstrated in a series of papers employing unitarized chiral perturbation theory that the phenomenology of the scalar open-charm state, the D ∗ 0 ( 2300 ) , can be understood as the interplay of two poles, corresponding to two scalar-isospin doublet states with different SU(3) ﬂavor content. Within this formalism the lightest open charm positive parity states emerge as being dynamically generated from the scattering of the Goldstone-boson octet off D mesons, a picture that at the same time solves various problems that the experimental observations posed. However, in recent lattice studies of D π scattering at different pion masses only one pole was reported in the D ∗ 0 channel, while it was not possible to extract reliable parameters of a second pole from the lattice data.


Introduction
The discovery of the charm-strange mesons D * s0 (2317) [1] and D s1 (2460) [2] with masses significantly lower than a e-mail: a.asokan@fz-juelich.de(corresponding author) b e-mail: tangmengna@mail.itp.ac.cn c e-mail: fkguo@itp.ac.cn d e-mail: c.hanhart@fz-juelich.de e e-mail: kamiya@hiskp.uni-bonn.def e-mail: meissner@hiskp.uni-bonn.dethe predictions for the lowest-lying scalar and axial-vector cs mesons from the quark model (see, e.g., Ref. [3]) in 2003 lead to intensive discussions on their nature.Closely related to these two hadrons, there were observations of broad bumps in the Dπ and D * π invariant mass distributions in B decays by BaBar, Belle and LHCb Collaborations [4][5][6][7].The bumps were fitted using a Breit-Wigner (BW) parametrization with energy-dependent widths, assuming the existence of one broad scalar (axial-vector) resonance coupled to D ( * ) π ; accordingly, such fits led to the D * 0 (2300) and D 1 (2430) entries that are listed in the Review of Particle Physics (RPP) [8]: where all numbers are given in units of MeV.However, the use of a BW form is not justified in these cases as constraints from chiral symmetry and coupled channel effects are not taken into account, see e.g.[9].Those are automatically built into unitarized chiral perturbation theory (UChPT), where all calculations find two D * 0 mesons and two D 1 mesons in the same energy region as the D * 0 (2300) and D 1 (2430), respectively [10][11][12][13][14][15][16][17]. 1 All these works tell a qualitatively coherent story, although e.g. the role of left-hands 1 What is meant here is that there are two states coupling to the π D ( * ) channel predominantly in S-wave, instead of only one D * 0 (2300)/D 1 (2430).Clearly, in addition to these states there is also the narrow D 1 (2420), decaying into π D * predominantly in D-wave (up to heavy quark spin symmetry violating contributions) which has a width of about 30 MeV. 0123456789().: V,-vol 123 cuts still needs be agreed upon [17,18].For instance, the parameters in the UChPT amplitude used in Refs.[14,15] are fixed from fitting to the results of a set of S-wave charmed-meson-light-pseudoscalar-meson (D ) scattering lengths computed using lattice quantum chromodynamics (QCD) [19].The two D * 0 poles in the UChPT amplitude of Refs.[14,15,19] are located at 2105 +6 −8 − i102 +10 −11 MeV and 2451 +35 −26 − i134 +7 −8 MeV.And it was demonstrated in Refs.[9,15] that the amplitudes are consistent with the LHCb data of the angular moment distributions from three-body B meson decays: [22], and B 0 → D0 π − K + [7].For a review on two-pole structures in QCD, see [23].
In seeming disagreement to these findings, the lattice QCD analysis of the Dπ -Dη-D s K coupled channel system by the Hadron Spectrum Collaboration (HadSpec) in Ref. [24] reported only one D * 0 state just below the Dπ threshold, with the pion mass of about 391 MeV.In this paper, we will discuss whether the higher D * 0 pole is consistent with the lattice data, and propose a K -matrix formalism constrained with the SU(3) flavor symmetry that can be used in analyzing coupled-channel lattice data.

Analysis of the amplitude from the lattice study
In Ref. [24] lattice data for the strangeness zero, isospin-1/2 channel at a pion mass of about 391 MeV were presented and analyzed with a sizable set of K -matrix parametrizations of the kind where i and j label the different reaction channels and m, g i and γ (n) i j are real parameters to be determined in the fit to the lattice data.From this, the T -matrix for the S-wave coupled-channel (Dπ -Dη-D s K ) scattering is given by with T K (s) defined as where the second term on the right-hand side contains the Chew-Mandelstam function, subtracted at the K -matrix pole parameter m.It is given by with where m (i) 1 and m (i) 2 are the masses of the two particles in channel i and s is the centre-of-mass (c.m.) energy squared.The imaginary part of T −1 K (s) i j is then given by the phasespace factor 2 , which automatically ensures the unitarity of the S-matrix.
The nine parametrizations presented in Ref. [24] differed by the set of parameters that was allowed to vary in the course of the fit.The parameters present in the different amplitudes along with their reduced χ 2 values from energy level fits performed in Ref. [24] are given in Table 1.

Pole search
The T -matrix is analytic over the whole complex energy plane except for poles and branch cuts along the real axis due to kinematic (right-hand cuts) and dynamic singularities (left-hand cuts).Dynamic singularities (left-hand cuts) are associated with the interactions in the crossed channels.Since those are usually distant, one assumes that their effect can be captured by polynomial terms allowed in the parametrization of the K -matrix used.Right-hand cuts start from branch points that appear whenever a channel opens.Accordingly, at each threshold the number of Riemann sheets of the complex energy (or s) plane gets doubled.Thus, the three-channel case studied here leads to eight Riemann sheets.The sheets are labeled as shown in Table 2, where the thresholds are arranged with increasing energies 1 = Dπ , 2 = Dη and 3 = D s K .For illustration we show in Fig. 1 the analogous labeling for two channels.See Fig. 3 of Ref. [25] for the three-channel case.
The poles correspond to bound states or resonances depending on their location on the Riemann sheets.Bound states correspond to poles on the physical sheet below the lowest threshold energy and resonances are poles in the complex plane of the unphysical sheets (in addition there are virtual state poles, located on the real axis of unphysical sheets, but those do not play a role in this work).The poles on the sheets closest to the physical sheet have the strongest influence on the scattering amplitude.In the current notation sheets RS211, RS221, and RS222 would be directly connected to the physical sheet, i.e., RS111, above the respective thresholds (c.f.Fig. 1).The poles of the T -matrix are given by the zeroes of the determinant of the matrix in Eq. (4), i.e., det Table 1 The parametrizations used in the analysis in Ref. [24].The check mark denotes a free parameter and "-" implies the parameter is fixed to zero.The channels are denoted with increasing threshold energies, with 1 = Dπ , 2 = Dη and 3   The unphysical sheets can be accessed by adding the discontinuity across the branch cut to Eq. (4).Via the Schwarz reflection principle the discontinuity across the branch cut is related to the imaginary part of the amplitude by where Im[T K (s+i )] needs to be understood as the analytic continuation of the imaginary part of the amplitude on the real axis above threshold.
Crossing from the physical sheet (RS111) to any sheet can be done by where the subscript X stands for the sheet number and Disc is a 3 × 3 matrix containing the relevant discontinuities needed for the sheet transition, e.g. for the transition from RS111 to RS211 we employ and for RS111 to RS221 This prescription is straightforwardly generalized to arbitrary transitions between sheets.
Table 3 The pole locations from amplitude parametrizations of Ref. [24], in units of MeV.Empty slots denote that a pole was not found within the search range for the particular parametrization on the cor-responding sheet.In the last line the results for the UChPT amplitude employed in Ref. [14]  At a pion mass of about 391 MeV, the lowest pole in the studied channel turns out to be a bound state, accordingly located on sheet RS111 [24]; the same conclusion was reached in UChPT in Ref. [14].This pole was found in the fits of all 9 parametrizations employed by the Hadron Spectrum Collaboration [24].At the same time, additional poles were found on sheets RS211, RS221, and RS222.These additional poles were found for almost all amplitude paramterizations employed in Ref. [24], which were, however, not reported in the publication since they not only scatter very much, but also are in parts located outside the energy region where the fit was performed.Table 3 shows the pole values found from the search with the corresponding sheets from the different amplitude parametrizations.The 1σ uncertainties of the pole values were calculated by the bootstrap method.
Graphically the poles on RS221 are displayed in Fig. 2. In the following we focus the discussion on this sheet, since this is the one where the UChPT amplitude has its most prominent higher D * 0 pole at physical [15] as well as the unphysical meson masses employed in the lattice study [14].The plots of the pole locations of the higher pole for the different parametrizations on the other Riemann sheets that connect closely to the physical axis (RS211 and RS222) are shown in the Appendix.Table 4 gives the location of the corresponding two particle thresholds.
Figures 2 and 13 in the Appendix and Table 3 clearly show two important features of the poles extracted from different parametrizations: (i) There is a significant correlation between real part and imaginary part of the poles, and the location of the pole extracted from the UChPT analysis is in line with that correlation.(ii) All poles are located on hidden sheets, which are the sheets that are not directly connected to the physical sheet.For example, the RS221 poles are well Fig. 2 The location of poles on sheet RS221 on the complex energy plane.The x-axis and y-axis show the real and imaginary part of energy, respectively.The poles from the amplitude parametrizations employed in Ref. [24] are shown in yellow.The pole from the UChPT amplitude [19] is shown in green [14].The vertical green and blue dashed lines represent the Dη and D s K thresholds, respectively.The error bars show the 1σ statistical uncertainty above the D s K threshold.Thus they are all shielded by the RS222 sheet and their effect on the amplitude can hardly be seen above the D s K threshold.As we discuss in the following, both features together guide one to an understanding that there indeed needs to be a second pole in an amplitude that describes the lattice data and that it is natural that the original analysis performed on the lattice data lead to badly constrained pole locations.The mechanism underlying this is that the distance from the threshold is overcome by an enhanced residue.This mechanism, also reported e.g. for the case of the f 0 (980) and a 0 (980), was observed before as a general feature of Flatté amplitudes [26].
Please note that the parameters M and , derived from the pole location, agree to those found e.g. in the BW fits only for narrow, isolated resonances-for details see the review on resonances in Ref. [8].Equally fundamental resonance properties are provided by the pole residues.A pole-residue quantifies the couplings of the resonance to the various channels.The residues of a pole located at s = s p are defined as The residues can be easily obtained using the L'Hôpital rule to compute the limit: Since the residues factorize according to R 2 i j = R ii R j j one can define an effective coupling via which has dimension [mass].The index r is meant to distinguish the residues from the parameters g i that appear in the K -matrix in Eq. ( 2).The couplings g r i characterize the transition strengths of the resonance to the channel.Those residues can also be extracted from production reactions and are independent of how the resonance was produced.
Since the poles of interest here are hidden, their effect on the physical axis is visible only at the thresholds irrespective of their exact pole locations.Moreover, the visible effect in the amplitude on the physical axis from a pole on a hidden sheet close to the threshold with a small residue is in fact hardly distinguishable from a faraway pole with a large residue.We regard this ambiguity as the most natural explanation for the large spread in the pole locations found in the analysis of the lattice data reported above.
To test this hypothesis, we now study the strengths of the residues as functions of the distance of the poles to the threshold.Clearly, there is some ambiguity in how to quantify the distance from the threshold.Since the channel couplings also drive the size of the imaginary part of the pole location and we want to avoid counting the effect of those couplings twice, we choose instead of (M − M thr. ) 2 + ( /2) 2 , which might appear more natural on the first glance, as a measure for the distance of the pole to the threshold.The pole mass M that appears above was introduced in Eq. ( 13) and M thr.denotes the threshold location relevant for the given sheet, e.g. in case of RS221 we have M thr.= M K + M D s .Table 5 shows the distance of the RS221 pole from the D s K threshold along with the square root of the absolute value of the residue to the three channels.The graphical representation for the Dπ -Dπ channel is shown in Fig. 3

(left).
A straight line is fit to the data to extrapolate the values at the threshold.The fitting was done using the MINUIT algorithm [27] from the iminuit interface [28,29].The uncertainties of the fit parameters quoted are from the MIGRAD routine of MINUIT.The 1σ statistical uncertainty of the fitted line was calculated using the bootstrap technique.From the straight Table 5 Distance of the RS221 pole from the D s K threshold and the square root of its residue to the respective channels in MeV for the amplitude parametrizations obtained in Ref. [24].In the last line the results for the UChPT amplitude employed in Ref. [14] are given for comparison

Amplitudes
Dist. from line fit the y-intercept was found to be at (5.8 ± 0.9) GeV.
The graphical representation for corresponding distances and residues for the Dη-Dη and the D s K -D s K channel is given in Fig. 3 (middle) and (right).The fit to the Dη (D s K ) residues provides an intercept of (2.4 ± 0.8) GeV ((4.6 ± 0.9) GeV).The corresponding results for the poles on sheets RS211 and RS222 are shown in the Appendix.While the linear fits shown do not provide an excellent representation of the extracted data for the different channels, they illustrate nicely that there is indeed a significant correlation between the distance of the poles from the threshold and the residues.In addition, the effect of the poles at threshold, encoded in the y intercepts deduced from the fits, is rather well constrained by the fits.We interpret this observation such that the lattice data require not only one bound state pole but also a higher pole as was also found in the various studies employing UChPT [10][11][12][13][14][15]17].
An interesting question is, if it is possible to come up with a parametrization to be used in the K -matrix fit that constrains better the pole location of the higher pole, with inputs of approximate symmetries of QCD.This will be the focus of the next section.

SU(3) symmetry
The parametrization dependence of the higher pole location calls for a stronger constrained amplitude.In the following we present a prescription of the K -matrix consistent with the SU(3) flavor symmetry.In the resulting scattering matrix SU(3) breaking comes only from the Chew-Mandelstam functions introduced in Eq. ( 5).Clearly, for the physical pion mass and low energies such a treatment is not justified, since the leading order chiral interaction scales with the energies of pion and kaons for Dπ and D s K scattering, respectively, which induces a sizeable SU(3) flavor breaking into the scattering potential.However, in this study we work at higher 4 Weight diagrams of the [15], [6] and [ 3] representations pion masses which leads to a much smaller pion-kaon mass difference.Moreover, we are mainly interested in the higher mass range, where the second pole is located.Under such circumstances the leading SU(3) breaking effect is induced by the loop functions which bring the cut structure to the amplitudes.
The flavor structure of the D interaction can be written as a direct product of an anti-triplet for the charmed mesons and an octet for the light pseudoscalar mesons.The direct product can be decomposed into a direct sum of the [ 3], [6] and [15] irreducible representations.Figure 4 shows the multiplet structure.
The SU(3) flavor basis and the isospin-symmetric particle basis are related via ⎛ where From the rotation matrix U , we may read off the following expressions for the SU(3) symmetric coupling structure of the Dπ , Dη and D s K (I = 1 2 , S = 0) coupled-channel system: The form of the K -matrix assuming the existence of two bare poles, in contrast to the one used in Ref. [24] which contains only one bare pole, reads Here the two bare poles are assumed to be in the two SU(3) multiplets with S-wave attractions from the leading order chiral dynamics [14].We have seven free parameters in total, g α , c α and m α .The overall factors in Eq. ( 22) are absorbed into the parameters of g α , c α .If there was no SU(3) constraint, a K -matrix with the same number of bare poles would contain 3 more parameters (a constant K matrix is symmetric and thus contains 6 parameters, instead of 3 c α 's here).With Eqs. ( 4) and ( 5), the T -matrix T K (s) can be calculated in the same way as in Sect. 2. The subtraction point for the Chew-Mandelstam function of the fits is chosen identical to the parameter m 3.

Fitting to lattice energy levels
Here we employ the flavor SU(3) constrained K -matrix to fit the lattice energy levels in the Dπ c.m. frame obtained in Ref. [24].To this end, we need to relate the T -matrix defined with the K -matrix in the continuum system and the energy levels in the finite volume system.In this study, we employ the scheme based on the effective field theory framework developed in Ref. [30].Below we briefly summarize this scheme.
With the Lippmann-Schwinger equation, the T -matrix in the continuum can be written as where V (s) is the interaction matrix and G(s) is the diagonal matrix of the scalar two-meson loop functions [31].With the momentum cut-off regularisation G(s) is given by where q max is the cut-off momentum and m 1,2 are the masses of the two particles in channel i. Similarly to Eq. ( 24), the T -matrix in a finite volume system T satisfies where Ṽ (s) and G(s) are the interaction and loop function in the finite volume system, respectively.With the spatial extension L of the cubic box, G(s) is given as Since V (s) is equal to Ṽ (s) up to exponentially suppressed corrections, the T -matrix in the finite volume system T is related to the T -matrix in the infinite volume T by with The lattice energy levels, which we need to fit, correspond to the zeroes of the determinant of T −1 provided in Eq. (30).
Table 6 The best fit values arrived in the different fits to the lattice energy levels of Ref. [24] based on Eq. ( 23), along with their χ 2 /dof.The symbol '-' is used for parameters set to zero (or absent) in the particular fit Fig. 5 Comparison of the energy levels from data and the fits.The energy levels from Ref. [24] are shown in yellow.The energy levels used as input in the fits are shown as circles and those not used is shown as crosses.The energy levels from Fit 1_4L, Fit 2_4L, Fit 3_4L and Fit 4_4L are shown in red, blue, black and light blue, respectively.The solid(dashed) red, green and blue lines show the Dπ , Dη and D s K non-interacting energy levels(thresholds) respectively

Results for the fits employing SU(3) constraints
We performed four different fits to the rest frame lattice energy levels in Ref. [24] using the MINUIT algorithm [27] with the Julia interface to the iminuit package [28,29]: -in Fit 1_4L all the parameters in Eq. ( 23) are included; -in Fit 2_4L we fix c 3 = 0 and c 6 = 0; -in Fit 3_4L we fix c 6 = 0; -in Fit 4_4L we fix g 6 = 0 to omit the explicit pole term of [6] (and thus m 6 is absent).
From every volume we use the lowest four energy levels (thus the addition _4L to the fit names) of the [000] A + 1 irreducible representation, 2 where the S-wave component gives the dominant contribution.In the next subsection we discuss the fit results for Fit 4_All that was performed including 2 We did not implement the discretization of our amplitude for moving frame data, since this is technically a lot more demanding and the usefulness of the SU(3) symmetry constraint can already be demonstrated with the rest frame fits.all the rest frame lattice levels.The obtained parameters for the different fits as well as the χ 2 values found are listed in Table 6. Figure 5 shows the energy levels obtained from the fits together with the data points in the lattice rest frame.The uncertainties of the fit parameters quoted are from the MIGRAD routine of MINUIT.Further, using the parameters from the fit, the poles of the T -matrix in the continuum on the different Riemann sheets are extracted.The resulting pole positions can be found in Table 7.It turns out that, contrary to the pole extraction employing Eq. ( 2), now in all fits the higher mass pole has a mass of about 2.5 GeV and is thus located close the Dη and D s K thresholds (see Fig. 6) A more detailed comparison of the performance of the fits shows that for [ 3] both a pole term and a constant term are needed to obtain an acceptable fit.Moreover, the uncertainties for the pole parameters that emerged from Fit 1_4L are a lot larger than for the other fits (in addition the fit even allows for an additional level very close to the fitting range).We therefore exclude both Fit 1_4L and Fit 2_4L from further discussions.On the other hand, fits of comparable quality emerge, if either the constant term in the [6] (Fit 3_4L) or the pole term in the [6] (Fit 4_4L) is abandoned.In the latter case the higher pole is generated via the unitarization.Note that in our case the number of parameters connected to the bound state pole is in any case 2 (one coupling constant and a mass), while in case of the fits performed by the Hadron Spectrum Collaboration this number is 4, for there an individual coupling is needed for each channel.The difference in the number of parameters needed for the [ 3] and the [6] channels can be understood straightforwardly from the observation that the pole in the [ 3] is a bound state and further parameters are needed to obtain a decent fit of the additional energy levels.The pole originated in the [6], on the other hand, sits rather high up in the spectrum, having a larger imaginary part, and thus naturally controls all energy levels that have a sizable contribution from this representation; either the bare pole or the constant term in [6] provides a seed for the [6] pole.
Table 7 also shows that in all fits poles appear on RS222 above the D s K threshold, however, with large uncertainties especially on the mass parameters.Since in this energy range RS222 connects directly to the physical sheet, these poles Fig. 6 Locations of the RS221 poles from the different fits based on Eq. ( 23), together with the pole reported in UChPT amplitude [14] in green, and the various extractions from the amplitudes extracted in Ref. [24] in yellow.The pole locations from Fit 1_4L, Fit 2_4L, Fit 3_4L and Fit 4_4L are shown in red, blue, black and light blue respectively.The green and blue vertical dashed lines represent the Dη and D s K thresholds, respectively show up as peaks in the amplitudes at high energies.Note that analogous poles were also present in the fits performed in the course of the analysis of Ref. [24] and in the UChPT amplitude (see Table 3), however, they appeared at significantly higher energies.In Fit 1_4L this pole can appear quite close to the energy region of interest, given the large uncertainty in the mass parameter.We interpret this phenomenon as reflecting a too large number of parameters in the fit.In the two best fits, namely Fit 3_4L and Fit 4_4L, on the other hand, the poles on RS222 are typically located deep inside the complex plane or rather high above the threshold, respectively, although within uncertainties it can appear rather close to threshold also for Fit 3_4L, with leads to the strong rise very close to the higher thresholds.The large spread in the amplitudes above the D s K threshold visible in Fig. 8 reflects the bad determination of the highest pole from the lattice data included in the fits.
A comparison of the RS111 pole locations from the SU(3) fits just reported and that the UChPT amplitude [14] and Ref. [24] is shown in Fig. 7.As expected the location of this bound state pole is consistent amongst all extractions.Fig. 7 Locations of the RS111 poles from the different fits together with the RS111 pole in the UChPT amplitude [14] and from HadSpec [24] Note that in all fits the constant term in the [15] representation turns out to be repulsive, in line with the expectations from leading order chiral perturbation theory.This is a nice and in fact non-trivial confirmation of the hypothesis that, even at pion masses as high as 391 MeV, already leading order chiral perturbation theory provides valuable guidance for the physics that leads to the emergence or not appearance of hadronic molecules.
We also tested if we can fit the lattice data when replacing the pole term in the [6] representation by a pole in the [15].Those fits, however, did not converge and are therefore not reported in the figures and tables.The amplitudes arrived from the fits are shown in Fig. 8, there, however, at  17.8 +0.8 −0.9 significantly higher energies.The appearance of these poles in a direct consequence of the K -matrix parametrization employed.Moreover, in Fit 1_4L and Fit 2_4L those poles are rather close to the threshold To complete the discussion, in Table 8 we show the square root of the absolute values of the RS221 pole residues to the respective channels derived from the fits.Further, Fig. 9 shows the distance of the RS221 pole to the D s K thresh-old, as defined in Eq. ( 17), versus the effective coupling of the pole to the Dπ , Dη and D s K channels, respectively.Besides Fit 4_4L, all values are statistically consistent with the intercepts arrived at in Sect.2.2 within errors.Our fit amplitudes are formulated as a momentum expansion.Accordingly, to find our main results, we performed fits with including only the lowest 4 lattice levels at each volume.However, to check for stability of our findings, we also performed fits with Fit 3 and Fit 4 to all rest frame levels ([000] A + 1 irrep.)-theresulting fits are labeled as Fit 3_ALL and Fit 4_All.It turns out that the former parametrization does not allow for a decent fit (the best χ 2 /dof we can achieve is 5.5).It calls for introducing additional bare poles or/and momentum-dependent contact terms.
In the left panel of Fig. 10 we show a comparison of the original fit, Fit 4_4L, with the full spectrum, in the right panel the fit results for Fit 4_All, where the higher levels are included in the fit.The parameters arrived from the new Fig. 12 Comparison of locations of the RS221 poles.The pole for Fit 4_All is shown in red, together with the pole reported in UChPT amplitude [14] in green, and the various extractions from the amplitudes extracted in Ref. [24] in yellow.The pole location from Fit 4_4L, when including only the lower energy levels, is shown in light blue.The green and blue vertical dashed lines represent the Dη and D s K thresholds, respectively fit are also shown in Table 6.Clearly, the fit is not excellent, however, when compared to only the rest frame levels, amplitude 4 and amplitude 6 from Ref. [24] reach χ 2 values of 36 and 24, respectively, and are thus of similar quality.The amplitude plots arrived from Fit 4_All are shown in Fig. 11.Though the bound state pole does not change, there is a change in the location of the higher pole in sheet RS221, which however has a real part similar to that from Fit 4_4L.The previous pole location (from Fit 4_4L) and the new pole location found from fitting including the higher energy levels are shown in Fig. 12.We therefore conclude that it is a stable result from our analysis that, as soon as SU(3) constraints are included in the fits, there are always poles close to the D s K and Dη threshold-we do not find anymore the large scatter of the original amplitudes.

Summary and discussion
We investigated the pole content of the nine K -matrix parametrizations provided in Ref. [24] in an analysis of lattice data for open charm states in the (S = 0, I = 1 2 ) channel.In addition to the bound state pole reported in Ref. [24], in every amplitude additional poles were found on unphysical Riemann sheets, however, their locations vary strongly between the different parametrizations.On the other hand various investigations employing UChPT find that the structure observed in various experiments in the channel (S = 0, I = 1 2 ) with open charm originates from the interplay of two D * 0 poles.In this paper we explain the origin of this seeming contradiction.In particular it is shown that also in the lattice analysis two poles are needed and that, although the poles scatter so dramatically in location, their effects on the amplitudes were comparable in all parametrizations.This is possible since all poles are located on an hidden sheet, such that their effect on the scattering amplitude becomes visible at the threshold.In such a situation the distance from the threshold can be overcome by an enhanced residue.This mechanism was observed before as a general feature of Flatté amplitudes [26].
To discuss the location of the higher pole directly from the lattice data, we propose to use an amplitude constrained by SU(3) flavor symmetry.The flavor constrained amplitude well reproduces the energy levels and produces a pole in the RS221 sheet close to Dη and D s K thresholds consistent to that of the UChPT amplitude.Such an SU(3) symmetric construction of the K matrix may be used in analyzing other lattice data and also experimental data where multiple channels are involved.

Fig. 1
Fig. 1 Illustration for the sheet labeling in the case of two channels

Fig. 3
Fig. 3 The distance of the real part of the pole on sheet RS221 from the D s K threshold versus the effective coupling of the pole to the Dπ channel (left), Dη channel (middle) and D s K channel (right).The red line shows the straight line fit.The red band encloses the 1σ uncertainty of the fit

Fig. 8 Fig. 9
Fig. 8 The resulting amplitudes from the various fits in the form of ρ 2 |T | 2 .The Dπ -Dπ, Dη-Dη and D s K -D s K amplitudes are shown in red, green and blue, respectively.The vertical red, green and blue

Fig. 10
Fig. 10 The left panel shows the finite volume energy levels arrived at for Fit 4_4L (only the first four lowest energy levels of each volume included in the fit), while the right panel shows them for Fit 4_All.The black crosses are the lattice energy level data.The dark red circles

Fig. 11
Fig. 11 The resulting amplitudes of Fig 4_ALL in the form of ρ 2 |T | 2 .The Dπ -Dπ, Dη-Dη and D s K -D s K amplitudes are shown in red, green and blue respectively.The vertical red, green and blue lines show the Dπ , Dη and the D s K thresholds, respectively.The error bands cover the 1σ statistical uncertainties

Fig. 14 Fig. 15
Fig. 14 The distance of real part of the pole on sheet RS211 from Dη threshold versus the effective coupling of the pole to the Dπ channel (left), Dη channel (center) and D s K channel (right).The red line shows the straight line fit.The red band encloses the 1σ uncertainty of the fit

Table 4
The two particle thresholds in MeV for the pion mass of 391 MeV dof

Table 7
The pole locations from the different fits

Table 8
The absolute value of the square root of the RS221 pole residues obtained by the SU(3) flavor constrained K -matrix to the respective channels