Charged charmoniumlike structures in the $e^+ e^- \to \psi(3686) \pi^+ \pi^-$ process based on the ISPE mechanism

TL;DR

This paper demonstrates that the initial single pion emission (ISPE) mechanism can accurately reproduce the charged charmonium-like structure observed in the $e^+ e^- o \psi(3686) \\pi^+ \\pi^-$ process, providing a plausible explanation for this phenomenon.

Contribution

It offers the first detailed analysis showing the ISPE mechanism can explain the charged charmonium-like structure in this process.

Findings

The charged structure can be reproduced by the ISPE mechanism.

The dipion invariant mass spectrum matches experimental data.

The cross section is consistent with observations.

Abstract

In 2017, the BESIII Collaboration announced the observation of a charged charmonium-like structure in the $\psi(3686)\pi^\pm$ invariant mass spectrum of the $e^+ e^- \to \psi(3686) \pi^+ \pi^-$ process at different energy points, which enables us to perform a precise study of this process based on the initial single pion emission (ISPE) mechanism. In this work, we perform a combined fit to the experimental data of the cross section of $e^+ e^- \to \psi(3686) \pi^+ \pi^-$, and the corresponding $\psi(3686)\pi^\pm$ and dipion invariant mass spectra. Our result shows that the observed charged charmonium-like structure in $e^+ e^- \to \psi(3686) \pi^+ \pi^-$ can be well reproduced based on the ISPE mechanism, and that the corresponding dipion invariant mass spectrum and cross section can be depicted with the same parameters. In fact, it provides strong evidence that the ISPE mechanism can be an underlying mechanism resulting in such novel a phenomenon.

Tables

Table 1. Table 1: The resonance parameters and dilepton decay widths of the involved states Tanabashi:2018oca . The dilepton decay widths of ψ ​ ( 4220 ) 𝜓 4220 \psi(4220) and ψ ​ ( 4380 ) 𝜓 4380 \psi(4380) are taken from Ref. Wang:2019

Meson	Mass (GeV)	Width (MeV)	${\Gamma }_{e^{+}{e}^{-}}$ (keV)
$\psi \left(4160\right)$	4.191	70	0.48
$\psi \left(4415\right)$	4.421	62	0.58
$\psi \left(4220\right)$	4.218	59	0.29
$\psi \left(4380\right)$	4.384	84	0.26

Table 2. Table 2: The fitted values of the parameters in three-charmonium scenario.

Prameters	$\psi \left(4160\right)$	$\psi \left(4220\right)$	$\psi \left(4415\right)$
$a_R$ (GeV-1)	6.0	5.0	4.0
${\varphi }_{\mathrm{Dir}}$	-1.631	-1.506	1.762
$F_{\mathrm{Dir}}$	-3.477	5.881	0.982
${\kappa }_{Dir}$	-1.829	-0.025	-0.918
${\varphi }_{\sigma }$	1.130	2.009	2.241
$f_{\sigma }$ (GeV2)	7.562	5.785	-3.971
$g_{\sigma }$	77.498	21.803	-1.920
${\phi }_{\sigma }$	1.186	-2.426	-1.923
${\varphi }_{\mathrm{ISPE}}$	-1.320	2.980	-2.813
$g_{RD\overline{D}\pi }$ (GeV-3)	1.717	1.574	0.002
$g_{R{D}^{\ast }\overline{D}\pi }$	-1.834	-0.670	0.025
$g_{R{D}^{\ast }{\overline{D}}^{\ast }\pi }$ (GeV-1)	0.120	-0.096	0.105
$g_{\mathrm{NoR}}$=18.940 GeV	$a_{\mathrm{NoR}}$=4.702 GeV-2
${\chi }^2/d.o.f.=3.704$

Table 3. Table 3: Branching ratios of ψ R → π + ​ π − ​ ψ ​ ( 3686 ) → subscript 𝜓 𝑅 superscript 𝜋 superscript 𝜋 𝜓 3686 \psi_{R}\to\pi^{+}\pi^{-}\psi(3686) obtained from the present fit in the unit of 10 − 3 superscript 10 3 10^{-3} .

	$\psi \left(4160\right)$	$\psi \left(4220\right)$	$\psi \left(4415\right)$	$\psi \left(4380\right)$
3R Fit	15.711	8.101	9.65	-
4R Fit	4.389	7.370	5.368	17.285

Table 4. Table 4: The same as Table 2 but for the four-charmonium scenario.

Parameter	$\psi \left(4160\right)$	$\psi \left(4220\right)$	$\psi \left(4415\right)$	$\psi \left(4380\right)$
$a_R$ (GeV-1)	5.5	4.0	3.3	3.0
${\varphi }_{\mathrm{Dir}}$	-0.146	1.661	-0.423	-0.345
$F_{\mathrm{Dir}}$	-12.243	3.361	0.504	1.523
${\kappa }_{\mathrm{Dir}}$	-0.151	-0.215	-1.190	-0.463
${\varphi }_{\sigma }$	0.517	-0.950	2.910	-0.328
$f_{\sigma }$ (GeV2)	11.953	1.813	1.832	-5.940
${\phi }_{\sigma }$	-1.680	0.287	-0.838	-1.455
$g_{\sigma }$	58.717	-6.612	-7.864	-25.906
${\varphi }_{\mathrm{ISPE}}$	-1.759	2.875	1.396	-5.475
$g_{RD\overline{D}\pi }$ (GeV-3)	-2.223	-1.736	-0.257	0.276
$g_{R{D}^{\ast }\overline{D}\pi }$	1.646	0.353	-0.470	1.249
$g_{R{D}^{\ast }{\overline{D}}^{\ast }\pi }$ (GeV-1)	0.323	0.131	0.133	-0.022
$g_{\mathrm{NoR}}$=-4.847 GeV		$a_{\mathrm{NoR}}$=2.058 GeV-2
${\chi }^2/d.o.f.=2.767$

Equations

$${\mathcal{M}}_{\mathrm{NoR}}^{(1)}=\bar{v}(k_2){\gamma }^{\mu }u(k_1)g_{\mathrm{NoR}}e\frac{-g_{\mu \nu }}{s}{ϵ}_{{\psi }^{\prime }}^{\ast \nu }{\mathcal{F}}_{\mathrm{NoR}}(s),$$

$${\mathcal{F}}_{\mathrm{NoR}}(s)=e^{-a_{\mathrm{NoR}}(\sqrt{s}-\Sigma m_f{)}^2},$$

$${\mathcal{M}}_R^{(i)}$$

$${\mathcal{F}}_R(s)=e^{-a_R\left(\sqrt{s}-m_R\right)}.$$

$${\mathcal{A}}_{\mu \nu }^{(2)}$$

$${\mathcal{A}}_{\mu \nu }^{(3)}$$

$${\Gamma }_{\sigma }(m_{\pi \pi })={\Gamma }_{\sigma }\frac{m_{\sigma }}{m_{\pi \pi }}\frac{\left|\vec{P}(m_{\pi \pi })\right|}{\left|\vec{P}(m_{\sigma })\right|}$$

$${\mathcal{A}}_{\mu \nu }^{(4)}$$

$${\mathcal{M}}^{\mathrm{Tot}}={\mathcal{M}}_{\mathrm{NoR}}+\sum _{R,i=2}^{i=4}{e}^{i{\varphi }_R^{(i)}}{\mathcal{M}}_R^{(i)}$$

$$d{\sigma }_{e^{+}{e}^{-}\to \psi (3686){\pi }^{+}{\pi }^{-}}=\frac{(2\pi {)}^4\sum {\left|\overline{{\mathcal{M}}^{\mathrm{Tot}}}\right|}^2}{4(k_1\cdot k_2)}d{\Phi }_3,$$

$$d{\Phi }_3=\frac{1}{32(2\pi {)}^{8}s}d{m}_{12}^{2}d{m}_{23}^{2}d(\cos \theta )d\eta ,$$

$${\mathcal{A}}_{\mu \zeta }$$

$${\mathcal{L}}_{\psi D^{(\ast )}{D}^{(\ast )}\pi }$$

$$=-i{g}_{RDD\pi }{\epsilon }^{\mu \nu \alpha \beta }{\psi }_{\mu }{\partial }_{\nu }D{\partial }_{\alpha }\pi {\partial }_{\beta }\bar{D}+g_{R{D}^{\ast }D\pi }{\psi }^{\mu }(D\pi {\bar{D}}_{\mu }^{\ast }+D_{\mu }^{\ast }\pi \bar{D})$$

$$-i{g}_{R{D}^{\ast }{D}^{\ast }\pi }{\epsilon }^{\mu \nu \alpha \beta }{\psi }_{\mu }{D}_{\nu }^{\ast }{\partial }_{\alpha }\pi {\bar{D}}_{\beta }^{\ast }-i{h}_{R{D}^{\ast }{D}^{\ast }\pi }{\epsilon }^{\mu \nu \alpha \beta }{\partial }_{\mu }{\psi }_{\nu }{D}_{\alpha }^{\ast }\pi {\bar{D}}_{\beta }^{\ast }.$$

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

$$=i{g}_{D^{\ast }D\pi }(D_{\mu }^{\ast }{\partial }^{\mu }\pi \bar{D}-D{\partial }^{\mu }\pi {\bar{D}}_{\mu }^{\ast })-g_{D^{\ast }{D}^{\ast }\pi }{\epsilon }^{\mu \nu \alpha \beta }{\partial }_{\mu }{D}_{\nu }^{\ast }\pi {\partial }_{\alpha }{\bar{D}}_{\beta }^{\ast }.$$

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

$$=i{g}_{\psi DD}{\psi }_{\mu }({\partial }^{\mu }D\bar{D}-D{\partial }^{\mu }\bar{D})-g_{\psi D^{\ast }D}{\epsilon }^{\mu \nu \alpha \beta }{\partial }_{\mu }{\psi }_{\nu }({\partial }_{\alpha }{D}_{\beta }^{\ast }\bar{D}$$

$$\displaystyle\quad+D\partial_{\alpha}\bar{D}^{\ast}_{\beta})-ig_{\psi D^{\ast}D^{\ast}}\big{\{}\psi^{\mu}(\partial_{\mu}D^{\ast\nu}\bar{D}^{\ast}_{\nu}-D^{\ast\nu}\partial_{\mu}\bar{D}^{\ast}_{\nu})$$

$$\displaystyle\quad+(\partial_{\mu}\psi_{\nu}D^{\ast\nu}-\psi_{\nu}\partial_{\mu}D^{\ast\nu})\bar{D}^{\ast\mu}+D^{\ast\mu}(\psi^{\nu}\partial_{\mu}\bar{D}^{\ast}_{\nu}-\partial_{\mu}\psi^{\nu}\bar{D}^{\ast}_{\nu})\big{\}}.$$

$$\frac{2g_{\pi }}{f_{\pi }}=g_{D^{\ast }{D}^{\ast }\pi }=\frac{g_{D^{\ast }D\pi }}{\sqrt{m_D{m}_{D^{\ast }}}},$$

$$g_{\psi DD}=g_{\psi D^{\ast }{D}^{\ast }}\frac{m_D}{m_{D^{\ast }}}=g_{\psi D^{\ast }D}\sqrt{m_D{m}_{D^{\ast }}}=\frac{m_{\psi }}{f_{\psi }}.$$

$${\mathcal{A}}_{\mu \zeta }^{(a)}$$

$${\mathcal{A}}_{\mu \zeta }^{(b)}$$

$${\mathcal{A}}_{\mu \zeta }^{(c)}$$

$${\mathcal{A}}_{\mu \zeta }^{(d)}$$

$${\mathcal{A}}_{\mu \zeta }^{(e)}$$

$${\mathcal{A}}_{\mu \zeta }^{(f)}$$

Full text

I Introduction

Since 2003, study on the charmoniumlike $XYZ$ states has become a hot spot of hadron physics and particle physics. Since it has the close relation to non-perturbative behavior of quantum chromodynamics (QCD), the relevant investigation of the charmoniumlike $XYZ$ states is helpful to deepen our understanding of strong interaction. In the past sixteen years, there have been extensive discussions about this issue (see reviews Liu:2013waa ; Chen:2016qju ; Liu:2019zoy ; Guo:2017jvc for more details).

To solve the puzzling phenomena existing in the hidden-bottom dipion decays of $\Upsilon(10860)$, the initial single pion emission (ISPE) mechanism was proposed in Refs. Chen:2011pv ; Liu:2011ky ; Chen:2012yr , where the charged $Z_{b}(10610)$ and $Z_{b}(10650)$ structures can be qualitatively reproduced. In 2012, considering the similarity between the bottomonium and charmonium families, the ISPE mechanism was applied to study the hidden-charm dipion decays of higher charmonia and charmoniumlike state $Y(4260)$, where the predicted charged charmoniumlike structures near $D\bar{D}^{*}$ or $D^{*}\bar{D}^{*}$ threshold may exist in the corresponding $J/\psi\pi^{+}$, $\psi(3686)\pi^{+}$, and $h_{c}\pi^{+}$ invariant mass spectra Chen:2011xk . These predictions given by the ISPE mechanism have provided valuable information to search for such charged charmoniumlike structures.

In 2013, the BESIII and Belle collaborations indeed discovered a charged charmonium-like $Z_{c}(3900)$ in the process $e^{+}e^{-}\to J/\psi\pi^{+}\pi^{-}$ around $\sqrt{s}=4.26\ \mathrm{GeV}$ Ablikim:2013mio ; Liu:2013dau , which was quickly and further confirmed by CLEO-c in the same process but at $\sqrt{s}=4.17$ GeV Xiao:2013iha . Later, another charged charmoniumlike structure named $Z_{c}(4020)$ was observed in the $h_{c}\pi^{+}$ invariant mass spectrum of the $e^{+}e^{-}\to\pi^{+}\pi^{-}h_{c}$ process Ablikim:2013wzq . Additionally, at the center-of-mass energy of 4.26 GeV, the BESIII Collaboration reported two charged charmoniumlike structures, $Z_{c}(3885)$ in the $(D^{\ast}\bar{D})^{\pm}$ invariant mass spectrum of $e^{+}e^{-}\to\pi^{\mp}(D^{\ast}\bar{D})^{\pm}$ process Ablikim:2013xfr and $Z_{c}(4025)$ in the $(D^{\ast}\bar{D}^{\ast})^{\pm}$ invariant mass spectrum of the $e^{+}e^{-}\to\pi^{\mp}(D^{\ast}\bar{D}^{\ast})^{\pm}$ process Ablikim:2013emm .

The story of finding charged charmoniumlike structures is continued. In 2017, the BESIII experiment made a progress on observing charged charmoniumlike structure in the $\psi(3686)\pi^{+}$ invariant mass spectrum by analyzing the data of cross section of the $e^{+}e^{-}\to\psi(3686)\pi^{+}\pi^{-}$ process at different energy points $\sqrt{s}=4.226,\,4.258,\,4.358,\,4.387,\,4.416,$ and $4.600$ GeV Ablikim:2017oaf . In fact, this observation deepens our confidence again. As shown in Fig. 1, the predicted charged charmoniumlike structures in the $J/\psi\pi^{+}$, $\psi(3686)\pi^{+}$, and $h_{c}\pi^{+}$ invariant mass spectrum have been reported in experiments until now.

The BESIII’s experimental measurement of the $e^{+}e^{-}\to\psi(3686)\pi^{+}\pi^{-}$ process forces us to carry out a further study of this process based on the ISPE mechanism. Before the present work, in Ref. Chen:2013coa , the similar idea was once adopted to study the $e^{+}e^{-}\to J/\psi\pi^{+}\pi^{-}$ process combined with the BESIII data of the $J/\psi\pi^{\pm}$ invariant mass spectrum and the corresponding dipion invariant mass spectrum of $e^{+}e^{-}\to J/\psi\pi^{+}\pi^{-}$ at one energy point $\sqrt{s}=4.26$ GeV, which shows that the $Z_{c}(3900)$ structure can be reproduced well by the ISPE mechanism. Different from Ref. Chen:2013coa , the present work will perform a combined fit to the released data of the $\psi(3686)\pi^{\pm}$ and $\pi^{+}\pi^{-}$ invariant mass spectrum from the $e^{+}e^{-}\to\psi(3686)\pi^{+}\pi^{-}$ process at different energy points. Comparing with the former work Chen:2013coa , the present work is obviously a tough task and is full of challenge. In this work, we will reproduce the structures observed in the $e^{+}e^{-}\to\psi(3686)\pi^{+}\pi^{-}$ process step by step based on the ISPE mechanism. It also reveals the relation between the observed charged charmoniumlike structures and the ISPE mechanism, which will be helpful to deepen our understanding on such a novel phenomenon.

This work is organized as follows. After introduction, we perform an analysis of the mechanisms working in the $e^{+}e^{-}\to\psi(3686)\pi^{+}\pi^{-}$, and then fit the present model with the experimental data. In Sec. III, we present our fit results to the cross sections for $e^{+}e^{-}\to\psi(3686)\pi^{+}\pi^{-}$ and the mass spectra of $\psi(3686)\pi^{\pm}$ and dipion invariant mass spectra. Sec. IV is devoted to conclusions and discussion.

II Mechanisms occurring in $e^{+}e^{-}\to\psi(3686)\pi^{+}\pi^{-}$

In Ref. Chen:2013coa , we once reproduced the $Z_{c}(3900)$ structure in the $e^{+}e^{-}\to J/\psi\pi^{+}\pi^{-}$ process based on the ISPE mechanism. In the present work, we will apply the same idea to discuss $e^{+}e^{-}\to\psi(3686)\pi^{+}\pi^{-}$. There exist serval different decay mechanisms occurring in the process $e^{+}e^{-}\to\psi(3686)\pi^{+}\pi^{-}$ shown in Fig. 2. Fig. 2 (1) represents the nonresonance contribution, where the final states $\psi(3686)\pi^{+}\pi^{-}$ directly couple to the photon produced by the electron-positron annihilation. Fig. 2 (2)-(4) correspond to the processes with an intermediate charmonium resonance, where the virtual photon first couples to a vector charmonium $\psi_{R}$, which decays into $\psi(3686)\pi^{+}\pi^{-}$. There are three different kinds of decay mechanisms in the higher charmonium dipion decays: i) As shown in Fig. 2 (2), the first is the one in which the dipion is produced by the gluon hadronization with gluons emitted by a charm quark. ii) In the second decay mechanism, the dipion is produced through a scalar meson decay, while the scalar meson and $\psi(3686)$ are connected to the initial $\psi_{R}$ by a charmed meson loop. This is presented in Fig. 2 (3). iii) The last decay mechanism is the ISPE mechanism, which is presented in Fig. 2 (4). We will list resonaces $R$ inserted into Fig. 2 later in Table 1.

In the following, we will present the amplitudes corresponding to different mechanisms. For Fig. 2 (1), its amplitude is expressed as

[TABLE]

where $g_{\rm NoR}$ is a coupling constant, $e=\sqrt{4\pi\alpha}$ denotes an electric charge with $\alpha=1/137$ being the fine structure constant, and $\mathcal{F}_{\rm NoR}(s)$ is a form factor. In the present work, we adopt the following form factor Chen:2013coa

[TABLE]

where $a_{\rm NoR}$ is a phenomenological parameter, which should be fixed when fitting the data. $\Sigma m_{f}=m_{\psi^{\prime}}+2m_{\pi^{\pm}}$ is the sum of the masses of the involved particles appearing in the final states.

For Figs. 2 (2)-(4), a general amplitude, with a help of the vector meson dominance (VMD) ansatz Lin:2013mka ; Xu:2015qqa , can be written as

[TABLE]

where $f_{R}$ is a decay constant of an intermediate $R$ and $p_{0}$ denotes the momentum of the intermediate $R$, which satisfies $p_{0}^{2}=(k_{1}+k_{2})^{2}=s$. $\mathcal{A}_{\mu\nu}$ indicates the decay amplitude of $\psi_{R}\to\psi(3686)\pi^{+}\pi^{-}$. Here, we also introduce the form factor $\mathcal{F_{R}}$ to balance the otherwise over-increased decay width with an increased phase space, which isgiven by,

[TABLE]

where $a_{R}$ is a parameter obtained when fitting the data.

For Fig. 2 (2), the dipion is produced directly by the gluon emission. Hereafter, such a process is named as a direct production process and the decay amplitude is parameterized as Novikov:1980fa ,

[TABLE]

where the terms in the first and second square brackets indicate that the orbital momenta of the dipion are $S$-wave and $D$-wave, respectively. $\Delta M=m_{R}-m_{\psi^{\prime}}$ is the mass splitting of the initial $\psi_{R}$ and final $\psi(3686)$. $F_{\rm Dir}$ and $\kappa_{\rm Dir}$ are parameters necessary for fitting and $\alpha$ is the angle between $\psi_{R}$ and $\pi^{-}$ in the $\pi^{+}\pi^{-}$ rest frame.

As shown in Fig. 2 (3), the dipion is produced through a scalar meson, and the scalar meson as well as $\psi(3686)$ are connected to $\psi_{R}$ via a charmed meson loop. Such a kind of a meson loop has been proven to be important in the dipion transitions between heavy quarkonia Chen:2011qx ; Chen:2011jp . Due to the kinetic limit of the involved process, we take the $\sigma$ meson into consideration. Moreover, in the present work, we mainly aim at the ISPE mechanism. Thus, we will not directly estimate the triangle diagrams but parameterize their contributions in the form Chen:2013coa ,

[TABLE]

where $f_{\sigma}$ and $g_{\sigma}$ are the $S$ and $D$ wave coupling constants between $\psi_{R}$ and $\psi(3686)\mathcal{S}$, respectively, and $\varphi_{\sigma}$ is the phase angle between the $S$ and $D$ wave couplings. We should notice that the mass and width of the $\sigma$ meson are comparable, and thus, we adopt a momentum dependent width, which is,

[TABLE]

where $\left|\vec{P}(m_{\pi\pi})\right|=\sqrt{m_{\pi\pi}^{2}/4-m_{\pi}^{2}}$ is the pion momentum in the dipion rest frame, while $|\vec{P}(m_{\sigma})|$ is the pion momentum with the on-shell $\sigma$ meson.

As for the ISPE mechanism, we directly estimate the triangle diagrams. To fit with the experimental data, we reduce the amplitudes of the ISPE mechanism in the form,

[TABLE]

where the coefficients $C_{ij},\ \{(i,j)=(0,1,2,3)\}$ are obtained by the loop integrals of the amplitudes in Appendix A.

With above preparations, we get the total amplitude, which could be the coherent sum of the amplitudes corresponding to different mechanisms, and reads as

[TABLE]

with $\phi_{R}^{(i)}$ being the phase angles in front of different amplitudes for each resonance. For convenience, we use $\phi_{\rm Dir}$, $\phi_{\sigma}$ and $\phi_{\rm ISPE}$ to refer to these phase angles in the amplitudes of direct production process, $\sigma$ production process and ISPE processes, respectively. The differential cross section is Xie:2015zga

[TABLE]

where $\sum\left|\overline{\mathcal{M}^{\rm Tot}}\right|^{2}$ indicates an average over the spin of initial electron and positron and a sum over the spin of the final state $\psi(3686)$. $d\Phi_{3}$ is the phase space factor, which is expressed as

[TABLE]

where $m_{23}$ corresponds to $m_{\psi^{\prime}\pi^{+}}$, $m_{12}$ is $m_{\pi^{+}\pi^{-}}$, and the definitions of angles $\theta$ and $\eta$ are shown in Fig. 3.

III Numerical Results and discussion

In the differential cross section in Eq. (10), there are two free parameters, $g_{\rm NoR}$ and $a_{\rm NoR}$, in the nonresonance amplitude. As for the intermediate charmonium contributions, there are 12 additional free parameters for each resonance in Table 1, which are

[TABLE]

In the following, we need to fit our results with the experimental data of $e^{+}e^{-}\to\psi(3686)\pi^{+}\pi^{-}$, which include the measured total cross section, the $\psi(3686)\pi^{\pm}$, and $\pi^{+}\pi^{-}$ invariant mass spectra released in Ref. Ablikim:2017oaf . For the $e^{+}e^{-}\to\psi(3686)\pi^{+}\pi^{-}$ process, there exist contributions from intermediate vector charmonia as shown in Fig. 2 (2)-(4). Thus, it is a key point how to select intermediate vector charmonia. In the present work, we adopt two scenarios, which will be addressed in Sec. III.1 and Sec. III.2.

III.1 Three-charmonium scenario to $e^{+}e^{-}\to\psi(3686)\pi^{+}\pi^{-}$

When checking the vector charmonium states located in 4.1-4.6 GeV listed in particle data book Tanabashi:2018oca , there exist only two observed charmonia, $\psi(4160)$ and $\psi(4415)$. Besides, we must mention that there is another well-established charmonium-like state $Y(4220)$ in this energy range. This state has been observed experimentally in the processes of $e^{+}e^{-}\to\chi_{c0}\omega$ Ablikim:2014qwy , $e^{+}e^{-}\to\pi^{+}\pi^{-}h_{c}$ BESIII:2016adj , $e^{+}e^{-}\to\pi^{+}\pi^{-}J/\psi$ Ablikim:2016qzw , $e^{+}e^{-}\to\pi^{+}D^{0}D^{\ast-}$ Ablikim:2018vxx , and $e^{+}e^{-}\to\pi^{+}\pi^{-}\psi(3686)$ Ablikim:2017oaf . The mass of this state is consistent with the predictions of a missing higher charmonium in the $J/\psi$ family (see the discussion in Refs. He:2014xna ; Wang:2019 ; Chen:2014sra ). Thus, $Y(4220)$ is replaced by the renamed $\psi(4220)$ in the following discussion.

In the present section, we first try to fit the experimental data by considering $\psi(4160)$, $\psi(4220)$, and $\psi(4415)$ as intermediate state contributions to $e^{+}e^{-}\to\pi^{+}\pi^{-}\psi(3686)$. Their resonance parameters are collected in Table 1.

After integrating over $m_{\pi^{\pm}\psi(3686)}$ or $m_{\pi^{+}\pi^{-}}$ in Eq. (10) at each fixed $\sqrt{s}$, one gets the differential cross sections of $e^{+}e^{-}\to\pi^{+}\pi^{-}\psi(3686)$ for each invariant mass squared distribution. The cross section depending on $\sqrt{s}$ can be obtained by integrating over both $m_{\psi(3686)\pi^{\pm}}$ and $m_{\pi^{+}\pi^{-}}$ in Eq. (10). Here, we perform a combined fit of the invariant mass squared dependent cross sections to the experimental data of the differential cross sections at a fixed $\sqrt{s}$. In the three-resonance scenario, there are 38 parameters, which should be determined by the present fit. Our fitted values for all the parameters are presented in Table 2. With these parameters, $\chi^{2}/{\rm d.o.f}$ is estimated to be $3.704$.

The fitted differential cross sections for $e^{+}e^{-}\to\pi^{+}\pi^{-}\psi(3686)$ at a fixed $\sqrt{s}$ in three-resonance scenario are presented in Fig. 4 (black dashed curves), where $\sqrt{s}$ is fixed at $4.226$, $4.258$, $4.358$, $4.387$, $4.416$, and $4.600$ GeV, respectively. Most of the differential cross sections can be reproduced. In particular, both $\pi^{\pm}\psi(3686)$ and $\pi^{+}\pi^{-}$ invariant mass spectra are well reproduced for $\sqrt{s}=4.226$ GeV. As for $\sqrt{s}=4.258$ GeV, the $\pi^{\pm}\psi(3686)$ invariant mass spectrum is a bit different from the experimental data. In the vicinity of $m^{2}_{\pi^{\pm}\psi(3686)}=15\ {\rm GeV}^{2}$, the experimental data show a jump, but such a feature is not reproduced in three-charmonium scenario. On the other hand, the fitted curve of $m_{\pi^{+}\pi^{-}}$ distribution is slightly larger than the experimental data in the low and high $m_{\pi^{+}\pi^{-}}$ energy range. Similar to the case for the fit to the data at 4.358 GeV, the fitted curve of the $m^{2}_{\pi^{+}\psi(3686)}$ distribution is a bit lower on the whole than the experimental data, while for $m^{2}_{\pi^{+}\pi^{-}}$ distribution, the fitted curve is a little bit larger than the experimental data in the lower $m_{\pi^{+}\pi^{-}}$ energy range. For the cases of $\sqrt{s}=4.387$ GeV and $\sqrt{s}=4.416$ GeV, the fitted curves for both $m^{2}_{\pi^{\pm}\psi(3686)}$ and $m^{2}_{\pi^{+}\pi^{-}}$ well reproduce the experimental data. As for the data at $4.600$ GeV, the fitted curve obtained under three-charmonium scenario is lower than the experimental data for the $m^{2}_{\pi^{\pm}\psi(3686)}$ invariant mass distribution, while for $m^{2}_{\pi^{+}\pi^{-}}$ distribution, since the errors of experimental data are very large, thus, the fitted curve is roughly consistent with the data.

The fitted cross section for $e^{+}e^{-}\to\pi^{+}\pi^{-}\psi(3686)$ is presented in Fig. 5 (black dashed curve). The structure in the vicinity of 4.2 GeV and the broad enhancement near 4.4 GeV are well reproduced. However, one can find that the fitted curve is above the experimental data in the vicinity of 4.3 GeV. Such a phenomenon indicates that there may exist other intermediate charmonium contributions to the $e^{+}e^{-}\to\pi^{+}\pi^{-}\psi(3686)$ process, which inspires us to introduce a four-charmonium scenario to study the $e^{+}e^{-}\to\pi^{+}\pi^{-}\psi(3686)$ process, which will be discussed in the next subsection. In Fig. 6, we present the individual contributions from the nonresonance background and three charmonia. each of which is the coherent sum of direct contributions, the scalar productions, and ISPE processes for each resonance. One can find the contribution from $\psi(4160)$ could reach up to nearly 100 pb, while the contributions from $\psi(4220)$ and $\psi(4415)$ are about 70 pb. From Fig. 6, one can find that resonances are still in a Breit-Wigner form. Thus, we simulate each lineshape in Fig. 6 with a Breit-Wigner resonance, which is in the form $12\pi\Gamma^{e^{+}e^{-}}\mathcal{B}(R\to\pi^{+}\pi^{-}\psi(3686))\Gamma_{R}/[(s-m_{R}^{2})^{2}+m_{R}^{2}\Gamma_{R}^{2}]$. With the dilepton decay widths listed in Table 1, one can roughly estimate the center values of the branching ratios for $\psi_{R}\to\pi^{+}\pi^{-}\psi(3686)$, where $\psi_{R}=(\psi(4160),\ \psi(4220),\ \psi(4415))$. As shown in Table 3, one can find the center values of branching ratios are of the order of $\sim 10^{-3}$ for $\psi(4220)$ and $\psi(4415)$ and $\sim 10^{-2}$ for $\psi(4160)$.

III.2 Four-charmonium scenario to $e^{+}e^{-}\to\psi(3686)\pi^{+}\pi^{-}$

As discussed in the three-charmonium scenario in Fig. 5, the experimental data in the vicinity of 4.3 GeV cannot be well reproduced. In Ref. Wang:2019 , the mass spectra and decay properties of higher charmonia have been investigated, and a charmonium $\psi(4380)$ as a partner of $\psi(4220)$ was predicted, which has the resonance parameters $m=4384$ MeV and $\Gamma=84$ MeV. This progress stimulates us to introduce the four-charmonium scenario to study $e^{+}e^{-}\to\psi(3686)\pi^{+}\pi^{-}$, namely, four charmonia $\psi(4160)$, $\psi(4415)$, $\psi(4220)$, and $\psi(4380)$ are considered when fitting the data. All the resonance parameters are listed in Table 1. Indeed, we find that the details of the cross section of $e^{+}e^{-}\to\psi(3686)\pi^{+}\pi^{-}$ around $\sqrt{s}=4.3$ GeV can be depicted.

The fitted parameters are presented in Table 4. With the center values of these parameters, one can get the differential cross sections for each invariant mass squared distribution and the cross section depending on $\sqrt{s}$, which are presented in Fig. 4 (red curves) and Fig. 5 (red curve), respectively. After including $\psi(4380)$, one finds that almost all the differential cross sections can be quantitatively fitted better than the three resonance-scenario. Especially, the jump near at 15 GeV2 in $m^{2}_{\pi^{\pm}\psi(3686)}$ invariant mass spectrum at $\sqrt{s}=4.258$ GeV can be well reproduced and $m^{2}_{\pi^{+}\pi^{-}}$ invariant mass spectrum at this energy point is well fitted. The fits to differential cross sections at $\sqrt{s}=4.358$ GeV are also much better than the one in the three-resonance scenario. As for the $\sqrt{s}$ dependent cross section, the four-resonance scenario can well reproduce experimental data, especially the data near at 4.3 GeV. Quantitatively, the $\chi^{2}/{\rm d.o.f}$ is reduced to 2.77 in the four-charmonium scenario. Similar to the case of the three-charmonium scenario, we present the individual contributions of the non-resonance background and resonances in Fig. 7. One can find that the resonance contributions are smaller than the corresponding one in the three-charmonium scenario. The fitted branching ratios of $\psi_{R}\to\pi^{+}\pi^{-}\psi(3686)$ are presented in Table 3, where the central values of these branching ratios for $\psi(4160)$, $\psi(4220)$, and $\psi(4415)$ are of the order of $\sim 10^{-3}$, while the one for $\psi(4380)$ is of the order of $\sim 10^{-2}$.

IV Summary

In the past years, the BESIII and Belle Collaborations have made a great progress on finding charged charmonium-ike $XYZ$ states. A series of charged charmonium-like structures, $Z_{c}(3900)$ Ablikim:2013mio , $Z_{c}(3885)$ Ablikim:2013xfr , $Z_{c}(4020)$ Ablikim:2013wzq , and $Z_{c}(4025)$ Ablikim:2013emm , were reported one by one. These experimental results have stimulated a great interest of both theorists and experimentalists since the charged property of these charmonium-like states indicates these states should contain at least four constituent quarks, thus they could be good candidates of tetraquark states. The estimates in frameworks of the QCD sum rule Qiao:2013dda ; Wang:2013vex ; Wang:2013exa ; Wang:2013llv ; Agaev:2017tzv ; Zhao:2014qva and potential model Zhu:2016arf ; Deng:2014gqa ; Anwar:2018sol ; Patel:2014vua ; Liu:2017mrh indicate that both $Z_{c}(3900)/Z_{c}(3885)$ and $Z_{c}(4020)/Z_{c}(4025)$ could be explained as tetraquark states with different configurations. Furthermore, the observed masses of $Z_{c}(3900)/Z_{c}(3885)$ and $Z_{c}(4020)/Z_{c}(4025)$ are very close to the thresholds of $D^{\ast}\bar{D}$ and $D^{\ast}\bar{D}^{\ast}$, thus, these charmonium-like states have been extensively investigated in deutron-like hadronic molecular scenarios. Based on the molecular assumption, mass spectrum, decay properties, and productions of these charmonium-like states have been investigated in different models, such as QCD sum rule Cui:2013yva ; Zhang:2013aoa ; Cui:2013vfa ; Wang:2013daa , potential model Liu:2008tn ; Liu:2008fh ; Sun:2012zzd ; He:2013nwa , and some other phenomenological Lagrangian approaches Gutsche:2014zda ; Chen:2015igx ; Li:2013xia ; Li:2014pfa ; Xiao:2018kfx ; Chen:2016byt .

However, before definitely identifying these charged charmomium-like structures as exotic tetraquark states, we need to exhaust all the possibilities of explaining them in the conventional theoretical framework. Here, we must mention our previous prediction of charged charmonium-like structures in our theoretical work Chen:2011xk , where these predicted charged charmonium-like structures by the ISPE mechanism may exist in the $J/\psi\pi^{\pm}$, $h_{c}\pi^{\pm}$, and $\psi(3686)\pi^{\pm}$ invariant mass spectra of higher charmonia, and in the charmonium-like state $Y(4260)$ decays into $J/\psi\pi^{+}\pi^{-}$, $h_{c}\pi^{+}\pi^{-}$ and $\psi(3686)\pi^{+}\pi^{-}$ Chen:2011xk . After two years of our paper, the BESIII and Belle collaborations discovered a charged charmonium-like $Z_{c}(3900)$ in the process $e^{+}e^{-}\to J/\psi\pi^{+}\pi^{-}$ around $\sqrt{s}=4.26\ \mathrm{GeV}$ Ablikim:2013mio ; Liu:2013dau , and CLEO-c confirmed it in the same process but at $\sqrt{s}=4.17$ GeV Xiao:2013iha . In 2013, BESIII observed another charged charmoniumlike structure named as $Z_{c}(4020)$ in the $h_{c}\pi^{+}$ invariant mass spectrum of the $e^{+}e^{-}\to\pi^{+}\pi^{-}h_{c}$ process Ablikim:2013wzq . These experimental observations make us believe that the ISPE mechanism indeed plays an important role in producing these novel phenomena. What is more important is that these measured experimental data make us possible to study them with higher precision. In Ref. Chen:2013coa , we fitted the experimental result of the charged $Z_{c}(3900)$ observed in the $J/\psi\pi^{\pm}$ invariant mass spectrum Ablikim:2013mio ; Liu:2013dau and indicated that the charged $Z_{c}(3900)$ structure can be well established based on the ISPE mechanism. This study further enforces theorist’s confidence to explain these experimental results in the conventional theoretical framework. Similar studies were proposed in Refs. Swanson:2014tra ; Szczepaniak:2015eza ; Swanson:2015bsa ; Bugg:2011jr ; Liu:2015taa . Later, the Lattice QCD calculation Ikeda:2017mee ; Ikeda:2016zwx also supports such a scenario.

In 2017, the BESIII Collaboration again observed charged charmonium-like structures in the $\psi(3686)\pi^{\pm}$ invariant mass spectrum of the $e^{+}e^{-}\to\psi(3686)\pi^{+}\pi^{-}$ process at different energy points $\sqrt{s}=4.226,\,4.258,\,4.358,\,4.387,\,4.416,\,4.600$ GeV Ablikim:2017oaf . This new observation inspires our interest in further testing the ISPE mechanism combined with the present precise experimental data. A crucial task of this work is to examine whether the cross section of the $e^{+}e^{-}\to\psi(3686)\pi^{+}\pi^{-}$ process, and the corresponding $\psi(3686)\pi^{\pm}$ and $\pi^{+}\pi^{-}$ invariant mass spectra at different $\sqrt{s}$ values can be reproduced by the ISPE mechanism. As illustrated by our calculation, we find that the reported charged charmonium-like structures in the $\psi(3686)\pi^{\pm}$ invariant mass spectrum of the $e^{+}e^{-}\to\psi(3686)\pi^{+}\pi^{-}$ process at different energy points can be depicted in a unified framework, which sheds light on the properties of the observed charged charmonium-like structures. We need to emphasize that our study also further supports the existence of a new higher charmonium $\psi(4380)$ which was predicted in Ref. Wang:2019 , since the fitting result under the four-charmonium scenario is better than that under the three-charmonium scenario. We strongly suggest the BESIII and BelleII collaborations to search for this predicted missing charmonium $\psi(4380)$, especially by analyzing the $e^{+}e^{-}\to\psi(3686)\pi^{+}\pi^{-}$ process in near future.

In the following years, studies on charmonium-like $XYZ$ states will still be an interesting research topic at the typical facilities like BESIII, LHCb, and BelleII. Since the observed charmonium-like $XYZ$ states can be a good candidate of an exotic tetraquark matter, it is a main task for both experimentalists and theorists how to identify them as a genuine exotic multiquark state. To achieve this aim, we need to check whether the charmonium-like $XYZ$ states can be assigned into a conventional charmonium family or can be settled down in a conventional theoretical framework. The present work is the effort along this line. Besides theoretical investigation, the Lattice QCD study about charmonium-like $XYZ$ will provide valuable information to shed light on the underlying properties of the $XYZ$ states. As emphasized in a recent review article Liu:2019zoy , a joint effort from phenomenological method, experimental analysis, and Lattice QCD calculation should be paid more attentions on, which must promote our knowledge of how to form charmonium-like $XYZ$ states.

Acknowledgements

This work is partly supported by the China National Funds for Distinguished Young Scientists under Grant No. 11825503, the National Natural Science Foundation of China under Grant No. 11775050, National Program for Support of Top-notch Young Professionals, and the Fundamental Re- search Funds for the Central Universities.

Appendix A Amplitudes of the ISPE mechanism

In this Appendix, we will present some details of the ISPE mechanism. Under the ISPE mechanism, $\psi_{R}$ transits into $D^{(\ast)}$ and $\bar{D}^{(\ast)}$ pair associated with a single pion emission. Since the emitted pion has continuous energy distribution, $D^{(\ast)}$ and $\bar{D}^{(\ast)}$ mesons with the low momentum can easily interact with each other and further transit into $\psi(3686)\pi^{\pm}$ by exchanging a $D^{(\ast)}$ meson. In Fig. 8, we present some typical diagrams of the ISPE mechanism with a $\pi^{-}$ emission at the hadronic level. For convenience, in these diagrams, we denote $\psi(3686)$ as $\psi^{\prime}$. Although there are only 6 diagrams in Fig. 8, to get the whole contribution of the ISPE mechanism, we must consider the diagrams caused by isospin symmetry and initial $\pi^{+}$ emission. For isospin symmetry, the diagrams can be obtained by the replacement $\{D^{(\ast)+},D^{(\ast)0},\bar{D}^{(\ast)0}\}\to\{\bar{D}^{(\ast)0},D^{(\ast)-},D^{(\ast)+}\}$, which are the same contributions as for a $\pi^{-}$ emission. For contributions of an initial $\pi^{+}$ emission, the typical diagrams of the ISPE mechanism can be obtained by applying charge conjugation on each Fig. 8 (i) (i=a$\sim$f). Thus, the total decay amplitude can be expressed as

[TABLE]

where the factor 2 reflects the isospin symmetry mentioned above, and $\mathcal{A}_{\mu\zeta[{p_{2}\rightleftharpoons p_{3}}]}^{(i)}$ means that the contributions of an initial $\pi^{+}$ emission can be obtained by applying $p_{2}\rightleftharpoons p_{3}$ transformation on $\mathcal{A}_{\mu\zeta}^{(i)}$.

To write the decay amplitudes $\mathcal{A}_{\mu\zeta}^{(i)}$ in Eq. (12), the effective Lagrangian approach is adopted. For the four particles interaction, the effective Lagrangian is Chen:2011pv ; Oh:2000qr

[TABLE]

We want to note here that, to reduce the number of free parameters in our fit, we assume $h_{RD^{\ast}D^{\ast}\pi}=g_{RD^{\ast}D^{\ast}\pi}$ in this work.

The effective Lagrangian that describes the interaction among a pair of charmed mesons and a pion is Chen:2011pv ; Casalbuoni:1996pg ; Colangelo:2003sa

[TABLE]

Finally, the interaction between $D^{(\ast)}$, $\bar{D}^{(\ast)}$, and $\psi$ can be described as Chen:2011pv ; Casalbuoni:1996pg ; Colangelo:2003sa

[TABLE]

For the coupling constants $g_{D^{(\ast)}D^{(\ast)}\pi}$ in Eq. (A) and $g_{\psi D^{(\ast)}D^{(\ast)}}$ in Eq. (A), heavy quark symmetry gives the relations between them separately, which can be written as

[TABLE]

In Eq. (17), ${g_{\pi}=0.59}$, ${f_{\pi}=0.132}$ GeV, and the decay constant of $\psi(3686)$ is $f_{\psi}=0.416$.

With above effective Lagrangians, the expressions for decay amplitudes $\mathcal{A}_{\mu\zeta}^{(i)}$ can be written as

[TABLE]

[TABLE]

[TABLE]

In $\mathcal{A}_{\mu\zeta}^{(i)}$, a monopole form factor $\mathcal{F}(q^{2})=\frac{q^{2}-\Lambda^{2}}{m_{E}^{2}-\Lambda^{2}}$ is introduced to reflect the structure effect at each vertex and the off-shell effect of the exchanged $D^{(\ast)}$ meson. In this form factor, $m_{E}$ is the mass of the exchanged $D^{(\ast)}$ meson, $\Lambda$ is parameterized as $\Lambda=m_{E}+\alpha_{\Lambda}\Lambda_{QCD}$ with $\Lambda_{QCD}=220$ MeV. Since Ref. Chen:2011xk has shown that the ISPE mechanism is weakly dependent on $\alpha_{\Lambda}$, thus in this work we also set $\alpha_{\Lambda}=1$ to present our results.