Light axial-vector and vector resonances $X(2100)$ and $X(2239)$

We study features of the resonances $X(2100)$ and $X(2239)$ by treating them as the axial-vector and vector tetraquarks with the quark content $ss \overline{s}\overline{s}$, respectively. The spectroscopic parameters of these exotic mesons are calculated in the framework of the QCD two-point sum rule method. Obtained prediction for the mass $m=(2067 \pm 84)~\mathrm{MeV}$ of the axial-vector state is in excellent agreement with the mass of the structure $X(2100)$ recently observed by the BESIII Collaboration in the decay $J/\psi \to \phi \eta \eta^{\prime}$ as the resonance in the $\phi \eta^{\prime}$ mass spectrum. We explore also the $S$-wave decays $X(2100) \to \phi\eta^{\prime}$ and $X(2100) \to \phi \eta$ using the QCD light-cone sum rule approach and technical methods of the soft-meson approximation. The width of the axial-vector tetraquark, $\Gamma=(130.2\pm 30.1)~\mathrm{MeV}$, saturated by these two decays is comparable with the measured full width of the resonance $X(2100)$. Our prediction for the vector $ss\overline{s}\overline{s }$ tetraquark's mass $\widetilde{m}=(2283\pm 114)~\mathrm{MeV}$ is consistent with the experimental result $2239.2 \pm 7.1 \pm 11.3~\mathrm{MeV}$ of the BESIII Collaboration for the mass of the resonance $X(2239)$.


I. INTRODUCTION
Hadrons with exotic structures and/or quantum numbers, which differ them from the conventionalqq mesons and qq ′ q ′′ baryons were and remain in agenda of the High Energy Physics community. Properties of the ordinary hadrons, i.e. their spectroscopic parameters as well as their strong, semileptonic and radiative transitions have been investigated in the framework of Quantum Chromodynamics (QCD), and successfully confronted with available experimental data. In the nonperturbative regime of momentum transfers, the relevant theoretical results have been obtained using methods and phenomenological models which use either the first principles of QCD or invoke additional assumptions about the internal structure and dynamics of hadrons.
At the same time, the QCD allows existence of not only the ordinary hadrons but also particles built of four, five, or more quarks, quark-gluon hybrids, and glueballs. The idea about the multi-quark nature of some observed particles was first applied to explain the unusual features of the light scalar mesons with masses m < 1 GeV [1]. The reason is that the nonet of scalar particles in the standard model of mesons should be realized as 1 3 P 0 quarkantiquark states. But masses of these scalars, in accordance with various model computations, are higher than 1 GeV. Moreover, the standard model could not correctly describe the mass hierarchy of the mesons inside the nonet. These problems can be evaded by assuming that the light scalars are four-quark exotic mesons, or at least contain substantial four-quark component. In the context of this scheme low masses of the scalar mesons, as well as the hierarchy inside of the nonet receive natural explanations. A recent model of the both light and heavy scalar nonets is based on suggestion about diquarkantidiquark structure of these particles which are mixtures of the spin-0 diquarks from (3 c , 3 f ) representation with spin-1 diquarks from (6 c , 3 f ) representation of the color-flavor group [2]. The spectroscopic parameters and width of the light scalar mesons f 0 (500) and f 0 (980) calculated by considering them as admixtures of the SU f (3) flavor octet and singlet tetraquarks are in a reasonable agreement with experimental data [3,4]. Other members of the light scalar nonet were also successfully explained as scalar particles with relevant diquark-antidiquark contents [5].
However, light quarks may not form stable tetraquarks: Theoretical studies proved that only tetraquarks composed of heavy and light diquarks may be stable against the strong decays. Thus, four-quark systems QQQQ and QQqq were studied in Refs. [6][7][8] by employing the conventional potential model with additive pairwise interaction of color-octet exchange type. Within this approach it was demonstrated that states QQqq may form the stable composites provided that the ratio m Q /m q is large enough. Experimental information on possible tetraquark candidates is also connected with the heavy resonances observed in various processes. Starting from discovery of the charmonium-like resonance X(3872) by Belle Collaboration [9], the exotic mesons are the objects of rapidly growing studies. Valuable experimental data collected during years passed from observation of the X(3872) resonance, as well as important theoretical achievements form now the physics of the exotic hadrons [10][11][12][13][14].
There are only few resonances seen in the experiments which may be considered as four-quark systems containing only the light quarks. One of such states is the famous structure Y (2175) discovered by the BaBar Collaboration in the process e + e − → γ ISR φf 0 (980) as a resonance in the φf 0 (980) invariant mass spectrum [15]. Existence of the Y (2175) later was confirmed by the BESII, Belle, and BESIII collaborations as well [16][17][18]. The mass and width of this state with spin-parities J P C = 1 −− is m = 2175 ± 10 ± 15 MeV and Γ = 58 ± 16 ± 20 MeV, respectively.
Other resonances which may be interpreted as light exotic mesons were observed recently by the BESIII Collaboration. Thus, the X(2239) was seen in the process e + e − → K + K − as a resonant structure in the cross section shape line [19]. The mass and width of this state were found equal to m = 2239.2 ± 7.1 ± 11.3 MeV and Γ = 139.8 ± 12.3 ± 20.6 MeV, respectively. The X(2100) was fixed in the process J/ψ → φηη ′ as a resonance in the φη ′ mass spectrum [20]. The collaboration studied the angular distribution of J/ψ → X(2100)η , but due to limited statistics could not clearly distinguish 1 + or 1 − assumption for the spin-parity J P of the X(2100). Therefore, the spectroscopic parameters of this resonance were determined using both of these assumptions. In the case J P = 1 − the mass and width of the X(2100) were measured to be m = 2002.1 ± 27.5 ± 21.4 MeV and Γ = 129 ± 17 ± 9 MeV. Alternatively, the assumption J P = 1 + led to the results m = 2062.8 ± 13.1 ± 7.2 MeV and Γ = 177 ± 36 ± 35 MeV.
Theoretical interpretations of these light resonances which may be considered as candidates for tetraquarks, as usua,l comprise all possible models and approaches available in high energy physics. Because the Y (2175) was discovered more than ten years ago, there are numerous and diverse articles in the literature devoted to its investigation. There are quite natural attempts to interpret it as an 2 3 D 1 excitation of the conventional ss meson [21,22]. Another traditional approach is to treat such states as dynamically generated resonances. As a dynamically generated state in the φKK system, the Y (2175) was examined in Ref. [23]. The similar dynamical picture may appear due to self-interaction between φ and f 0 (980) mesons as well [24]. Alternative explanations of the Y (2175) resonance's structure include a hybrid meson ssg, or a baryon-antibaryon qqsqqs state that couples strongly to the ΛΛ channel (for relevant references and other models, see Ref. [19]).
The resonance Y (2175) as a vector tetraquark with ssss or ssss content was explored in Refs. [25] and [26,27], respectively. In these works the authors used the QCD sum rule method and evaluated spectroscopic parameters of these states. The newly found structures X(2100) and X(2239) (hereafter X 1 and X 2 , respectively) were also analyzed as vector or axial-vector tetraquarks. Thus, in Ref. [28] the mass spectrum of the ssss tetraquark states was investigated within the relativized quark model. The authors concluded that the resonance X 2 can be assigned as a P -wave 1 −− ssss tetraquark. In the framework of the QCD sum rule method the X 1 resonance was studied in Refs. [29,30]. Predictions obtained there allowed the authors to interpret it as the axial-vector ssss tetraquark with the quantum numbers J P C = 1 +− . In accordance with Ref. [31], the X 1 may be identified as the second radial excitation of the conventional meson h 1 (1380).
As is seen, theoretical interpretations of observed light resonances are numerous and sometimes contradict to each other. There is a necessity to consider this prob-lem in a more detailed form and analyze not only spectroscopic parameters of the light resonances, but also to explore their decay channels and widths. In the present work we study the axial-vector and vector light tetraquarks ssss and compute their masses and couplings. By confronting theoretical predictions and experimental data we identify the observed resonances Y (2175), X 1 and X 2 with these tetraquark structures. It turns out that the resonance X 1 can be interpreted as a axial-vector tetraquark state. We calculate the width of the decays X 1 → φη ′ and X 1 → φη which are essential for our interpretation of the X 1 . Among the vector resonances Y (2175) and X 2 , parameters of the latter is closer to our result.
Calculations in the present paper are performed in the context of the QCD sum rule method, which is one of the powerful nonperturbative approaches in high energy physics [32,33]. The masses and couplings of the fourquark systems are evaluated using two-point QCD sum rules with an accuracy higher than in existing samples. To find the width of the decays X 1 → φη ′ and X 1 → φη we employ sum rules on the light cone and technical tools of the soft-meson approximation [34,35].
This paper is structured in the following form: In Sections II and III we analyze the spectroscopic parameters of the axial-vector and vector tetraquarks ssss and provide details of relevant sum rule calculations. In Sec. IV the strong couplings g X1φη ′ and g X1φη corresponding to the vertices X 1 φη ′ and X 1 φη are found using the QCD light-cone sum rule method. These coupling are required to evaluate the width of the decays X 1 → φη ′ and X 1 → φη, respectively. Section V contains summary of the obtained results and our conclusions.

II. MASS AND COUPLING OF THE AXIAL-VECTOR TETRAQUARK ssss
In this section we compute the mass and coupling of the axial-vector tetraquark T AV = ssss. As it has been emphasized above, to this end we use the QCD sum rules method which is based on first principles of QCD and allows one, via a quark-hadron duality assumption, to express physical parameters of hadrons in terms of the universal nonperturbative quantities, i.e. vacuum expectation values of local quark, gluon, and mixed operators. This method was successfully applied to explore parameters not only of conventional hadrons, but also to study various multi-quark systems [36].
To derive the required sum rules we consider the twopoint correlation function Π µν (p), which is defined by the formula where J µ (x) is the interpolating current for the axialvector tetraquark ssss. The choice of J µ (x) in one of the main operations in the sum rule computations. The tetraquark with content ssss and spin-parities J P C = 1 +− can be interpolated using different currents. The current that leads to a reliable prediction for the mass and coupling of the axial-vector state has the following form [29] Here a and b are the color indices and C is the charge conjugation operator. The sum rules necessary to calculate the mass m and coupling f of the T AV can be derived in accordance with prescriptions of the method, which require first to express the correlation function Π µν (p) using the tetraquark's physical parameters . We consider T AV as a ground-state particle, and after isolating the first term in Π Phys µν (p) get Equation (3) is obtained by saturating the correlation function with a complete set of J P = 1 +− states and carrying out the integration over x. Effects of higher resonances and continuum states are denoted above by dots.
To simplify further the correlator Π Phys µν (p), it is convenient to introduce the matrix element where ǫ µ is the polarization vector of the T AV state. Then the correlation function Π Phys µν (p) takes the simple form Equation (5) determines the physical or phenomenological side of the sum rules. The correlation function Π µν (p) calculated by employing the quark propagators constitutes the QCD side of the sum rules. It is given by the expression where S ab (x) is the s-quark propagator and In calculations we employ the x-space light-quark prop- where γ E ≃ 0.577 is the Euler constant, and Λ is the QCD scale parameter. We use also the notation G µν ab ≡ G µν A t A ab , A = 1, 2, . . . 8, and t A = λ A /2, with λ A being the Gell-Mann matrices.
The propagator (8) contains various light quark, gluon and mixed condensates of different dimensions. The term sg s σGs written down in Eq. (8) as well as other ones proportional to ss 2 , and ss g 2 s G 2 are obtained using the factorization hypothesis of the higher dimensional condensates. It is known, however, that the factorization assumption is not precise and violates is the case of higher dimensional condensates [37]. Thus, for the condensates of dimension 10 even an order of magnitude of such a violation is unclear. But, contributions to sum rules arising from higher dimensional condensates are very small, therefore, in what follows, we ignore uncertainties generated by this violation.
At the next stage we calculate the resultant fourx Fourier integrals in Π OPE µν (p). The correlation function Π OPE µν (p) obtained by this way contains two Lorentz structures which may be chosen to derive the sum rules. For our purposes terms ∼ g µν both in Π Phys µν (p) and Π OPE µν (p) are convenient, because scalar particles do not contribute to these terms. Afterwards we equate the corresponding invariant amplitudes Π Phys (p 2 ) and Π OPE (p 2 ), and find an expression in momentum space which, after some manipulations, can be used to derive the desired sum rules. Indeed, to suppress contributions of the higher resonances and continuum states we apply to both sides of the obtained equality the Borel transformation. The last operation to be carried out is continuum subtraction, which is achieved by invoking assumption on quark-hadron duality. After these manipulations the equality depends on auxiliary parameters of the sum rules M 2 and s 0 : M 2 is the Borel parameter appeared due to corresponding transformation, s 0 is the continuum subtraction parameter that separates the ground-state and higher resonances from each another.
To find the sum rules for m and f we need an additional expression which can be obtained by acting d/d −1/M 2 to the first equality. The sum rules for m and f have the perturbative and nonperturbative components. The nonperturbative components contain the quark, gluon, and mixed vacuum condensates, which appears after sandwiching relevant terms in Π OPE (p) between vacuum states. Our analytical results contain the nonperturbative terms up to dimension-20. We keep all of them in numerical computations bearing in mind that higher dimensional terms appear due to the factorization hypothesis as product of basic condensates, and do not encompass all dimension-20 contributions.
A very important problem of calculations is a proper choice for the Borel M 2 and continuum threshold s 0 parameters. These parameters are not arbitrary, but should meet some known requirements: At maximum of the Borel parameter the pole contribution (PC) has to constitute a fixed part of the correlation function, whereas at minimum of M 2 it must be a dominant contribution. We define PC in the form where Π(M 2 , s 0 ) is the Borel transformed and subtracted invariant amplitude Π OPE (p 2 ). The minimum of M 2 is fixed from convergence of the sum rules, i.e. at M 2 min contribution of the last term (or a sum of last few terms) cannot exceed, for example, 0.01 part of the whole result.
In the case of multi-quark hadrons at M 2 max one, as usual, requires PC > 0.2. There is an another restriction on the lower limit M 2 min : at M 2 min the perturbative contribution has to prevail over the nonperturbative one.
The sum rule predictions should not depend on the parameters M 2 and s 0 . But in real calculations m and f demonstrate sensitiveness to the choice of M 2 and s 0 . Hence, the parameters M 2 and s 0 have to be fixed in such a manner that to reduce this effect to a minimum. Performed analysis allows us to find the working regions which obey all the aforementioned constraints. In Fig. 1 we depict the pole contribution as functions of M 2 and s 0 : at M 2 = 1.4 the pole contribution is 0.68, whereas at M 2 = 2 it becomes equal to 0.39. The prediction for the mass m is plotted in Fig. 2, where one can see its weak dependence on the parameters M 2 and s 0 . The results for the spectroscopic parameters of the tetraquark T AV read: m = (2067 ± 84) MeV, f = (0.89 ± 0.11) × 10 −2 GeV 4 .
Theoretical errors in the sum rule computations appear due to different sources. The auxiliary parameters M 2 The result obtained for the mass of the axial-vector tetraquark T AV is in excellent agreement with the mass of the structure X 1 reported by the BESIII Collaboration. Therefore, it is possible to identify T AV with the resonance X 1 . Our conclusion is also in accord with previous theoretical predictions obtained by means of the QCD sum rules method. Thus, the mass of the resonance X 1 was estimated in Refs. [29,30] m = 2000 +100 −90 MeV, m = (2080 ± 120) MeV, (12) respectively. As is seen, all these calculations support the assumption on the axial-vector tetraquark nature of the structure X 1 . But one needs to explore its decay channels X 1 → φη ′ and X 1 → φη, and find width of this resonance: only after successful comparison with experimental data it is legitimate to make more strong conclusion about X 1 .
We are going to address this problem in Sec. IV.

III. SPECTROSCOPIC PARAMETERS OF THE VECTOR TETRAQUARK ssss
In the previous section we have explored the axialvector tetraquark T AV and identified it as a candidate for the resonance X 1 . But there are two other light states which should be classified within the four-quark picture. In the present section we are going to analyze the vector tetraquark T V = ssss with the quantum numbers J P C = 1 −− and compare the obtained result for its mass with the experimental information of BaBar and BESIII collaborations.
Calculations of the T V tetraquark's mass m and coupling f do not differ considerably from ones fulfilled in the previous section. There are only some qualitative differences on which we want to concentrate. First of all, the interpolating current for the vector state is defined by the expression [27] The physical side of the sum rule is given by Eq. (5) with evident replacements. The correlation function Π OPE µν (p) that determines the QCD side of the sum rule has the following expression The remaining operations have been explained above. Therefore we present only final results of performed analysis. The working windows for the Borel and continuum threshold parameters in the case of the vector tetraquark T V are determined by the intervals It is seen that these regions differ from ones presented in Eq. (10) by only small shift of the parameter s 0 . The windows (10) comply all constraints necessary in the sum rule computations. In fact, at M 2 = 1.4 the pole contribution is 60%, whereas at M 2 = 2 it is equal to 30% of the whole result. Convergence of the sum rules is also satisfied. The mass and coupling of the vector tetraquark T V are: Fig. 3 we plot the spectroscopic parameters m and f as functions of M 2 and s 0 .
Comparing the mass of the vector state T V and experimental information on the resonances Y (2175) and X 2 , one can see that it can be identified with the X 2 . In fact, difference between the masses of T V and X 2 is approximately 60 MeV smaller than between T V and Y (2175). The similar conclusion was drawn also in Ref. [28]. The mass m X2 = 2227 MeV of the four-quark vector system ssss found there is consistent with BESIII data.
The mass of the vector tetraquark ssss was computed using the QCD sum rule method in Refs. [30] and [27] as well. The prediction for the mass of this four-quark meson m = (3080 ± 110) MeV made in Refs. [30] disfavors classifying it as the resonance Y (2175). Comparing this result with recent measurements of the BESIII Collaboration, we see that it also cannot be assigned to be the resonance X 2 . To study vector tetraquarks with the ssss content, in Ref. [27] the authors constructed two independent interpolating currents which couple to J P C = 1 −− states. These currents led to slightly different predictions m 1 = (2410 ± 250) MeV and m 2 = (2340±170) MeV. In accordance with [27] the first state might correspond to a structure in the φf 0 (980) invariant mass spectrum at around 2.4 GeV. The second one was interpreted in Ref. [27] as the resonance Y (2175) but, from our point of view, it is closer to the structure X 2 .
Within the framework of the QCD sum rule method the decay X 1 → φη ′ [and X 1 → φη] can be investigated by means of different approaches. In fact, a key quantity to calculate the width of this decay is the coupling g X1φη ′ describing the strong interaction in the vertex X 1 φη ′ . The coupling g X1φη ′ can be evaluated using, for example, the QCD three-point sum rule method. Alternatively, one can extract it from the relevant QCD light-cone sum rule (LCSR), which has some advantages when calculating tetraquark-meson-meson vertices containing light mesons. The reason is that the LCSRs for tetraquark-meson-meson vertices differ from ones involving only conventional mesons. Thus, the LCSR for vertices of conventional mesons depends on various distribution amplitudes (DAs) of one of the final mesons, which encode all information about nonperturbative dynamical properties of the meson. In the case of the tetraquarkmeson-meson vertices due to four-quark nature of the tetraquark, after contracting relevant quark fields instead of DAs of a the final meson the sum rule contains only local matrix elements of this meson. Then to satisfy the four-momentum conservation at vertices the momentum of a final light meson should be set q = 0. This leads to crucial changes in the calculational scheme, because now one has to accompany the LCSR method with technical tools of the soft-meson approximation [35,39].
Let us consider the dominant process X 1 → φη ′ in a detailed form. The second decay mode X 1 → φη , as we shall see below, can be analyzed in the same manner. The starting point to explore the decay X 1 → φη ′ is the correlation function where J φ µ (x) is the interpolating current of the φ meson Following the standard recipes, we write down Π µν (p, q) in terms of the physical parameters of the particles X 1 , φ and η ′ Π Phys where p ′ and p, q are momenta of the initial and final particles, respectively. In Eq. (19) contributions of excited resonances and continuum states are indicated by dots. By utilizing the matrix elements one can considerably simplify Π Phys µν (p, q). The matrix element 0|J φ µ (x)|φ(p) is expressed in terms of φ meson's mass m φ , decay constant f φ and polarization vector ε µ . The matrix element of the vertex X 1 φη ′ is written down using the strong coupling g X1φη ′ which has to be evaluated from the sum rule. In the soft limit q → 0 we get p ′ = p, as a result instead of two-variable Borel transformation we have to perform one-variable Borel transformation, which yields where m 2 = (m 2 φ + m 2 )/2. In Eq. (21) we still keep p ν = p ′ µ to make clear the Lorentz structure of the obtained expression. To derive the LCSR for the strong coupling g X1φη ′ we will employ the structure ∼ g µν .
In the soft approximation the physical side of the sum rule has more complicated structure than in the case of full LCSR method. The complications are connected with behavior of contributions arising from higher resonances and continuum states in the soft limit. The problem is that in the soft limit some of these contributions even after the Borel transformation remain unsuppressed and appear as contaminations in the physical side [35]. Therefore, before performing the continuum subtraction in the final sum rule they should be removed by means of some operations. This problem is solved by acting on the physical side of sum rule by the operator [35,40] that singles out the ground-state term. It is natural that the same operator P(M 2 , m 2 ) should be applied also to the QCD side of the sum rule. But before these manipulations the correlation function Π OPE µν (p, q) has to be calculated in the soft-meson approximation and expressed in terms of the η ′ meson's local matrix elements.
In the soft limit Π OPE µν (p) is given by the formula where α and β are the spinor indices. It is seen that Π OPE µν (p) really depends on local matrix elements of the η ′ meson. But these matrix elements should be converted to forms suitable to express them in terms of the standard matrix elements of the η ′ meson. To this end, we continue calculations by employing the expansion where Γ j is the full set of Dirac matrices Then operators s(0)Γ j s(0) , as well as ones appeared due to G µν insertions from propagators S(±x), generate standard local matrix elements of the η ′ meson. Substituting Eq. (23) into the expression of the correlation function and carrying out the summation over color indices in accordance with rules described in a detailed form in Ref.
[39], we find local matrix elements of the η ′ meson that contribute to Π QCD (p). Performed analysis demonstrates that in the softmeson approximation only twist-3 matrix element η ′ |siγ 5 s|0 gives non-zero contribution to the correlation function Π OPE µν (p). The matrix elements of the η and η ′ differ from ones of other pseudoscalar mesons: This is connected with mixing phenomena in the η − η ′ system. Thus, due to the mixing both the η ′ and η mesons have ss components. Of course, ss is dominant for the η ′ meson, whereas it plays a subdominant role in the η meson's quark content. Nevertheless, through the strange components both of these mesons can appear in the final state of the decays X 1 → φη ′ and X 1 → φη.
The mixing in the η − η ′ system can be described in different basis: For our purposes, the quark-flavor basis is more convenient than the octet-singlet basis of the flavor SU f (3) group. The quark-flavor basis was used in our previous papers to study different exclusive processes with η ′ and η mesons [41][42][43]. In the quark-flavor basis the twist-3 matrix element η ′ |siγ 5 s|0 can be written down in the following form where the parameter h s η ′ is defined by the equality In Eq. (25) m η ′ and f s η ′ are the mass and s-component of the η ′ meson decay constant. Here the A η ′ is the matrix element which appear due to U (1) axial-anomaly. The parameter h s η ′ may be computed by employing Eqs. (24) and (25), but we use its phenomenological value extracted from analysis of relevant exclusive processes. Thus, we have where ϕ = 39 • .3 ± 1 • .0 is the mixing angle in the quarkflavor basis.
Our result for the Borel transform of the invariant function Π OPE (p 2 ) corresponding to the structure ∼ g µν reads where It is worth noting that the spectral density ρ pert. (s) is computed as the imaginary part of the relevant term in the correlation function. The Borel transform of nonperturbative terms are found directly from Π OPE (p 2 ) and includes terms up to dimension six. After acting the operator P(M 2 , m 2 ) to Π OPE (M 2 ) one can perform the continuum subtraction. This implies replacement ∞ → s 0 in the first term, whereas terms ∼ (M 2 ) 0 and ∼ 1/M 2 should be left in their original forms [35]. The width of the decay X 1 → φη ′ is determined by the formula where In numerical computations, the parameters M 2 and s 0 are varied within the limits The mass of the final-state mesons φ and η ′ are borrowed from Ref. [38] m φ = (1019.461 ± 0.019) MeV, m η ′ = (957.78 ± 0.06) MeV, f φ = (215 ± 5) MeV.
The X 1 → φη ′ is the main S-wave decay channel of the tetraquark X 1 . The partial width of the second process X 1 → φη can be easily evaluated by employing expressions obtained in the present section. The differences between two decays stem from the twist-3 matrix element, which for this decay is given by the formula 2m s η|siγ 5 s|0 = −h s sin ϕ, Let us note that |g X1φη | has been extracted from the sum rule at s 0 ∈ [5.8, 6.8] GeV 2 .
Saturating the full width of the X 1 resonance by these two decays we get: This estimate does not coincide with full width of the resonance X 1 , but is comparable with it.

V. SUMMARY AND CONCLUSIONS
In the present work we have studied the axial-vector and vector tetraquarks with the quark content ssss. The mass m = (2067 ± 84) MeV of the axial-vector state obtained in the present work is in excellent agreement with measurements of the BESIII Collaboration. The width of this state Γ = (130.2 ± 30.1) MeV within both theoretical and experimental errors is consistent with the data. These facts have allowed us to interpret the resonance X(2100) discovered recently the BESIII Collaboration as an axial-vector state with quark content ssss.
The vector ssss tetraquark with the mass m = (2283 ± 114) MeV can be identified with the structure X(2239) rather than with the resonance Y (2175). There is still the light resonance Y (2175) which in the present scheme may be considered as a conventional vector meson, because its mass is small to classify it as a vector tetraquark. One should take into account also a possible structure in the φf 0 (980) invariant mass spectrum at 2.4 GeV. In our present work we have tried to answer questions on nature of two light resonances. It is evident that the whole family of such structures deserves further detailed investigations.