Several experiments at the electron-proton ($ep$) collider HERA have collected high precision data on the spectrum of leading-proton and leading-neutron carrying a large fraction of the proton's energy. In this paper, we have analyzed recent experimental data on the production of leading-nucleon in deep inelastic scattering (DIS) processes at HERA in the framework of a perturbative QCD (pQCD). An approach based on the fractures functions framework has been used, and the nucleon fracture functions (nucleon FFs) ${\cal M}_2^{(n/p)} (x, Q^2, x_L)$ have been extracted from global QCD analysis of DIS data measured by ZEUS collaboration at HERA. We show that the approach of fracture functions formalisem allows us phenomenologically parametrize the nucleon FFs at the input scale, $Q_0^2$.
Considering leading-nucleon production data in the DIS processes, we present the results for the separate parton distributions for all parton species.
The extracted results from the $t$-integrated leading-baryon fracture functions, $F_2^{\rm LB(3)} (x, Q^2, x_L)$ are in good agreement with all DIS data analyzed, for a wide range of longitudinal momentum fraction $x_L$ as well as scaled fractional momentum variable $x$.