The S- and P-wave fully charmed tetraquark states and their radial excitations

Inspired by recent progresses in observations of the fully charmed tetraquark states by LHCb, CMS, and ATLAS collaborations, we perform a systematic study of the ground states and the first radial excitations of the S- and P-wave ccc¯c¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textrm{cc}\bar{\textrm{c}}\bar{\textrm{c}}$$\end{document} system. Their mass spectra, root mean square (r.m.s.) radii and radial density distributions are studied with the relativized quark model. The calculations show that there is no stable bound states for the fully charmed tetraquark states, and the r.m.s. radii of these tetraquark states are smaller than 1 fm. Our results support assigning X(6600) structure, MX(6600)=6552±10±12\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_{X(6600)}=6552\pm 10\pm 12$$\end{document} MeV, as one of the 0++\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0^{++}$$\end{document}(1S) and 2++\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$2^{++}$$\end{document}(1S) states or their mixtures. Another structure also named as X(6600) by CMS Collaboration, MX(6600)=6.62±0.03-0.01+0.02\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_{X(6600)}=6.62\pm 0.03^{+0.02}_{-0.01}$$\end{document} GeV, may arise from the lowest 1P states with JPC=0-+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J^{PC}=0^{-+}$$\end{document}, 1-+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1^{-+}$$\end{document}, and 2-+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$2^{-+}$$\end{document}. The possible assignments for X(6900) include the 0++\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0^{++}$$\end{document}(2S), 2++\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$2^{++}$$\end{document}(2S) states, and the highest 1P state with JPC=0-+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J^{PC}=0^{-+}$$\end{document}. As for X(7200), it can be interpreted as one of the highest 2P states with JPC=0-+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J^{PC}=0^{-+}$$\end{document}, 1-+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1^{-+}$$\end{document}, and 2-+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$2^{-+}$$\end{document}, and the 3S states can not be completely excluded from the candidates.


Introduction
Since the first observation of X(3872) by Belle in 2003 [1], a series of exotic hadrons named as XYZ states emerged like bamboo shoots after a spring rain.These discoveries in experiments have motivated theorists to devote a great deal of energy in studying the inner structure of these new states.In general, these new observed states can not be categorized as the conventional mesons or baryons, and they are commonly explained as compact tetraquark states, hadronic molecular states or admixture of these two states.Among these new hadrons, the states composed of fully heavy quark components are especially interesting.It was supposed that the interactions between heavy quarks(antiquarks) are dominated by the short range one-gluon-exchange(OGE) potential rather than the long-range potential resulted from light meson exchanges.Thus, the configurations comprised of four heavy quarks are easier to form compact tetraquark states instead of the molecules.
In 2020, the LHCb collaboration studied the J/ψJ/ψ invariant mass spectrum using the pp collision data at center-of-mass energies of √ s=7, 8 and 13 TeV [2].They observed a broad structure above the J/ψJ/ψ threshold ranging from 6.2 to 6.8 GeV and a narrow structure around 6.9 GeV/c 2 with the significance of larger than 5σ.At the ICHEP 2022 conference, the ATLAS collaboration reported their observation of several fully charmed tetraquark excesses decaying into a pair of charmonium states in the four µ final states [3].They not only confirmed the existence of the X(6900) structure reported earlier by LHCb Collaboration but also observed a broad structure at lower mass and a structure labelled as X(6600).In addition, the CMS Collaboration also reported their studies on the J/ψJ/ψ mass spectrum, which were carried out by using proton-proton data at center-of-mass energies of 13 TeV [4].The existence of X(6900) was confirmed with significance larger than 9.4σ.Besides X(6900), the CMS collaboration also observed two structures in the J/ψJ/ψ channel, which were labeled as X(6600) and X(7200).All of the experimental data about the fully charmed tetraquark states are collected in Table I.Actually, the full-heavy tetraquark states were already studied in the literatures  before the experimental observations.The recent breakthroughs in experiments in searching for the fully charmed tetraquark states have inspired again the intensive discussions.People employed many methods/models to carry out their studies on the structure, the production, and decay property of these new states .The most popular interpretations about these new discoveries are compact tetraquark states [70][71][72][73][74][75][76][77][78][79].In Ref. [34], Q. F. Lü et al. analyzed the mass spectrum of the S-wave full-heavy tetraquark states with a extended relativistic quark model.In their studies, the X(6900) was categorized as the first radial excitation of cccc system.Using the method of QCD sum rules, Z.
In our previous work, we used the relativistic quark model to study the mass spectra, r.m.s.radii and the radial density distributions of the singly and doubly charmed baryons.At present, we extend our previous method to analyze the fully charmed tetraquark states.We hope this study can help to shed more light on the nature of these new exotic states.The paper is organized as follows.After the introduction, we briefly describe the phenomenological method adopted in this work in Sec.II.
In Sec.III we present our numerical results and discussions about the full-charmed tetraquark states.
And Sec IV is reserved for our conclusions.
2 Phenomenological method adopted in this work The relativistic quark model has been successfully extended to study the mass spectrum of the S-wave tetraquark state [34].In the following, we give a brief introduction to the Hamiltonian of relativized quark model.The Hamiltonian for a four-body system is composed of the relativistic kinetic energy term, the confining potentials and one-gluon exchange potentials [87,88], where the relativistic kinetic energy term is, The confining potential V conf ij is written as, with In Eq.( 3), F i • F j stands for the color matrix and F reads for quarks, for antiquarks (5) with The one-gluon exchange potential All of the one-gluon exchange potentials can be expressed in terms of a smeared one-gluon exchange propagator G which can be written as where and The hyperfine interaction V hyp includes tensor interaction and the contact interaction, with For the spin-orbit interaction, it can also be divided into two parts which can be written as, with and In Eqs.(10), (11),( 13) and ( 14) and V so(s) ii are achieved from the G(r ij ) and confining potential V conf ij (r ij ) by introducing momentum-dependent factors, with , where p ij is the magnitude of the momentum of either of the quarks in the ij center-of-mass frame.
The internal motions of the quarks in a four-body system can be expressed by three sets of Jacobi coordinates as shown in Fig. 1.As for the Jacobi coordinates in Fig. 1(a), they can be defined as, According to the Jacobi transformation, the other two sets of coordinates can also be expressed in respectively.In the flavor space, the fully charmed tetraquark state where the square bracket denotes the flavor symmetry.For a double quark(untiquark) system in the fully charmed tetraquark state, its spin wave function is antisymmetric singlet or symmetric triplet and they can be expressed as In this work, the Gaussian basis function is employed to construct the spatial wave function of the tetrquark state, and the Gaussian basis reads with In Eq.( 25), n max is the maximum number of the Gaussian basis functions.The total wave function of the spin and spatial parts with the angular momentum (J,M ) can be written as, where κ is the quantum numbers {n 12 ,l 12 ,s 12 ,j a ,n 34 ,l 34 ,s 34 ,j b ,n,l,j} of the basis.Theoretically, a tetraquark state with total angular momentum J is the superposition of all the bases.The basis can be classified by the total orbital angular momentum L=l 12 +l 34 +l.The tetraquark states with different orbital angular momentum L but with same J P C will mix with each other, for example, the 0 ++ state of L=0 will couple with L=2, and L=4 states through spin-orbital and tensor potentials.
This interaction will influence the mass spectrum slightly, thus we neglect this mixing mechanism in this work.
According to the Pauli exclusion principle, the total wave function of a tetraquark should be antisymmetric, and all possible configurations for fully charmed tetraquark are presented in the third column of Table IV.For a S-wave tetraquark state, its possible spin and parity quantum numbers are J P C =0 ++ , 1 +− , and 2 ++ .For the P -wave states, there are two orbital excited modes, the ρ-mode with the orbital excitation in the diquark or antidiquark, i.e., (l 12 ,l 34 ,l)=(1,0,0) or (0,1,0), and the λ-mode with the orbital excitation between the two clusters, i.e., (l 12 ,l 34 ,l)=(0,0,1).Besides of the conventional quantum numbers, i.e., J P C = 0 −+ , 1 −− , 2 −+ , 3 −− , the P -wave can also access exotic quantum numbers, i.e., J P C = 0 −− , 1 −+ , 2 −− .For simplicity, each configuration in Table IV can be expressed as |c C ρ/λ ; 2S+1 L J , where c=3 or 6 stand for the color configuration, C and S are the C-parity and the total spin angular momentum of the configuration.
For a four-body system, the calculations of the Hamiltonian matrix elements become laborious even with Gaussian basis functions.This process can be simplified by introducing the ISG basis functions.
These new sets of basis functions can be written as [89], where ε is the shifted distance of the Gaussian basis.Taking the limit ε → 0 is to be carried out after the matrix elements have been calculated analytically.For more details about the calculations of the Hamiltonian matrix elements, one can consult the Refs.[89,90].
After all of the matrix elements are evaluated, the mass spectra can be obtained by solving the generalized eigenvalue problem, Here, H ij denotes the matrix element in the total color-flavor-spin-spatial base, E is the eigenvalue, C j stands for the corresponding eigenvector, and N ij is the overlap matrix elements of the Gaussian functions, which arises from the nonorthogonality of the bases and can be expressed as, 3 Numerical results and discussions 3.1 S-and P-wave tetraquark states The results of the relativized quark model depend on the input parameters such as the constituent quark mass and the parameters in the Hamiltonian.Up to now, there has been no solid experimental data for the spectrum of pure tetraquark states, it is impossible for us to fix these parameters by fitting present experimental data.In most cases, these parameters were determined by fitting them to the experimental masses of the mesons or baryons [30,33,35,[91][92][93].In the present work, all input parameters are listed in Table II and taken from the original Refs.[87,88] where they successfully reproduced the excremental data of mesons and baryons.In our previous work, predictions for heavy baryons with these parameters were indeed consistent well with the experimental data [90].Recently, these parameters were already extended to study the tetraquark states in a uniform frame.As for the accuracy of the relativized quark model, it depends on the quenched approximation and relativistic corrections.Considering these two effects in Ref. [87], they claimed that the average accuracies are 25 MeV for light and heavy-light mesons and 10 MeV for heavy mesons, respectively.We expect that the uncertainties of predicted masses of the teatraquark states are limited in a reasonable range.
In order to investigate the convergence and stability of the numerical results, we plot the masses of the configurations |3 − ; 3 S 1 and |3 + λ ; 3 P 0 in Fig. 2(a) and their r.m.s.radii r 2 12 in Fig. 2(b).In these figures, the wave function is expanded with Gaussian basis n 3 max =2 3 , 3 3 , 4 3 , 5 3 , 6 3 , 7 3 , and 8 3 , respectively.We can see that the results decrease with the basis number and converge to a stable value when n 3 max > 6 3 .In this work, the calculations are carried out with n 3 max = 8 3 Gaussian basis to ensure the stability of the final results.In general, the physical state of a tetraquark are the mixtures of different configurations with the same quantum number of J P C .Thus, the calculations are carried out in two stages.In the first stage, the masses of different configurations shown in the third column in Table IV are obtained by solving the Schrödinger equation with the variational method.The masses and the r.m.s.radii for the ground and the first radially excited states are also presented in Table IV.In the second stage, the mixing effect is considered and the masses of physical states are obtained by diagonalizing the mass matrix in the basis of eigenstates obtained in the first stage.The mass matrix, the eigenvalues, and the eigenvectors for the ground states and the first radial excitations are summarized in Tables V-VI, respectively.
In the OGE model, the interaction between the two quarks within a color-sextet diquark is repulsive, while that in the color-anti-triplet one is attractive.On the other hand, the interaction between the diquark and antidiquark of |(Q 1 Q 2 ) 6 ( Q3 Q4 )6 configuration is attractive and is much stronger than that of the |(Q 1 Q 2 )3( Q3 Q4 ) 3 one.From Table IV, one can see that |6 + ; 1 S 0 configuration is located lower than the |3 + ; 1 S 0 .This can be explained by the stronger attractive potential between diquark and antidiquark in |6 + configuration.This also applies to the ρ-mode excitations, i.e., 6 ρ <3 ρ .If the λ-mode excitations are also considered, the relationship is slightly complicated, where the relationship between different color configurations is 6 ρ <3 λ <3 ρ <6 λ , which can be seen in Table IV.This is the result of the competition of various interactions among quarks.
From Table IV, it is shown that the r.m.s.radii r 2 12/34 of |3 + and |3 ± ρ states are smaller than those of the |6 + and |6 ± ρ , respectively.This is also due to the attractive interactions between the two quarks within a color-antitriplet diquark and repulsive ones in a color-sextet diquark.On the other hand, the stronger attraction between diquark and antidiquark in where Ω 12/34 and Ω are the solid angles spanned by vectors r 12/34 and r, respectively.Some of the results are shown in Figs.3-8.It can be seen that the radial density distributions are in the range of 1 fm, which indicates that the four quarks are confined into a compact state.Second, the peaks for the first radial excitations are located more outward compared with their ground states.Third, for the λ-mode excitation, e.g.|3 + λ ; 3 P 0 , it is shown in Fig. 6 that the r 2 ω(r) peak is located more outward than that of r 2 12/34 ω(r 12/34 ), while the situation is opposite for ρ-mode excitation, e.g.|3 + ρ ; 3 P 0 in Fig. 7.
To investigate the inner structure of the fully charmed tetraquark, we also obtain its proportions in the meson-meson configuration and the r.m.s radii of these states, which are shown in the last four columns of Tables VII-VIII.For 0 ++ (6450) as an example, it is composed of 37  [94,95] and a lattice-QCD inspired quark model [48].In Ref. [48], the author analyzed a complete set of four-body configurations, including diquark-antidiquark, meson-meson, K-type configurations and their couplings.Here, we only list the result of 8 c ⊗ 8 c (meson-meson) configuration.In Ref. [70], they studied the 3 c ⊗ 3 c (diquark-antidiquark) and molecular configurations using the inverse Laplace transform sum rule.Another description for the tetraquark are the butterfly(short for but) and flipflop configurations adopted in Ref. [96], where their results are larger than those of others.The results in the last column were obtained using a two-body Bethe-Salpeter equation by considering the mesonmeson components plus diquark-antidiquark components.It can be seen that our results are roughly compatible with those in Refs.[34,48,70,94,95].Our results are consistent well with those of Ref. [34] where they adopted the same method in their calculations except not considering the spin-dependent interactions in the Hamiltonian.

The newly observed tetraquark states
In order to exhibit the mass spectrum of the fully charmed tetraquark states more obviously, the results are also displayed in Fig. 9.The model predictions for the 1S-wave tetraquarks are 0 ++ (6450), 0 ++ (6534), 1 +− (6517), and 2 ++ (6544), which are higher than the J/ψJ/ψ threshold.This indicates that there may not exist bound tetraquark states cccc in the scheme of relativistic quark model.From Fig. 9, one can see that the mass gaps between the ground states and the first radial excitations are about 300∼400 MeV.This behavior is very similar with that of the doubly charmed baryon spectra.
Another important feature about the mass spectrum of the tetraquark states is that the predicted tetraquark states in theory are much richer than the experimental data.One possible explanation about this feature is that some higher resonances have large decay widths and their signals annihilate in the background signals.Even if these resonances are observed in experiments, they may appear as a broad structure in the invariant spectrum.Finally, some states are located very near with each other, e.g.0 ++ (6450), 0 ++ (6534) and 2 ++ (6544).If the mass splitting of these states is smaller than their decay widths, they will overlap with each other and contribute to a broad structure in the invariant spectrum.range of this structure, thus these 1S states are good candidates for this structure.For the 1P -wave tetraquarks, the low-lying 0 −+ (6666), 1 −+ (6624), and 2 −+ (6647) are all located in a structure which was also named as X(6600) by ATLAS Collaboration but with a measured mass to be 6.62 ± 0.03 +0.02 −0.01 GeV.This structure may be interpreted as a mixture of these low-lying 1P states.
As for the X(6900), the reported results by LHCb, CMS, and ATLAS Collaborations do not consist well with each other.It is shown in Fig. 10 that the lowest 2S state 0 ++ (6861) is consistent with the experimental data 6886 ± 11 ± 1 MeV and 6870 ± 30 +60 −10 MeV.On the other hand, the 2S states 0 ++ (6955) and 2 ++ (6952), and a P -wave 0 −+ (6939) state are compatible with the results 6927 ± 9 ± 5 MeV.If X(6900) was treated as a S-wave radial excitation, it can decay into the lowest twocharmonium channels η c η c and J/ψJ/ψ via S-wave transition.The decay widths of these decay channels will be broad because of the large phase space.As a P -wave tetraquark state, X(6900) can also decay into J/ψJ/ψ channel with a P -wave orbital excitation, but its width is kinetically suppressed.
The CMS Collaboration also reported a tetraquark state named as X(7200) with a measuring mass to be 7287 ± 19 ± 5. It can be seen in Fig. 10 that the high-lying 2P state 0 −+ (7298) is located in the range of this structure.This means the 0 −+ (7298) is a possible candidate for X(7200) in the framework of relativistic quark model.Considering the decay width of X(7200) plus the uncertainties of the quark model, the 1 −+ (7338) and 2 −+ (7321) can not be completely excluded.In our present work, we do not include the predicted results of the 3S states.According to the mass gaps between the ground states and the first radial excitations, we can roughly estimate the order of magnitude of 3S states to be 7200∼7300 MeV.This implies the 3S excitations might also be the candidates of X(7200) state.

Conclusions
In this work, we have systematically investigate the mass spectra, the r.m.s.radii and the radial density distributions of the S-and P -wave cccc system.The calculation is carried out in the frame work of relativistic quark model, where the Coulomb term, confining potential, tensor potential, contact interaction, and the spin-orbital interaction are all considered.We first analyze the masses, r.m.s.radii, and radial density distributions of different color configurations which are categorized by 3c ⊗ 3 c and 6 c ⊗ 6c .Then, we obtain the mass spectrum and r.m.s.radii of the physical states by considering the mixing effect.
The results show that the r.m.s.radii of the cccc system are less than 1 fm, and the radial density distributions are also in this range.This implies that the full-charmed systems are compact tetraquark states.For the S-wave states, the predicted masses are all higher than the thresholds of heavy charmonium pairs, which means no stable bound states exist for cccc system.The results also show that the relationship of different color configurations among the mass spectrum is 6 ρ <3 λ <3 ρ <6 λ for P -wave full-charmed tetraquarks.

FIG. 2 :
FIG. 2:The dependence of the mass and r.m.s.radius on the number of Gaussian basis.

TABLE I :
The experimental data for the fully charmed tetraquark states.

TABLE II :
Relevant parameters of the relativized quark model fm and 0.309 fm for the ground state of |3 + configuration, while the results are 0.456 fm and 0.235 fm for |6 + one.To further understand the inner structures of different configurations, we also plot the radial density distributions which are obtained by the wave functions from quark model.The distribution functions are defined as, ω(r 12/34 ) = |Ψ(r 12 , r 34 , r)| 2 drdr 34/12 dΩ 12/34 ω(r) = |Ψ(r 12 , r 34 , r)| 2 dr 12 dr 34 dΩ .6% |3 + ; 1 S 0 (3 c ⊗ 3 c ) and 62.4% |6 + ; 1 S 0 (6 c ⊗ 6 c ) components.This state contain 54.1% 1 c ⊗ 1 c configuration and 45.9% 8 c ⊗ 8 c ones.In addition, predicted masses of other collaborations for the ground states of the cccc system are listed in Table III.The 8 c ⊗ 8 c configurations in the table are studied in different QCD sum rules

TABLE III :
Predicted masses of the ground states of tetraquark state in different configurations.All results are in units of GeV

TABLE V :
The mass matrix (MeV), the eigenvalue (MeV) and the eigenvector for the ground state by diagonalizing the mass matrix.

TABLE VI :
The mass matrix (MeV), the eigenvalue (MeV) and the eigenvector for the first radial excited state by diagonalizing the mass matrix.