PHENIX Results on $J/\psi$ production in small systems

<!--HTML-->The production of heavy quarkonia states in p+A collisions is sensitive to a range of initial- and final-state effects that must be quantified in order to fully understand the dynamics of heavy quarks in the QCD medium. One way to disentangle the various effects is by measuring charmonia production in different rapidity and momentum regions, using various beam species. Here we present finalized PHENIX results on $J/\psi$ modification from the RHIC system-size scan, which include p+Al, p+Au, and $^{3}$He+Au measurements at forward and backward rapidity. Comparisons with state-of-the-art theory calculations and LHC data will be discussed.