Masses of fully heavy tetraquarks $QQ \bar Q \bar Q$ in an extended relativized quark model

Inspired by recent measurement of possible fully charmed tetraquarks in LHCb Collaboration, we investigate the mass spectra of fully heavy tetraquarks $QQ \bar Q \bar Q$ in an extended relativized quark model. Our estimations indicate that the broad structure around 6.4 GeV should contain one or more ground $cc \bar c \bar c$ tetraquark states, while the narrow structure near 6.9 GeV can be categorized as the first radial excitation of $cc \bar c \bar c$ system. Moreover, with the wave functions of the tetraquarks and mesons, the strong decays of tetraquarks into heavy quarkonium pair are qualitatively discussed, which can be further checked by the LHCb and CMS Collaborations.


I. INTRODUCTION
Since the observation of X(3872) in 2003 [1], the searching for hadrons beyond the conventional mesons and baryons have become one of intriguing topics in the past decades. On the experimental side, a growing number of new hadron states have been observed experimentally. Some of these states cannot be accommodated into the traditional mesons or baryons, which can be good candidates of molecular or tetraquark states. Recent experimental and theoretical status can be found in the literature reviews [2][3][4][5][6][7][8][9][10][11].
Among the observed new hadron states, those with heavy quark components are particularly interesting, since the spectroscopy of traditional mesons and baryons with heavy quarks are much clear than the light hadrons. Moreover, the interactions involved heavy quarks are supposed to be dominated by the short range one gluon exchange potential rather than the long range potential resulted from light meson exchanges. Thus, the new hadron states composed by four heavy quarks should be good candidates of compact tetraquark states rather than deuteron-like molecular states.
In 2017, the CMS Collaboration reported their measurement of exotic structure in four lepton channel and found an excess in 18.4 ± 0.1(stat.) ± 0.2(syst.)GeV/c 2 with a global significance of 3.6 σ [12]. This structure indicates a possible bbbb tetraquark state [12][13][14]. It should be noticed that this structure is below the threshold of bottomonium meson pair, which demonstrates that the decays into bottomonium meson pair through quark rearrangement should be hindered. Later, the LHCb and CMS Collaborations analyzed the invariant mass distributions of Υ(1S )µ + µ − , but no evident structure was observed [15,16].
On the theoretical side, the compact tetraquark states composed of bbbb have been investigated extensively, but the conclusions are model dependent. In Refs. [17][18][19][20][21][22][23][24][25], the lowest * Electronic address: lvqifang@hunnu.edu.cn † Electronic address: chendy@seu.edu.cn ‡ Electronic address: dongyb@ihep.ac.cn bbbb tetraquark state is estimated to be below the threshold of bottomonium meson pair, while in Refs. [26][27][28][29][30][31][32][33][34], all the bbbb tetraquark states are above the threshold. To further distinguish different model and reveal the underlying dynamics of fully heavy tetraquark states, more efforts are needed, especially from the experimental side. Very recently, the LHCb Collaboration reported their measurement of the J/ψ pair invariant mass spectrum and a structure near 6.9 GeV/c 2 was observed with the significance greater than 5σ [35]. The resonance parameters are fitted to be in a no-interference scenario, or m = 6886 ± 11(stat.) ± 11(syst.) MeV/c 2 Γ = 168 ± 33(stat.) ± 69(syst.) MeV/c 2 , in an interference scenario. Besides the structure near 6.9 GeV, the experimental data also indicated another two structures in the vicinity of 6.4 GeV and 7.2 GeV, respectively [35]. These structures may be the evidence of compact tetraquark state composed by cccc, which can be a criterion for different models.
After the observation of the LHCb Collaboration, the state around 6.9 GeV has been investigated in different models. In Ref. [36], this state was interpreted as a P−wave tetraquark state in a nonrelativistic quark model, while the QCD sum rule estimations indicated that it could be a second radial excited S −wave tetraquark state [37]. The results in Refs. [38,39] suggested that the resonances with J P = 0 + and 1 + are about 6.4 ∼ 6.6 GeV, while the 2 + state is about 7.0 GeV, which are consistent with the structures reported by LHCb Collaboration [35].
In Ref. [40], we extended the relativized quark model proposed by Godfrey and Isgur to investigate the doubly heavy tetraquarks with the same model parameters. With such an extension, the tetraquaks and conventional mesons can be described in a uniform frame. In the present work, we further study the full heavy tetraquarks QQQQ in the extended relativized quark model and give possible interpretation of the newly observed state around 6.9 GeV. Moreover, the newly observed structures are above the threshold of heavy quarkonium pair, thus, these states can decay into heavy quarkonium pair by quark rearrangement. For simplicity, the decay amplitude should be proportional to the overlap of wave functions of the initial and final states, thus, we can qualitatively discuss the decay behaviors of tetraquarks with the wave functions estimated from the relativized quark model. This work is organized as follows. In section II, we present a review of the extended relativized quark model used in the present work. The numerical results of the masses and decays for the tetraquarks are given in Section III. The last section is devoted to a brief summary.

II. EXTENDED RELATIVIZED QUARK MODEL
To investigate the masses of fully heavy tetraquarks Q 1 Q 2Q3Q4 , we employ an extended relativized quark model, which has been developed very recently for the tetraquark states [40]. It is an extension of the relativized quark model to deal with the four-body systems. The Hamiltonian for a Q 1 Q 2Q3Q4 state can be expressed as where is the relativistic kinetic energy, V oge i j is the one gluon exchange potential including the spin-spin interaction, and V conf i j stands for the confining part. The explicit formula and parameters of relativized potentials can be found in Refs. [40,41].
The wave function of a Q 1 Q 2Q3Q4 state is composed of color, flavor, spin, and spatial parts. In the color space, two types of colorless states with determinate permutation properties exist where the |33 and |66 are antisymmetric and symmetric under the exchange of Q 1 Q 2 orQ 3Q4 , respectively. In the flavor space, the combinations of {cc}, {cc}, {bb}, and {bb} are always symmetric, where the braces { } are adopted to stand for symmetric flavor wave functions.
For the spin part, the six spin bases can be written as, where (Q 1 Q 2 ) 0 and (Q 3Q4 ) 0 are antisymmetric and the (Q 1 Q 2 ) 1 and (Q 3Q4 ) 1 are symmetric for the two fermions under permutations. The matrix elements of the color and spin parts are same as the doubly heavy tetraquarks [40].
In the spatial space, the Jacobi coordinates are presented in Figure 1. For the cccc and bbbb systems, we can define and Then, other relevant coordinates of this system can be obtained in terms of r 12 , r 34 , and r. For a S −wave state, we adopt a set of Gaussian functions to approach its realistic spatial wave function [42] Ψ(r 12 , r 34 , r) = n 12 ,n 34 ,n C n 12 n 34 n ψ n 12 (r 12 )ψ n 34 (r 34 )ψ n (r), (17) where C n 12 n 34 n are the expansion coefficients. The ψ n 12 (r 12 )ψ n 34 (r 34 )ψ n (r) is the position representation of the basis |n 12 n 34 n , where It should be stressed that our final results are independent on geometric Gaussian size parameters r 1 , a, and N max when sufficiently large bases are chosen [42]. The ψ n 12 (r 12 ) and ψ n 34 (r 34 ) can be written in a similar way, and the momentum representation the basis |n 12 n 34 n can be obtained by the Fourier transformation.
According to the Pauli exclusion principle, the total wave function of a tetraquark should be antisymmetric, and possible configurations for cccc and bbbb systems are presented in Table I. With the full wave functions, all the matrix elements of the Hamiltonian can be worked out. Then, the masses can be obtained by solving the following generalized eigenvalue problem where the H i j are the matrix elements in the total bases, N i j is the overlap matrix elements of the Guassian functions arising form the nonorthogonality of bases, E stands for the mass, and C j are the eigenvector corresponding to the coefficients C n 12 n 34 n of spatial wave function. Moreover, for a given system, different configurations with same J PC can mix with each other. The mixing effects are taken into account by diagonalizing the mass matrix of these configurations.

III. RESULTS AND DISCUSSIONS
In present work, we adopt N 3 max = 10 3 Gaussian bases to estimate the S −wave QQQQ spectra. With these large bases, the numerical results are stable enough for our quark model calculations. The predicted masses of ground states for cccc and bbbb systems are presented in Table II. For the cccc system, the masses of four ground states lie in the range 6435 ∼ 6543MeV, which are higher than the J/ψJ/ψ threshold. Compared with the experimental data, we expect that these states should correspond to the broad structure in the vicinity of 6.4 GeV. This broad structure may be one state or an overlap of several states from current data, and more experimental information are needed to clarify its nature. For the bbbb system, the masses are also above the relevant ΥΥ thresholds. These results for the ground states are consistent with nonrelativistic quark model calculations where the pairwise potentials are adopted properly [31][32][33][34]43].
Besides the masses, we can also calculate the proportions of hidden color components and the root mean square radii. In addition to |33 and |66 , other sets of color representations can be defined as and Then, the relations among three sets of color representations can be expressed as follows, and In present work, we adopt the |11 and |88 representations to stand for the neutral color and hidden color components, respectively. The color proportions and root mean square radii of the calculated ground states are displayed in Table III. For the cccc and bbbb systems, the expectations satisfy the following relations From Table III, it can be seen that these states have significant hidden color components and small root mean square radii, and this phenomena indicates that all of them can be regarded as compact tetraquarks. The low-lying radial excitations of cccc and bbbb systems are also calculated in our approach and the results are presented in Table IV. Theoretically, there are two types of radial excitations, the QQ orQQ mode, and the one between QQ andQQ subsystems. The physical states should correspond to the mixture of these two modes. Our results show that the first excitations for the cccc states lie around 6900 MeV, which should correspond to the observed structure near 6.9 GeV by LHCb Collaboration. Given the J/ψJ/ψ decay mode, the J PC of this structure should equal to 0 ++ or 2 ++ . Current information is insufficient to determine its spin-parity, and more theoretical and experimental efforts are needed. Moreover, we find that another set of excitations are around 7050 MeV, which have not been observed by LHCb Collaboration. In addition, the predicted excitations for the bbbb systems are around 19600 and 19730 MeV, respectively, which can easily fall apart into the bottomonium pairs. It is claimed that there is no signal in the Υµ + µ − channels by previous LHCb and CMS experiments [15,16]. We expect that it is due to the small Υµ + µ − branching ratio or lower bbbb production rates relative to the charm sector. Finally, the S −wave spectra for cccc and bbbb systems are plotted in Figure 2 for reference.
The proportions of the hidden color components and root mean square radii for excited QQQQ states are listed in Table V. The small root mean square radii for these states suggest that they are compact tetraquarks, and the sketch of the cccc structure near 6.9 GeV is presented in Figure 3. For the excited states, besides the J/ψJ/ψ and ΥΥ final states, lots of fall-apart channels are also open. The possible decay modes via S −wave and P−wave are listed in Table VI. Since the structure near 6.9 GeV corresponds to a radial excited state in our calculation, the ψ(2S )J/ψ channel is expected to be significant. The ratios can be adopted to describe the relative magnitudes between J/ψJ/ψ and ψ(2S )J/ψ final states. For simplicity, one can assume the decay amplitudes are proportional to the overlap of initial and final states, and the proportional coefficient can be canceled in the final ratios. Here, the wave functions for initial tetraquarks have been obtained by solving the generalized eigenvalue problem, and the wave functions of J/ψ, ψ(2S ), Υ, and Υ(2S ) can be got within the relativized quark model as well. With these wave functions, the ratios for 0 ++ and 2 ++ states can be estimated to be R[cccc(6948)] = 0.075.
Combined with the branching ratios of J/ψ → µ + µ − and ψ(2S ) → µ + µ − , one can further define Then, the ratios R 4µ are predicted to be It can be found that the ψ(2S )J/ψ channel for the excited states is important even though the phase spaces are smaller. The similar situation occurs for the lower excited bbbb states, where the R 4µ of bbbb(19567), bbbb(19625), and bbbb(19633) states are 0.113, 0.111, and 0.084, respectively. These ratios indicate that the lower excited bbbb states can decay to µ + µ − µ + µ − final states through Υ(2S )Υ more easily than ΥΥ mode. Future experiments can search for these states in ψ(2S )J/ψ and Υ(2S )Υ final states.
Our estimations indicate that the broad structure around 6.4 GeV should contain one or more ground cccc tetraquark states, while the narrow structure near 6.9 GeV can be categorized as the first radial excitation of cccc system. The signifi-cant hidden color component and small root mean square radii demonstrate that these states are compact tetraquarks. For the radial excited states, the decay ratios between the J/ψJ/ψ and ψ(2S )J/ψ [or ΥΥ and Υ(2S )Υ] modes are also qualitatively discussed with the wave functions of the tetraquarks and mesons. Our results show that the ψ(2S )J/ψ or Υ(2S )Υ channel is significant for these excited tetraquarks. We hope our sophisticated calculations of the fully heavy tetraquarks may provide valuable information for future experimental searches.