Parton distribution and fragmentation functions with massive gluons

The correct description of the hadron’s structure requires understanding how quarks and gluons form the observable hadrons and how they are distributed within them. Two key nonperturbative quantities encapsulate this information: parton distribution functions (PDFs) and fragmentation functions (FFs). The former define a probabilistic light-front momentum distribution of partons within a hadron, whereas the latter describe the hadronization process of high-energy partons. Computing these functions analytically poses significant challenges, as it demands models that accurately incorporate the nonperturbative infrared dynamics of quantum chromodynamics (QCD). In this work, we compute the pion PDF and its elementary and full FFs using the Curci-Ferrari (CF) model. This model enables the exploration of nonperturbative QCD effects by introducing a gluon-mass scale within the Landau gauge QCD Lagrangian. The two-point quark and gluon correlation functions derived from the CF model agree well with lattice-QCD results and reproduce the pion decay constant in the chiral limit, consistent with chiral perturbation theory. The resulting pion PDF and FFs computed with the CF quark propagator and pion Bethe-Salpeter amplitude are in good qualitative and quantitative agreement with those obtained using the Qin-Chang model, a benchmark approach to nonperturbative QCD. These findings support the broader applicability of the Curci-Ferrari model in hadron phenomenology.