Zc(3900)/Zc(3885)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z_c(3900)/Z_c(3885)$$\end{document} as a virtual state from πJ/ψ-D¯∗D\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pi J/\psi -\bar{D}^*D$$\end{document} interaction

In this work, we study the πJ/ψ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\pi \!}J{\!/}\psi $$\end{document} and D¯∗D\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\bar{D}^*D$$\end{document} invariant mass spectra of the Y(4260) decay to find out the origin of the Zc(3900)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z_c(3900)$$\end{document} and Zc(3885)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z_c(3885)$$\end{document} structures. The πJ/ψ-D¯∗D\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\pi \!}J{\!/}\psi -\bar{D}^{*}D$$\end{document} interaction is studied in a coupled-channel quasipotential Bethe–Saltpeter equation approach, and embedded to the Y(4260) decay process to reproduce both π-J/ψ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\pi ^{-}\!}J{\!/}\psi $$\end{document} and D∗-D0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D^{*-}D^0$$\end{document} invariant mass spectra observed at BESIII simultaneously. It is found out that a virtual state at energy about 3870 MeV is produced from the interaction when both invariant mass spectra are comparable with the experiment. The results support that both Zc(3900)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z_c(3900)$$\end{document} and Zc(3885)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z_c(3885)$$\end{document} have the same origin, that is, a virtual state from πJ/ψ-D¯∗D\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\pi \!}J{\!/}\psi -\bar{D}^*{D}$$\end{document} interaction, in which the D¯∗D\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\bar{D}^*{D}$$\end{document} interaction is more important and the coupling between D¯∗D\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\bar{D}^*{D}$$\end{document} and πJ/ψ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\pi \!}J{\!/}\psi $$\end{document} channels plays a minor role.


Introduction
In recent years, many exotic resonance structures were observed near the threshold of two hadrons, which are difficult to put into the conventional quark model. The exotic resonance structures near theD * D threshold (in this work we will remark hidden charmed system with a vector D * /D * meson and a pseudoscalarD/D meson asD * D if the explicit is not necessary) are good examples of such phenomena. The first XY Z particle, X (3872) is almost on theD * D threshold, which was interpreted as aD * D hadronic molecular state immediately after its observation [1,2]. Later, an isovector resonant structure named Z ± c (3900) was observed in 2013 at BESIII and Bell in the π ± J/ψ invariant mass spectrum of e + e − → π + π − J/ψ at √ s = 4.26 GeV [3,4], and further confirmed at CLEO-c in the same channel at √ s = 4.17 GeV [5]. The observed Z c (3900) is also near theD * D threshold, thus, it is natural to explain it as an isovector partner of X (3872) in theD * D molecular scenario and expected to be observed in the (D * D) I =1 channel. It was confirmed by the a e-mail: junhe@njnu.edu.cn observation of Z c (3885) in the DD * invariant mass spectrum of Y (4260) decay in process e + e − → π ± (DD * ) ∓ [6]. Recently, the neutral partners of Z ± c (3900) and Z ± c (3885) were also observed at BESIII [7,8]. Besides, the spin parity of these state has been determined as J P = 1 + by a partial wave analysis [9].
Thanks to the experiments at BESIII, Belle and CLEOc, the isovector Z c (3900)/Z c (3885) has been established.
Since the Z ± c (3900)/Z ± c (3885) carries a charge, it cannot be explained as a cc state, which must be neutral. After the observation at BESIII, many interpretations of the origin of Z c (3900)/Z c (3885) have been proposed, which includes the hadronic molecular state [10][11][12][13][14], tetraquark state [15][16][17], initial-single-pion-emission mechanism [18,19], cusp effect from triangle singularity [20]. And its decays were also studied though intermediate meson loop [21][22][23]. Due to its closeness to theD * D threshold, the hadronic molecular state is an important picture to explain the Z c (3900)/Z c (3885) structure. In Refs. [24,25], theD * D interaction as well as theB * B interaction was studied in a one-boson-exchange model. No bound state was found with the light-meson exchange. In the chiral unitary approach, the contribution of exchange of heavy meson was included into theD * D interaction, which provides attraction strong enough to produce a bound state [26]. The importance of the heavy-meson exchange was confirmed by a further study in the one-boson-exchange model combined with a quasipotential Bethe-Salpeter equation approach [27].
The early studies in the hadronic molecular picture focus on how to produce a bound state corresponding to the Z c (3900)/Z c (3885) from theD * D interaction. It is interesting to study if the bound state obtained in those studies can reproduce the original experimental data of the invariant mass spectra. In Ref. [26], the invariant mass spectra were studied in a chiral unitary approach while the Breit-Wigner form with mass and width obtained from the interaction were adopted. In Ref. [19] the invariant mass spectra were studied in an initial-single-pion-emission mechanism, where the D * D rescattering were not considered. Besides, when this study was done, only π J/ψ invariant mass spectrum was available. In Refs. [28][29][30][31], the mass invariant mass spectra was explicitly studied and fitted, and the poles corresponding to the Z c (3900)/Z c (3885) were extracted, especially in Ref. [28] the analysis suggested the Z c (3885)/Z c (3900) maybe originate from a virtual state. However, in theses studies, all coupling constants of the interactions were chosen as free parameters. A recent lattice work suggests that the off-diagonal π J/ψ −D D * coupling is more important than theD D * interaction, and a semiphenomenological analysis was adopted to study the invariant mass spectrum [32]. In their comparison with data, two general free parameters were adopted to π J/ψ andD D * channel, respectively, which smeared an important experiment results about the relative magnitudes of decays in these two channels. BESIII reported that the decay width of Z c (3885) in DD * is still much larger than that of Z c (3900) in J/ψπ channel with a ratio 6.2±1.1±2.7, though which is much smaller than conventional charmonium states above the open charm threshold [33].
In this work, we try to reproduce both line shapes and relative magnitudes of the π J/ψ andD * D invariant mass spectra, simultaneously. It is performed by studying the Y (4260) decay with reacattering of π J/ψ −D * D, which is calculated in a quasipotential Bethe-Salpeter equation approach. In the calculation, the interaction is constructed with the Lagrangians from the heavy-quark effective theory. And in this work, we only consider the system with negative charge, the positive and neutral cases are analogous due to the SU(3) symmetry.
In the next section, the formalism adopted to calculate the three-body decay of the Y (4260) in the current work is presented. The interaction potential is constructed with an effective Lagrangian and the quasipotential Bethe-Salpeter equation will be introduced briefly. The numerical results are given in Sect. 4. A brief summary is given in the last section.

Formalism of three-body decay
The Z c (3900)/Z c (3885) resonance structure were observed at the e + e − → Y (4260) → π Z c → π(π J/ψ/D * D) process at BESIII. The internal structure of the Y (4260) is still in the debate. To avoid the complexity, we adopt a phenomenological vertex Y (4260) → π(D * D). The effect of the e − e − → Y (4260) is also absorbed into this vertex. Hence, to study the invariant mass spectrum, we consider the three-body-decay diagram in Fig. 1. In this work, we focus on the invariant mass spectrum near theD D * threshold and the πJ/ψ interaction is suppressed by the OZI rule, The explicit calculation in our model also suggests that the intermediate π J/ψ channel will be suppressed and D * − D 0 channels in the loop between direct decay vertex and the rescattering from beginning, which was also adopted in Ref. [31]. When calculating the rescattering amplitude T , With the decay amplitude M the invariant of the mass spectrum can be obtained from the differential decay width of Y (4260) as where the M is the mass of the Y (4260) and phase space can be written as Here cm means the center of mass frame of particles 2 and 3. The explicit deduction is given in Appendix A.
The key to study the decay amplitude is to write rescattering amplitude T of the J/ψ − − D − D * 0 − D * − D 0 interaction, which can be obtained with the help of the Bethe-Salpeter equation as shown in Fig. 2.
To avoid difficulty of solving a four-dimensional equation in the Minkowski space, with quasipotential approximation, the Bethe-Salpeter equation is often reduced to a three-dimensional equation, which can be further reduced to a one-dimensional equation by partial wave decomposition. In this work, the OBE interaction will be adopted. The off-shellness of two constituent hadrons should be kept to avoid the unphysical singularity below the threshold. The covariant spectator theory, in which the heavier constituent is put on shell [34][35][36][37][38], will be adopted in our study of the π − J/ψ − D − D * 0 − D * − D 0 interaction. Such treatment was explained explicitly in the appendices of Ref. [27] and has been applied to studied the X(3250), the Z c (3900) and the LHCb pentaquarks and its strange partners [13,[39][40][41]. The partial-wave Bethe-Salpeter equation with fixed spin parity J P of system is written as [27] iT J P with the reduced propagator written down in the center-ofmass frame with P = (M, 0) as . (4) Here the heavier particle (remarked with h) is put on shell, The p 0 l for the lighter particle (remarked as l) is then W − E h (p ). Here and hereafter we will adopt a definition p = | p|. And the momentum of particle 2 p 2 = − p and the momentum of where λ 32 = λ 3 − λ 2 and η = P P 2 P 3 (−1) J −J 2 −J 3 with J (2,3) and P (2,3) being the spin and parity of constituent 2 or 3. Without loss of generality the initial and final relative momenta are chosen as p = (0, 0, p) and p = (p sin θ, 0, p cos θ), and the d J λλ (θ ) is the Wigner d-matrix. In most cases, the integral in Eq. (3) is non-convergent. In this work an exponential regularization is introduced by a replacement of the propagator as We would like to remind that the regularization of heavier particle vanishes because it is put onshell in the quasipotential approximation adopted. With the regularization, the contributions at large momentum p is suppressed heavily at the energies higher than 2 GeV [42], which guarantees the convergence of the integral. if we multiply exponential factor on both sides of the Eq. (3), it can be found that the regularization factor can be seen as a form factor introduced due to the off-shell effect of particle 1 in a form of e −(k 2 −m 2 ) 2 /Λ 4 . The interested reader is referred to Ref. [27] for further information about the regularization.
The write the amplitude of three-body decay of the Y (4260), we adopt an effective Lagrangian for the Y → π D D * as, The partial-wave amplitudes with spin parity J P is With Lagrangian we adopted, only the J P = 0 + and 1 − partial wave survive as The total three-body decay with the rescattering is written as The distribution can be obtained as Here the explicit form of A J λ 2 ,λ 3 ;λ (p cm 3 , Ω 1 ) with J = 0, 1 in Eq. (9) is applied.
The distribution can be further rewritten with the partial wave amplitudes withe J P as witĥ where i and j denote the independent λ 2,3 and λ, and the factors f i=0 = 1/ √ 2 and f i =0 = 1 are inserted. In this work we introduce parameterized background contribution witĥ The parameters will be determined by comparing with experiment.
The couplings of heavy-light charmed mesons to J/ψ follow form, The three couplings in (18) are related to the single parameter g 2 as With above Lagrangians, the potential for the interactions can be constructed, which is presented explicitly in Appendix B.

The numerical results
The amplitude T for the π − J/ψ − D − D * 0 − D * − D 0 interaction can be obtained by discretizing the momenta p, p , and p in the integral equation (3) by the Gauss quadrature with a weight w(p i ). After such treatment, the integral equation can be transformed to a matrix equation [27] The propagator G is a diagonal matrix as with on-shell momentum The rescattering amplitude T can be solved as T = (1 − V G) −1 V . The pole of rescattering amplitude can be found at |1 − V G| = 0 after analytic continuation total energy W into the complex plane as z. The amplitude for the Y (4260) decay M can be written as M = A bk +T G A with the on-shell element being chosen.
In our model, the parameters in the Lagrangians are determined by the heavy quark symmetry. The free parameters are the cutoff Λ and the a, b and c for the background. The cutoffs in the regularization and in the form factor for the exchanged meson have the same value for simplification. In this work, we try to reproduce the line shapes and relative magnitudes of the π − J/ψ and D * − D 0 invariant mass spectra by varying the parameters. When comparing the theoretical results and the experimental data, we should be careful about the number of the events of two channels which were obtained with different efficiencies in experiments. Fortunately, in original report of BESIII [33], both the cross sections and corresponding numbers of the events for the Z c (3900) in π − J/ψ channel and Z c (3885) in D * − D channel were presented as 13.5 pb with 307 events and 83.5 pb with 502 events, respectively. Here we adopt the M max (π ± J/ψ) distribution as the π − J/ψ invariant mass spectrum to avoid the reflection peak because the background contribution is parametrized in this work. The theoretical results for events can be obtained by multiplied the efficiencies on theoretical decay distribution for π − J/ψ and D * − D 0 , respectively. Besides, the different bin sizes adopted in two channels are also considered in the calculation. After such treatment and a general normalization to the experimental data, the comparison between the theoretical and experimental results can be carried out. It is found that with a cutoff Λ = 1.85 GeV the invariant mass spectra can be reproduced as shown in Fig. 3, and the corresponding parameters for the background are (a, b, c) = (0.5, 1.2, 3.6) and (1.0, 0.05, 8.5) for the π − J/ψ and D * − D 0 invariant mass spectra, respectively.
At low energies, the π − J/ψ invariant mass spectrum is mainly from the background contribution, which decreases with increase of the energies near and higher than the D * − D 0 threshold. A sharp peak arises near the threshold due to the π J/ψ −D * D rescattering, which effect decreases a little slower at energies above the threshold than at energies below the threshold. The full model can reproduce the π − J/ψ invariant mass spectrum generally. The peak seems too sharp compared with the experiment, which may be from the contributions neglected in this work. As in our previous work in Ref. [27] where only D * D scattering were considered, the D * − D 0 invariant mass spectrum of the Y (4260) decay Though peaks can be produced in the invariant mass spectra, we still need to find out that the resonance structures are from a pole or just cusps. In the literatures [47,48], the category of the pole from the two-body interaction has been studied. In the order of the attraction of interaction from strong to weak, there exist four types of poles, bound state which is below threshold and usually called molecular state, virtual state which is also below the threshold but in the second Riemann surface, virtual state with width which is below the threshold but has an imaginary part, and resonance which is beyond the threshold and has an imaginary part. Hence, we adjust the cutoff, with which the strength of the interaction has positive correlation. The poles produced form the π − J/ψ − D * − D 0 − D − D * 0 interaction with typical cutoffs are listed in Table 1.
First, we list both the results for the π − J/ψ − D * − D 0 − D − D * 0 interaction and these after turning off the π − J/ψ channel. From the results, one can find that the interaction is dominant with the D * D interaction. If the cutoff larger than about 2.1 GeV, a pole below theD * D threshold is produced from the interaction. Compared with the results with- Table 1 The bound states from the π − J/ψ − D * − D 0 − D − D * 0 interaction with typical cutoffs Λ. The cutoff Λ and energy W are in units of GeV, and MeV, respectively  out π − J/ψ channel, the imaginary part of the pole is obviously from the coupled-channel effect. By varying the values if cutoff a little, the results in full model and these without π − J/ψ channel are almost same. Hence, it is a bound state mainly from theD * D interaction. With the decrease of the cutoff, the interaction becomes weaker and the pole is running to the threshold. When the cutoff is smaller than about 2.1 GeV, a pole will appear in the second Riemann surface of theD * D interaction. This pole is leaving the threshold with the decrease of the cutoff and will merge with the lower pole as shown in Fig. 3. If the cutoff decreases further, the pole dies away and no virtual state with width is produced, though a peak still can be produced as a cusp near theD * D threshold in the invariant mass spectra which is much wider than the experimental Z c (3885). Combined with results in Fig. 3 and in Table 1, one can find both Z c (3900) and Z c (3885) are from a virtual bound state mainly from theD * D interaction. The invariant mass spectra can not be reproduced with a bound state or cusp effect without pole.
The above results suggest that the π − J/ψ channel plays a minor role in the π − J/ψ − D * − D 0 − D − D * 0 interaction, which also leads to a relatively small decay width in π J/ψ channel compared with that inD * D channel reported at BESIII [33]. It is interesting to give the results only with coupling between π − J/ψ channel and D * − D 0 − D − D * 0 channel. From the results in Fig. 4, if we increase the cutoff to a value about 3 GeV, the peak in the π − J/ψ invariant mass spectrum can be reproduced. However, the peak in the D * − D 0 invariant mass spectrum is much wider than the experiment. No pole is produced from the interaction, and the peaks are from the cusp effect. Hence, in our model, with only the coupling of π − J/ψ and D * − D 0 − D − D * 0 channels, the D * − D 0 invariant mass spectrum can not be explained. This result supports two structures are from the virtual state mainly from theD * D interaction.

Summary
In this work, the π J/ψ andD * D invariant mass spectra of the Y (4260) decay is studied with rescattering of π J/ψ − D * D, which is calculated in a quasipotential Bethe-Salpeter equation approach. The theoretical invariant mass spectra are compared with BESIII experiment to determine the pole structure of π J/ψ −D * D interaction.
The peaks in both invariant mass spectra are reproduced from the π J/ψ −D * D rescattering in the Y (4260) decay. When the experimental data at BESIII is reproduced, the π J/ψ −D * D interaction produce a virtual state at energy of about 3870 MeV. TheD * D channel plays important role to produce the virtual state and the coupling between π J/ψ andD * D is relatively small, which is consistent with experimentally observed larger cross section of Z c (3855) inD * D channel than that of Z c (3900) in the π J/ψ channel. After turning off theD * D interaction and keeping only the coupling between π J/ψ andD * D, the cusp effect still can give peaks near theD * D threshold. However, the peak in theD * D invariant mass spectrum is quite broad, which conflicts with the BESIII experiment.
Here, we present the explicit form of the one-boson-exchange potential. The potentials for the D − D * 0 → D − D * 0 interaction with vector V, J/ψ and σ meson exchanges are (B.9) The potential for the D * − D 0 → D * − D 0 interaction can be obtained from these for the D − D * 0 → D − D * 0 interaction by alternating particle 1 and particle 2.
A form factor is introduced to compensate the off-shell effect of exchanged meson [50,51] f (q 2 ) = Λ 4 + q 2 t − m 2 2 /4 Λ 4 + q 2 − q 2 t + m 2 /2 2 , (B.12) where q 2 t denotes the value of q 2 at the kinematical threshold. The kinematical regime between the threshold and the on-shell point of the exchange particle is stressed and tchannel contributions at threshold are directly given by their couplings. The form factor is only function of the Lorentz invariant q 2 , pole free on the real q 2 axis, normalized to 1 for q 2 = m 2 and q 2 = q 2 t , but does not have its maximum at q 2 = m 2 . In the propagator of the meson exchange we make a replacement q 2 → −|q 2 | to remove the singularities as Ref. [52].