QCD sum rule studies on the $s s \bar s \bar s$ tetraquark states with $J^{PC} = 1^{+-}$

We apply the method of QCD sum rules to study the structure $X$ newly observed by the BESIII Collaboration in the $\phi \eta^\prime$ mass spectrum in 2.0-2.1 GeV region in the $J/\psi \rightarrow \phi \eta \eta^\prime$ decay. We construct all the $s s \bar s \bar s$ tetraquark currents with $J^{PC} = 1^{+-}$, and use them to perform QCD sum rule analyses. One current leads to reliable QCD sum rule results and the mass is extracted to be $2.00^{+0.10}_{-0.09}$ GeV, suggesting that the structure $X$ can be interpreted as an $s s \bar s \bar s$ tetraquark state with $J^{PC} = 1^{+-}$. The $Y(2175)$ can be interpreted as its $s s \bar s \bar s$ partner having $J^{PC} = 1^{--}$, and we propose to search for the other two partners, the $s s \bar s \bar s$ tetraquark states with $J^{PC} = 1^{++}$ and $1^{-+}$, in the $\eta^\prime f_0(980)$, $\eta^\prime K \bar K$, and $\eta^\prime K \bar K^*$ mass spectra.


I. INTRODUCTION
Recently, the BESIII Collaboration reported their observation of a new structure X in the φη ′ mass spectrum in 2.0-2.1 GeV region, when studying the J/ψ → φηη ′ decay [1]. This experiment gives two possibilities: 1. After assuming X to have the spin-parity quantum numbers J P = 1 − , its mass and decay width are determined to be 2. After assuming X to have the spin-parity quantum numbers J P = 1 + , its mass and decay width are determined to be Here, the first uncertainties are statistical and the second systematic. The significances are 5.3σ and 4.9σ, respectively, so these two assumptions can not be distinguished at BESIII. One possible theoretical explanation is to interpret it as an isoscalar axial-vector meson with I(J P ) = 0(1 + ), the second radial excitation of h 1 (1380) [2].
In our previous studies [18,19] we have applied the method of QCD sum rules to systematically study the ssss tetraquark states with J P C = 1 −− . There we found two independent ssss tetraquark currents with J P C = 1 −− , and the masses are evaluated to be 2.34 ± 0.17 GeV and 2.41 ± 0.25 GeV, not far from each other [19]. These two values are both significantly larger than the first mass value listed in Eq. (1), suggesting that the structure X is difficult to be interpreted as an ssss tetraquark state of J P C = 1 −− . Instead, the Y (2175) can be well interpreted as an ssss tetraquark state of J P C = 1 −− [18,19]. Moreover, the above two mass values are extracted from two diagonalized currents, which do not strongly correlate to each other and may couple to two different physical states: one is the Y (2175), and the other is around 2.4 GeV. There have been some evidences for the latter structure in the previous experiments [11,14,15,17], and we refer to Ref. [19] for detailed discussions.
In the present study we follow the same approach to study the ssss tetraquark states with J P C = 1 +− , and examine whether the structure X can be explained. Again, we shall find that there are two independent ssss tetraquark currents with J P C = 1 +− , which we shall use to perform QCD sum rule analyses. The internal structures of exotic hadrons are always complicated. For each internal structure we can construct the relevant interpolating current, and there are usually many interpolating currents when studying multiquark states. In this case, the only two independent currents make it possible to study their mixing. Note that we have done this in Ref. [19] when studying the ssss tetraquark states with J P C = 1 −− . By doing this we can carefully examine the relations between physical states and the relevant interpolating currents, and further understand the internal structures of exotic hadrons. This paper is organized as follows. In Sec. II, we systematically construct the ssss tetraquark currents with J P C = 1 +− , using both diquark/antidiquark fields and quark-antiquark pairs. These currents are then used to perform QCD sum rule analyses in Sec. III, and numerical analyses in Sec. IV. Their mixing are investigated in Sec. V. Sec. VI is a summary.

II. INTERPOLATING CURRENTS
The ssss tetraquark currents with the quantum numbers J P C = 1 −− have been systematically constructed in Ref. [18]. See also Refs. [20][21][22] where many other vector and axial-vector tetraquark currents are systematically constructed. In this section we follow the same approach to construct the ssss tetraquark currents with the quantum numbers J P C = 1 +− . We find two non-vanishing diquark-antidiquark currents: , where a and b are color indices; C = iγ 2 γ 0 is the chargeconjugation operator; the sum over repeated indices is taken. These two diquark-antidiquark currents are independent of each other. Recalling that the diquark fields s T a Cs b /s T a Cγ µ γ 5 s b /s T a Cγ µ s b /s T a Cσ µν γ 5 s b have the quantum numbers J P = 0 − /1 − /1 + /1 ± , respectively, the first current η 1µ only contains excited diquark and antidiquark fields, but the second one η 2µ contains (at least) one ground-state diquark/antidiquark field. Hence, η 2µ has a more stable internal structure and may lead to better sum rule results.
Besides the above diquark-antidiquark currents, we find that there are four mesonic-mesonic currents: The following relations can be verified by using the Fierz transformation, so the number of independent mesonicmesonic currents is also two: Moreover, we can use the Fierz transformation to relate the diquark-antidiquark and mesonic-mesonic currents: Therefore, these two constructions are equivalent.
In the following we shall use η 1µ and η 2µ to perform QCD sum rule analyses.

III. QCD SUM RULE ANALYSES
QCD sum rules [23][24][25], a powerful and successful nonperturbative method, have been widely applied to study various exotic hadrons [26][27][28][29][30][31][32][33][34][35][36][37]. In this method we calculate the two-point correlation function at both the hadron and quark-gluon levels: At the hadron level, we can express Π(q 2 ) in the form of the dispersion relation with a spectral function ρ(s): Then we adopt a parametrization of one pole dominance and a continuum contribution: where X is the ground state. At the quark-gluon level, we insert η 1µ and η 2µ into Eq. (9), which are then calculated using the method of operator product expansion (OPE). After performing the Borel transformation at both the hadron and quark-gluon levels, we obtain Then we approximate the continuum using the spectral density of OPE above a threshold value s 0 , and obtain the following sum rule equation (13) Finally, we can use this equation to calculate M X , the mass of the ground state X, through For the currents η 1µ and η 2µ , we have calculated the OPE up to dimension twelve.
Explicitly, we have calculated the perturbative term, the gluon condensate g 2 s GG , the quark condensate ss , the quark-gluon condensate g ss σGs , and their combinations g 2 s GG ss , g 2 s GG ss 2 , g 2 s GG g ss σGs , g 2 s GG ss g ss σGs , ss 2 , ss 3 , ss 4 , g ss σGs 2 , ss g ss σGs , and ss 2 g ss σGs . The results for η 1µ and η 2µ are shown in Eqs. (15) and (16), respectively. For completeness, we have also calculated the sum rules for the off-diagonal term: After performing the Borel transformation to Π η1η2 (q 2 ), we obtain Π η1η2 (M 2 B ) whose explicit expression is shown in Eq. (17).

IV. NUMERICAL ANALYSES
In this section we use the currents η 1µ and η 2µ to perform numerical analyses, for which we use the following g 2 s GG = (0.48 ± 0.14) GeV 4 , g sq σGq = −M 2 0 × qq , M 2 0 = (0.8 ± 0.2) GeV 2 . Different from Ref. [18] where there are two independent ssss tetraquark currents with J P C = 1 −− leading to similar QCD sum rule results, in the present study we find that the two ssss tetraquark currents with J P C = 1 +− , η 1µ and η 2µ , lead to totally different sum rule results. This can be clearly seen in Fig. 1, where we show the Borel transformed correlation functions Π η1η1 (M 2 B , s 0 ) and Π η2η2 (M 2 B , s 0 ) as functions of the threshold value s 0 . We find that Π η1η1 (M 2 B , s 0 ) is negative, and so non-physical, in the region s 0 < 10 GeV 2 . Hence, it can not strongly couple to any structure that is smaller than 3.0 GeV. The situation for η 2µ is different since Π η2η2 (M 2 B , s 0 ) is positive and well defined. This behavior seems to be reasonable because η 1µ only contains excited diquark and antidiquark fields, while η 2µ contains (at least) one ground-state diquark/antidiquark field and so more stable.
In the following we shall only use the current η 2µ to perform numerical analyses. After carefully investigating a) the OPE convergence, b) the pole contribution, and c) the mass dependence on the two free parameters M B and s 0 , we obtain reliable QCD sum rule results in the regions 1.6 GeV 2 < M 2 B < 2.0 GeV 2 and 5.5 GeV 2 < s 0 < 6.5 GeV 2 : • First we study the convergence of the operator product expansion. After taking s 0 to be ∞ and the integral subscript 16m 2 s to be zero, we obtain the numerical series of the OPE as a function of From this equation, we clearly see that the OPE convergence is quite good: the dimension 12 terms (∼ M −2 B ) are significantly smaller than the dimension 10 terms (∼ M 0 B ), which are again significantly smaller than the dimension 8 terms (∼ M 2 B ). Numerically, we show the ratio in Fig. 2 as a function of the Borel mass M B . We find it to be smaller than 5% in the regions 1.6 GeV 2 < M 2 B < 2.0 GeV 2 and 5.5 GeV 2 < s 0 < 6.5 GeV 2 .
• Then we study the pole contribution, defined as We show it as a function of the Borel mass M B in Fig. 3. We find it to be 30% < PC < 58%  in the regions 1.6 GeV 2 < M 2 B < 2.0 GeV 2 and 5.5 GeV 2 < s 0 < 6.5 GeV 2 . This amount of pole contribution is acceptable when one applies the method of QCD sum rules to study multiquark states. • Finally we study the mass dependence on the two free parameters, the Borel mass M B and the thresh-old value s 0 . To clearly see this, we show M η2 , the mass extracted from the current η 2µ , in Fig. 4 as a function of M B and s 0 .
In the left panel we show M η2 as a function of the Borel mass M B , and find it quite stable in the Borel window 1.6 GeV 2 < M 2 B < 2.0 GeV 2 . Comparing this figure with Fig. 2 and Fig. 3, we find that one can obtain a still larger pole contribution by choosing a smaller Borel mass (as shown in Fig. 3), but at the same time the convergence of OPE would become worse (as shown in Fig. 2) and the mass dependence on the Borel mass would become stronger (as shown in the left panel of Fig. 4). Considering all these behaviours, we find it suitable to fix the Borel window to be 1.6 GeV 2 < M 2 B < 2.0 GeV 2 . In the right panel we show M η2 as a function of the threshold value s 0 . We find that the mass curves moderately depend on the threshold value s 0 . Especially, we evaluate the mass to be 1.94 GeV< M η2 < 2.06 GeV in the region 5.5 GeV 2 < s 0 < 6.5 GeV 2 . This uncertainty is about 6%, quite typical in QCD sum rule studies.
Summarizing the above analyses, we have used the ssss tetraquark current η 2µ with J P C = 1 +− to perform QCD sum rule analyses. After carefully choosing the working regions to be 1.6 GeV 2 < M 2 B < 2.0 GeV 2 and 5.5 GeV 2 < s 0 < 6.5 GeV 2 , we extract the mass to be

V. MIXING OF CURRENTS
In the previous section we have used the two single ssss tetraquark currents with J P C = 1 +− , η 1µ and η 2µ , to perform QCD sum rule analyses. In this section we further study their mixing. We shall follow the procedures used in Ref. [19], where the mixing of two ssss tetraquark currents with J P C = 1 −− is carefully investigated.
To do this, first let us examine how large is the overlap between η 1µ and η 2µ . We show the off-diagonal term Π η1η2 (M 2 B ) in the left panel of Fig. 5 as a function of the Borel mass M B , compared with Π η1η1 (M 2 B ) and Π η2η2 (M 2 B ). This term has been defined in Eq. (18) and its explicit expression has been given in Eq. (17). From these figures, it is difficult to judge whether the off-diagonal term is important or not, because it is neither too large nor too small. Hence, we further diagonalize the following matrix at around M 2 B = 1.8 GeV 2 and s 0 = 6.0 GeV 2 Then we obtain the mixing angle θ = 2.7 o and two new currents J 1µ and J 2µ defined as: J 2µ = sin θ η 1µ + cos θ i η 2µ .
These two new currents do not strongly correlate to each other in the region 1.6 GeV 2 < M 2 B < 2.0 GeV 2 , as shown in the right panel of Fig. 5.
We use J 1µ and J 2µ to perform QCD sum rule analyses, and the results obtained are almost the same as those extracted from η 1µ and η 2µ : a) J 1µ does not lead to reliable QCD sum rule results because Π J1J1 (M 2 B , s 0 ) is negative in the region s 0 < 10 GeV 2 , and b) the mass extracted from J 2µ is about 2.00 GeV, the same as the one extracted from η 2µ .

VI. SUMMARY AND DISCUSSIONS
In this work we systematically construct all the ssss tetraquark currents with the quantum numbers J P C = 1 +− . We find there are two independent ones (η 1µ and η 2µ ), which are then used to perform QCD sum rule analyses. The sum rules extracted from η 1µ and η 2µ are much different from each other: a) η 1µ does not lead to reliable results because Π η1η1 (M 2 B , s 0 ) is negative, and so non-physical, in the region s 0 < 10 GeV 2 , and b) η 2µ leads to reliable results and the mass is extracted to be 2.00 +0. 10 −0.09 GeV, consistent with the second mass value listed in Eq. (2), 2062.8 ± 13.1 ± 4.2 MeV. The mixing between η 1µ and η 2µ has been taken into account, and the results are the same. Hence, our results suggest that the structure X observed at BESIII [1] has the spin-parity quantum numbers J P = 1 +− , and it can be interpreted as an ssss tetraquark state.
Recalling that in Refs. [18,19] we have systematically investigated the ssss tetraquark states with J P C = 1 −− . There we also found two independent ssss tetraquark currents with J P C = 1 −− , but they lead to similar sum rule results, i.e., the masses are extracted to be 2.34±0.17 GeV and 2.41 ± 0.25 GeV, not far from each other [19]. These two values are both larger than the first mass value listed in Eq. (1), 2002.1 ± 27.5 ± 15.0 MeV, suggesting that the structure X observed at BESIII [1] is difficult to be interpreted as an ssss tetraquark state of J P C = 1 −− .
Besides these isoscalar ssss tetraquark states, in Ref. [22] we have systematically constructed all the isovector tetraquark currents of I G J P C = 1 + 1 +− /1 + 1 −− /1 − 1 ++ /1 − 1 −+ , and found a one-to-one correspondence among them, i.e., for every tetraquark current of I G J P C = 1 + 1 +− one can construct a corresponding one of I G J P C = 1 + 1 −− , etc. These tetraquark   currents have been used to perform QCD sum rule analyses in Refs. [20,22], and the results are summarized in Table I, where q denotes an up or down quark, and s denotes a strange quark. Note that the sum rule results do not have the above one-to-one correspondence, for examples: a) there are four qqqq currents and four qsqs currents with I G J P C = 1 + 1 +− , and the masses extracted from these currents are all around 1.47-1.66 GeV; b) there are also four qqqq currents and four qsqs currents with I G J P C = 1 + 1 −− , but the masses extracted from the former four are around 1.60-1.73 GeV and the masses extracted from the latter four are around 1.91-2.13 GeV. This behaviour may relate to their internal structures, such as internal orbital excitations.
Similarly, there is a one-to-one correspondence among the ssss tetraquark currents with J P C = 1 +− /1 −− /1 ++ /1 −+ . Those with J P C = 1 +− and 1 −− have been used to perform QCD sum rule analyses in the present study as well as in Refs. [18,19]. The results are also summarized in Table I. From this table, we propose to search for the ssss tetraquark states with J P C = 1 ++ and 1 −+ in future experiments. We are now studying them following the same approach used in the present study. Their masses may also be around 2.0-2.4 GeV, and the possible decay channels to observe them are η ′ f 0 (980), η ′ KK, and η ′ KK * , etc.
When studying light tetraquark states, it is usually difficult to determine the experimental signal as a genuine four-quark state other than a conventionalqq meson, because the signal always has a quite large decay width. For example, besides the ssss tetraquark state of J P C = 1 +− , there are many other possible interpretations to explain the structure X, such as the second radial excitation of h 1 (1380) having I(J P ) = 0(1 + ) [2]. However, with the large amount of data collected at BE-SIII, this problem may be partly solved, and it is promising to continuously study light exotic hadrons. Together with those studies on charmonium-like XY Z states, our understudying on the nature of exotic hadrons can be significantly improved. I: Masses extracted from the vector and axial-vector tetraquark currents. Possible experimental candidates are listed for comparisons. We use q to denote an up or down quark, and s to denote a strange quark. The mass value 2.00 +0. 10 −0.09 GeV denoted by † is obtained in the present study.