A diffusion Monte Carlo calculation of fully heavy pentaquarks

Using a diffusion Monte Carlo algorithm, we calculated the spectra of all possible $S$-wave fully heavy pentaquarks within the framework of the quark model. Our aim was to compare the masses of different spin-color configurations corresponding to the same spatial wavefunction for each quark composition. In particular, we computed the masses of the configuration with the maximum number of undistinguishable quarks and those of all the possible baryon+meson splittings compatible with the total number of $b$'s and $c$'s for a given pentaquark. What we found was that, in all cases, compact baryon+meson configurations had lower masses than their more symmetric counterparts. Moreover, those "molecular" associations were unstable (with larger masses) that their infinitely separated non-interacting pieces except int the case of $(ccb)(c\bar{b})$ and $(bbc)(b\bar{c})$ pentaquarks. In any case, even for unstable arrangements, the excess masses are so small ($\sim$20-60 MeV) as to make those compact molecular structures serious candidates to be experimentally detected.