Is $X(7200)$ the heavy anti-quark diquark symmetry partner of $ X(3872)$?

The $D^{(\ast)}\Xi_{cc}^{(\ast)}$ system and $\bar{\Xi}_{cc}^{(\ast)}\Xi_{cc}^{(\ast)}$ system can be related to the $D^{(\ast)}\bar{D}^{(\ast)}$ system via heavy anti-quark di-quark symmetry (HADS). In this work, we employ a contact-range effective field theory to systematically investigate the likely existence of molecules in these systems in terms of the hypothesis that X(3872) is a $1^{++}$~$D\bar{D}^{\ast}$ bound state in the isospin symmetry limit, with some of the unknown low energy constants estimated using the light-meson saturation approximation. In the meson-meson system, a $J^{PC}=2^{++}$~$\bar{D}^{\ast}D^{\ast}$ molecule commonly referred to as $X(4013)$ is reproduced, which is the heavy quark spin partner of $X(3872)$. In the meson-baryon system, we predict two triply charmed pentaquark molecules, $J^{P}=1/2^{-}$~$D^{\ast}\Xi_{cc}$ and $J^{P}=5/2^{-}$~$D^{\ast}\Xi_{cc}^{\ast}$. In the baryon-baryon system, there exist seven di-baryon molecules, $J^{PC}=0^{-+}$~$\bar{\Xi}_{cc}\Xi_{cc}$, $J^{PC}=1^{--}$~$\bar{\Xi}_{cc}\Xi_{cc}$, $J^{PC}=1^{-+}$~$\bar{\Xi}_{cc}\Xi_{cc}^{\ast}$, $J^{PC}=1^{--}$~$\bar{\Xi}_{cc}\Xi_{cc}^{\ast}$, $J^{PC}=2^{-+}$~$\bar{\Xi}_{cc}\Xi_{cc}^{\ast}$, $J^{PC}=2^{-+}$~$\bar{\Xi}_{cc}^{\ast}\Xi_{cc}^{\ast}$ and $J^{PC}=3^{--}$~$\bar{\Xi}_{cc}^{\ast}\Xi_{cc}^{\ast}$. Among them, the $J^{PC}=0^{-+}$~$\bar{\Xi}_{cc}\Xi_{cc}$ and/or $J^{PC}=1^{--}$~$\bar{\Xi}_{cc}\Xi_{cc}$ molecules may contribute to the $X(7200)$ state recently observed by the LHCb Collaboration, which implies that $X(7200)$ can be related to $X(3872)$ via HADS. As a byproduct, with the heavy quark flavor symmetry we also study likely existence of molecular states in the $B^{(\ast)}\bar{B}^{(\ast)}$, $\bar{B}^{(\ast)}\Xi_{bb}^{(\ast)}$, and $\bar{\Xi}_{bb}^{(\ast)}\Xi_{bb}^{(\ast)}$ systems.


I. INTRODUCTION
Hadrons can be classified either as mesons made of a pair of quark and antiquark or as baryons made of three quarks in the conventional quark model, which has been rather successful at least up to 2003. According to the quark model [1], more complicated quark configurations, such as multiquark states, should also exist, which are also allowed by QCD. In recent years a lot of XY Z and P c states that can not easily fit into the conventional quark model were discovered by many collaborations all over the world, which not only challenges the conventional quark model but also opens a new era for hadron physics [2][3][4][5][6].
Being close to the mass threshold of DD * , it is naturally considered to be a hadronic molecule. In addition, another convincing evidence of the molecular nature of X(3872) is its isospin breaking decays into J/Ψππ and J/Ψπππ [15,16]. A series of theoretical works interpreted X(3872) as a J P C = 1 ++D D * molecule [17][18][19][20][21][22][23][24][25]. Although the molecular interpretation seems to be the most popular, its nature has not been completely unveiled. 1 Nonetheless, given current experimental and theoretical results, one may conclude that the molecule components of X(3872) must play a relevant (if not dominant) role. In our previous work [25], assuming that X(3872) is a J P C = 1 ++D D * molecule, we have used the one boson exchange model to predict the binding energies and scattering lengths of the heavy quark spin symmetry (HQSS) partners of X(3872).
HQSS plays an important role in describing the spectra of heavy hadrons including both conventional and exotic ones. The D and D * mesons as well as Σ c and Σ * c baryons can be regarded as doublets of HQSS. Applying HQSS to their bottom partners, the mass splittings of doublets become small, which indicates that HQSS works well in the heavy quark limit. The interactions between heavy hadrons are also constrained by HQSS, which can further decrease the number of unknown couplings and increase the predictive power of effective field theories (EFT) as well. In the meson-meson system, D * s0 (2317) and D s1 (2460) are assigned as DK and D * K hadronic molecules, respectively, belonging to a HQSS doublet [35][36][37][38][39][40][41][42]. A lot of theoretical works described Z c (3900) and Z c (4020) as isovector resonances generated by theDD * andD * D * interactions, respectively, which indicates that Z c (3900) can be related to Z c (4020) via HQSS [23,43,44]. Very recently Z cs (3985) was proposed to be the SU (3)-flavor partner of Z c (3900) [45][46][47]. In the heavy meson and baryon systems, P c (4312), P c (4440), and P c (4457) are nicely arranged into a multiplet ofD ( * ) Σ ( * ) c hadronic molecules [48][49][50][51][52][53][54][55], and in terms of HQSS there should exist other four states which need to be confirmed by future experiments.
Heavy quark flavor symmetry has already shown its power in constructing the spectra of charmed and bottom hadrons. In short, if there exists a charmed hadron there should also exist the corresponding bottom partner. Heavy anti-quark di-quark symmetry(HADS) dictates that in terms of color degree of freedom a pair of heavy quarks can be related to a heavy antiquark in the heavy quark limit [56]. From this symmetry we can derive the following mass relation: [57]. Although this relation has not been confirmed by experiments so far, it is consistent with lattice QCD calculations and phenomenological studies [58][59][60][61]. In our previous work, we have employed HADS to extend theD ( * ) Σ Along the line of Ref. [62], we can extend the D ( * )D( * ) system to the D ( * ) Ξ spectra of these two systems could not be determined. In this work we employ the light-meson saturation approach to estimate the ratio of the two couplings [65]. In our previous work [55], we found that the two couplings of theD ( * ) Σ ( * ) c system determined by the contact-range EFT and the light-meson saturation approach are quite similar. In Ref. [66] Peng. et al used the lightmeson saturation approach to estimate the contact-range potential betweenDΞ c andDΞ c , yielding results consistent with the local hidden gauge approach [67,68]. Therefore, in this work, to fully determine the mass spectra of these systems, we will combine the hypothesis that X(3872) is a Recently, the LHCb Collaboration reported one resonant state with a mass of 6.9 GeV in the J/ψ-J/ψ invariant mass spectrum, labelled as X(6900). In addition, there exist two more structures, a broad one from 6.2 to 6.8 GeV and a vague structure around 7.2 GeV [69]. These tetraquark states containing fully charmed quarks undoubtedly open another new era for hadron physics since the discovery of X(3872), which has motivated a lot of theoretical studies [70][71][72][73][74][75][76][77][78][79][80][81][82][83][84][85][86][87].
Two popular interpretation of these structures are resonances generated by several coupled channels of charmonium pairs and compact tetraquark states of diquark and anti-diquark. In this work we mainly focus on the structure near 7.2 GeV, namely X(7200), whose mass is close to the mass threshold ofΞ cc Ξ cc . Since theΞ ( * ) cc Ξ ( * ) cc system can be related to theD ( * ) D ( * ) system via HADS, we further explore whether there existΞ cc Ξ cc molecules from the hypothesis that X(3872) is a J P C = 1 ++ bound state, in other words, whether X(7200) can be related to X(3872) via HADS.
The manuscript is structured as follows. In Sec. II we present the contact-range potentials of

II. FORMALISM
Generically, the interaction between two heavy hadrons can be decomposed into a long-and a short-range piece within the EFT. The long-range piece corresponds to the one-pion-exchange potential, while the short-range piece can be described by a series of contact-range potentials with unknown couplings. In the present context, the short-range interaction between two heavy hadrons is constrained by heavy quark symmetry. In the following we explain how to derive the cc interactions. In the line of Refs. [62,88], their interactions can be determined in the EFT approach. One should note that we just consider the leading order contactrange potentials because our previous studies indicated that the pion exchanges are perturbative in the charm sector [89,90]. The short-range interaction of the meson-meson system can be written as where C a and C b are two unknown couplings and H is a non-relativistic superfield with H = 1 √ 2 [P + P * · σ], and P and P * denote pseuoscalar and vector charmed meson fields. The D ( * ) Ξ ( * ) cc system can be related to the D ( * )D( * ) system through HADS, resulting in the following Lagrangian where T = 1 √ 3 σΞ cc + Ξ * cc is a non-relativistic superfield of Ξ cc and Ξ * cc respecting HQSS. The Lagrangian of theΞ The term C a is independent of the spin, while C b accounts for the spin spin interaction. Using the same approach of Ref. [62], the contact potentials of the D ( * )D( * ) , cc systems can be easily constructed and are shown in Table I. To make concrete predictions for likely existence of molecular states we have to solve a nonrelativistic bound state equation with the contact-range potentials of Table I. In the momentum space, we can solve the Lippmann-Schinwinger equation to find bound state as where φ(k) is the vertex function, B the binding energy, and µ the reduced mass. To solve the equation we regularize the contact potential in the following way with Λ the cutoff, f (x) a regular function, and C(Λ) the running coupling constant. A typical choice of the cutoff is Λ = 0.5 − 1 GeV, while for the regulator we choose a Gaussian type baryon, and doubly heavy baryon-doubly heavy anti-baryon systems depending on two unknown coupling constants, C a and C b . These coupling constants can be determined from the sum of C a and C b by reproducing the mass of X(3872) and the ratio of C a and C b by the light-meson saturation approach. state In this work, to be consistent with the cutoff adopted in the light-meson saturation approach, we use a cutoff of Λ = 0.7 GeV. Considering only S-wave contact interactions, the integral equation simplifies to

III. NUMERICAL RESULTS AND DISCUSSIONS
In this section we study the likely existence of molecular states in the D ( * )D( * ) , D ( * ) Ξ ( * ) cc , and Ξ ( * ) ccΞ ( * ) cc systems. The contact-range potentiasl of these systems are only dependent on two couplings, C a and C b . To determine these couplings, the best approach is to fit to the experimental data. However, there exists only one molecular candidate X(3872) to determine the sum of C a and C b . We, therefore, resort to the light-meson saturation approach to estimate the ratio of C a   and C b .
The assumption of light-meson saturation implies that the values of the EFT couplings are saturated by light meson exchanges. Following the formalism of Ref. [65], we estimate the couplings of C a and C b by the sigma and vector meson exchange saturation, where the cutoff Λ ∼ m σ , m V implies that the saturation works at an EFT cutoff close to the masses of exchanged mesons, i.e., 0.6 ∼ 0.8 GeV. The values of saturated couplings are expected to be proportional to the potentials of light meson exchanges in the OBE model once we have removed the spurious Dirac-delta potential. Thus C a and C b can be written as where g σ 1 denotes the charm meson coupling to sigma meson, g v1 and f v1 denote electric-type and magnetic-type couplings between charm mesons and light vector mesons, and M is a mass scale to render f v1 dimensionless, see, e.g., Table II. The proportionality constant is unknown and depends on the details of the renormalization procedure. However, assuming that the constant is the same for C sat a and C sat b , we can calculate their ratio.
In the following, we explain how to determine the values of C a and C b in detail. Assuming that X(3872) is a J P C = 1 ++ D * D bound state with the contact potential C a + C b , we can determine the sum of C a and C b as 31.4 GeV −2 by Eq. (6), where the cutoff is fixed at 0.7 GeV. With this sum we can only predict the binding energies of states for which the contact potentials are C a +C b .
Fortunately, combining the ratio of C a and C b determined by the light-meson saturation, we can determine the values of C a and C b . The ratio is Comparing the ratio C b /C a ≈ 0.134 in the meson-baryon system of Ref. [55], we find that the spin-spin term plays a more important role in the meson-meson system, which could be the main reason why we can not obtain a complete multiplet of hadronic molecules in the heavy meson-meson system. Using the ratio determined by the light-meson saturation we finally obtain  Table III. For the D ( * )D( * ) system, we obtain a J P C = 2 ++ D * D * bound state with a binding energy of 5 MeV, namely X(4013), which has already been predicted by many theoretical studies [22,23,93]. However, its existence has not been confirmed by experiments. Except for X(3872) and X(4013) there are no other bound states in this system, which is consistent with the conclusion of the OBE model [25]. From the scattering length of J P C = 0 ++D D we find that the potential betweenD and D is attractive, but is too weak to form a bound state. Recently a Lattice QCD study obtained a shallow bound state located 4 MeV below the DD mass threshold [94].
If we use the OBE model to reproduce the 4 MeV binding energy, the cutoff needed is about Λ = 1.4 GeV, whose value is close to but larger than that obtained in our previous work to reproduce X(3872), i.e., Λ = 1.01 GeV. We note that several theoretical works have predicted   the existence of a bound state in the DD channel [21,38,[95][96][97].
With heavy quark flavor symmetry we can also study the B ( * )B( * ) system. We find 5 bound states, J P C = 0 ++B B, J P C = 1 ++B B * , J P C = 1 +−B B * , J P C = 1 +−B * B * , and J P C = 2 ++B * B * , consistent with the results of the OBE model [25]. Among them we can regard the J P C = 1 ++B B * and J P C = 2 ++B * B * states, referred to as X b and X b , as the hidden bottom partners of X(3872) and X(4013), respectively. We note that Törnqbist has already predicted X(4013) with a binding energy 45 MeV in the one pion exchange model [98], but its existence has not been confirmed by experiments.

B. Meson-baryon system
Using HADS we can extend the D ( * )D( * ) system to the D ( * ) Ξ ( * ) cc system as shown in Fig. 1, where a HQSS multiplet of meson-baryon molecules can be expected. The contact-range potentials of the D ( * ) Ξ ( * ) cc system are expressed as combinations of C a and C b , which can be related to the couplings arising in the D ( * )D( * ) system by HADS. Therefore, we take the same input to calculate the binding energies and scattering lengths of the D ( * ) Ξ ( * ) cc system, and the results are given in Table IV tems, J P = 1/2 − D * Ξ cc and J P = 5/2 − D * Ξ * cc , always remain bound. Thus in this work we predict the existence of two triply charmed pentaquark molecules. It is to be noted that in our previous study, we also predicted a J P = 1/2 − D * Ξ cc bound state in the OBE model [99], while in Ref. [64] a J P = 5/2 − D * Ξ * cc bound state is predicted using a contact-range potential. These two molecules can decay into Ω ccc ω via triangle diagrams, and we expect that they can be detected in future experiments.
With heavy quark flavor symmetry we also study theB ( * ) Ξ ( * ) bb system, and find that bound states exist in all the channels within uncertainties induced by heavy quark symmetry breaking, which indicates that we obtain a complete HQSS multiplet of hadronic molecules. Comparing the results with those of the meson-meson system, we find that the meson-baryon molecular states are more bound, in agreement with naive expectations because of the smaller reduced masses in the latter system.

C. Baryon-anti-baryon system
In this section, we further extend the D ( * ) Ξ ( * ) cc system to theΞ molecule   Table V. We find seven bound states, i.e., J P C = 0 −+Ξ cc Ξ cc , J P C = 1 −−Ξ cc Ξ cc , J P C = 1 −+Ξ cc Ξ * cc , J P C = 1 −−Ξ cc Ξ * cc , J P C = 2 −+Ξ cc Ξ * cc , J P C = 2 −+Ξ * cc Ξ * cc , and J P C = 3 −−Ξ * cc Ξ * cc . TheΞ cc Ξ cc system [100] and theΞ * cc Ξ * cc system [101] have been studied in the OBE model, which show that all of these systems are bound. Assuming heavy quark flavor symmetry, we further investigate their hidden bottom partnersΞ bb , and all the states are found to be bound as shown in Table V, suggesting the emergence of a complete multiplet of hadronic molecules.
Among these molecules, we note that there are two molecules, J P C = 0 −+Ξ cc Ξ cc and J P C = 1 −−Ξ cc Ξ cc , located close to 7200 MeV, which could contribute to the structure X(7200) observed by the LHCb Collaboration recently. These two molecules are predicted in terms of the hypothesis that X(3872) is a J P C = 1 ++D D * bound state and utilyzing the HADS. It is well known that the heavy meson-anti meson channel affects the spectra of conventional charmonia. According to FIG. 2. Left: the two charmed meson channel heavily affects the mass spectra of charmonia close to their mass thresholds. Right: HADS implies that the two doubly charmed baryon channel could also affect the mass spectra of fully charmed tetraquark states close to their mass thresholds. the relationship between X(7200) and X(3872) proposed by the present work it is rather likely that the two doubly charmed baryon-antibaryon channel would affect the spectra of fully charmed tetraquark states as shown in Fig. 2, which indicates that X(7200) might contain bothΞ cc Ξ cc molecular and charmed diquark-anti diquark tetraquark components. It is interesting to note that our study shows that X(7200) might be a superposition of two molecular states.

IV. SUMMARY
In this work we extended the D ( * )D( * ) system to the D ( * ) Ξ ( * ) cc system then to theΞ ( * ) cc Ξ ( * ) cc system utilyzing the heavy antiquark diquark symmetry (HADS). The contact-range potentials of these systems respecting heavy quark spin symmetry (HQSS) are constrained to two couplings at the leading order, C a and C b . Assuming that X(3872) is a J P C = 1 ++D D * bound state, we determined the sum of C a and C b . To fully determine the value of C a and C b , we used the lightmeson saturation approach to determine its ratio. With the so-determined C a and C b we calculated the binding energies (if bound states exist) and scattering lengths of all these systems as well as their heavy quark flavor partners.
In theD ( * ) D ( * ) system, we obtained a J P C = 2 ++D * D * molecule with a mass of 4013 MeV referred to as X(4013), which can be regarded as the heavy quark spin symmetry partner of X(3872). Moreover, in the bottom sector we predicted the existence of five molecules. We extended theD ( * ) D ( * ) to the D ( * ) Ξ ( * ) cc system by HADS. Two molecules, J P = 1/2 − D * Ξ cc and J P = 5/2 − D * Ξ * cc , are obtained. There exist a complete multiplet of hadronic molecules in the corresponding bottom sector. Using HADS again we studied theΞ ( * ) cc Ξ ( * ) cc system, and found 7 bound states, J P C = 0 −+Ξ cc Ξ cc , J P C = 1 −−Ξ cc Ξ cc , J P C = 1 −+Ξ cc Ξ * cc , J P C = 1 −−Ξ cc Ξ * cc , J P C = 2 −+Ξ cc Ξ * cc , J P C = 2 −+Ξ * cc Ξ * cc and J P C = 3 −−Ξ * cc Ξ * cc . There also exist a complete multiplet of hadronic molecules in the bottom sector. Among them, the J P C = 0 −+Ξ cc Ξ cc and J P C = 1 −−Ξ cc Ξ cc molecules may contribute to X(7200) recently observed by the LHCb Collaboration, which are related to X(3872) via HADS.