## Fully strange tetraquark resonant states as the cousins of X (6900)

Yao Ma , 1, ∗ Wei-Lin Wu , 2, † Lu Meng , 3, ‡ Yan-Ke Chen , 2, § and Shi-Lin Zhu 1, ¶

1 School of Physics and Center of High Energy Physics, Peking University, Beijing 100871, China 2 School of Physics, Peking University, Beijing 100871, China

3 Institut f¨r Theoretische Physik II, Ruhr-Universit¨t Bochum, D-44780 Bochum, Germany u a

We conduct systematic calculations of the S-wave fully strange systems with 'normal' ( J PC = 0 ++ , 1 + -, 2 ++ ) and 'exotic' ( J PC = 0 + -, 1 ++ , 2 + -) C-parities, which are the strange analogue of the fully charmed tetraquark state X (6900) . Within a constituent quark potential model, we employ the Gaussian expansion method to solve the four-body Schr¨ odinger equation and the complex scaling method to identify resonant states. We obtain a series of resonant states and zero-width states in the mass range of 2.7 to 3.3 GeV, with their widths ranging from less than 1 MeV to about 50 MeV. Their rms radii strongly indicate that they are compact tetraquark states. Among these states, the T 4 s, 2 ++ (2714) may be the most likely one to be observed experimentally. We urge the experimental exploration of the 2 ++ ssss ¯¯ state around 2.7 GeV in the ϕϕ channel. Since the lowest S-wave ssss ¯¯ state is around 2.7 GeV, the compact P-wave ssss ¯¯ states are expected to be heavier. Hence, ϕ (2170) and X (2370) are unlikely to be compact tetraquark states.

## I. INTRODUCTION

- · An evidence for a structure X (2100) in the ϕη ′ mass spectrum of the J/ψ → ϕηη ′ decay was reported in 2019 by the BESIII Collaboration [24].

Since 2003, the interest in multiquark states has been reignited by the experimental discovery of a series of candidates (for recent reviews, see Refs. [1-11]). These discoveries have significantly enriched the hadron spectrum. The study of their properties and structures enhances our understanding of the non-perturbative nature of the strong interaction.

Inspired by their many strangeness decays, there are several promising candidates for the fully strange tetraquark states.

- · The state ϕ (2170) (was Y (2175) ) was first reported in 2006 by the BaBar Collaboration in the process of e + e -→ ϕ (1020) f 0 (980) [12-15]. Later it was confirmed by the BES and BESIII Collaborations in the J/ψ → ηϕf 0 (980) process [16-18], and the Belle Collaboration in the e + e -→ ϕπ + π -and e + e -→ ϕf 0 (980) processes [19].
- · The state X (2370) was first observed in 2010 by the BESIII Collaboration in the π + π - ′ η invariant-mass spectrum [20]. Later it was confirmed in 2019 by the BESIII Collaboration in the KKη ¯ ′ invariant-mass distribution [21], which suggests that the X (2370) likely contains four valence s or ¯ s quarks.
- · The pseudoscalar state X (2500) was first observed in 2016 by the BES Collaboration in the process J/ψ → γϕϕ [22], together with the scalar state f 0 (2100) , and three tensor states, the f 2 (2010) , f 2 (2300) and f 2 (2340) .
- · The X (2239) was seen in 2018 by the BESIII Collaboration in the cross section of e + e -→ K K + -[23].

∗ yaoma@pku.edu.cn

† wlwu@pku.edu.cn

‡ lu.meng@rub.de

§ chenyanke@stu.pku.edu.cn

¶ zhusl@pku.edu.cn

- · Another structure in the ϕf 0 cross section around 2.4 GeV, labelled as X (2436) [25], was reported in 2007 by the BaBar Collaboration [13]. It is also seen by the BES Collaboration [16] and Belle Collaboration [19], and evidenced by the BESIII Collaboration [17]. However, the latest BESIII results indicate that its statistical significance is no more than 2 σ [26].

Among the aforementioned possible T ssss ¯¯ candidates, all the established ones have negative parity. The existence of the positive parity states remains uncertain. Additionally, the LHCb Collaboration recently discovered a fully charmed tetraquark candidate X (6900) [27], which was later confirmed by the CMS and ATLAS [28, 29]. The QCD Lagrangian leads to the same quark gluon interaction if one omits the quark mass difference. The discovery of X (6900) implies the possible existence of the strange analogue T ssss ¯¯ .

Inspired by these experimental observations, various theoretical investigations have been conducted on the fully strange tetraquark systems. Regarding the ϕ (2170) state with J PC = 1 --, some authors used the QCD sum rule [30-33], quark model [34] and diquark-antidiquark scheme [35] to interpret it as an ssss ¯¯ tetraquark state, although other interpretations were not excluded [36-42]. Similarly, the X (2370) state was also explained as an ssss ¯¯ tetraquark state with 0 -+ [43, 44]. The pseudoscalar state X (2500) was interpreted as a P -wave ssss ¯¯ tetraquark state of J PC = 0 -+ [43-45]. In Ref. [46], the X (2500) , X (2239) , f 2 (2300) were assigned as 0 -+ , 1 --, 2 ++ ssss ¯¯ tetraquark states, respectively. Moreover, the structure X (2100) was interpreted as a 1 + -ssss ¯¯ state [44, 47, 48]. In Ref. [49], the resonant states X (2100) and X (2239) were treated as the axial-vector and vector tetraquarks with the quark content ssss ¯¯ , respectively. The structure X (2436) was explained as a P -wave ssss ¯¯ tetraquark state of J PC = 1 --[33, 44]. The ssss ¯¯ tetraquark states with other quantum numbers were also expolred [44, 45, 50-54].

So far, there is no experimental conclusion regarding the

existence of positive parity ssss ¯¯ states, and theoretical investigations on these states are also limited. Although some theoretical calculations on the fully strange tetraquark system in quark models have been conducted [8, 46, 50, 52], they are relatively sparse and less refined compared to those on tetraquark systems containing heavy quarks. The studies in Refs. [46, 50] focus exclusively on the diquark-antidiquark picture, obtaining only the states above the dimeson thresholds. These states would become scattering states or resonant states once the dimeson configuration is incorporated. Only the resonant states have physical significance and should be identified using nontrivial approaches. In Refs. [8, 52], the ssss ¯¯ tetraquark resonant states were investigated via the real scaling method. This method detects resonant states by avoiding the cross pattern of energy levels across a parameter controlling the spatial extension of the wave function, which is usually quite rough.

To address these issues, we will conduct systematic calculations of the positive parity S-wave fully strange system. We will employ the Gaussian expansion method (GEM) [55] to solve the four-body Schr¨ odinger equation within a nonrelativistic constituent quark potential model. Additionally, we will adopt the complex scaling method (CSM) [56-58] to identify genuine resonant states from meson-meson scattering states. In our recent works, we have incorporated both the dimeson and diquark-antidiquark configurations to carry out benchmark calculations for tetraquark bound states [59-61] and resonant states [62-64]. We will continue to include both configurations for the fully strange tetraquark system. In this study, we will investigate both the 'normal' C-parity states ( J PC = 0 ++ , 1 + -, 2 ++ ) and the 'exotic' C-parity states ( J PC = 0 + -, 1 ++ , 2 + -), with the 'exotic' C-parity systems referring to those that cannot be composed of two S-wave ground quarkonia.

This paper is arranged as follows. In Sec. II, we introduce the theoretical framework, including the Hamiltonian in the constituent quark potential, the construction of the wave function, the complex scaling method, and the approach to analyzing the spatial structures. In Sec. III, we present the numerical results for the properties of the fully strange resonant states. Finally, we give a brief summary and discussion in Sec. IV.

## II. THEORETICAL FRAMEWORK

## A. Hamiltonian

The nonrelativistic Hamiltonian of a four-quark system reads

$$H = \sum _ { i } ^ { 4 } \left ( m _ { i } + \frac { p _ { i } ^ { 2 } } { 2 m _ { i } } \right ) + \sum _ { i < j = 1 } ^ { 4 } V _ { i j } \,, \quad \ \ ( 1 ) \quad \stackrel { \text{exch} } { ( 1 - } } _ { \text{anti} }$$

where m i and p i are the mass and momentum of quark i .

In this study, we will employ the AL1 constituent quark potential model [65, 66] containing the one-gluon-exchange (OGE) interaction and a linear confinement term. The AL1 potential reads

$$V _ { i j } = - \frac { 3 } { 1 6 } & \lambda _ { i } ^ { c } \cdot \lambda _ { j } ^ { c } \Big ( - \frac { \kappa } { r _ { i j } } + \lambda r _ { i j } - \Lambda \\ & + \frac { 8 \pi \kappa ^ { \prime } } { 3 m _ { i } m _ { j } } \frac { \exp ( - r _ { i j } ^ { 2 } / r _ { 0 } ^ { 2 } ) } { \pi ^ { 3 / 2 } r _ { 0 } ^ { 3 } } s _ { i } \cdot s _ { j } \Big ), \quad ( 2 )$$

where r 0 = A Ä 2 m m i j m i + m j ä -B . λ c represents the SU(3) color Gell-Mann matrix and s i is the spin operator of quark i . The parameters of the potential have been determined by fitting the meson spectra in Ref. [66] and are listed in Table I. The calculated meson masses in AL1 model are shown in Table II.

TABLE I. The parameters in the AL1 quark potential model.

|      κ |   λ [GeV 2 ] |   Λ[GeV] |    κ ′ |   m s [GeV] |   A [GeV B - 1 ] |      B |
|--------|--------------|----------|--------|-------------|------------------|--------|
| 0.5069 |       0.1653 |   0.8321 | 1.8609 |       0.577 |           1.6553 | 0.2204 |

TABLE II. The theoretical masses (in MeV) and rms radii (in fm) of the ss ¯ mesons in the AL1 model, compared with the experimental results taken from Ref. [67].

| J PC        | Meson m Exp .   | m Theo .       | r rms Theo .   |
|-------------|-----------------|----------------|----------------|
| 0 - + η ′ a | -               | 713 . 5 0 . 54 |                |
| η ′ (2 S )  | -               | 1565 . 2       | 1 . 17         |
| η ′ (3 S )  | -               | 2140 . 9       | 1 . 65         |
| 1 -- ϕ      | 1019 . 5        | 1021 . 0       | 0 . 70         |
| ϕ (2 S )    | 1680            | 1695 . 1       | 1 . 25         |
| ϕ (3 S )    | 2188            | 2231 . 6       | 1 . 70         |

a For simplicity, we assume that there is no mixing effects between the I = 0 η nn ( ¯) and η ′ ( ss ¯) .

## B. Wave function construction

The wave function bases of the ssss ¯¯ tetraquark states can be expressed as the direct product of spatial wave function ϕ , spin wave function χ s , and color wave function χ c .

$$\psi = \mathcal { A } \left ( \phi \otimes \chi _ { s } \otimes \chi _ { c } \right ),$$

where A is the antisymmetrization operator, representing the exchange of identical quarks. For the ssss ¯¯ system, A = (1 -P 12 ) (1 -P 34 ) , where P ij permutes the i th and j th (anti)quarks.

For the spatial wave function, the Gaussian expansion method (GEM) [55] is employed. Namely, the spatial wave

function is expanded using the following basis:

$$\phi _ { n l m } ( r ) = \sqrt { \frac { 2 ^ { l + 5 / 2 } } { \Gamma \left ( l + \frac { 3 } { 2 } \right ) r _ { n } ^ { 3 } } } \left ( \frac { r } { r _ { n } } \right ) ^ { l } e ^ { - \frac { r ^ { 2 } } { r _ { n } ^ { 2 } } } Y _ { l m } ( \hat { r } ), \quad ( 4 )$$

where the r n is taken in geometric progression, r n = r a 0 n -1 . Y lm is the spherical harmonics.

For a four-body system, one can work on the Jacobi coordinates to exclude the center-of-mass motion. In principle, one could construct a complete basis functions using any type of Jacobi coordinates. However, to get rid of the complexity from the angular momentum, we use only S-wave bases constructed in different types of Jacobi coordinates to partially compensate for higher partial waves (see Ref. [61] for details). In this work, we include two types of Jacobi coordinates, the diquark-antidiquark structure and the dimeson structures, as shown in Fig. 1.

FIG. 1. Two structures for the tetraquark system. (a) diquarkantidiquark structure. (b) dimeson structure.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00503v1-figure-2f4a14bb-34d0-48bf-945a-1bc865205c1c.jpg)

In our calculation, we choose n max = 12 bases for each coodinate, and the settings of the basis parameters are

$$\begin{cases} \ r _ { 0 } = 0. 4 \, f m, r _ { n _ { \max } } = 2. 0 \, f m \, \ s - s \, \text{or} \, \bar { s } - \bar { s } & \text{ferent} \quad & \text{in the} \\ \ r _ { 0 } = 0. 4 \, f m, r _ { n _ { \max } } = 2. 0 \, f m \, \ ( s s ) - ( \bar { s } \bar { s } ) & \text{,} \quad ( 5 ) \\ \ r _ { 0 } = 0. 4 \, f m, r _ { n _ { \max } } = 1. 3 \, f m \, \ s - \bar { s } \\ \ r _ { 0 } = 0. 5 \, f m, r _ { n _ { \max } } = 4. 5 \, f m \, \ ( s \bar { s } ) - ( s \bar { s } ) & \text{Th} \quad & \text{to obt} \end{cases}$$

where we use relatively large r n max for the coordinate between two ss ¯ clusters to better describe the meson-meson scattering states and possible molecular states.

For the spin wave function, the bases for total spin S = 0 1 2 , , are

$$& \text{are} \\ & S = 0 \colon \begin{cases} \left [ ( s _ { 1 } s _ { 2 } ) _ { 0 } \left ( \bar { s } _ { 3 } \bar { s } _ { 4 } \right ) _ { 0 } \right ] _ { 0 } & \text{Unc} \\ \left [ ( s _ { 1 } s _ { 2 } ) _ { 1 } \left ( \bar { s } _ { 3 } \bar { s } _ { 4 } \right ) _ { 1 } \right ] _ { 0 } & \text{H} \\ & S = 1 \colon \begin{cases} \left [ ( s _ { 1 } s _ { 2 } ) _ { 0 } \left ( \bar { s } _ { 3 } \bar { s } _ { 4 } \right ) _ { 1 } \right ] _ { 1 } & \text{(6)} \\ \left [ ( s _ { 1 } s _ { 2 } ) _ { 1 } \left ( \bar { s } _ { 3 } \bar { s } _ { 4 } \right ) _ { 0 } \right ] _ { 1 }, & \text{Me} \\ \left [ ( s _ { 1 } s _ { 2 } ) _ { 1 } \left ( \bar { s } _ { 3 } \bar { s } _ { 4 } \right ) _ { 1 } \right ] _ { 1 } & \text{non} \\ & S = 2 \colon \end{cases} \text{$[(s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s=s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s+s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,-s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s.s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,.s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s:s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,s,...} \end{cases}$$

$$S = 2 \colon \ \left [ ( s _ { 1 } s _ { 2 } ) _ { 1 } \left ( \bar { s } _ { 3 } \bar { s } _ { 4 } \right ) _ { 1 } \right ] _ { 2 }.$$

Two optional color wave function bases sets are

$$\text{color-I: } \begin{cases} \left [ ( s _ { 1 } s _ { 2 } ) _ { \bar { 3 } } \left ( \bar { s } _ { 3 } \bar { s } _ { 4 } \right ) _ { 3 } \right ] _ { 1 } \\ \left [ ( s _ { 1 } s _ { 2 } ) _ { 6 } \left ( \bar { s } _ { 3 } \bar { s } _ { 4 } \right ) _ { \bar { 6 } } \right ] _ { 1 } \end{cases}, \quad \text{(7)}$$

$$\text{color-II: } \begin{cases} \left [ ( s _ { 1 } \bar { s } _ { 3 } ) _ { 1 } ( s _ { 2 } \bar { s } _ { 4 } ) _ { 1 } \right ] _ { 1 } \\ \left [ ( s _ { 1 } \bar { s } _ { 3 } ) _ { 8 } ( s _ { 2 } \bar { s } _ { 4 } ) _ { 8 } \right ] _ { 1 } \end{cases}. \quad \quad ( 8 )$$

Since both sets are orthogonal and complete, using either set will yield exactly the same results [61].

Besides the antisymmetrization operation in Eq. (3), since the ssss ¯¯ system carries definite C-parity, it is necessary to add (or subtract) the C-transformation of each term in Eq. (3) to ensure that the overall wave function has the desired positive (or negative) C-parity. The specific C-transformation is

$$& \mathfrak { P } \\ & \quad \left [ ( q _ { 1 } q _ { 2 } ^ { \prime } ) _ { \bar { c } } ^ { s _ { 1 } } \left ( \bar { q } _ { 3 } \bar { q } _ { 4 } ^ { \prime } \right ) _ { c } ^ { s _ { 2 } } \right ] _ { 1 _ { c } } ^ { S } \phi \left ( r _ { 1 }, r _ { 2 }, r _ { 3 }, r _ { 4 } \right ) \\ & \quad \stackrel { C } { \to } ( - 1 ) ^ { S - s _ { 1 } - s _ { 2 } } \left [ ( q _ { 1 } q _ { 2 } ^ { \prime } ) \right ) _ { \bar { c } } ^ { s _ { 2 } } \left ( \bar { q } _ { 3 } \bar { q } _ { 4 } ^ { \prime } \right ) _ { c } ^ { s _ { 1 } } \right ] _ { 1 _ { c } } ^ { S } \phi \left ( r _ { 3 }, r _ { 4 }, r _ { 1 }, r _ { 2 } \right )$$

for the bases with color-I set, and

$$\begin{array} { c } \left [ \left ( q _ { 1 } \bar { q } _ { 3 } \right ) _ { c } ^ { s _ { 1 } } \left ( q _ { 2 } ^ { \prime } \bar { q } _ { 4 } ^ { \prime } \right ) _ { c } ^ { s _ { 2 } } \right ] _ { 1 _ { c } } ^ { S } \phi \left ( r _ { 1 }, r _ { 2 }, r _ { 3 }, r _ { 4 } \right ) \\ \overset { C } { \longrightarrow } \left ( - 1 \right ) ^ { s _ { 1 } + s _ { 2 } } \left [ \left ( q _ { 1 } \bar { q } _ { 3 } \right ) _ { c } ^ { s _ { 1 } } \left ( q _ { 2 } ^ { \prime } \bar { q } _ { 4 } ^ { \prime } \right ) _ { c } ^ { s _ { 2 } } \right ] _ { 1 _ { c } } ^ { S } \phi \left ( r _ { 3 }, r _ { 4 }, r _ { 1 }, r _ { 2 } \right ) \\ \end{array}$$

for the bases with color-II set. The particle q and q ′ can be different quarks, although they both represent the same s quark in the ssss ¯¯ system.

## C. Complex scaling method

The complex scaling method (CSM) is a direct approach to obtain the energies and the decay widths of resonant states in a many-body system by performing an analytical continuation of the Schr¨dinger equation [56-58]. o This is achieved by carrying out a complex rotation on the coordinate r and momentum p , given by

$$U ( \theta ) r = r e ^ { i \theta }, \ \ U ( \theta ) p = p e ^ { - i \theta }. \quad \quad ( 1 0 )$$

Under the rotation, the Hamiltonian in Eq. (1) becomes

$$H ( \theta ) = \sum _ { i = 1 } ^ { 4 } \left ( m _ { i } + \frac { p _ { i } ^ { 2 } e ^ { - 2 i \theta } } { 2 m _ { i } } \right ) + \sum _ { i < j = 1 } ^ { 4 } V _ { i j } \left ( r _ { i j } e ^ { i \theta } \right ). \ ( 1 1 )$$

Meanwhile, for the resonant states with pole positions within the range of the rotated angle, their wave functions become normalizable by integration, thereby solvable through localized Gaussian bases in the same way as bound states. As a result, solving the complex-scaled Schr¨dinger equation will o

simultaneously yield the eigenenergies of bound states and resonant states within the rotated angle.

A typical pattern of the solved eigenenergies in the complex energy plane is: The bound states lie on the negative real axis of the energy plane. The continuum states align along beams originating from thresholds with Arg( E ) = -2 θ . The resonant states with mass M R and width Γ R are located at E R = M R -i Γ R / 2 , and only those within | Arg( E R ) | &lt; 2 θ can be solved. The positions of the bound and resonant states remain unchanged with the variation of the rotation angle. One can find more details in Refs. [62, 68, 69].

## D. Spatial structure

The root-mean-square (rms) radius is a good physical quantity for reflecting the spatial structure of the tetraquark states, and a widely used criterion for distinguishing between compact and molecular tetraquark states.

The conventional definition of the rms radius under CSM is

$$r _ { i j } ^ { r m s, C } \equiv \text{Re} \left [ \sqrt { \frac { \left ( \Psi ( \theta ) \left | r _ { i j } ^ { 2 } e ^ { 2 i \theta } \right | \Psi ( \theta ) \right ) } { \left ( \Psi ( \theta ) \left | \Psi ( \theta ) \right ) } } \right ], \quad ( 1 2 ) \text{ \quad initio}$$

where the Ψ( ) θ is the obtained complex wave function of the tetraquark state. The round bra-ket represents the so-called c-product [70] defined as

$$( \phi _ { n } \ | \ \phi _ { m } ) \equiv \int \phi _ { n } ( r ) \phi _ { m } ( r ) \mathrm d ^ { 3 } r, \quad \ \ ( 1 3 ) \quad \ \ T$$

without taking complex conjugate of the bra-state. This procedure ensures the function inside the integral is analytic, thereby the expectation value of the physical quantity remains stable as the rotation angle changes. The rms radius calculated from the c-product is generally not real; however, its real part can still reflect the internal quark clustering behavior if the resonant state is not too broad, as discussed in Ref. [71].

Nevertheless, for tetraquark states containing identical quarks, this definition may fail to identify the molecular nature of their spatial structure due to the antisymmetrization of identical quarks, as discussed in Ref. [62-64]. To eliminate this ambiguity, we have proposed a new approach to calculate the rms radii of the fully heavy QQQQ ¯ ¯ system in Ref. [63], which can also be applied to the ssss ¯¯ system. In this new approach, the total wave function is non-orthogonally decom- posed into

$$\mathfrak { s } \text{ and } \quad \text{posed into} \\ \text{com-} e \text{ real } \Psi ( \theta ) = \sum _ { \substack { s _ { 1 } \leq s _ { 2 } \\ \text{along} } } \left [ [ ( q _ { 1 } \bar { q _ { 3 } } ^ { \prime } ) _ { 1 _ { c } } ^ { s _ { 1 } } ( q _ { 2 } \bar { q _ { 4 } } ^ { \prime } ) _ { 1 _ { c } } ^ { s _ { 2 } } ] _ { 1 _ { c } } ^ { S } \phi ( r _ { 1 }, r _ { 2 }, r _ { 3 }, r _ { 4 } ; \theta ) \\ \text{$\mathfrak { s } $ real } \text{ along } \quad \text{$-} \left [ ( q _ { 2 } \bar { q _ { 3 } } ^ { \prime } ) _ { 1 _ { c } } ^ { s _ { 1 } } ( q _ { 1 } \bar { q _ { 4 } } ^ { \prime } ) _ { 1 _ { c } } ^ { s _ { 2 } } ] _ { 1 _ { c } } ^ { S } \phi ( r _ { 2 }, r _ { 1 }, r _ { 3 }, r _ { 4 } ; \theta ) \\ < 2 \theta \\ \text{states } \quad \text{$+} \left [ ( q _ { 2 } \bar { q _ { 4 } } ^ { \prime } ) _ { 1 _ { c } } ^ { s _ { 1 } } ( q _ { 1 } \bar { q _ { 3 } } ^ { \prime } ) _ { 1 _ { c } } ^ { s _ { 2 } } ] _ { 1 _ { c } } ^ { S } \phi ( r _ { 2 }, r _ { 1 }, r _ { 4 }, r _ { 3 } ; \theta ) \right ] \\ \text{angle.} \\ = \mathcal { A } \left [ \sum _ { \substack { s _ { 1 } \leq s _ { 2 } \\ \text{along} } } \left [ ( q _ { 1 } \bar { q _ { 3 } } ^ { \prime } ) _ { 1 _ { c } } ^ { s _ { 1 } } ( q _ { 2 } \bar { q _ { 4 } } ^ { \prime } ) _ { 1 _ { c } } ^ { s _ { 2 } } ] _ { 1 _ { c } } ^ { S } \phi ( r _ { 1 }, r _ { 2 }, r _ { 3 }, r _ { 4 } ; \theta ) \right ] \\ \equiv \mathcal { A } \Psi _ { 1 3, 2 4 } ( \theta ), \\ \text{where } s _ { 1 }, s _ { 2 } \text{ sum over all spin configurations with total spin} \\ \sim \text{$-} \text{$-} \text{$-} \text{$-} \text{$-} \text{$-} \text{$-} \text{$-} \text{$-} \text{$-} \text{$-} \text{$-} \text{$-} \text{$-} \text{$-} \text{$-} \text{$-} \text{$-} \text{$-} \text{$-} \text{$-} \text{$-} \text{$-} \text{$-} \text{$-} \text{$-} \text{$-} \text{$-} \text{$-} \text{$-} \text{$-} \text{$-} \text{$-} \text{$-} \end{$$

where s , s 1 2 sum over all spin configurations with total spin S . The particle q and q ′ can be different quarks, although they both represent s quark in the ssss ¯¯ system. We use the decomposed wave function Ψ 13 24 , ( θ ) to define the new rms radius

$$\text{a-} \\ \text{ite}, \quad \text{$r_{ij}$} \equiv \text{Re} \left [ \sqrt { \frac { \left ( \Psi _ { 1 3, 2 4 } ( \theta ) \left | r _ { i j } ^ { 2 } e ^ { 2 i \theta } \right | \Psi _ { 1 3, 2 4 } ( \theta ) \right ) } { \left ( \Psi _ { 1 3, 2 4 } ( \theta ) \left | \Psi _ { 1 3, 2 4 } ( \theta ) \right ) } } \right ]. \ ( 1 5 )$$

This newly defined rms radius can better describe the spatial structure of the molecular-type tetraquark states containing identical quarks. We will provide the results for both definitions of the rms radius in our calculation.

## III. NUMERICAL RESULTS

## A. States with 'normal' C-parity

The 'normal' C-parity ssss ¯¯ systems refer to those that can be composed of two S-wave ground quarkonia, namely, the systems with J PC = 0 ++ , 1 + -, 2 ++ . The complex eigenenergies yielded are shown in the upper panel of Fig. 2. We choose different complex scaling angles θ = 9 ◦ , 12 ◦ , 15 ◦ to distinguish bound states and resonant states, which remain stationary as the angle changes, from the continuum states. The markers that align along beams originating from mesonmeson thresholds with Arg( E ) = -2 θ are the continuum state eigenenergies. The results show that no bound state is found below any of the thresholds, while a series of resonant states are obtained and marked out by black circles. Their complex energies, proportions of different color configurations, and two kinds of rms radii are summarized in Table III. For convenience, we label the tetraquark resonant states obtained in our calculations as T 4 s,J PC ( M ) , where M is the mass of the state.

The masses of the resonant states are in the region (2 7 . , 3 3) . GeV, and their widths range from a few MeV to several tens of

FIG. 2. The complex energy eigenvalues of the ssss ¯¯ states in the AL1 potential with varying θ in the CSM. The solid lines represent the continuum lines rotating along Arg( E ) = -2 θ . The resonant states do not shift as θ changes and are marked out by the black circles.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00503v1-figure-b061610a-4e67-439d-b2f2-4b5503118a53.jpg)

MeV. For the resonant states above the 1 S S 2 di-meson threshold, our calculation actually underestimates their widths because the decay width of the 2 S meson has not been accounted for. The rms radii of all the obtained states approximately match the sizes of the 2S mesons. Furthermore, the newly defined r rms M , ij does not exhibit any molecular characteristics, indicating that all these states are compact tetraquark states.

In tetraquark systems, there are abundant dimeson thresholds. Each of these thresholds represents a branch cut connecting two different Riemann sheets. In our calculations, we must choose appropriate rotation angle to make the pole of interest emerge. Meanwhile, due to numerical reasons, the pole close to the continuum solutions will be affected by them, as illustrated in Refs. [64, 71]. Therefore, to reduce uncertainties, the rotation angle should be chosen such that the continuum states are not too close to the resonant state. In practices, we adjust the angle to make each pole and its rms radii stable.

Comparing the fully strange systems with cccc ¯¯ using the AL1 model in our previous work [63], we observe a similar overall pattern. However, the fully strange systems exhibit an additional quasi-bound state with 2 ++ below the ϕϕ (2 S ) threshold. This difference may stem from the mass-dependent color-magnetic term in the potential.

For the 0 ++ states, the five resonant states can all decay into η η ′ ′ , ϕϕ , η η ′ ′ (2 S ) , ϕϕ (2 S ) final states, although there may be some effects if the mixing between the I = 0 η nn ( ¯) and η ′ ( ss ¯) is taken into account. The T 4 s, 0 ++ (2917) and T 4 s, 0 ++ (3133) have an additional decay channel into η η ′ ′ (3 S ) . The higher T 4 s, 0 ++ (3175) and T 4 s, 0 ++ (3248) have two additional channel into η η ′ ′ (3 S ) and η ′ (2 S η ) ′ (2 S ) . For the 1 + -states, the five resonant states can all decay into η ϕ ′ , η ϕ ′ (2 S ) and η ′ (2 S ϕ ) . The T 4 s, 1 + -(3142) has an additional decay channel into η ϕ ′ (3 S ) . And the higher T 4 s, 1 + -(3228) and T 4 s, 1 + -(3237) have two extra η ϕ ′ (3 S ) and η ′ (3 S ϕ ) decay channels. For the 2 ++ states, the lowest T 4 s, 2 ++ (2714) below the ϕϕ (2 S ) threshold have only one decay channel ϕϕ . Other two 2 ++ resonant states T 4 s, 2 ++ (2993) and T 4 s, 2 ++ (3164) have an additional decay channel into ϕϕ (2 S ) . The T 4 s, 2 ++ (3266) has an extra decay channel to ϕϕ (3 S ) .

Among these states, the state T 4 s, 2 ++ (2714) is the most likely one to be observed experimentally. Its mass is the lowest, and its decay channel ϕϕ predominantly decays into K K K K + -+ -, which are all charged particles that can be detected directly by detectors, leading to high detection efficiency and resolution.

## B. States with 'exotic' C-parity

The 'exotic' C-parity systems refer to those that cannot be composed of two S-wave ground quarkonia, namely, the system with J PC = 0 + -, 1 ++ , 2 + -. The complex eigenenergies yielded are shown in the lower panel of Fig. 2. We obtain a series of resonant states and zero-width states, whose energies, proportions of different color configurations, and two kinds of rms radii are summarized in Table IV. The rms radii approximately match the sizes of the 2S mesons. The newly defined r rms M , ij does not exhibit any molecular characteristics, indicating that all these states are compact tetraquark states.

For the J PC = 0 + -and 2 + -systems, there is no S-wave

TABLE III. The complex energies E = M -i Γ 2 / (in MeV), proportions of different color configurations and rms radii (in fm) of the ssss ¯¯ resonant states with 'normal' C-parity in the AL1 potential. The r rms C , ij and r rms M , ij denote the conventional rms radius and the newly defined rms radius based on the decomposed wave function, respectively. The last column shows the spatial structure of the states, where C. and M. represent the compact and molecular structure, respectively.

| J PC   | M - i Γ / 2   | χ 3 c ⊗ 3 c   | χ 6 c ⊗ 6 c   | χ 1 c ⊗ 1 c   | χ 8 c ⊗ 8 c   | r ij     | r ij       | rms , C         | rms , C         | rms , C                           | rms , C                           | structure   |
|--------|---------------|---------------|---------------|---------------|---------------|----------|------------|-----------------|-----------------|-----------------------------------|-----------------------------------|-------------|
| J PC   | M - i Γ / 2   | χ 3 c ⊗ 3 c   | χ 6 c ⊗ 6 c   | χ 1 c ⊗ 1 c   | χ 8 c ⊗ 8 c   | r rms ss | r rms ss ¯ | r rms s 1 ¯ s 3 | r rms s 2 ¯ s 4 | r rms s 1 s 2 = r rms ¯ s 3 ¯ s 4 | r rms s 1 ¯ s 4 = r rms s 2 ¯ s 3 | structure   |
| 0 ++   | 2852 - 40 i   | 86%           | 14%           | 38%           | 62%           | 0.95     | 1.16       | 1.20            | 1.20            | 0.91                              | 1.19                              | C.          |
| 0 ++   | 2917 - 9 i    | 40%           | 60%           | 53%           | 47%           | 1.23     | 1.21       | 1.12            | 1.12            | 1.18                              | 1.35                              | C.          |
| 0 ++   | 3133 - 7 i    | 58%           | 42%           | 47%           | 53%           | 1.51     | 1.44       | 1.27            | 1.27            | 1.48                              | 1.66                              | C.          |
| 0 ++   | 3175 - 4 i    | 46%           | 54%           | 51%           | 49%           | 1.30     | 1.27       | 1.14            | 1.14            | 1.28                              | 1.39                              | C.          |
| 0 ++   | 3248 - 10 i   | 35%           | 65%           | 55%           | 45%           | 1.37     | 1.36       | 1.31            | 1.31            | 1.32                              | 1.49                              | C.          |
| 1 + -  | 2819 - 3 i    | 63%           | 37%           | 46%           | 54%           | 1.01     | 1.11       | 1.04            | 1.05            | 1.00                              | 1.18                              | C.          |
| 1 + -  | 2940 - 46 i   | 87%           | 13%           | 38%           | 62%           | 1.02     | 1.12       | 1.18            | 1.13            | 1.03                              | 1.15                              | C.          |
| 1 + -  | 3142 - 12 i   | 77%           | 23%           | 41%           | 59%           | 1.15     | 1.43       | 1.51            | 1.28            | 1.10                              | 1.49                              | C.          |
| 1 + -  | 3228 - 2 i    | 66%           | 34%           | 45%           | 55%           | 1.22     | 1.37       | 1.27            | 1.29            | 1.22                              | 1.47                              | C.          |
| 1 + -  | 3237 - 4 i    | 64%           | 36%           | 45%           | 55%           | 1.22     | 1.33       | 1.18            | 1.22            | 1.22                              | 1.44                              | C.          |
| 2 ++   | 2714 - 6 i    | 75%           | 25%           | 42%           | 58%           | 1.06     | 1.09       | 1.11            | 1.11            | 0.98                              | 1.14                              | C.          |
| 2 ++   | 2993 - 48 i   | 85%           | 15%           | 38%           | 62%           | 1.00     | 1.03       | 1.14            | 1.14            | 1.03                              | 1.03                              | C.          |
| 2 ++   | 3164 - 3 i    | 92%           | 8%            | 36%           | 64%           | 0.94     | 1.47       | 1.47            | 1.47            | 0.93                              | 1.50                              | C.          |
| 2 ++   | 3266 - 2 i    | 66%           | 34%           | 45%           | 55%           | 1.29     | 1.38       | 1.22            | 1.22            | 1.28                              | 1.50                              | C.          |

di-quarkonium threshold. So the calculated states appear as bound states on the complex energy plane. However, there may actually be higher-wave thresholds below these states, which are not included in the calculation. If the higher-wave mesons are included, the properties of these states may change and their widths may no longer be zero.

parity J PC = 0 ++ , 1 + -, 2 ++ and 'exoitic' C-parity J PC = 0 + -, 1 ++ , 2 + -using the AL1 quark potential model. We employ the GEM to solve the four-body Schr¨ odinger equation, and the CSM to identify genuine resonant states from meson-meson scattering states.

For the J PC = 1 ++ system, the masses of the resonant states range in (2 7 . , 3 3) . GeV, and their widths are small, with the largest being only a few MeV. Generally, one would expect the widths of these resonant states to arise from their decays into lower-energy meson-meson channels. However, their narrow widths suggest that such decays are suppressed. The reason for the suppression has been discussed in detail in Ref. [63]. Simply put, when the rms radii between different i, j quarks are roughly equal, the transition operator of the color-electric term becomes approximately proportional to the identity operator. Therefore, the color mixing transition is suppressed. In the ssss ¯¯ system, however, the contribution from the color-magnetic term is not as significantly suppressed by the quark mass as it is in the heavy quark sector. Therefore, the transition can still occur via spin mixing, which, although limited, can lead to slightly larger widths of these resonant states. All three J PC = 1 ++ resonant states can decay via ϕϕ (2 S ) , though this process is suppressed.

## IV. SUMMARY AND DISCUSSION

We calculate the mass spectrum of the S-wave fully strange tetraquark systems with 'normal' C-

We obtain a series of resonant and zero-width states in the mass region (2 7 . , 3 3) , GeV, with widths ranging from less than 1 MeV to around 50 MeV. Apart from the state T 4 s, 2 ++ (2714) , all the others lie above the M (1 S M ) ′ (2 S ) dimeson thresholds.

We calculate two kinds of rms radius for all states, the conventional rms radius r rms C , ij and the newly defined rms radius r rms M , ij based on the decomposed wave function. For all states, both kinds of rms radii approximately match the sizes of the 2S mesons, and the newly defined r rms M , ij does not exhibit any molecular characteristics, indicating that all these states are compact tetraquark states.

Among these states, the state T 4 s, 2 ++ (2714) is the most likely to be observed experimentally due to its lowest mass and its decay into ϕϕ , which predominantly produces directly detectable charged K K + -pairs, ensuring high detection efficiency and resolution.

Recalling the experimental candidates for ssss ¯¯ tetraquark states, their masses are all in range of 2.1 to 2.5 GeV , whereas the lowest S-wave resonant state in our calculations is around 2.7 GeV. This suggests the states like ϕ (2170) and X (2370) are unlikely to be compact ssss ¯¯ tetraquark states. Since the lowest S-wave state is already as high as 2.7 GeV, it's unlikely that a P-wave compact state would be lower than the S-wave ones. However, they could still be P-wave molecular-

TABLE IV. The complex energies E = M -i Γ 2 / (in MeV), proportions of different color configurations and rms radii (in fm) of the ssss ¯¯ resonant states with 'exotic' C-parity in the AL1 potential. The r rms C , ij and r rms M , ij denote the conventional rms radius and the newly defined rms radius based on the decomposed wave function, respectively. The last column shows the spatial structure of the states, where C. and M. represent the compact and molecular structure, respectively.

| J PC   | M - i Γ / 2    | χ 3 c ⊗ 3 c   | χ 6 c ⊗ 6 c   | χ 1 c ⊗ 1 c   | χ 8 c ⊗ 8 c   | r rms ij   | r rms ij   | r rms , M ij    | r rms , M ij    | r rms , M ij                      | r rms , M ij                      | structure   |
|--------|----------------|---------------|---------------|---------------|---------------|------------|------------|-----------------|-----------------|-----------------------------------|-----------------------------------|-------------|
| J PC   | M - i Γ / 2    | χ 3 c ⊗ 3 c   | χ 6 c ⊗ 6 c   | χ 1 c ⊗ 1 c   | χ 8 c ⊗ 8 c   | r rms ss   | r rms ss ¯ | r rms s 1 ¯ s 3 | r rms s 2 ¯ s 4 | r rms s 1 s 2 = r rms ¯ s 3 ¯ s 4 | r rms s 1 ¯ s 4 = r rms s 2 ¯ s 3 | structure   |
| 0 + -  | 2725           | 33%           | 67%           | 56%           | 44%           | 1.13       | 0.96       | 0.96            | 0.96            | 1.10                              | 0.96                              | C.          |
|        | 2873           | 65%           | 35%           | 45%           | 55%           | 1.17       | 1.03       | 1.02            | 1.02            | 1.13                              | 1.02                              | C.          |
|        | 3148           | 21%           | 79%           | 60%           | 40%           | 1.44       | 1.20       | 1.20            | 1.20            | 1.41                              | 1.20                              | C.          |
|        | 3285           | 78%           | 22%           | 41%           | 59%           | 1.30       | 1.28       | 1.28            | 1.28            | 1.24                              | 1.28                              | C.          |
| 1 ++   | 2723 - 0 . 5 i | 59%           | 41%           | 47%           | 53%           | 1.14       | 1.09       | 0.90            | 1.03            | 0.99                              | 1.13                              | C.          |
|        | 2863 - 4 i     | 99%           | 1%            | 34%           | 66%           | 1.07       | 1.01       | 1.05            | 0.95            | 1.07                              | 1.01                              | C.          |
|        | 3151 - 0 . 1 i | 66%           | 34%           | 45%           | 55%           | 1.17       | 1.31       | 1.14            | 1.23            | 1.18                              | 1.42                              | C.          |
| 2 + -  | 2896           | 100%          | 0%            | 33%           | 67%           | 1.10       | 1.02       | 1.02            | 1.02            | 1.10                              | 1.02                              | C.          |

type resonant states. In this work, a state T 4 s, 2 ++ (2714) is found near the threshold of ground ϕ combined with a radially excited ϕ (2 S ) , although it does not exhibit a molecular spatial structure. Introducing P-wave thresholds might yield molecular-type resonant states near the threshold of an S-wave meson combined with a P-wave meson, which could be lower in mass. This possibility awaits careful investigation in the future.

We emphasize that the fully charmed and fully strange tetraquark resonant states share many common characteristics. Their internal structures and underlying dynamics are nearly the same. All of them are compact states. The discovery of X (6900) implies the existence of its fully strange cousins, which can be searched for in the ϕϕ or ϕϕ (2 S ) channels. While X (6900) can only be produced at LHC now, the fully strange tetraquark resonances may be observed at BE-

- [1] H.-X. Chen, W. Chen, X. Liu, and S.-L. Zhu, The hidden-charm pentaquark and tetraquark states, Phys. Rept. 639 , 1 (2016), arXiv:1601.02092 [hep-ph].
- [2] A. Hosaka, T. Iijima, K. Miyabayashi, Y. Sakai, and S. Yasui, Exotic hadrons with heavy flavors: X, Y, Z, and related states, PTEP 2016 , 062C01 (2016), arXiv:1603.09229 [hep-ph].
- [3] R. F. Lebed, R. E. Mitchell, and E. S. Swanson, HeavyQuark QCD Exotica, Prog. Part. Nucl. Phys. 93 , 143 (2017), arXiv:1610.04528 [hep-ph].
- [4] A. Esposito, A. Pilloni, and A. D. Polosa, Multiquark Resonances, Phys. Rept. 668 , 1 (2017), arXiv:1611.07920 [hep-ph].
- [5] A. Ali, J. S. Lange, and S. Stone, Exotics: Heavy Pentaquarks and Tetraquarks, Prog. Part. Nucl. Phys. 97 , 123 (2017), arXiv:1706.00610 [hep-ph].
- [6] F.-K. Guo, C. Hanhart, U.-G. Meißner, Q. Wang, Q. Zhao, and B.-S. Zou, Hadronic molecules, Rev. Mod. Phys. 90 , 015004 (2018), [Erratum: Rev.Mod.Phys. 94, 029901 (2022)], arXiv:1705.00141 [hep-ph].
- [7] N. Brambilla, S. Eidelman, C. Hanhart, A. Nefediev, C.-P. Shen, C. E. Thomas, A. Vairo, and C.-Z. Yuan, The XYZ

SIII, Belle II and LHC. The future discovery of these states shall shed light on both the charmed and strange compact tetraquark states.

## ACKNOWLEDGMENTS

We are grateful to Zi-Yang Lin, Jun-Zhang Wang and Liang-Zhen Wen for the helpful discussions. This project was supported by the National Natural Science Foundation of China (11975033 and 12070131001). This project was also funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation, Project ID 196253076-TRR 110). The computational resources were supported by Highperformance Computing Platform of Peking University.

- states: experimental and theoretical status and perspectives, Phys. Rept. 873 , 1 (2020), arXiv:1907.07583 [hep-ex].
- [8] Y.-R. Liu, H.-X. Chen, W. Chen, X. Liu, and S.-L. Zhu, Pentaquark and Tetraquark states, Prog. Part. Nucl. Phys. 107 , 237 (2019), arXiv:1903.11976 [hep-ph].
- [9] H.-X. Chen, W. Chen, X. Liu, Y.-R. Liu, and S.-L. Zhu, An updated review of the new hadron states, Rept. Prog. Phys. 86 , 026201 (2023), arXiv:2204.02649 [hep-ph].
- [10] L. Meng, B. Wang, G.-J. Wang, and S.-L. Zhu, Chiral perturbation theory for heavy hadrons and chiral effective field theory for heavy hadronic molecules, Phys. Rept. 1019 , 1 (2023), arXiv:2204.08716 [hep-ph].
- [11] M. Mai, U.-G. Meißner, and C. Urbach, Towards a theory of hadron resonances, Phys. Rept. 1001 , 1 (2023), arXiv:2206.01477 [hep-ph].
- [12] B. Aubert et al. (BaBar), A Structure at 2175-MeV in e + e -→ ϕ f0(980) Observed via Initial-State Radiation, Phys. Rev. D 74 , 091103 (2006), arXiv:hep-ex/0610018.
- [13] B. Aubert et al. (BaBar), The e+ e--&gt; K+ K- pi+ pi-, K+ K- pi0 pi0 and K+ K- K+ K- cross-sections measured

- with initial-state radiation, Phys. Rev. D 76 , 012008 (2007), arXiv:0704.0630 [hep-ex].
- [14] B. Aubert et al. (BaBar), Measurements of e + e -→ K K η K + -, + K π -0 and K K 0 s ± π ∓ cross- sections using initial state radiation events, Phys. Rev. D 77 , 092002 (2008), arXiv:0710.4451 [hep-ex].
- [15] J. P. Lees et al. (BaBar), Cross Sections for the Reactions e+e--&gt; K+ K- pi+pi-, K+ K- pi0pi0, and K+ K- K+ K- Measured Using Initial-State Radiation Events, Phys. Rev. D 86 , 012008 (2012), arXiv:1103.3001 [hep-ex].
- [16] M. Ablikim et al. (BES), Observation of Y(2175) in J / psi -&gt; eta phi f(0)(980), Phys. Rev. Lett. 100 , 102003 (2008), arXiv:0712.1143 [hep-ex].
- [17] M. Ablikim et al. (BESIII), Study of J/ψ → ηϕπ + π -at BESIII, Phys. Rev. D 91 , 052017 (2015), arXiv:1412.5258 [hepex].
- [18] M. Ablikim et al. (BESIII), Observation of e + e -→ ηY (2175) at center-of-mass energies above 3.7 GeV, Phys. Rev. D 99 , 012014 (2019), arXiv:1709.04323 [hep-ex].
- [19] C. P. Shen et al. (Belle), Observation of the phi(1680) and the Y(2175) in e+e- -&gt; phi pi+ pi-, Phys. Rev. D 80 , 031101 (2009), arXiv:0808.0006 [hep-ex].
- [20] M. Ablikim et al. (BESIII), Confirmation of the X (1835) and observation of the resonances X (2120) and X (2370) in J/ψ → γπ + π - ′ η , Phys. Rev. Lett. 106 , 072002 (2011), arXiv:1012.3510 [hep-ex].
- [21] M. Ablikim et al. (BESIII), Observation of X (2370) and search for X(2120) in J/ψ → γKKη ¯ ′ , Eur. Phys. J. C 80 , 746 (2020), arXiv:1912.11253 [hep-ex].
- [22] M. Ablikim et al. (BESIII), Observation of pseudoscalar and tensor resonances in J/ψ → γϕϕ , Phys. Rev. D 93 , 112011 (2016), arXiv:1602.01523 [hep-ex].
- [23] M. Ablikim et al. (BESIII), Measurement of e + e -→ K K + -cross section at √ s = 2 00 . -3 08 . GeV, Phys. Rev. D 99 , 032001 (2019), arXiv:1811.08742 [hep-ex].
- [24] M. Ablikim et al. (BESIII), Observation and study of the decay J/ψ → ϕηη ′ , Phys. Rev. D 99 , 112008 (2019), arXiv:1901.00085 [hep-ex].
- [25] C. P. Shen and C. Z. Yuan, Combined fit to BaBar and Belle Data on e+ e- -&gt; phi pi+ pi- and phi f0(980), Chin. Phys. C 34 , 1045 (2010), arXiv:0911.1591 [hep-ex].
- [26] M. Ablikim et al. (BESIII), Measurement of e + e -→ ϕπ + π -cross sections at center-of-mass energies from 2.00 to 3.08 GeV, (2021), arXiv:2112.13219 [hep-ex].
- [27] R. Aaij et al. (LHCb), Observation of structure in the J/ψ -pair mass spectrum, Sci. Bull. 65 , 1983 (2020), arXiv:2006.16957 [hep-ex].
- [28] A. Hayrapetyan et al. (CMS), New Structures in the J/ ψ J/ ψ Mass Spectrum in Proton-Proton Collisions at s=13 TeV, Phys. Rev. Lett. 132 , 111901 (2024), arXiv:2306.07164 [hep-ex].
- [29] G. Aad et al. (ATLAS), Observation of an Excess of Dicharmonium Events in the Four-Muon Final State with the ATLAS Detector, Phys. Rev. Lett. 131 , 151902 (2023), arXiv:2304.08962 [hep-ex].
- [30] Z.-G. Wang, Analysis of the Y(2175) as a tetraquark state with QCD sum rules, Nucl. Phys. A 791 , 106 (2007), arXiv:hepph/0610171.
- [31] H.-X. Chen, X. Liu, A. Hosaka, and S.-L. Zhu, The Y(2175) State in the QCD Sum Rule, Phys. Rev. D 78 , 034012 (2008), arXiv:0801.4603 [hep-ph].
- [32] H.-X. Chen, C.-P. Shen, and S.-L. Zhu, A possible partner state of the Y (2175) , Phys. Rev. D 98 , 014011 (2018), arXiv:1805.06100 [hep-ph].
- [33] Y.-W. Jiang, W.-H. Tan, H.-X. Chen, and E.-L. Cui, Strong de-
- cays of the ϕ (2170) as a fully-strange tetraquark state, (2023), arXiv:2310.16612 [hep-ph].
- [34] C. Deng, J. Ping, F. Wang, and T. Goldman, Tetraquark state and multibody interaction, Phys. Rev. D 82 , 074001 (2010).
- [35] N. V. Drenska, R. Faccini, and A. D. Polosa, Higher Tetraquark Particles, Phys. Lett. B 669 , 160 (2008), arXiv:0807.0593 [hepph].
- [36] G.-J. Ding and M.-L. Yan, A Candidate for 1- strangeonium hybrid, Phys. Lett. B 650 , 390 (2007), arXiv:hep-ph/0611319.
- [37] G.-J. Ding and M.-L. Yan, Y(2175): Distinguish Hybrid State from Higher Quarkonium, Phys. Lett. B 657 , 49 (2007), arXiv:hep-ph/0701047.
- [38] A. Martinez Torres, K. P. Khemchandani, L. S. Geng, M. Napsuciale, and E. Oset, The X(2175) as a resonant state of the phi K anti-K system, Phys. Rev. D 78 , 074031 (2008), arXiv:0801.3635 [nucl-th].
- [39] L. Alvarez-Ruso, J. A. Oller, and J. M. Alarcon, On the phi(1020) f0(980) S-wave scattering and the Y(2175) resonance, Phys. Rev. D 80 , 054011 (2009), arXiv:0906.0222 [hepph].
- [40] Y. Dong, A. Faessler, T. Gutsche, Q. L¨, and V. E. Lyubovitskij, u Selected strong decays of η (2225) and ϕ (2170) as Λ ¯ Λ bound states, Phys. Rev. D 96 , 074027 (2017), arXiv:1705.09631 [hep-ph].
- [41] J. Ho, R. Berg, T. G. Steele, W. Chen, and D. Harnett, Is the Y (2175) a Strangeonium Hybrid Meson?, Phys. Rev. D 100 , 034012 (2019), arXiv:1905.12779 [hep-ph].
- [42] S. S. Agaev, K. Azizi, and H. Sundu, Nature of the vector resonance Y (2175) , Phys. Rev. D 101 , 074012 (2020), arXiv:1911.09743 [hep-ph].
- [43] R.-R. Dong, N. Su, H.-X. Chen, E.-L. Cui, and Z.-Y. Zhou, QCD sum rule studies on the ssss ¯¯ tetraquark states of J PC = 0 -+ , Eur. Phys. J. C 80 , 749 (2020), arXiv:2003.07670 [hepph].
- [44] N. Su and H.-X. Chen, S- and P-wave fully strange tetraquark states from QCD sum rules, Phys. Rev. D 106 , 014023 (2022), arXiv:2204.13959 [hep-ph].
- [45] F.-X. Liu, M.-S. Liu, X.-H. Zhong, and Q. Zhao, Fully-strange tetraquark ssss ¯¯ spectrum and possible experimental evidence, Phys. Rev. D 103 , 016016 (2021), arXiv:2008.01372 [hep-ph].
- [46] Q.-F. L¨ u, K.-L. Wang, and Y.-B. Dong, The ssss ¯¯ tetraquark states and the newly observed structure X (2239) by BESIII Collaboration, Chin. Phys. C 44 , 024101 (2020), arXiv:1903.05007 [hep-ph].
- [47] E.-L. Cui, H.-M. Yang, H.-X. Chen, W. Chen, and C.-P. Shen, QCD sum rule studies of ssss ¯¯ tetraquark states with J PC = 1 + -, Eur. Phys. J. C 79 , 232 (2019), arXiv:1901.01724 [hepph].
- [48] Z.-G. Wang, Light tetraquark state candidates, Adv. High Energy Phys. 2020 , 6438730 (2020), arXiv:1901.04815 [hep-ph].
- [49] K. Azizi, S. S. Agaev, and H. Sundu, Light axial-vector and vector resonances X (2100) and X (2239) , Nucl. Phys. B 948 , 114789 (2019), arXiv:1906.04061 [hep-ph].
- [50] D. Ebert, R. N. Faustov, and V. O. Galkin, Masses of light tetraquarks and scalar mesons in the relativistic quark model, Eur. Phys. J. C 60 , 273 (2009), arXiv:0812.2116 [hep-ph].
- [51] R.-R. Dong, N. Su, and H.-X. Chen, Highly excited and exotic fully-strange tetraquark states, Eur. Phys. J. C 82 , 983 (2022), arXiv:2206.09517 [hep-ph].
- [52] Y. Tan, Y. Wu, H. Huang, and J. Ping, Investigating Fully Strange Tetraquark System with Positive Parity in a Chiral Quark Model, Universe 10 , 17 (2024).
- [53] H.-Z. Xi, Y.-W. Jiang, H.-X. Chen, A. Hosaka, and N. Su, Fully-strange tetraquark states with the exotic quantum num-

- bers J PC = 0 + -and 2 + -, Phys. Rev. D 108 , 094019 (2023), arXiv:2307.07819 [hep-ph].
- [54] Q.-N. Wang, D.-K. Lian, and W. Chen, Light tetraquark states with exotic quantum numbers J PC = 2 + -, (2024), arXiv:2404.18547 [hep-ph].
- [55] E. Hiyama, Y. Kino, and M. Kamimura, Gaussian expansion method for few-body systems, Prog. Part. Nucl. Phys. 51 , 223 (2003).
- [56] J. Aguilar and J. M. Combes, A class of analytic perturbations for one-body schroedinger hamiltonians, Commun. Math. Phys. 22 , 269 (1971).
- [57] E. Balslev and J. M. Combes, Spectral properties of manybody schroedinger operators with dilatation-analytic interactions, Commun. Math. Phys. 22 , 280 (1971).
- [58] S. Aoyama, T. Myo, K. Kat¯ o, and K. Ikeda, The complex scaling method for many-body resonances and its applications to three-body resonances, Progress of theoretical physics 116 , 1 (2006).
- [59] Y. Ma, L. Meng, Y.-K. Chen, and S.-L. Zhu, Ground state baryons in the flux-tube three-body confinement model using diffusion Monte Carlo, Phys. Rev. D 107 , 054035 (2023), arXiv:2211.09021 [hep-ph].
- [60] Y. Ma, L. Meng, Y.-K. Chen, and S.-L. Zhu, Doubly heavy tetraquark states in the constituent quark model using diffusion Monte Carlo method, Phys. Rev. D 109 , 074001 (2024), arXiv:2309.17068 [hep-ph].
- [61] L. Meng, Y.-K. Chen, Y. Ma, and S.-L. Zhu, Tetraquark bound states in constituent quark models: Benchmark test calculations, Phys. Rev. D 108 , 114016 (2023), arXiv:2310.13354 [hep-ph].
- [62] Y.-K. Chen, W.-L. Wu, L. Meng, and S.-L. Zhu, Unified description of the Qsq ¯ ¯ q molecular bound states, molecular resonances, and compact tetraquark states in the quark potential model, Phys. Rev. D 109 , 014010 (2024), arXiv:2310.14597 [hep-ph].
- [63] W.-L. Wu, Y.-K. Chen, L. Meng, and S.-L. Zhu, Benchmark calculations of fully heavy compact and molecular tetraquark states, Phys. Rev. D 109 , 054034 (2024), arXiv:2401.14899 [hep-ph].
- [64] W.-L. Wu, Y. Ma, Y.-K. Chen, L. Meng, and S.-L. Zhu, Fully heavy tetraquark resonant states with different flavors, (2024), arXiv:2406.17824 [hep-ph].
- [65] C. Semay and B. Silvestre-Brac, Diquonia and potential models, Z. Phys. C 61 , 271 (1994).
- [66] B. Silvestre-Brac, Spectrum and static properties of heavy baryons, Few Body Syst. 20 , 1 (1996).
- [67] R. L. Workman et al. (Particle Data Group), Review of Particle Physics, PTEP 2022 , 083C01 (2022).
- [68] Z.-Y. Lin, J.-B. Cheng, and S.-L. Zhu, T + cc and X (3872) with the complex scaling method and DD D π ( ¯ ) three-body effect, (2022), arXiv:2205.14628 [hep-ph].
- [69] Y.-K. Chen, L. Meng, Z.-Y. Lin, and S.-L. Zhu, Virtual states in the coupled-channel problems with an improved complex scaling method, Phys. Rev. D 109 , 034006 (2024), arXiv:2308.12424 [hep-ph].
- [70] W. J. Romo, Inner product for resonant states and shell-model applications, Nucl. Phys. A 116 , 617 (1968).
- [71] M. Homma, T. Myo, and K. Kat¯ o, Matrix elements of physical quantities associated with resonance states, Progress of theoretical physics 97 , 561 (1997).
Yao Ma Wei-Lin Wu Lu Meng Yan-Ke Chen Shi-Lin Zhu
Fully strange tetraquark resonant states as the cousins of $X(6900)$
We conduct systematic calculations of the S-wave fully strange systems with
``normal" $\left(J^{P C}=0^{++}, 1^{+-}, 2^{++}\right)$ and ``exotic"
$\left(J^{P C}=0^{+-}, 1^{++}, 2^{+-}\right)$ C-parities, which are the strange
analogue of the fully charmed tetraquark state $X(6900)$. Within a constituent
quark potential model, we employ the Gaussian expansion method to solve the
four-body Schr\"odinger equation and the complex scaling method to identify
resonant states. We obtain a series of resonant states and zero-width states in
the mass range of 2.7 to 3.3 GeV, with their widths ranging from less than 1
MeV to about 50 MeV. Their rms radii strongly indicate that they are compact
tetraquark states. Among these states, the $T_{4s,2^{++}}(2714)$ may be the
most likely one to be observed experimentally. We urge the experimental
exploration of the $2^{++}$ $s s \bar{s} \bar{s}$ state around 2.7 GeV in the
$\phi\phi$ channel. Since the lowest S-wave $s s \bar{s} \bar{s}$ state is
around 2.7 GeV, the compact P-wave $s s \bar{s} \bar{s}$ states are expected to
be heavier. Hence, $\phi(2170)$ and $X(2370)$ are unlikely to be compact
tetraquark states.