We calculate the energy spectrum of the S-wave doubly heavy tetraquark systems, including the QQ(′)q¯q¯, QQ(′)s¯q¯, and QQ(′)s¯s¯ (Q(′)=b, c and q=u, d) systems within the constituent quark model. We use the complex scaling method to obtain bound states and resonant states simultaneously, and the Gaussian expansion method to solve the complex-scaled four-body Schrödinger equation. With a novel definition of the root-mean-square radii, we are able to distinguish between meson molecules and compact tetraquark states. The compact tetraquarks are further classified into three different types with distinct spatial configurations: compact even tetraquarks, compact diquark-antidiquark tetraquarks, and compact diquark-centered tetraquarks. In the I(JP)=0(1+) QQq¯q¯ system, there exists the D*D molecular bound state with a binding energy of −14 MeV, which is the candidate for Tcc(3875)+. The shallow B¯*B¯ molecular bound state is the bottom analog of Tcc(3875)+. Moreover, we identify two resonant states near the D*D* and B¯*B¯* thresholds. In the JP=1+ bbq¯q¯(I=0) and bbs¯q¯ systems, we obtain deeply bound states with a compact diquark-centered tetraquark configuration and a dominant χ3¯c⊗3c component, along with resonant states with similar configurations as their radial excitations. These states are the QCD analog of the helium atom. We also obtain some other bound states and resonant states with “QCD hydrogen molecule” configurations. Moreover, we investigate the heavy quark mass dependence of the I(JP)=0(1+) QQq¯q¯ bound states. We strongly urge the experimental search for the predicted states.
Published by the American Physical Society
2024