Production of loosely bound hadron molecules from bottomonium decays

Abstract

We present multiple results on the production of loosely bound molecules in bottomonium annihilations and ${e}^{+}{e}^{\ensuremath{-}}$ collisions at $\sqrt{s}=10.58\text{ }\text{ }\mathrm{GeV}$. We perform the first comprehensive test of several models for deuteron production against all the existing data in this energy region. We fit the free parameters of the models to reproduce the observed cross sections, and we predict the deuteron spectrum and production and the cross section for the ${e}^{+}{e}^{\ensuremath{-}}\ensuremath{\rightarrow}d\overline{d}+X$ process both at the $\mathrm{\ensuremath{\Upsilon}}(1,2,3S)$ resonances and at $\sqrt{s}=10.58\text{ }\text{ }\mathrm{GeV}$. The predicted spectra show differences but are all compatible with the uncertainties of the existing data. These differences could be addressed if larger datasets are collected by the Belle II experiment. Fixing the source size parameter to reproduce the deuteron data, we then predict the production rates for $H$ dibaryon and hypertriton in this energy region using a simple coalescence model. Our prediction on the $H$ dibaryon production rate is below the limits set by the direct search at the Belle experiment, but in the range accessible to the Belle II experiment. The systematic effect due to the MC modeling of quarks and gluon fragmentation into baryons is reduced, deriving a new tuning of the pythia 8 Monte Carlo generator using the available measurement of single- and double-particle spectra in $\mathrm{\ensuremath{\Upsilon}}$ decays.