Nature of the D_0 meson in the D pi scattering with chiral symmetry

TL;DR

This paper investigates the scalar D_0 meson using chiral symmetry, modeling it as a resonance in D pi scattering, and finds its properties are robust against changes in chiral condensate, emphasizing its molecular nature.

Contribution

It constructs a chiral-symmetric D pi scattering amplitude showing the D_0 as a chiral partner and highlights the meson's molecular dressing effects.

Findings

D_0 appears as an s-wave resonance in D pi scattering.

Properties of D_0 are insensitive to chiral condensate reduction.

D_0's spectral function and pole are influenced by D pi molecular components.

Abstract

We study the nature of the scalar D_0 meson from the viewpoint of chiral symmetry. With the linear representation of chiral symmetry, we construct the D pi scattering amplitude satisfying the chiral low-energy theorem, in which the D_0 meson appears as an s-wave resonance. We show that the properties of the D_0 meson can be successfully reproduced as the chiral partner of the D meson coupled with the D pi scattering states. At the same time, we find that the spectral function and the pole position of D_0 are not very sensitive to the reduction of the chiral condensate, indicating the importance of the dressing of the bare state by the D pi molecular component.

Equations

$$H_{L,R}\to S{H}_{L,R}{g}_{L,R}^{†},$$

$$M\to g_{L}M{g}_R^{†},$$

$$\mathcal{L}$$

$$-\frac{{\Delta }_m}{2f_{\pi }}\text{Tr }[H_{L}M{\bar{H}}_R+H_R{M}^{†}{\bar{H}}_L],$$

$$H$$

$$G$$

$$\mathcal{L}$$

$$+\frac{m{\Delta }_m}{2f_{\pi }}[D(M+M^{†})D^{†}-D_0(M+M^{†})D_0^{†}$$

$$-D_0(M-M^{†})D^{†}+D(M-M^{†})D_0^{†}]+\cdots ,$$

$$\mathcal{L}$$

$$\displaystyle\quad-\frac{\lambda}{4}(MM^{{\dagger}})^{2}+\epsilon(M+M^{{\dagger}})\Bigr{]},$$

$$m_{\pi }^2$$

$$m_{\sigma }^2$$

$$M_D$$

$$M_{D_0}$$

$$\mathcal{L}$$

$$V_{ij}(s,t,u)$$

$$+\left(\begin{array}{cc}1& \sqrt{2}\\ \sqrt{2}& 0\end{array}\right)\frac{m^2{\Delta }_m^2}{f_{\pi }^2}\frac{1}{u-M_{D_0}^2}$$

$$+\left(\begin{array}{cc}-1& 0\\ 0& -1\end{array}\right)\frac{m{\Delta }_m(m_{\sigma }^2-m_{\pi }^2)}{f_{\pi }\bar{\sigma }}\frac{1}{t-m_{\sigma }^2},$$

$$V^{1/2}(s,t,u)$$

$$-\frac{(M_{D_0}^2-M_D^2{)}^2}{4{\bar{\sigma }}^2}\frac{1}{u-M_{D_0}^2}$$

$$-\frac{(M_{D_0}^2-M_D^2)(m_{\sigma }^2-m_{\pi }^2)}{2{\bar{\sigma }}^2}\frac{1}{t-m_{\sigma }^2},$$

$$V^{3/2}(s,t,u)$$

$$-\frac{(M_{D_0}^2-M_D^2)(m_{\sigma }^2-m_{\pi }^2)}{2{\bar{\sigma }}^2}\frac{1}{t-m_{\sigma }^2},$$

$$\left(\begin{array}{c}|D\pi (I=1/2)⟩\\ |D\pi (I=3/2)⟩\end{array}\right)$$

$$V^I$$

$$=\frac{1}{{\bar{\sigma }}^2}\times \{\begin{array}{ll}-2& (I=1/2)\\ 1& (I=3/2)\end{array}.$$

$$s$$

$$t$$

$$u$$

$$V^{1/2}$$

$$-\frac{1}{4}\frac{M_{D_0}^2-M_D^2}{{\bar{\sigma }}^2}\left[-1+\frac{2M_D\omega }{M_{D_0}^2-M_D^2}\right]$$

$$+\frac{2}{4}\frac{M_{D_0}^2-M_D^2}{{\bar{\sigma }}^2}+\mathcal{O}(Q^2)$$

$$=-\frac{2M_D\omega }{{\bar{\sigma }}^2}+\mathcal{O}(Q^2),$$

$$V^{3/2}$$

$$V(\sqrt{s})=\frac{1}{2}{\int }_{-1}^{1}d\cos \theta V(\sqrt{s},\cos \theta ).$$

$$p$$

$$V(\sqrt{s})$$

$$-\frac{(M_{D_0}^2-M_D^2{)}^2}{4{\bar{\sigma }}^2}\frac{1}{4p^2}\ln \frac{M_{D_0}^2-\frac{(m_{\pi }^2-M_D^2{)}^2}{s}}{M_{D_0}^2-\frac{(m_{\pi }^2-M_D^2{)}^2}{s}+4p^2}$$

$$-\frac{(M_{D_0}^2-M_D^2)(m_{\sigma }^2-m_{\pi }^2)}{2{\bar{\sigma }}^2}\frac{1}{4p^2}\ln \frac{m_{\sigma }^2}{m_{\sigma }^2+4p^2}.$$

$$T(\sqrt{s})$$

Full text

I Introduction

There have been a lot of newly observed hadrons Tanabashi:2018oca thanks to the recent developments of high-energy experiments. The findings of various exotic states suggest possible configurations of hadrons beyond the simple $\bar{q}q$ and $qqq$ classifications in conventional quark models. The suggestions for exotic configurations include, for instance, multiquarks, hadronic molecules, gluon hybrids, and so on. Investigations of understanding the nature of the exotic hadrons are intensively performed both experimentally and theoretically Brambilla:2010cs ; Hosaka:2016pey ; Guo:2017jvc .

Among many hadrons, charmed mesons composed of $c\bar{u}$ or $c\bar{d}$ are of particular interest from the viewpoint of QCD symmetries. In the massless limit of the light up ($u$) and down ($d$) quarks, QCD exhibits chiral SU(2)${}_{L}\times$SU(2)R symmetry. Important low-energy properties of hadrons are constrained by chiral symmetry and its spontaneous breaking Hosaka:2001ux ; Scherer:2012xha . In the limit of infinitely heavy mass of the charm ($c$) quark, hadrons with a single charm quark follows the heavy quark symmetry Manohar:2000dt , which also governs the dynamics of those hadrons. Thus, the heavy-light mesons provide a unique playground where two distinct symmetries in QCD interplay with each other. These symmetries are also important to study of the property change of the charmed mesons in nuclear matter Hosaka:2016ypm .

In this work, we focus on the scalar $D_{0}$ meson, which was firstly observed as a broad peak in the $D\pi$ invariant mass distribution in 2004 Link:2003bd . Because $D_{0}$ has an opposite parity to the ground state pseudoscalar $D$ meson, it can be regarded as a chiral partner of the $D$ meson Nowak:1992um ; Bardeen:1993ae ; Bardeen:2003kt ; Suenaga:2017deu ; Suenaga:2018kta . On the other hand, its broad decay width to the $D\pi$ state implies the picture of the hadronic molecular state which is dynamically generated by the meson-meson interaction Kolomeitsev:2003ac ; Hofmann:2003je ; Guo:2006fu ; Gamermann:2006nm ; Guo:2009ct ; Altenbuchinger:2013vwa ; Du:2017ttu ; Du:2017zvv ; Guo:2018tjx ; Du:2019oki . $D_{0}$ is also regarded as a flavor partner of the $D_{s0}(2317)$ state in the strangeness sector, which is known as a candidate of the exotic hadron.

The aim of this paper is to combine two pictures on the $D_{0}$ state, the chiral partner of $D$ and the hadronic molecule. To this end, we utilize the linear representation of chiral symmetry to incorporate the chiral partner structure, and describe the $D_{0}$ meson as a resonance in the scattering amplitude through the nonperturbative resummation. This framework is analogous to the analysis of the $\sigma$ meson in the $\pi\pi$ scattering in Refs. Basdevant:1969sz ; Basdevant:1970nu ; Jido:2000bw ; Hyodo:2010jp . The nature of the $D_{0}$ resonance is then studied by the property change with respect to the chiral condensate.

The paper is organized as follows. In Sec. II, we start from the effective Lagrangian of charmed mesons, respecting chiral and heavy quark symmetries in QCD. The $D\pi$ scattering amplitude is then constructed, with the special emphasis of the chiral low-energy theorem. In Sec. III, we numerically study the properties of the $D_{0}$ meson which is realized as a resonance in the $D\pi$ scattering amplitude. The last section is devoted to a summary of this work.

II Formulation

II.1 Effective Lagrangian

To include the chiral partner structure of charmed mesons, we employ the effective Lagrangian based on the linear representation of chiral symmetry Nowak:1992um ; Bardeen:1993ae ; Bardeen:2003kt ; Suenaga:2017deu . We first introduce the heavy-light meson fields $H_{L,R}\sim Q\bar{q}_{L,R}$, which transform as

[TABLE]

under chiral symmetry $(g_{L},g_{R})\in$ SU(2)${}_{L}\times$SU(2)R and heavy quark spin symmetry $S\in$ SU(2)S. The light mesons $\sigma$ and $\pi^{a}$ $(a=1,\ldots,3)$ are collected in the matrix field in the isospin space $M=\sigma+i\pi^{a}\tau^{a}\sim q_{L}\bar{q}_{R}$, which transforms as

[TABLE]

under chiral SU(2)${}_{L}\times$SU(2)R symmetry. The effective Lagrangian with no derivatives of the light meson fields, which is invariant under chiral and heavy quark spin symmetries, can be constructed as

[TABLE]

where $\bar{H}_{L,R}\equiv\gamma^{0}H_{L,R}^{{\dagger}}\gamma^{0}$, $v_{\mu}$ is the vector to specify the reference frame of the heavy quark, and the trace is taken over the Dirac space.111Note that $H_{L,R}$ ($\bar{H}_{L,R}$) is a row (column) vector in the isospin space, because it contains a light antiquark (quark). In the interaction term, the coupling constant $\Delta_{m}$ is normalized by the pion decay constant $f_{\pi}$ for later convenience. To establish the connection to physical heavy-light mesons, we first introduce the negative (positive) parity heavy-light meson multiplet $H$ ($G$) as $H_{L,R}=[G\pm iH\gamma_{5}]/\sqrt{2}$. In the charm quark sector, $H$ ($G$) contains pseudoscalar $D$ and vector $D^{*}$ (scalar $D_{0}$ and axial vector $D_{1}$) as the heavy quark spin doublet, which are collected as

[TABLE]

Substituting these into Eq. (3), we obtain the effective Lagrangian of heavy-light mesons. The relevant terms for the present study are given by

[TABLE]

where $D=(D^{0},D^{+})$, $D^{{\dagger}}=(\bar{D}^{0},D^{-})^{t}$ and we have adopted the relativistic notation as in Ref. Suenaga:2017deu with $m$ being the chiral invariant mass of the charmed mesons. The first line of Eq. (14) denotes the free fields of $D$ and $D_{0}$, having the common mass $m$ due to the chiral partner structure, and the second and third lines represent the interactions with $\sigma$ and $\pi$, respectively. It is a consequence of chiral symmetry that the couplings to $\sigma$ and $\pi$ are dictated by a common coupling constant $\Delta_{m}$.

The light meson part of the effective Lagrangian is

[TABLE]

where the trace is taken over the isospin space. For a negative $m_{0}^{2}$, chiral symmetry is spontaneously broken at the mean field level. Denoting the mean field value of the $\sigma$ field as $\bar{\sigma}$ and redefining the $\sigma$ field as the fluctuation around $\bar{\sigma}$, we can read off the hadron masses from Eqs. (14) and (15) as

[TABLE]

In the light meson sector, three parameters in the Lagrangian $m_{0}$, $\lambda$, and $\epsilon$ are determined by the physical values of $m_{\pi}$, $m_{\sigma}$, and the chiral condensate in vacuum, which corresponds to the pion decay constant $\bar{\sigma}=f_{\pi}$. Therefore, the mass difference between $D$ and $D_{0}$ in vacuum is given by $\Delta_{m}=M_{D_{0}}-M_{D}$.

As mentioned above, the parameter $\Delta_{m}$ also determines the $D_{0}D\pi$ vertex in the third line of Eq. (14). Thus, for a given $M_{D_{0}}$, we can predict the decay width $\Gamma_{D_{0}}$ in the perturbative calculation. However, if we use the central value adopted by the Particle Data Group (PDG) Tanabashi:2018oca $M_{D_{0}}=2318$ MeV to determine $\Delta_{m}$, we obtain $\Gamma_{D_{0}}=1115$ MeV, which largely deviates from the PDG value $\Gamma_{D_{0}}=267\pm 40$ MeV. In other words, $D_{0}$ as a pure chiral partner state is not consistent with the experimental data.222See also the recent discussion on the PDG value in Ref. Du:2019oki . In fact, huge decay width obtained by the perturbative estimate indicates the strong coupling, for which the nonperturbative treatment of the coupling is needed. In the following, we regard $D_{0}$ field in the Lagrangian as a bare state, and describe the physical $D_{0}$ state as a resonance in the nonperturbative $D\pi$ scattering amplitude.

II.2 Tree-level $D\pi$ interaction

Let us construct the $D\pi$ scattering amplitude at tree level. From Eq. (14), we have the $DD\sigma$ and $D_{0}D\pi$ vertices

[TABLE]

and from Eq. (15), we have $\sigma\pi\pi$ vertex

[TABLE]

From these vertices, the $D\pi$ scattering amplitude at tree level can be constructed by the $s$-channel $D_{0}$ exchange, $u$-channel $D_{0}$ exchange, and $t$-channel $\sigma$ exchange, as shown in Fig. 1.333In this work, we consider the single-channel $s$-wave $D\pi$ scattering, and the couplings to $s$-wave $D_{0}\sigma$ channel and $p$-wave $D^{*}\sigma$, $D_{1}\pi$ channels are not included. These channels are in higher energies and do not contribute to the chiral low-energy theorem. Straightforward calculation leads to the tree-level amplitude in the particle basis. For instance, in the total charge $Q=+1$ sector, we obtain

[TABLE]

with $\pi^{0}D^{+}$ and $\pi^{+}D^{0}$ for $i=1$ and $2$, respectively. This coupled-channel interaction can be diagonalized by taking the isospin basis where the interaction $V^{I}$ for isospin $I$ reads

[TABLE]

with

[TABLE]

In Eqs. (23) and (24), we replace the coupling strength $m\Delta_{m}$ by the masses of the $D$ mesons using Eqs. (18) and (19). Because of the isospin symmetry of the effective Lagrangians, there is no transition between $I=1/2$ and 3/2 sectors. Note also that the $s$-channel pole term only exists in the $I=1/2$ sector, reflecting the isospin of the $D_{0}$ state.

II.3 Chiral low-energy theorem

Before proceeding to the calculation of the nonperturbative scattering amplitude, it is instructive to examine the results (23) and (24) from the viewpoint of chiral low-energy theorem. Regarding pions as Nambu-Goldstone (NG) bosons associated with the spontaneous breaking of chiral symmetry, we can impose several conditions on the low-energy dynamics of pions. In particular, low-energy scattering of pions with any target hadron is determined only by the pion decay constant and the total isospin of the system, as a consequence of the Weinberg-Tomozawa theorem Weinberg:1966kf ; Tomozawa:1966jm . In the case of the $D\pi$ scattering, the low-energy scattering should behave as

[TABLE]

This is an important constraint on the hadron scattering amplitude in the following sense. First, from Eq. (26) the sign (either attractive or repulsive) and strength of the interaction are uniquely determined, according to the total isospin of the system. Second, the relation is insensitive to the internal structure of the target hadron, except for the isospin of the state. This means that, for instance, the $\pi N$ system also follows the same relation (26), which was confirmed experimentally to a good accuracy by the $\pi N$ scattering lengths Scherer:2012xha . In chiral perturbation theory, the Weinberg-Tomozawa theorem is automatically satisfied by the leading order calculation, while a careful examination is needed in the linear representation as in the present study.

To study the low-energy (near-threshold) behaviors of Eqs. (23) and (24), we assign the meson momenta as $D(k)\pi(p)\to D(k^{\prime})\pi(p^{\prime})$. We regard the three-momentum of the mesons in the center-of-mass frame $|\bm{p}|=|\bm{k}|=|\bm{p}^{\prime}|=|\bm{k}^{\prime}|$ as a small quantity. Denoting it as $\mathcal{O}(Q)$, we expand the interaction $V^{I}$ in powers of $Q$. Masses of non-NG bosons, $m_{\sigma}$, $M_{D}$, and $M_{D_{0}}$ are counted as $\mathcal{O}(Q^{0})$, while the pion mass is regarded as $m_{\pi}\sim\mathcal{O}(Q^{1})$ according to the chiral counting. In this case, the energy of pions is counted as $\mathcal{O}(Q^{1})$ that we define $p_{0}=p_{0}^{\prime}\equiv\omega$. The energy of the $D$ mesons $k_{0}=k_{0}^{\prime}$ is counted as $M_{D}+\mathcal{O}(Q^{2})$. Then the Mandelstam variables are expanded as

[TABLE]

We finally obtain the low-energy behaviors of Eqs. (23) and (24) as

[TABLE]

We confirm that the results satisfy the relation (26) by the Weinberg-Tomozawa theorem. Furthermore, the energy dependence agrees with the prediction by chiral perturbation theory $V\propto(s-u)/4=M_{D}\omega+\mathcal{O}(Q^{2})$ Guo:2006fu .

It is worth noting that the contribution from each channel has an $\mathcal{O}(Q^{0})$ term, in contradiction to the low-energy theorem which requires $\mathcal{O}(Q^{1})$ as the leading order. In other words, the use of the $s$-channel pole term with energy independent coupling is not consistent with the chiral low-energy theorem. By summing up all the channels required by the chiral invariant Lagrangian, the $\mathcal{O}(Q^{0})$ terms cancel out with each other, and we obtain $V^{I}\sim\mathcal{O}(Q^{1})$ in agreement with the low-energy theorem. This is analogous to the $\pi N$ scattering in the linear sigma model, where the sum of the $s$- and $u$-channel Born terms with pseudoscalar coupling leads to $\mathcal{O}(Q^{0})$ contribution, and the inclusion of the $t$-channel $\sigma$ exchange recovers the result of the Weinberg-Tomozawa theorem.

II.4 Nonperturbative $D\pi$ scattering amplitude

To describe the $D_{0}$ meson as a resonance in the $D\pi$ scattering, we project the interaction $V(s,t,u)$ to $s$ wave, and perform nonperturbative resummation to satisfy the unitarity. From now on, we concentrate on the $I=1/2$ channel and suppress the isospin index. In the center of mass frame, the $t$ and $u$ variables can be specified by the total energy $\sqrt{s}$ and the scattering angle $\cos\theta=\bm{p}\cdot\bm{p}^{\prime}/|\bm{p}\cdot\bm{p}^{\prime}|$ with the momentum assignment $D(k)\pi(p)\to D(k^{\prime})\pi(p^{\prime})$. The $s$-wave part $V(\sqrt{s})$ is obtained by averaging over the scattering angle as

[TABLE]

Defining the magnitude of the three momentum as

[TABLE]

with the Källén function $\lambda(x,y,z)=x^{2}+y^{2}+z^{2}-2xy-2yz-2zx$, we obtain the expression for the $s$-wave projected interaction:

[TABLE]

Because the tree-level interaction does not satisfy the unitarity condition of the $S$ matrix, various resummation (unitarization) schemes have been proposed. Here we employ the N/D method Oller:1998zr ; Hyodo:2010jp . By neglecting the contribution from the left-hand cut, we obtain the general expression of the scattering amplitude $T(\sqrt{s})$ from the tree-level interaction $V(\sqrt{s})$ as

[TABLE]

with

[TABLE]

where $a(\mu)$ is the subtraction constant at the subtraction scale $\mu$. This $G(\sqrt{s})$ function is obtained by the once-subtracted dispersion integral on the unitarity cut, and is essentially the same with the $D\pi$ loop function with the dimensional regularization. With this identification, Eq. (32) corresponds to the resummation of the loop diagrams in the Bethe-Salpeter equation.

A resonance in the $D\pi$ scattering is expressed as the pole of the analytically continued scattering amplitude in the second Riemann sheet of the complex energy plane. When the amplitude has a pole at $\sqrt{s}=z\in\mathbb{C}$, the mass $M_{R}$ and width $\Gamma_{R}$ of the resonance can be identified by

[TABLE]

Because of the nonperturbative resummation, the bare mass $M_{D_{0}}$ in the Lagrangian is shifted to $M_{R}-i\Gamma_{R}/2$, including the decay width to the $D\pi$ channel. Therefore, this resummation procedure can also be regarded as the renormalization of the $D_{0}$ propagator with the Dyson equation. We however note that the self-energy should include not only the simple $D\pi$ loop, as was done in Ref. Suenaga:2017deu , but also the crossed diagrams, in order to satisfy the chiral low-energy theorem.

As we discussed in Sec. II.3, the low-energy limit of our interaction $V(\sqrt{s})$ reduces to the Weinberg-Tomozawa term, which can dynamically generate a $D\pi$ molecule state without the explicit bare state Kolomeitsev:2003ac ; Hofmann:2003je ; Guo:2006fu ; Gamermann:2006nm ; Guo:2009ct ; Altenbuchinger:2013vwa ; Du:2017ttu ; Du:2017zvv ; Guo:2018tjx ; Du:2019oki . On top of that, our interaction contains the explicit coupling to the bare $D_{0}$ state which is introduced as a chiral partner of $D$. Thus, our model is capable of describing both the chiral partner state and the $D\pi$ molecule state. In the following, we calibrate the model by empirical data and study the nature of the $D_{0}$ resonance in the $D\pi$ scattering amplitude.

III Numerical analysis

III.1 $D_{0}$ resonance

Now we numerically study the property of the $D_{0}$ resonance. We first set the parameters of the model to the empirical values,

[TABLE]

In addition, we need to specify one degree of freedom to determine the ultraviolet cutoff in Eq. (33). Here we set $\mu=1000$ MeV and use the subtraction constant $a$ at this scale as a free parameter for the cutoff degree of freedom. The bare mass $M_{D_{0}}$ in the Lagrangian is also used as a free parameter.

As a physical input to determine the free parameters, we use the mass and width of $D_{0}$ in PDG Tanabashi:2018oca . According to Eq. (34), we should have a pole of the $D\pi$ scattering amplitude at

[TABLE]

Note however that the PDG values are obtained by the Breit-Wigner parametrization which is recently challenged by Ref. Du:2019oki . At present, we simply adopt Eq. (37) as a representative value, but it is also possible to use the updated pole positions as suggested in Ref. Du:2019oki . Such analysis is reserved for a future work.

The best fit values of the bare mass and the subtraction constant are found to be $M_{D_{0}}=2024$ MeV and $a=2.5$, with which we obtain the $D_{0}$ pole at

[TABLE]

in good agreement with Eq. (37). We also find a virtual state pole at $\sqrt{s}=1795$ MeV. Although the location of the pole is far below the threshold (2005 MeV), this may correspond to the one found in the sigma propagator in Refs. Hidaka:2003mm ; Hidaka:2004ht . We show the reduced $D\pi$ cross section $\Theta(\sqrt{s}-m_{\pi}-M_{D})|T(\sqrt{s})|^{2}/s$ by solid line in Fig. 2. This corresponds to the spectrum of the scalar channel with isospin $I=1/2$, and the peak of the $D_{0}$ resonance is seen. We also find that the pole position of Eq. (37) cannot be reproduced with some values of the bare mass, even if we vary the subtraction constant freely. Namely, the bare mass should lie in the region $2001\text{ MeV}\leq M_{D_{0}}\leq 2045$ MeV in order to satisfy Eq. (37) within the uncertainty.

III.2 Reduction of chiral condensate

While the chiral partner state is included in the model as a bare state, the bare mass of 2024 MeV is more than 300 MeV away from the actual pole position. This may indicate that the bare pole contribution is not very relevant to the structure of the physical state. To check this qualitatively, we vary the value of the chiral condensate $\bar{\sigma}$ from its vacuum value, and study the response of the spectrum and the pole position. This is partly motivated by the reduction of the chiral condensate in nuclear medium Hayano:2008vn , but the quantitative discussion on the in-medium spectrum will require the actual many-body calculation as was done in Refs. Suenaga:2017deu ; Suenaga:2018kta .

Here we introduce a reduction parameter $\Phi$ as

[TABLE]

and vary it within $0.7\leq\Phi\leq 1$. The masses of the $D$ mesons are given by Eqs. (18) and (19), with fixed $m$ and $\Delta_{m}$. Namely, the mass of $D$ (bare mass of $D_{0}$) increases (decreases) as we decrease $\Phi$ from unity. For the light mesons, we assume that the pion mass is independent of $\Phi$ and the mass of sigma varies as $m_{\sigma}^{2}=m_{\pi}^{2}+2\lambda\bar{\sigma}^{2}$ with $\lambda$ fixed, as suggested in Ref. Hyodo:2010jp . Namely, the mass of sigma decreases along with the reduction of $\Phi$.

The modification of the $D_{0}$ spectrum along with the reduction of $\bar{\sigma}$ is shown in Fig. 2 by the dashed lines. We find that the spectrum is barely changed, even if the chiral condensate is reduced by 30 %. This shows that the resonance is not sensitive to the change of the chiral condensate. This can also be seen by the trajectory of the pole shown in Fig. 3. The change of the pole position is a few tens of MeV with the $30$ % reduction of $\bar{\sigma}$. Moreover, the real part of the pole increases when the chiral condensate is reduced, although the bare mass is moved to the lower energies. We also note that the virtual state pole at 1795 MeV moves toward the threshold. This is in line with the findings in Refs. Hidaka:2003mm ; Hidaka:2004ht , although the pole locates still far away from the threshold ($\sqrt{s}=1864$ MeV at $\Phi=0.7$).

As found in Ref. Hyodo:2010jp , the pole of the chiral partner state is naively expected to move toward the threshold with reducing the imaginary part due to the decrease of the phase space. We confirm this behavior by artificially tuning the subtraction constant ($a=0$) to have a narrow $D_{0}$ resonance (pole at $\sqrt{s}=2069-32i$ MeV with $\Phi=1$). This indicates that, in the physical $D_{0}$ resonance (with $a=2.5$), the dressing by the $D\pi$ molecule component is so large that it washes out the chiral partner nature of the bare state. In this way, we conclude that the $D\pi$ molecule component is important for the physical $D_{0}$ meson, on top of the bare state which is the chiral partner of the $D$ meson.

IV Summary

We have studied the properties of the $D_{0}$ resonance in the $D\pi$ scattering, utilizing the chiral effective Lagrangian. We introduce the chiral partner of the $D$ meson as an explicit bare state in the linear representation of chiral symmetry. At the same time, we respect the chiral low-energy theorem which leads to the dynamical generation of the $D\pi$ molecule state. Because the empirical data can be explained by this mode, we demonstrate that the physical $D_{0}$ resonance can be described as a superposition of the chiral partner of $D$ and the $D\pi$ molecular state. By studying the response to the modification of the chiral condensate, we find that the bare $D_{0}$ state is highly dressed by the $D\pi$ molecular component to realize the physical $D_{0}$ resonance. Thus, we conclude that the $D\pi$ molecular component dominates the internal structure of the $D_{0}$ meson.

As future prospects, we can improve the empirical information of the $D_{0}$ state as suggested by Ref. Du:2019oki . It is possible to use the updated pole position as an input, or to directly analyze the $B$ meson decay data. The generalization to three-flavor case is also important, because the thresholds of $D\eta$ and $DK$ channels are located around 2400 MeV, whose contribution may also affect to the $D_{0}$ resonance. In this case, the $D_{s0}(2317)$ can also be studied in the same framework. Another direction would be to explore the modification of the $D_{0}$ property in nuclear medium, which can be pursued with the techniques developed in Refs. Suenaga:2017deu ; Suenaga:2018kta .

Acknowledgments

The authors thank Masayasu Harada, Atsushi Hosaka, Feng-Kun Guo, Daisuke Jido, Daiki Suenaga, and Shigehiro Yasui for fruitful discussions. This work is supported in part by JSPS KAKENHI Grant No. 24740152, and by the Yukawa International Program for Quark-Hadron Sciences (YIPQS).