Francesco Giovanni Celiberto
We present the fragmentation functions for the production of all-heavy (fully heavy) S -wave tetraquarks ( T 4 Q ) with scalar ( 0 + + ), axial-vector ( 1 + − ), and tensor ( 2 + + ) quantum numbers in high-energy hadronic collisions. This work extends the previous framework by incorporating nonconstituent heavy-quark contributions and introducing a replica-based uncertainty-quantification strategy derived from multiscale variations (MHOUs). The construction follows a nonrelativistic quantum chromodynamics (QCD) factorization approach, combining gluon- and heavy-quark-initiated fragmentation channels at leading power. Initial-scale inputs are modeled through updated potential-inspired wave functions, while the subsequent Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution is performed via the threshold-aware scheme. A comprehensive systematic analysis of uncertainties is carried out, with contributions from color-composite long-distance matrix elements (LDMEs) and perturbative multiscale inputs. The resulting grids, publicly released in format, provide the first complete phenomenological set for all-heavy exotics, enabling precise studies of all-charm tetraquark production and jet-associated observables within the environment. This article completes the high-energy resummation-driven generation of the program and establishes a definitive baseline for future collider-oriented analyses of all-heavy multiquark dynamics.