Diffusion Monte Carlo calculations of fully-heavy (multiquark) bound states
Abstract
We use a diffusion Monte Carlo method to solve the many-body Schr\"odinger equation describing fully-heavy tetraquark systems. This approach allows to reduce the uncertainty of the numerical calculation at the percent level, accounts for multi-particle correlations in the physical observables, and avoids the usual quark-clustering assumed in other theoretical techniques applied to the same problem. The interaction between particles was modeled by the most general and accepted potential, i.e. a pairwise interaction including Coulomb, linear-confining and hyperfine spin-spin terms. This means that, in principle, our analysis should provide some rigorous statements about the mass location of the all-heavy tetraquark ground states, which is particularly timely due to the very recent observation made by the LHCb collaboration of some enhancements in the invariant mass spectra of J/-pairs. Our main results are: (i) the cc c c, cc b b (bb c c) and bb b b lowest-lying states are located well above their corresponding meson-meson thresholds; (ii) the JPC=0++ cc c c ground state with preferred quark-antiquark pair configurations is compatible with the enhancement(s) observed by the LHCb collaboration; (iii) our results for the cc c b and bb c b sectors seem to indicate that the 0+ and 1+ ground states are almost degenerate with the 2+ located around 100\,MeV above them; (iv) smaller mass splittings for the cb c b system are predicted, with absolute mass values in reasonable agreement with other theoretical works; (v) the 1++ cb c b tetraquark ground state lies at its lowest S-wave meson-meson threshold and it is compatible with a molecular configuration.
Turn this paper into a lesson
ArcXiv compiles a structured reading guide from this paper's metadata: plain-English importance, contributions, prerequisite concepts, which sections to read first, flashcards, and a quiz. Grounded in the abstract, never invented.