Exotic vector charmonium and its leptonic decay width

We propose a novel type of interpolating field operators, which manifests the hybrid-like configuration that the charm quark-antiquark pair recoils against gluonic degrees of freedom. A heavy vector charmonium-like state with a mass of $4.33(2)\,{\rm GeV}$ is disentangled from the conventional charmonium states in the quenched approximation. This state has affinity for the hybrid-like operators but couples less to the relevant quark bilinear operator. We also try to extract its leptonic decay constant and give a tentative upper limit that it is less than one tenth of that of $J/\psi$, which corresponds to a leptonic decay width about dozens of eV. The connection of this state with $X(4260)$ is also discussed.


I. INTRODUCTION
X(4260) was observed by many experiments as a ππJ/ψ resonance structure in the initial state radiation (ISR) process e + e − → γ ISR J/ψππ [1][2][3]. Its resonance parameters are determined now to be M X = 4251 (9) MeV and Γ X = 120 (12) MeV [4]. According to its production mode, X(4260) must have the quantum number J P C = 1 −− . In addition, the ratio of X(4260) → π + π − J/ψ and X(4260) → π 0 π 0 J/ψ events observed by the CLEO collaboration is consistent with X(4260) being an isoscalar. In other words, X(4260) has the same quantum number as that of vector charmonia J/ψ, ψ ′ , etc. However, in contrast to ψ states, X(4260) has not been observed directly in the e + e − annihilation yet. On the other hand, its mass is well above the DD threshold, but it has been observed only in the J/ψπ + π − system instead of DD ones. These facts may imply that X(4260) has a large branch fraction for the J/ψπ + π − decay mode. Thus, the small combined width Γ(X → e + e − )Br(X → J/ψππ) = 9.2±1.0 eV can be understood as that X(4260) has a very small e + e − width. These features motivate conjectures that X(4260) might be an exotic state, for example, a hybrid charmonium [5][6][7]. Anyway, more theoretical information for X(4260) is needed in order to unravel its nature, among which the leptonic decay width of X(4260) is an important quantity.
As far as the hybrid charmonium is concerned, extensive lattice QCD studies have been devoted to J P C = 1 −+ channel. In the constituent quark model picture, this quantum number cannot appear in the qq system, therefore it is usually conjectured that additional degrees of freedom should be involved and the minimum * cheny@ihep.ac.cn configuration can be qqg where g is a constituent gluon. The corresponding qqg interpolation fields are used in lattice calculations which predict the mass of the 1 −+ charmonum-like state to be around 4.3 GeV [8][9][10][11][12][13][14][15][16][17]. Similar studies have been also extended to investigate possible hybrids with the conventional quantum number, but the challenging task is to distinguish these states from conventional mesons. The state-of-art approach for this goal is the variational method based on large enough operator sets built through sophisticated methods. In the vector channel, there is a state observed with a mass round 4.4 GeV [14], which couples weakly to the quark bilinear operator but seems intimate with quark-antiquark-gluon operators.
Generally speaking, the appearance of an interpolating field operator does not necessarily reflect the inner structure of a hadron state. However, for heavy quark systems where the non-relativistic picture may be available to some extent, the coupling of the operator to a specific state may bear some useful information of its status. Taking a mesonic hybrid for example, even though it is an ambiguous concept from the point of quantum chromodynamics (QCD) , it is always thought of a hadron state made up of a quark-antiquark pair plus a constituent gluonic component in the constituent quark model picture. Of course one can also relax the definition of a hybrid to an exotic object which has additional degree of freedom apart from the constituent quarks. This kind of additional degree of freedom can be a fluctuating flux tube in the flux tube model, the color bag of the MIT bag model, etc. Anyway, the essence of the exotic nature of a meson state is that the constituent qq pair acquires a center-of-mass motion by recoiling against the additional degrees of freedom, which is distinct from the conventional hadron states. This is our starting point to build a novel type of hybrid-like operators. We split the charm quark-antiquark pair component and the gluon field in qqg operator into two parts with different spatial separa-tions. In the momentum space, this manifests the centerof-mass motion of qq pair in the rest frame of the state. We calculate the correlation functions of these operators, from which we try to extract the possible exotic charmonium state. Since the operators with different spatial separation provide different correlation functions, we fit them simultaneously along with the correlation function involving the electromagnetic current to obtain the decay constants of the states which contribute significantly. This paper is organized as follows: Section 2 contains the description of the construction of the new lattice interpolation operators for the hybrid-like vector meson. The lattice parameters and the numerical techniques are presented in Sec. 3. We discuss our results and their connections with X(4260) in Sec. 4. The conclusions and a summary can be found in Sec. 5.

II. NEW INTERPOLATION FIELD FOR EXOTIC VECTOR CHARMONIUM
In this work, we will focus on the exotic vector charmonium (J P C = 1 −− ) by assuming a hybrid-like configuration ccg. A simple and straightforward local operator possibly reflecting this constituent configuration is O , where a, b are color indices, i the spatial index, and B ab i (x) = 1 2 ǫ ijk F ab jk the chromomagnetic field tensor. This kind of operator can be compared with the commonly used quark bilinear operator for the vector O (M) i =cγ i c(x). In order to find the nonrelativistic form of these interpolation operators, we use the Foldy-Wouthuysen-Tani transformation [18] to decompose the charm quark and antiquark fields (Dirac spinor) in terms of the Pauli spinors φ/φ † which annihilates/creates a charm quark, and χ/χ † which creates/annihilates a charm antiquark. The explicit expressions of the operators O to the lowest order of the nonrelativistic approximation can be written as where one can see that the block operator is a spin singlet and color octet, while that of O . In order to resemble the center-of-mass motion of the cc recoiling against an additional degree of freedom, we split the operator O (H) i into two spatial parts,c a γ 5 c b and B ab i , separated by an explicit spatial displacement r. In a fixed gauge (the Coulomb gauge in this work), we get a set of It is expected that the coupling of this type of operator to the conventional charmonia (without the center-ofmass motion of charm quark-antiquark pair in the nonrelativistic picture) would be suppressed, while the coupling to the exotic state can be enhanced.

III. NUMERICAL DETAILS
We use the tadpole-improved gauge action [19][20][21] to generate gauge configurations on anisotropic lattices with the temporal lattice spacing much finer than the spatial one. The aspect ratio takes ξ = a s /a t = 5, where a s and a t are the spatial and temporal lattice spacing, respectively. Two lattices L 3 × T = 8 3 × 96(β = 2.4) and 12 3 × 144(β = 2.8) with different lattice spacings are used to check the discretization artifacts and the relevant input parameters are listed in Table I, where a s values are determined from r −1 0 = 410 (20) MeV. We use the tadpole-improved clover action to calculate the quark propagators. The relevant parameters in the fermion action are tuned carefully by requiring that the physical dispersion relations of vector and pseudoscalar mesons are correctly reproduced at each bare quark mass [22,23]. The bare charm quark masses at different β are determined from the physical mass of J/ψ, m J/ψ = 3.097 GeV. The spatial extension of both lattices is ∼ 1.7 fm, which is tested to be large enough for charmonium states. The ground state masses of 1S and 1P charmonia calculated on these two lattices (see Fig. 2 and Table II of Ref. [24] for the details) show that the finite a s effects are small. Since the spatial extended interpolation operators O (H) i discussed above are gauge variant, we carry out the calculation of the quark propagators and correlation functions after transforming each configuration to the Coulomb gauge.

A. Data analysis strategy
Our first task is to verify the existence of the exotic vector charmonium. We use the following source operator to calculate the correlation functions, where τ refers to the source time slice. For the sink operator O (H) i , the two-point functions we calculate are where S F (x, y) stands for the charm quark propagator. Accordingly, there are two types of wall-source quark propagators to be calculated. One of them uses the usual wall source by setting the source element to unity at each spatial site of the source time slice. The other one uses the source by multiplying the chromomagnetic field tensor B i (z, τ ) to each site of the plain wall source.
In order to increase the statistics additionally, for each configuration we calculate T charm quark propagators S F ( x, t; 0, τ ) by setting the corresponding source vectors on each time slice τ . This permits us to average over the temporal direction when calculating the two-point functions.
In practice, the two-point functions C (H) (r, t; τ ) with the same r = |r| are averaged, such that the quantum number is kept to be right J P C = 1 −− . After averaging over the time direction, the practical two-point functions we calculate are where N r is the degenerate degree of r = |r|. In the data analysis stage, we perform a simulitaneous multi-exponential fit to C (H) (r, t)'s by using a correlated minimal-χ 2 fit method with the jackknife covariance matrix (we use three mass terms throughout this work).

B. Masses and Bethe-Salpeter amplitudes of conventional vector charmonia
In order to test the reliability of the fitting strategy mentioned above, we first carry out a similar analysis to the correlation functions C (M) (r, t) of the spatially extended version of operator O (M) on the β = 2.4 lattice. The procedure is detailed as follows. The spatially ex- whose correlation function with the corresponding wall source operator , say, C ′ (M) (r, t + τ ; τ ), is defined similarly as in Eq. (4). After averaging over the temporal direction, we have where Φ ′ i (r) is the r-dependent spectral weight of the i-th state.
In the fitting procedure, we fix a maximal t (denoted by t max and varying the lower bound t min of the fit window, then we obtain the masses of the lowest three states keep constant to some extent for a series of t min , as shown in the upper panel of Figure 1. We average the masses in this range with each value weighted by its error and get the values m 1 = 3.097(1) GeV, m 2 = 3.679 (19) GeV, and m 3 = 4.007(57) GeV, respectively. These three states may correspond to J/ψ, ψ ′ (3686), and ψ(4040) and we can almost reproduce their experimental spectrum. We also plot the Φ ′ i (r/a s )'s (normalized as Φ ′ i (0) = 1) in the right panel of Figure 1 where one can find that, there is no radial node for Φ ′ 1 (r/a), one radial node for Φ ′ 2 (r/a s ), and there are two radial nodes for Φ ′ 3 (r/a s ). Given the quark model assignments n 3 S 1 state for J/ψ, ψ(3686), and ψ(4040) with n = 1, 2, 3, respectively, this is actually not surprising, since Φ ′ i (r) is proportional to the Coulomb Bethe-Salpeter amplitude of the i-th S-wave charmonium, which, at the leading order of the non-relativistic approximaiton, corresponds to the radial wave functions in the quark model [25,26].

C. Existence of an exotic vector charmonia and its mass
From above one can see that our data analysis strategy is robust for the conventional vector charmonium states, therefore we perform the similar study for the correlation functions described in Eq. (4) and (5). Figure 2 shows the plots of Φ i (r) with respect to r (in physical units) through a three-mass-term fit (the upper panel is for β = 2.4 at t min = 12a t , and the lower panel for β = 2.8), whose masses are fitted to be 3.100 (7) GeV, 3.58(9) GeV, and 4.6(2) GeV for β = 2.4, and 3.090(6) GeV, 3.54(5) GeV, and 4.6(1) GeV for β = 2.8 at t min = 16a t . Φ 1 (r) and Φ 2 (r) damp more rapidly and are close to zero near r ∼ 0.3 fm while Φ 3 (r) is still relatively large. The lowest two states correspond very possibly to the conventional vector charmonia J/ψ and ψ ′ according to their masses. In contrast, the third state, with a much higher mass, still dominates the two point functions with r larger than 0.3 fm. This may signal the exotic nature of this state that is reflected by the spatially extended sink operator O (H) i . Of course, the higher conventional vector charmonia, such as ψ(4040) and ψ(4415), should also contribute to the two-point function C (H) (r, t), however in our data analysis procedure, C (H) (r, t) cannot accommodate more statistically meaningful states. The presence of the higher conventional charmonia may result in the small shift of the masses of the fitted states, for example, the mass of the second state deviates from that of the would-be ψ ′ state.
It is seen from Fig. 2 the r-behaviors of the first (J/ψ) and the second state(ψ ′ ) are similar up to an overall factor. If this is the case for all the conventional charmonium states since the r behavior may depict the center-of-mass motion of the cc component and the conventional charmonia are free of this in the nonrelativistic approximation, we can conjecture that the spectral weights of the i-th conventional charmonia Φ i (r) can be factorized into Φ(r)W i where Φ(r) is approximately uniform and insensitive to the different conventional charmonia, such that the two-point functions C (H) (r, t) with different r can be linearly combined to eliminate the contribution from the conventional charmonium states. In practice, we combine linearly the correlation functions C (H) (r, t) at two specific r 1 and r 2 as where ω is a tunable parameter. For each lattice, an optimal ω can be obtained by the requirement that the effective mass plateau of C(ω, t) is as long as possible. To be specific, for β = 2.4, we use r 1 = 0 and r 2 = a s and set ω = 2.576, which is very close to the central value of the ratio Φ 1 (0)/Φ 1 (a s ) = 2.583 through the three-mass fit illustrated in Fig 2. Using this omega and the fit results in Fig 2, the spectral weights of the three states in C(ω, t) are roughly 0.002, -0.1, and -1.1, respectively. This implies that the relative contribution of the third state is strongly enhanced by this subtraction scheme. Similarly, for β = 2.8 we use r 1 = a s and r 2 = √ 3a s and set ω = 3.658 which is also close to the central value of the ratio Φ 1 (a s )/Φ 1 ( √ 3a s ) = 3.681. The spectral weights of the three states are roughly 0.007, -0.2 and -1.9, respectively. In Fig. 3 we plot the effective mass plateaus of C(ω, t)'s for β = 2.4 and β = 2.8. The t and the masses are expressed in the physical units according to the lattice spacings listed in Table I. One can see that both plateaus are fairly good and lie on each other. Since both the sink operator O (H) (r) and the source operator O (W ) (whose correlation functions are C (H) (r, t)) are expected to couple strongly to hybrid-like states and the contribution from the conventional charmonia is subtracted largely by the above scheme, we take the state reflected by the observed plateau as the exotic vector charmonium and name it as X in the rest part of this work (we keep the name of the experimental state X(4260)). The horizontal line shows the fitted mass M X = 4.33(2) GeV through a one-exponential fit in the time range from 0.3 to 0.9 fm.
To this end, we claim that a vector charmonium-like state with a mass of 4.33(2) GeV has been unambiguously singled out, whose exotic nature may be reflected by its distinct coupling to the special interpolation field O (H) (r, t) in comparison with the conventional vector charmonia. It should be noted that this state has also been observed by previous lattice studies using variational methods based on lattice operator sets [14], however, the spatially extended operators O (H) (r, t) we use give more clear picture of its inner structure.

D. Leptonic decay constant of the exotic vector charmonium
Since this hybrid-like charmonium can be disentangle from the conventional charmonia with the prescription above, its leptonic decay constant can be investigated accordingly. The leptonic decay constant f V of a vector meson state V is defined by where J (em) µ (0) is the electromagnetic current and ǫ µ ( p, r) is the polarization vector of V at momentum p. For vector charmonium states, J (em) µ (0) can be approximated bycγ µ c(0) if the contribution from other quark flavors through annihilation diagrams is neglected. Since the vector current J em µ (x) defined in the continuum limit is no longer conserved on the lattice, we perform a nonperturbative renormalization procedure [27] to extract the multiplicative renormalization constant Z V of the current. The renormalization constant of the spatial componets of J em µ is determined to be Z V . In order to obtain this matrix elements, we also calculate other two categories of correlation functions in addition to C (H) (r, t), where averaging over the temporal direction is also taken implicitly in the above expressions. After the intermediate states insertion to C (J) (t), C H) (r, t), and C (W ) (t), we have where m n is the mass of the n-th state and the parameter Z (K) n with K refering to H or W is defined as Accordingly the leptonic decay constant f Vn can be derived from Z (J) n from the definition Eq. (10) as where C is a overall constant prefactor owing to the redefinition of our quark fields and the anisotropic lattices we are using. The time dependence of the correlation functions is usually observed from their effective mass plots. The effective masses of different correlation functions C (X) (t) are defined as where K stands for J, H, and W and illustrated in Fig. 4, where the mass of the ground state (J/ψ) is also plotted as a horizontal dashed line to guide eyes. The left panel of Fig. 4 is for the β = 2.4 lattice, and the right panel is for β = 2.8. In the plots, we present in the third row the effective mass of C (H) (r, t) at a specific r (r = √ 6a s for β = 2.4 and r = 3a s , respectively), where C (H) (r, t) is dominated by the third state in the short time range, changes the sign and finally saturated by the ground state when t increases. This is manifested in the effective mass plot by the phenomenon that M eff (t) shows a meta-stable plateau (roughly 4.4 GeV) higher than the ground state before the discontinuity time, and then converges to the mass of the ground state. This phenomenon also implies that the spectral weight of the third state is much larger than that of the lower states and the signs of the two spectral weights are different. In the fitting procedure, we fix the maximal time t max of the fitting window (t max = 30a t for β = 2.4 and t max = 40a t for β = 2.8), and let the minimal time vary in a range. The fitted spectral weights of C (J) (t) and C (W ) (t) are listed in Table II (raw data), from which the decay constants can be derived. The fitted masses and the decay constants are converted to the values in the physical unit (GeV) and are presented in Table III. The spectral weights of C (H) (r, t) are less relevant and omitted here to save space (one can refer to Fig. 2 to see their relative magnitudes for the three states). All the errors are statistical and obtained through a jackknife analysis. In order to illustrate the fit quality, we also show the effective masses (in blue bands) of the functions in Eq. (12) using the fitted parameters at t min = 12a t for β = 2.4 and t min = 16a t for β = 2.8. It is seen that the fit functions describe the measured data very well. Actually for all the fits, the χ 2 /N d.o.f 's are around 1 even though the number of the degree-of-freedom N d.o.f are always several hundreds. Table II and Table III, the fitted parameters are almost stable and insensitive to t min . The masses of the first and the second states are consistent with the those of J/ψ and ψ ′ , while the mass of the third state is a little higher than the hybrid-like state we obtain before. This can be attributed to the contamination from even higher states to some extent. In the right part of Table II, the spectral weights |Z (W ) 3 | 2 are an order of magnitude larger than those of the lowest two states (corresponding mostly to J/ψ and ψ ′ ). This is not strange since Z (W ) n is the coupling of the hybrid-like wall-source operator O (W ) to the n-th state and is expected to be enhanced when The spectral weights of the three states at different tmin on both lattices through a correlated three-mass-term fit. The errors are obtained through a jackknife analysis. The experimental results are also listed for comparison.    Table I. The errors are obtained through a jackknife analysis. The experimental results are also listed for comparison.
3.097 0.407(5) 3.686 0.290 (2) coupling to a hybrid-like state. In contrast, the spectral weights |Z  Table II) of the lowest two states, even suppressed by Z (W ) n , are much larger than that of the third state which are close to zero with errors. This may imply that the decay constant of the third state is very small, since Z n , from which the decay constant of the n-th state can be derived using Eq. (14). The derived decay constants of the three states are also listed in Table III, where the last row lists the experimental values for comparison. For J/ψ, we get its decay constant to be roughly 0.46(2) GeV at β = 2.4 and 0.43(2)GeV at β = 2.8, which are close to the experimental value although 5%-10% larger. The deviation can be attributed to the artifact of the finite lattice spacings (and also the uncertainty owing to the quenched approximation). The derived decay constant of ψ ′ , f ψ ′ , seems compatible with the experimental value, but with huge errors which come mainly from the uncertainty of |Z The most striking observation is that the decay constant f 3 of the third state is consistent with zero within the error. Superficially, it seems that the exotic vector charmoium has a nearly zero decay constant. However, there is a possibility that this is a mixing effect of two nearby states (for example, the would-be exotic state and ψ(4415)), whose contribution to C (J) (t) cancels to some extent, because we perform the simultaneous fit using only three mass terms. We have addressed that an exotic vector charmoninum state does exist and contributes substantially to the correlation functions C (H) (r, t) in previous context. This is also the case for C (W ) (t) since |Z (W ) 3 | 2 is one magnitude or even more larger than |Z (W ) 1 | 2 and |Z (W ) 2 | 2 . Based on these facts and considering the possible admixture of a conventional charmonium state to the exotic state, we try to estimate the upper limit of the decay constant of the exotic vector chrmonium. If the third state is actually contributed from the would-be exotic state X and the adjacent vector charmonium state ψ(4415), then the spectral weight Z According to Eq. (14), this is equivalent to Since |Z where we take α QED = 1/134 at the charm quark mass scale. Therefore, f X can be roughly estimated to be So we can assign a safer upper limit of f X as which gives the upper limit of the leptonic decay width of the exotic vector charmonium, Γ(X → e + e − ) < 40 eV.

IV. DISCUSSION
The suprisingly small e + e − decay width of the exotic vector charmonium X is in sharply contrast to that of conventional vector charmonia, which are usually of the order of keV. In other words, if this hybrid-like vector charmonium does exist in the real world, its contribution to the inclusive cross section of e + e − annihilation is rather small. Actually the R value scan versus the invariant energy √ s of e + e − collision does not show any indication of an extra vector charmonium-like state around √ s = 4.3 GeV (there is however a small dip in this energy range). The BEPCII/BESIII in Beijing is now accumulating the data of e + e− collision in this energy range and will hopefully give more precise line shape of R-value here to clarify the situation. On the other hand, as mentioned above, the vector charmonium state X(4260) was observed by several experiments in the initial state radiation of e + e − annihilation into J/ψπ + π − . The combined decay width of X(4260) is Γ(X(4260) → e + e − )Br(X(4260) → J/ψππ) = 9.2±1.0eV.
If X(4260) is tentatively assigned to the X state investigated in this study, combining the above value with the leptonic decay width of X, we can give an estimate of the branch ratio of X(4260) decaying into J/ψπ + π − Br(X(4260) → J/ψππ) > 20%, (24) which means that J/ψππ is one of the most important decay mode of X(4260). This can naturally explain why X(4260) was only observed in this channel till now. Furthermore, given the likely hybrid nature of the X state, it can be expected that the spin singlet cc component of X prefers a hadronic transition into spin singlet charmonium, such as h c . So X(4260) → h c ππ can be also an important decay mode of X(4260). Recently the BESIII Collaboration studied the e + e − → π + π − h c process at the center-of-mass energies from 3.90 GeV to 4.42 GeV. They found that the cross sections are of the same order of magnitude as, but have different line shape from those of e + e− → π + π − J/ψ [29]. Anyway, it is highly desired to investigate whether they are from the same resonance structure or not. The reason for the large branching ratio of X(4260) → J/ψππ can be depicted as follows. In the e + e − annihilation, the charm quark-antiquark pair cc is produced in the short range through the virtual photon. During the hadronization procedure, the charm quark and the charm antiquark emit soft gluons continuously, which form a colored gluon halo around the gradually localized color octet cc (in a relative sense). Finally a meta-stable state is formed as the X(4260) particle. Obviously, the color octet cc kernel can be readily neutralized into a color singlet charmonium by absorbing (emitting) soft gluons from (to) the halo, and the gluon halo thereby becomes color neutral and are emitted as light hadrons, for example, the ππ pair. If this is actually the case, the color flux between the charm quark and antiquark have less chance to be excited to a high enough energy to break, and thus the possibility of the DD decay modes are suppressed. There is a little similarity between this 'halo charmonium' picture and the so-called 'hadro-charmonium' picture, where the relatively localized color neutral cc kernel is surrounded by a light hadron cloud [30]. However, the advantage of the former resides in that the direct color interaction between the halo and the kernel provides an obvious binding mechanism, while in the 'hadroncharmonium picture, there need more phenomenological assumptions to describe the interaction between the meson cloud and the charmonium kernel.

V. CONCLUSION
To summarize, we use a new type of spatially extended hybrid-like operator to investigate the possible existence of exotic vector charmonia. In the non-relativistic approximation of these operators, the localized color octet charm quark-antiquark component is in the spin singlet state and separates from the chromo-magnetic field strength with a spatial distance. These operators couple preferably to a higher vector state X with a mass of 4.33(2) GeV when the distance increases. This observation indicates that the charm quark-antiquark pair of X may acquire a center-of-mass motion by recoiling against additional degrees of freedom depicted by the chromomagnetic field strength operator, which are necessarily gluonic in the quenched approximation. In this meaning, the state X can be taken as a hybrid-like vector charmonium. In addition, through a simultaneous multi-state fit to different correlation functions built from the vector current operator and the the hybrid operator mentioned above, the leptonic decay constant of X is tentatively determined to be roughly one order of magnitude smaller than f J/ψ , say, f X < 40 MeV, which gives a very small leptonic decay width Γ(X → e + e − ) < 40 eV. This is a very important characteristic parameter for X to be identified from experiments. Obviously the mass and the leptonic decay width of X are consistent with the production and decay properties of X(4260), which escapes so far from the direct measure in the e + e − annihilation. Based on the combined width Γ ee Γ J/ψπ + π − /Γ tot = 9.2 ± 1.0 eV of X(4260), if it can be assigned to the X state in this study, its decay branch fraction of J/ψππ mode can be larger than 20%, which also naturally explains why X(4260) is dominantly observed in J/ψππ. By virtue of the inner structure of X, X(4260) should also be observed in the h c ππ channel.