Multimodal fragmentation of all-heavy pentaquarks: Uncertainty-aware predictions for hadron colliders

We present an uncertainty-aware description of leading-power fragmentation for all-charm pentaquark states ( S -wave | c c c c c ¯ ⟩ ) at hadron colliders. We construct a multimodal set of collinear fragmentation functions, , incorporating both perturbative and nonperturbative uncertainties. Perturbative effects are estimated via missing higher-order variations, while the nonperturbative wave function is modeled through controlled modifications of its transverse-momentum structure, consistently combined within a replicalike framework. The initial-scale input for constituent charm fragmentation is refined to describe both compact multiquark and diquark-driven production mechanisms. We employ the (sym) interface to study NLL / NLO + (next-to-leading logarithmic/next-to-leading order) semi-inclusive pentaquark-plus-jet production at the HL-LHC and future Future Circular Collider. The bottom sector is left to future dedicated studies due to its enhanced sensitivity to nonperturbative modeling. Our results provide a flexible framework for uncertainty-controlled predictions, bridging exotic-hadron structure, heavy-flavor fragmentation, and high-energy QCD.