-
Being closely related to color confinement, hadron spectroscopy is one of the most important aspects of the nonperturbative quantum chromodynamics (QCD). Many new hadron resonances or resonance-like structures have been observed since 2003 at worldwide experiments, including Belle, BaBar, BESIII, CDF, LHCb and so on. In particular, many of them were observed in the heavy-quarkonium mass region, and have properties difficult to be understood from the quark model point of view. Thus, they are called
$ XYZ $ states, and have spurred plenty of experimental and phenomenological investigations, as well as studies using lattice QCD. For recent reviews, we refer to Refs. [1–13].In order to understand the physics behind the messy spectrum of the
$ XYZ $ structures, we need to be careful about interpreting experimental observations. Most of these structures were discovered by observing a peaking structure in the invariant mass distribution of two or three particles in the final state. Peaking structures, in particular the narrow ones, are often due to singularities of the$ S $ -matrix, which have different kinds including poles and branch points. Resonances are poles of the$ S $ -matrix, while the branch points arise from unitarity and are due to the on-shellness of intermediate particles. The simplest one of the latter is the two-body threshold cusp. It is a square-root branch point and shows up exactly at all$ S $ -wave thresholds coupled to the measured energy distributions. The strength of the cusp depends on the masses of the involved particles and the interaction strength of the rescattering from the intermediate two particles to the final states. Triangle singularity is more complicated. It is due to three on-shell intermediate particles, see Fig. 1, and happens on the physical boundary when the interactions at all the three vertices happen as classical processes in spacetime[14]. The triangle singularity is a logarithmic branch point, and thus can lead to drastic observable effects if it is located close to the physical region. The threshold cusp and triangle singularity are just two examples of the more general Landau singularities[15]. For a recent review of threshold cusps and TSs in hadronic reactions, we refer to Ref. [13].The triangle singularities related to the
$ XYZ $ structures have been discussed in Refs. [16–36]. Here we focus on those related to the charged charmonium-like structures observed in the$ e^+e^-\to J/\psi\pi^+\pi^- $ and$ \psi(3686)\pi^+\pi^- $ reactions. Notice that the triangle singularity mechanism is not supposed to be the only reason behind the relevant observed peaks. In particular, there can be$ Z_c $ resonances in addition, and the production of near-threshold resonances can get enhanced by the mechanism as emphasized in Ref. [13]. Thus, the formalism considering the triangle singularity effects should also have the freedom to include$ Z_c $ contributions through the$ \psi\pi $ ,$ D\bar D^* $ and$ D^*\bar D^* $ coupled-channel$ T $ -matrix. We shall present formulae that will be useful for an analysis of the data taking into account triangle singularities for such processes. In Sec. 2, we discuss the relevant triangle diagrams for the reactions of interest, point out that both$ S $ - and$ D $ -wave couplings of the$ D_1(2420) $ to the$ D^*\pi $ need to be taken into account, and give the expressions which can be used to account for the triangle singularity effects in the analysis of the$ e^+e^-\to \psi\pi^+\pi^- $ data. A brief summary is given in Sec. 3. The scalar triangle loop integral is evaluated in the appendix. -
The location of the peak induced by a triangle singularity is normally not far from the corresponding two-body threshold. Thus, for the
$ Z_c(3900) $ [37–41] and$ Z_c(4020) $ [42-43], it is important to discuss the triangle diagrams with intermediate$ D\bar D^*+c.c. $ and$ D^*\bar D^* $ , respectively, coupled to the final states where the$ Z_c $ structures were observed. It has been pointed out in Refs. [16-17, 23, 27] that the$ D_1\bar D D^*+c.c. $ loops are important for the understanding of the$ Z_c(3900) $ .①In order to see the possible impact of the
$ D_1\bar D D^* $ triangle diagram on the structures in both the final state$ \psi\pi $ and initial state$ \psi\pi\pi $ line shapes, a 3D plot for the$ |I(D_1, \bar D, D^*)|^2 $ is shown in Fig. 2, where$ I $ refers to the loop integral$ I(m_1^2, m_2^2, m_3^2, p_{12}^2, p_{13}^2, p_{23}^2) $ defined in Eq. (13) in the appendix and we use the particle names to represent the corresponding$ m_i^2 $ and have neglected$ p_{ij}^2 $ for simplicity. The$ \psi\pi $ pair comes from the$ \bar D D^* $ rescattering (see Fig. 3(b) below), and the$ D_1 $ width is taken into account by using a complex mass of the form$ m_1-{\rm i}\Gamma_1/2 $ . It is clear that the triangle loop integral is able to produce a peak in both the final state$ \psi\pi $ invariant mass and the initial energy$ \sqrt{s} $ distributions.Figure 2. (color online) The absolute value squared of the scalar three-point loop integral for the
$ D_1\bar D D^*$ intermediate particles as a function of the$ \psi\pi$ invariant mass and the initial energy$ \!\sqrt{s}$ .The sharp peak in the 3D plot is due to the presence of a triangle singularity which would be on the physical boundary if the width of the
$ D_1 $ is neglected. One sees that the peak is around the$ \bar D D^* $ threshold and the$ D_1\bar D $ threshold in the$ m_{J/\psi\pi} $ and$ \sqrt{s} $ distributions, respectively. In addition, there is a clear cusp at the$ \bar D D^* $ threshold in the$ m_{J/\psi\pi} $ distribution, which is the manifestation of the two-body threshold as a subleading singularity of the triangle diagram. Because of the finite width of the$ D_1 $ there is no such an evident cusp in the$ \sqrt{s} $ distribution. Because of the singular behavior, it is thus crucial to have the triangle diagrams included in the analysis of the data in order to extract the resonance parameters of the$ Z_c $ or even to conclude whether it is necessary to introduce a$ Z_c $ . This is the point of Ref. [16] which concluded the necessity of the$ Z_c(3900) $ , an opinion shared in Refs. [23, 27], and suggested the importance of the$ D_1\bar DD^* $ triangle singularity for the first time.The BESIII data for both the
$ J/\psi\pi $ and$ D\bar D^* $ invariant mass distributions were later on reanalyzed considering such triangle diagrams in Refs. [23, 26], while Ref. [23] concluded that there should exist a$ Z_c(3900) $ as either a resonance pole above the$ D\bar D^* $ threshold or a virtual state below it, Ref. [26] concluded that the data could be fitted comparably well without introducing the$ Z_c(3900) $ . One notable difference in the treatment of the triangle diagrams in these two references is that the$ D_1D^*\pi $ coupling was treated as$ D $ -wave in Ref. [23] and as$ S $ -wave in Ref. [26]. -
The smallness of the
$ D_1(2420) $ width,$ (31.7\pm2.5) $ MeV[44], suggests that it is approximately a charmed meson with$ j_\ell^P\!=\!\frac32^+ $ , where$ j_\ell $ is the angular momentum of the light degrees of freedom in the$ D_1 $ , including the light quark spin and the orbital angular momentum. Because$ j_\ell^P\!=\!\frac12^- $ for the ground state$ D^{(*)} $ , a$ \frac32^+ $ meson decays into the$ D^{(*)}\pi $ in a$ D $ wave. Thus, the$ D $ -wave$ D_1D^*\pi $ coupling was used in Refs. [16, 23]. However, it turns out that the$ D $ -wave decay can only account for about half of the$ D_1(2420) $ decay width as we discuss now. There is another$ \frac32^+ $ charmed meson$ D_2(2460) $ which is the spin partner of the$ D_1(2420) $ . Its decays into the$ D^{(*)}\pi $ are purely$ D $ -wave. Thus, one can fix the$ D $ -wave decay coupling constant$ h_D $ defined in the following Lagrangian (see Refs. [45-46] for the Lagrangian using the four-component notation)$$ {\cal{L}}_D = \frac{h_D}{2F_\pi} {\rm{Tr}} \left[T_b^i\sigma^j H_a^{\dagger } \right] \partial^i\partial^j\pi_{ba}, $$ (1) which satisfies heavy quark spin symmetry (HQSS), where
$$ \begin{split} H_a =& {{D}}^*_a\cdot {{\sigma}} + D_a, \\ T_a^i =& D_{2a}^{ij} \sigma^j + \sqrt{\frac23}\, D_{1a}^i + {\rm i} \sqrt{\frac16}\, \epsilon_{ijk} D_{1a}^j \sigma^k, \end{split} $$ (2) represent the
$ j_\ell\!=\!\frac12^- $ and$ \frac32^+ $ spin multiplets, respectively,$ {{\sigma}} $ is the Pauli matrices in the spinor space, Tr$ [\cdot] $ denotes the trace in the spinor space, and$ \pi_{ba} $ represents the pion fields with the subindices$ a,b $ the indices in the light flavor space:$$ \pi = \left( {\begin{array}{*{20}{c}} {{\pi ^0}/\sqrt 2 }&{{\pi ^ + }}\\ {{\pi ^ - }}&{ - {\pi ^0}/\sqrt 2 } \end{array}} \right). $$ (3) From reproducing the central value of the
$ D_2 $ width,$ (47.5\pm1.1)\, {\rm MeV}$ , one gets$ |h_D|\!=\!1.17\, {\rm{GeV}}^{-1} $ . Using this value, one gets the$ D $ -wave contribution to the$ D_1(2420)\to D^*\pi $ width as 15.2 MeV, which is only about half of the$ D_1(2420) $ width. Assuming that the$ D^*\pi $ (and the sequential decay to$ D\pi\pi $ ) modes dominate the$ D_1 $ width, the rest of the$ D_1 $ width, about 16.5 MeV, should come from$ S $ -wave decays. The$ S $ -wave$ D_1D^*\pi $ coupling can be written as$$ {\cal{L}}_S = {\rm{i}}\frac{h_S}{\sqrt{6}F_\pi} {{D}}_{1b} \cdot {{D}}^{*{\dagger }}_a \partial^0 \pi_{ba} , $$ (4) with the coupling constant
$ |h_S|\!=\!0.57 $ . Here we require the$ S $ -wave coupling to be proportional to the pion energy to satisfy the Goldstone theorem because the pions are the pseudo-Goldstone bosons of the spontaneous breaking of chiral symmetry in QCD. -
Let us now construct the amplitude of the triangle diagrams for the reaction
$ e^+e^-\to \psi\pi^+\pi^-$ , and we consider all the possible THH-type diagrams with$ T $ and$ H $ representing the$ j_\ell^P\!=\!\frac32^+ $ and$ \frac12^- $ charmed mesons, respectively.The relevant THH-type triangle diagrams include
$ D_1\bar D D^*+c.c. $ ,$ D_1\bar D^* D^*+c.c. $ ,$ D_2\bar D^* D+c.c. $ and$ D_2\bar D^* D^*+c.c. $ The$ D_1\bar D D^*+c.c. $ diagrams for the charged$ D_1\bar D $ are shown in Fig. 3, and the analogous diagrams for the neutral$ D_1\bar D $ are not shown. The diagrams for the other mentioned THH triangles are similar. They were considered in Ref. [19]. Here we only consider the narrow$ D_{1(2)} $ mesons. However, one should keep in mind that their production together with a$ \bar D^{(*)} $ is suppressed in the heavy quark limit[47], and the broad$ D_{0}(D_1') $ might also play a role here. But the$ D_{0}(D_1') $ properties are still under discussion and the values listed in the Review of Particle Physics (RPP)[44], where were extracted from fitting to the$ D^{(*)}\pi $ invariant mass distributions using the Breit-Wigner parametrization, are not trustworthy (for detailed discussion, see Refs. [48-49]).Figure 3. The
$D_1\bar D D^*+c.c.$ triangle diagrams with charged$D_1\bar D$ contributing to the process$Y\to \psi\pi^+\pi^-$ with the initial$Y$ coming from$e^+e^-.$ For the narrow
$ D_1(2420) $ and$ D_2(2460) $ decays into$ D^{(*)}\pi $ , we use the coupling discussed in Sec. 2.2. In principal, the rescattering from the intermediate$ D^{(*)}\bar D^{(*)} $ into$ \psi\pi $ needs to be described through a coupled-channel$ T $ -matrix, see the treatment in Ref. [23]. For simplicity, one may also approximate the$ T $ -matrix by that from a$ Z_c $ exchange with the$ Z_c $ parametrized using a Flatté form to account for the$ D^{(*)}\bar D^* $ threshold.With the
$ Y, Z_c $ and$ \pi $ as the external particles, the amplitude for a triangle loop mentioned above shown as Fig. 3(a) with a$ D $ -wave TH$\! \pi $ coupling is proportional to$$ \begin{split}& {\cal{A}}_{D,a} \propto I(m_1^2, m_2^2, m_3^2, s, M_\pi^2, p_Z^2)\times \\ &\;\;\;\;\;\;\;\;\left(3 {{p}}_3\cdot {{\varepsilon}}_{Y} {{p}}_3\cdot {{\varepsilon}}_{Z} - {{p}}_3^{\,2} \,{{\varepsilon}}_{Y}\cdot {{\varepsilon}}_{Z} \right), \end{split} $$ (5) and that with an
$ S $ -wave$ D_1D^*\pi $ coupling is proportional to$$ {\cal{A}}_{S,a} \propto I(m_1^2, m_2^2, m_3^2, s, M_\pi^2, p_Z^2) E_3{{\varepsilon}}_{\psi}\cdot {{\varepsilon}}_{Z} , $$ (6) where
$ s $ is the c.m. energy squared for the$ e^+e^- $ ,$ p_Z^2 $ is the invariant mass squared of the meson pair coupled to the$ Z_c $ , and$ {{\varepsilon}}_{Y(Z)} $ are the spatial components of the polarization vectors for the$ Y(Z_c) $ . The polarization sum for each of them is$ \sum_{\lambda}\varepsilon_{(\lambda)}^i \varepsilon_{(\lambda)}^j\!=\!\delta^{ij} $ . Notice that the intermediate charmed mesons are treated nonrelativistically, so that the Lorentz boost effect from the$ D_1 $ rest frame to the$ e^+e^- $ c.m. frame is of higher order. Thus, the partial waves in the$ D_1D^*\pi $ coupling lend directly to the partial waves between the pion and$ Z_c $ (or the pair of the$ \psi $ and the other pion). The intermediate particles in the diagram shown in Fig. 3(b) is charge conjugated to those in the one in Fig. 3(a). The amplitude can be obtained by changing$ p_3 $ to$ p_2 $ in the above expressions.Here let us give expressions and relations for some kinematic variables entering into the analysis. We define the following variables:
$$ m_{12}^2 = (p_1+p_2)^2, \; \; m_{13}^2 = (p_1+p_3)^2, \; \; m_{23}^2 = (p_2+p_3)^2 . $$ (7) They satisfy
$$ m_{12}^2 + m_{13}^2 + m_{23}^2 = s + \sum\limits_{i = 1,2,3} M_i^2, $$ (8) where
$ M_i $ 's are the masses of the external particles in the final state. In terms of these variables,$ p_Z^2 $ in Eqs. (5) and (6) is$ m_{12}^2 $ , and that for Fig. 3(b) is$ m_{13}^2 $ .In the rest frame of the initial state, we express all momenta and energies in terms of
$ m_{12} $ and$ m_{23} $ :$$ \begin{split} |{{p}}_2| =&\frac1{2 \sqrt{s}} \sqrt{\lambda(s, M_2^2, m_{13}^2)}, \\ |{{p}}_3| =& \frac1{2 \sqrt{s}} \sqrt{\lambda(s, M_3^2, m_{12}^2)}, \\ E_2 =& \sqrt{M_2^2 + {{p}}_2^{\;2}}, \\ E_3 =& \frac1{2\sqrt{s}}\left( s+M_3^2-m_{12}^2\right), \\ {{p}}_2\cdot {{p}}_3 =& \frac12 \left(M_2^2+M_3^2 + 2E_2E_3 -m_{23}^2 \right). \end{split} $$ (9) Before we move on, let us first show the impact of considering different partial waves for the
$ D_1D^*\pi $ coupling in the structure induced by the$ D_1\bar D D^*+c.c. $ triangle diagrams for the$ e^+e^-\to J/\psi\pi^+\pi^- $ reaction in Fig. 4. The sharp peak in the$ J/\psi\pi^+ $ energy distribution (the right band in the Dalitz plot) is due to the diagram in Fig. 3(a), and the broader peak (the left band in the Dalitz plot) is the kinematical reflection due to the diagram in Fig. 3(b). One sees that the pion-momentum dependence at the$ D_1D^*\pi $ vertex has a significant impact on the distributions. The reason is that the magnitude of the three-momentum of a pion can be as large as about 1 GeV, and a larger background in the$ J/\psi\pi^+ $ distribution can be caused by the$ D $ -wave coupling than that by the$ S $ -wave one. This could be the main reason for the different conclusions reached in Refs. [23, 26]. One also sees that the$ D $ -wave coupling also leads to a double-bump structure in the$ \pi^+\pi^- $ invariant mass distribution.Figure 4. (color online) The Dalitz plot and its projections to the
$J/\psi\pi^+$ and$\pi^+\pi^-$ energy distributions for the$e^+e^-\to J/\psi\pi^+\pi^-$ reaction induced by the$D_1\bar D D^*$ triangle diagrams. The$e^+e^-$ c.m. energy is taken at 4.26 GeV. Left:$D$ -wave$D_1D^*\pi$ coupling; middle:$S$ -wave$D_1D^*\pi$ coupling; Right: the coupling contains both the$S$ -wave and$D$ -wave parts as given in Sec. 2.2.Next, let us consider all the THH-type diagrams. For accounting for the whole set of the THH-type triangle diagrams, one needs to decide on the relative couplings between
$ D_1\bar D+c.c. $ ,$ D_1D^*+c.c. $ and$ D_2\bar D^*+c.c. $ with the initial state. In principal, for an analysis of the experimental data, one may assume them to be independent. This is because at different$ e^+e^- $ c.m. energies, the charmonium or charmonium-like state that is important in that energy region could be different internal structures. This makes difficult the use of HQSS to relate these couplings. Then, we can write the amplitude for Fig. 3(a) considering all theTHH-type diagrams as:$$ \begin{split} {\cal{A}}_a = &{\rm{i}} \Bigg\{ 10h_S\, E_3 {{\varepsilon}}_Y\cdot {{\varepsilon}}_\psi\times \\ &\left[ c_1 I\left(D_1,\bar{D},D^{{*}},s,M_3^2, m_{12}^2\right)T_*(m_{12}^2) \right.+\\ & \left. c_2 I\left(D_1,\bar{D}^{{*}},D^{{*}}, s,M_3^2, m_{12}^2\right) T_{**}(m_{12}^2)\right]+ \\ & h \left(3 {{p}}_3\cdot \varepsilon_Y {{p}}_3\cdot {{\varepsilon}}_\psi-{{p}}_3^{\;2} {{\varepsilon}}_Y\cdot {{\varepsilon}}_\psi\right)\times \\ &\left[10 c_1 I\left(D_1,\bar{D},D^{{*}},s,M_3^2, m_{12}^2\right) T_*(m_{12}^2)\right.- \\ & \left. 5 c_2 I\left(D_1,\bar{D}^{{*}},D^{{*}},s,M_3^2, m_{12}^2\right)T_{**}(m_{12}^2) \right.- \\ &3 c_3 I\left(D_2,\bar{D}^{{*}},D^{{*}},s,M_3^2, m_{12}^2\right)T_{**}(m_{12}^2)- \\ & \left.2 c_3 I\left(D_2,\bar{D}^{{*}},D,s,M_3^2, m_{12}^2\right)T_{*}(m_{12}^2)\right] \Bigg\}\\ \equiv & B_{S,a} {{\varepsilon}}_Y\cdot {{\varepsilon}}_\psi + B_{D,a} \left(3 {{p}}_3\cdot {{\varepsilon}}_Y {{p}}_3\cdot {{\varepsilon}}_\psi-{{p}}_3^{\;2} {{\varepsilon}}_Y\cdot {{\varepsilon}}_\psi\right) , \end{split} $$ (10) with
$ M_3\!=\!M_2\!=\!M_\pi $ . Here, we use$ c_{1,2} $ and$ c_3 $ to account for the couplings of the$ D_1\bar D $ ,$ D_1 D^* $ and$ D_2\bar D^* $ pairs to the initial state, respectively, and use the intermediate particles to represent the corresponding$ m_i^2 $ .$ T_* $ and$ T_{**} $ represent the$ T $ -matrix elements (with the polarization vectors amputated) for the rescattering processes$ \bar D D^*\to\psi\pi^+ $ and$ \bar D^* D^*\to\psi\pi^+ $ , respectively. The amplitude for the diagrams with charge-conjugated intermediate particles,$ {\cal{A}}_b $ , is obtained by replacing$ M_3 $ by$ M_2 $ ,$ {{p}}_3 $ by$ {{p}}_2 $ ,$ E_3 $ by$ E_2 $ , and$ m_{12} $ by$ m_{13} $ :$$ {\cal{A}}_b = B_{S,b} {{\varepsilon}}_Y\cdot {{\varepsilon}}_\psi + B_{D,b} \left(3 {{p}}_2\cdot {{\varepsilon}}_Y {{p}}_2\cdot {{\varepsilon}}_\psi-{{p}}_2^{\;2} {{\varepsilon}}_Y\cdot {{\varepsilon}}_\psi\right). $$ (11) The sum
$ {\cal{A}}_a + {\cal{A}}_b $ , after having parameterized$ T_* $ and$ T_{**} $ for the$ D^{(*)}\bar D^{(*)}\to \psi\pi $ rescattering which allows for the existence of a$ Z_c $ pole, may be used to account for the triangle singularity effects in the analysis of the data for the$ e^+e^-\to \psi\pi^+\pi^- $ reaction.In the following, we choose specific ratios among
$ c_{1,2,3} $ to show the triangle-singularity induced structures. Here, we assume that the initial vector source couples to the${\rm T\bar{H}}+c.c. $ pairs like a$ D $ -wave ($ j_\ell\!=\!2 $ ) charmonium. That is to assume that the HQSS breaking happens at the charmonium level, instead of at the charmed meson level. The reason for this choice is that the$ D_1 $ width is small so that the$ D_{1(2)} $ mesons are well approximated by$ j_\ell^P\!=\!\frac32^+ $ mesons and the$ S $ -$ D $ mixing for the$ D_1 $ mesons should be small (the mixing angle was determined to be$ -0.10\pm0.04 $ rad by Belle[50], noting however that the broad$ D_1 $ extracted from that reference is problematic as discussed in Refs. [48-49]), while higher charmonia typically have a large$ S $ -$ D $ mixing as can be seen from the fact that their dileptonic decay widths do not differ much[44] (see the discussion in Ref. [51]; the$ S $ -$ D $ mixing angle for vector charmonia above 4 GeV could be as large as about 34$ ^\circ\!=\!0.59 $ rad[52]). This assumption amounts to take$ c_1\!=\!c_2\!=\!c_3 $ in Eq. (10). Then, the absolute value squared of the amplitude reads as$$ \begin{split} |{\cal{A}}_a \!+\!{\cal{A}}_b|^2 \!=\!& \, 6 \Bigg\{ \frac12\left(|B_{S,a}|^2 \!+\! |B_{S,b}|^2\right) \!+\! {\rm{Re}} \left(B_{S,a} B_{S,b}^*\right)\!+ \\ & |B_{D,a}|^2{{p}}_3^{\;4} + |B_{D,b}|^2{{p}}_2^{\;4} + {\rm{Re}} \left(B_{D,a} B_{D,b}^*\right)\times \\ &\left[ 3({{p}}_2\cdot {{p}}_3)^2 - {{p}}_2^{\;2} {{p}}_3^{\;2}\right] \Bigg\}.\\[-13pt] \end{split} $$ (12) In order to see clearly the triangle singularity effects, we also switch off any nontrivial structure in the rescattering matrix elements
$ T_* $ and$ T_{**} $ by setting them to the same constant. The resulting Dalitz plot distributions for$ e^+e^-\to \psi(3686)\pi^+\pi^- $ at four different c.m. energy values,$ \sqrt{s}\!=\!4.226 $ , 4.258, 4.358 and 4.416 GeV, are shown in Fig. 5. We also show the BESIII data of the$ \psi(3686)\pi $ invariant mass distributions reported in Ref. [53] for a rough comparison. Notice that we did not perform a fit to the data. The sensitivity of the triangle-singularity induced structure on the energy is evident, as already observed for this reaction in Ref. [19]. For a fit to the data, which is beyond the scope of this special issue, one needs to parameterize the rescattering$ T $ -matrix, which contains possible$ Z_c $ states, and treat$ c_{1,2,3} $ as free parameters for each energy point. In addition, one also needs to include contributions other than the triangle diagrams such as the direct production of$ \psi(3686)\pi^+\pi^- $ and the two-body rescattering (including the$ \pi\pi $ final state interaction and the$ T_* $ and$ T_{**} $ as those in Eq. (10)). The BESIII data for this reaction reported in Ref. [53] was analyzed in Ref. [54] without considering the triangle singularities.Figure 5. (color on line) The Dalitz plot and its projections to the
$\psi(3686)\pi^+$ and$\pi^+\pi^-$ energy distributions for the$e^+e^-\to \psi(3686)\pi^+\pi^-$ reaction induced by the THH-type triangle diagrams. The amplitude is given by the sum of Eqs. (10) and (11) with$T_*$ and$T_{**}$ set to the same constant and$c_1\!=\!c_2\!=\!c_3$ . The BESIII data of the$\psi(3686)\pi$ invariant mass distributions[53] are shown for a rough comparison. -
The triangle singularity effects need to be properly taken into account in order to establish the exotic hadron spectrum and extract the resonance parameters more reliably. In this paper, we present the formalism of considering the triangle singularities that can be used in the experimental analysis of the
$ e^+e^-\to \psi\pi^+\pi^- $ data. Notice that in a complete analysis, the triangle diagrams are not supposed to provide the only mechanism for the reactions. For instance, there can be a direct production of the$ J/\psi\pi\pi $ and$ D^{(*)}\bar D^*\pi $ , followed by$ \pi\pi $ and$ D^{(*)}\bar D^* $ -$ J/\psi\pi $ final state interactions. While the$ D^{(*)}\bar D^* $ -$ J/\psi\pi $ final state interactions are the same as$ T_* $ and$ T_{**} $ in Eq. (10), the$ \pi\pi $ final state interaction also needs to be taken into account, which may be done by using the Omnès dispersive formalism as that used in Refs. [54-55].
-
摘要: 对强子谱的研究有助于加深量子色动力学中色禁闭现象的理解。强子谱中存在诸多疑难问题,特别是2003年以来各高能物理实验上发现了很多难以理解的新的强子共振结构,其中大量是在重夸克偶素质量区间发现的,它们被统称为
$XYZ$ 结构。强子共振态对应于$S$ 矩阵的极点。除此之外,$S$ 矩阵还具有其它类型的奇点,比如由于三个中间态粒子同时在 壳导致的三角奇点。本文中,我们将讨论三角奇点在$XYZ$ 物理中的实现,注重于构造可以用来分析 其对$e^+e^-\to\psi\pi^+\pi^-$ 反应的$\psi\pi$ 不变质量谱中的带电$Z_c$ 结构的影响的振幅。-
关键词:
- Exotic hadrons /
- XYZ states /
- triangle singularity
Abstract: The spectrum of hadrons is important for understanding the confinement of quantum chromodynamics. Many new puzzles arose since 2003 due to the abundance of experimental discoveries with the$XYZ$ structures in the heavy quarkonium mass region being the outstanding examples. Hadronic resonances correspond to poles of the$S$ -matrix, which has other types of singularities such as the triangle singularity due to the simultaneous on-shellness of three intermediate particles. Here we briefly discuss a few possible manifestations of triangle singularities in the$XYZ$ physics, paying particular attention to the formalism that can be used to analyze the data for charged$Z_c$ structures in the$\psi\pi$ distributions of the reaction$e^+e^-\to \psi\pi^+\pi^-$ .-
Key words:
- Exotic hadrons /
- XYZ states /
- triangle singularity
$ D_0^*(2400) $ was considered in Ref. [21] instead of the narrow$ D_1(2420) $ .
-
Figure 4. (color online) The Dalitz plot and its projections to the
$J/\psi\pi^+$ and$\pi^+\pi^-$ energy distributions for the$e^+e^-\to J/\psi\pi^+\pi^-$ reaction induced by the$D_1\bar D D^*$ triangle diagrams. The$e^+e^-$ c.m. energy is taken at 4.26 GeV. Left:$D$ -wave$D_1D^*\pi$ coupling; middle:$S$ -wave$D_1D^*\pi$ coupling; Right: the coupling contains both the$S$ -wave and$D$ -wave parts as given in Sec. 2.2.Figure 5. (color on line) The Dalitz plot and its projections to the
$\psi(3686)\pi^+$ and$\pi^+\pi^-$ energy distributions for the$e^+e^-\to \psi(3686)\pi^+\pi^-$ reaction induced by the THH-type triangle diagrams. The amplitude is given by the sum of Eqs. (10) and (11) with$T_*$ and$T_{**}$ set to the same constant and$c_1\!=\!c_2\!=\!c_3$ . The BESIII data of the$\psi(3686)\pi$ invariant mass distributions[53] are shown for a rough comparison. -
[1] CHEN H X, CHEN W, LIU X, et al. Phys Rept, 2016, 639: 1. doi: 10.1016/j.physrep.2016.05.004 [2] ESPOSITO A, PILLONI A, POLOSA A D. Phys Rept, 2017, 668: 1. doi: 10.1016/j.physrep.2016.11.002 [3] HOSAKA A, ⅡJIMA T, MIYABAYASHI K, et al. PTEP, 2016, 2016: 062C01. [4] RICHARD J M. Few Body Syst, 2016, 57: 1185. doi: 10.1007/s00601-016-1159-0 [5] LEBED R F, MITCHELL R E, SWANSON E S. Prog Part Nucl Phys, 2017, 93: 143. doi: 10.1016/j.ppnp.2016.11.003 [6] GUO F K, HANHART C, MEISSNER U G, et al. Rev Mod Phys, 2018, 90: 015004. doi: 10.1103/RevModPhys.90.015004 [7] ALI A, LANGE J S, STONE S. Prog Part Nucl Phys, 2017, 97: 123. doi: 10.1016/j.ppnp.2017.08.003 [8] OLSEN S L, SKWARNICKI T, ZIEMINSKA D. Rev Mod Phys, 2018, 90: 015003. doi: 10.1103/RevModPhys.90.015003 [9] KOU E, URQUIJO P , ALTMANNSHOFER W, et al(Belle-II Collaboration). PTEP, 2019, 2019: 123C01. Erratum: PTEP, 2020, 2020: 029201. [10] CERRI A, GLIGOROV V V, MALVEZZI S, et al. CERN Yellow Rep Monogr, 2019, 7: 867. doi: 10.23731/CYRM-2019-007.867 [11] LIU Y R, CHEN H X, CHEN W, et al. Prog Part Nucl Phys, 2019, 107: 237. doi: 10.1016/j.ppnp.2019.04.003 [12] BRAMBILLA N, EIDELMAN S, HANHART C, et al. Phys Rept, 2020, 873: 1. doi: 10.1016/j.physrep.2020.05.001 [13] GUO F K, LIU X H, SAKAI S. Prog Part Nucl Phys, 2020, 112: 103757. doi: 10.1016/j.ppnp.2020.103757 [14] COLEMAN S, NORTON R E. Nuovo Cim, 1965, 38: 438. doi: 10.1007/BF02750472 [15] LANDAU L D. Nucl Phys, 1959, 13: 181. doi: 10.1016/0029-5582(59)90154-3 [16] WANG Q, HANHART C, ZHAO Q. Phys Rev Lett, 2013, 111: 132003. doi: 10.1103/PhysRevLett.111.132003 [17] WANG Q, HANHART C, ZHAO Q. Phys Lett B, 2013, 725: 106. doi: 10.1016/j.physletb.2013.06.049 [18] LIU X H, LI G. Phys Rev D, 2013, 88: 014013. doi: 10.1103/PhysRevD.88.014013 [19] LIU X H. Phys Rev D, 2014, 90: 074004. doi: 10.1103/PhysRevD.90.074004 [20] LIU X H, OKA M. Phys Rev D, 2016, 93: 054032. doi: 10.1103/PhysRevD.93.054032 [21] SZCZEPANIAK A P. Phys Lett B, 2015, 747: 410. doi: 10.1016/j.physletb.2015.06.029 [22] SZCZEPANIAK A P. Phys Lett B, 2016, 757: 61. doi: 10.1016/j.physletb.2016.03.064 [23] ALBALADEJO M, GUO F K, HIDALGO-DUQUE C, et al. Phys Lett B, 2016, 755: 337. doi: 10.1016/j.physletb.2016.02.025 [24] LIU X H. Phys Lett B, 2017, 766: 117. doi: 10.1016/j.physletb.2017.01.008 [25] BONDAR A E, VOLOSHIN M B. Phys Rev D, 2016, 93: 094008. doi: 10.1103/PhysRevD.93.094008 [26] PILLONI A, FERNÁNDEZ-RAMÍREZB C, JACKURAC A, et al( JPAC Collaboration). Phys Lett B, 2017, 772: 200. doi: 10.1016/j.physletb.2017.06.030 [27] GONG Q R, PANG J L, WANG Y F, et al. Eur Phys J C, 2018, 78: 276. doi: 10.1140/epjc/s10052-018-5690-7 [28] GUO F K. Phys Rev Lett, 2019, 122: 202002. doi: 10.1103/PhysRevLett.122.202002 [29] BRAATEN E, HE L P, INGLES K. Phys Rev D, 2019, 100: 031501. doi: 10.1103/PhysRevD.100.031501 [30] BRAATEN E, HE L P, INGLES K. Phys Rev D, 2020, 101: 014021. doi: 10.1103/PhysRevD.101.014021 [31] BRAATEN E, HE L P, INGLES K. Phys Rev D, 2019, 100: 074028. doi: 10.1103/PhysRevD.100.074028 [32] BRAATEN E, HE L P, INGLES K. Phys Rev D, 2019, 100: 094006. doi: 10.1103/PhysRevD.100.094006 [33] ACHASOV N N, SHESTAKOV G N. Phys Rev D, 2019, 99: 116023. doi: 10.1103/PhysRevD.99.116023 [34] NAKAMURA S X, TSUSHIMA K. Phys Rev D, 2019, 100: 051502. doi: 10.1103/PhysRevD.100.051502 [35] NAKAMURA S X. Phys Rev D, 2019, 100: 011504. doi: 10.1103/PhysRevD.100.011504 [36] NAKAMURA S X. arXiv: 1912.11830[hep-ph]. [37] ABLIKIM M, ACHASOV M N, AI X C, et al(BESIII Collaboration). Phys Rev Lett, 2013, 110: 252001. doi: 10.1103/PhysRevLett.110.252001 [38] LIU Z Q, SHEN C P, YUAN C Z, et al(Belle Collaboration). Phys Rev Lett, 2013, 110: 252002. Erratum: [Phys Rev Lett, 2013, 111: 019901]. [39] ABLIKIM M, ACHASOV M N,ALBAYRAK O, et al(BESIII Collaboration). Phys Rev Lett, 2014, 112: 022001. doi: 10.1103/PhysRevLett.112.022001 [40] ABLIKIM M, ACHASOV M N,ALBAYRAK O, et al(BESIII Collaboration). Phys Rev D, 2015, 92: 092006. doi: 10.1103/PhysRevD.92.092006 [41] ABLIKIM M, ACHASOV M N, AI X C, et al(BESIII Collaboration). Phys Rev Lett, 2017, 119: 072001. doi: 10.1103/PhysRevLett.119.072001 [42] ABLIKIM M, ACHASOV M N, ALBAYRAK O, et al(BESIII Collaboration). Phys Rev Lett, 2014, 112: 132001. doi: 10.1103/PhysRevLett.112.132001 [43] ABLIKIM M, ACHASOV M N, ALBAYRAK O, et al(BESIII Collaboration). Phys Rev Lett, 2013, 111: 242001. doi: 10.1103/PhysRevLett.111.242001 [44] TANABASHI M. Phys Rev D, 2018, 98: 030001. doi: 10.1103/PhysRevD.98.030001 [45] CASALBUONI R, DEANDREA A, BARTOLOMEO N D, et al. Phys Rept, 1997, 281: 145. doi: 10.1016/S0370-1573(96)00027-0 [46] COLANGELO P, DE FAZIO F, FERRANDES R. Phys Lett B, 2006, 634: 235. doi: 10.1016/j.physletb.2006.01.021 [47] LI X, VOLOSHIN M B. Phys Rev D, 2013, 88: 034012. doi: 10.1103/PhysRevD.88.034012 [48] DU M L, ALBALADEJO M, FERNÁNDEZ-SOLER P, et al. Phys Rev D, 2018, 98: 094018. doi: 10.1103/PhysRevD.98.094018 [49] DU M L, GUO F K, MEISSNER U G. Phys Rev D, 2019, 99: 114002. doi: 10.1103/PhysRevD.99.114002 [50] ABE K, ABE T,ADACHI I, et al(Belle Collaboration). Phys Rev D, 2004, 69: 112002. doi: 10.1103/PhysRevD.69.112002 [51] WANG Q, CLEVEN M, GUO F K, et al. Phys Rev D, 2014, 89: 034001. doi: 10.1103/PhysRevD.89.034001 [52] BADALIAN A M, BAKKER B L G, DANILKIN I V. Phys Atom Nucl, 2009, 72: 638. doi: 10.1134/S1063778809040085 [53] ABLIKIM M, ACHASOV M N, AHMED S, et al(BESIII Collaboration). Phys Rev D, 2017, 96: 032004. Erratum: [Phys Rev D, 2019, 99: 019903]. [54] MOLNAR D A S, DANILKIN I, VANDERHAEGHEN M. Phys Lett B, 2019, 797: 134851. doi: 10.1016/j.physletb.2019.134851 [55] CHEN Y H, DAI L Y, GUO F K, et al. Phys Rev D, 2019, 99: 074016. doi: 10.1103/PhysRevD.99.074016 -