Full bibliography

We report measurements of ϒ(1S), ϒ(2S) and ϒ(3S) production in p + p collisions at √s = 500 GeV ffiffi by the STAR experiment in year 2011, corresponding to an integrated luminosity Lint = 13 pb−1. The results provide precise cross sections, transverse momentum (pT) and rapidity (y) spectra, as well as cross section ratios for pT < 10 GeV=c and |y| < 1. The dependence of the ϒ yield on charged particle multiplicity has also been measured, offering new insights into the mechanisms of quarkonium production. The data are compared to various theoretical models: the color evaporation model (CEM) accurately describes the ϒ(1S) production, while the color glass condensate + nonrelativistic quantum chromodynamics (CGC + NRQCD) model overestimates the data, particularly at low pT. Conversely, the color singlet model (CSM) underestimates the rapidity dependence. These discrepancies highlight the need for further development in understanding the production dynamics of heavy quarkonia in high-energy hadronic collisions. The trend in the multiplicity dependence is consistent with CGC/saturation and string percolation models or ϒ production happening in multiple parton interactions modeled by PYTHIA8. © 2025 American Physical Society