This study addresses the central question that how can a precise measurement of the gluon Parton Distribution Function (PDF) be attained, thereby enhancing our understanding of the intricate quark-gluon composition within nucleons? Through a dual approach involving Quantum Chromodynamics (QCD) and high-energy particle collisions, the research combines computational tools and experimental techniques, focusing on protons and neutrons. Scattering experiments and advanced detectors are utilized to probe nucleon structure. Theoretical frameworks like QCD equations and the Momentum Sum Rule investigate Parton Distribution Functions at specific momentum fractions and energy scales. The primary goal is to achieve an accurate measurement of the gluon PDF, crucial for unraveling strong force complexities and internal nucleon dynamics. The study's theoretical plots and equations illustrate hypothetical behaviors of gluon distribution, quark distribution, and splitting functions, contributing to a deeper understanding of fundamental particles.

J/ψ, a charmonium bound state made of a charm and an anti-charm quark, has been discovered in the 1970s and confirmed the quark model. Because the mass of charms quarks is significantly above the Quantum Chromodynamics (QCD) scale ΛQCD, charmonia are considered as excellent probes to test perturbative Quantum Chromodynamics (pQCD) calculations. Over the past decades, they have been studied extensively at different high energy colliders. However, their production mechanisms, which involves multiple scales, are still not very well understood. Recently, in high multiplicity p+p collisions at RHIC and at the LHC, a significant enhancement of J/ψ production yield has been observed, which suggests a strong contribution from the Multi-Parton Interaction (MPI). This is different from the traditional pQCD picture where charm quark pairs are produced from a single hard scattering between partons in p+p collisions. In this work, we will report the J/ψ normalized production yield as a function of normalized charged particle multiplicity over a board range of rapidity and event multiplicity in the J/ψ→μ+μ− channel with PHENIX Run 15 p+p data at \sqrt{s}=200 GeV. The results are compared with PYTHIA 8 simulations with the MPI option on and off. Finally, the outlooks of J/ψ in p+Au and Au+p collisions along with Color Glass Condensate (CGC) predictions and of multiplicity dependent ψ(2S)/J/ψ ratio in p+p data will be briefly discussed.