11institutetext: Department of Physics Education, Farhangian University, P.O.Box 14665-889, Tehran, Iran. 22institutetext: School of Particles and Accelerators, Institute for Research in Fundamental Sciences (IPM), P.O.Box 19395-5531, Tehran, Iran 33institutetext: Department of Physics, Faculty of Nano and Bio Science and Technology, Persian Gulf University, 75169 Bushehr, Iran. 44institutetext: AGH University, Faculty of Physics and Applied Computer Science, Al. Mickiewicza 30, 30-055 Kraków, Poland. 55institutetext: Department of Physics, University of Science and Technology of Mazandaran, P.O.Box 48518-78195, Behshahr, Iran.

AAK24: Global QCD analysis on polarized parton distribution in the presence of A2subscript𝐴2A_{2}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT asymmetry measurements

Fatemeh Arbabifar [email protected]11    Shahin Atashbar Tehrani [email protected]2233    Hamzeh Khanpour [email protected]442255
(Received: April 25, 2024 / Revised version: August 29, 2024)
Abstract

This article introduces AAK24, a Next-to-Leading Order (NLO) QCD analysis of polarized data from both polarized Deep Inelastic Scattering (DIS) and Semi-Inclusive Deep Inelastic Scattering (SIDIS) experiments on the nucleon. The AAK24 QCD analysis incorporates SU(2) and SU(3) symmetry breaking, specifically δu¯δd¯δs¯𝛿¯𝑢𝛿¯𝑑𝛿¯𝑠\delta\bar{u}\neq\delta\bar{d}\neq\delta\bar{s}italic_δ over¯ start_ARG italic_u end_ARG ≠ italic_δ over¯ start_ARG italic_d end_ARG ≠ italic_δ over¯ start_ARG italic_s end_ARG, while assuming δs¯𝛿¯𝑠\delta\bar{s}italic_δ over¯ start_ARG italic_s end_ARG and δs𝛿𝑠\delta sitalic_δ italic_s are equal. Emphasizing the significance of the semi-inclusive data, the study explores the determination of polarized sea quark distributions. Recent experimental data from JLAB17, COMPASS16, and COMPASS17, including the A2subscript𝐴2A_{2}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT asymmetry measurements along with SIDIS observables, are thoroughly examined for their impact on the central values of polarized PDFs, their uncertainties, and overall fit quality. Additionally, we include the nonperturbative target mass corrections (TMC) as well as higher-twist terms (HT) which are particularly important. In this work, the uncertainties are quantified using the standard Hessian method. The main results and findings of the AAK24 QCD analysis show overall good agreement with the analyzed experimental data, aligning well with other polarized PDF determinations, particularly DSSV14, LSS10, JAM17, and AKS14, all considering SU(2) and SU(3) symmetry breaking.

pacs:
PACS-keydiscribing text of that key and PACS-keydiscribing text of that key

1 Introduction

In recent years, advancements in the determination of the nucleon partonic distributions and spin projections, particularly from high-energy experimental data, have significantly enhanced our understanding in this area Deur:2018roz ; Aidala:2012mv ; Ethier:2020way ; AbdulKhalek:2021gbh ; Anderle:2021dpv ; COMPASS:2015mhb . Extracting the polarized and unpolarized partonic distributions from QCD analysis is pivotal for studying the phenomenology of hard scattering processes Aidala:2012mv . The determination of polarized parton distribution functions (polarized PDFs) and along their uncertainties has rapidly expanded, driven by polarized deep inelastic scattering (DIS) experiments at CERN, SLAC, DESY, and JLAB DeFlorian:2019xxt ; Khan:2022vot ; Adamiak:2021ppq ; Khanpour:2017cha ; Arbabifar:2023hok ; Taghavi-Shahri:2016idw ; Hirai:2008aj ; Blumlein:2010rn ; Sato:2016tuz ; Leader:2014uua ; deFlorian:2009vb , along with the inclusion of semi-inclusive deep inelastic scattering (SIDIS) data by various phenomenological groups Adamiak:2023yhz ; deFlorian:2014yva ; Arbabifar:2013tma ; Leader:2010rb .

In this work, we utilize comprehensive proton and deuteron SIDIS asymmetry data sets from SMC and COMPASS group at CERN, including the semi-inclusive asymmetries for charged pion (π±superscript𝜋plus-or-minus\pi^{\pm}italic_π start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT), kaon (k±superscript𝑘plus-or-minusk^{\pm}italic_k start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT) and light charged hadron (h±superscriptplus-or-minush^{\pm}italic_h start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT) production from a polarized proton target. Using the SIDIS asymmetry data sets, we take into account the flavor SU(2) and SU(3) symmetry breaking as well. The recent experimental data from JLAB17 CLAS:2017qga , COMPASS16 COMPASS:2016jwv , and COMPASS17 COMPASS:2015mhb experiments have been incorporated into our data samples. More importantly, this analysis marks the first simultaneous incorporation of SIDIS data and the A2subscript𝐴2A_{2}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT asymmetry measurements for proton, neutron and deuteron from SLAC/E143, SMC, HERMES, and Jlab-Hall A experiments as well. We have thoroughly examined the impact of these new data sets on both the central values and extracted uncertainties, as well as their effect on the overall fit quality. Additionally, we calculate the uncertainty of polarized PDFs using the standard Hessian method.

This paper is organized as follows: Sec. 2 outlines the basic theoretical formalism utilized in the analysis of polarized PDFs within the context of AAK24. Sec. 3 discusses the details of the polarized DIS and SIDIS datasets for proton, neutron, and deuteron targets. In Sec. 4, we present and discuss the AAK24 parametrization of polarized PDFs and the assumptions that have considered. Sec. 5 introduces the AAK24 methodology employed for extracting the polarized PDFs, the fitting procedure, and estimating their uncertainties, specifically through the standard Hessian method. Sec. 6 presents the full results and findings of this work, comparisons with other models, and discussion of data/theory agreement as well. We also discuss the effects arising from the inclusion of the new data sets added to this analysis. Finally, Sec. 7 provides a comprehensive summary and conclusion of the work.

2 Theoretical framework

In the following, we present a brief review of the basic theoretical formalism employed in the analysis of polarized PDFs within the context of AAK24. Our discussion in this section primarily focuses on two key aspects: the spin-dependent structure function within the framework of QCD, and the unpolarized structure function as well. The inclusion of the unpolarized structure function is essential as it plays a pivotal role in extracting the polarized PDFs from the measured asymmetries. Furthermore, we also present and discuss the observables used in this work to extract the AAK24 polarized PDFs.

In perturbative QCD, the spin structure function g1subscript𝑔1g_{1}italic_g start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT for Q2Λ2much-greater-thansuperscript𝑄2superscriptΛ2Q^{2}\gg\Lambda^{2}italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ≫ roman_Λ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT (without specifying the nucleon target label N𝑁Nitalic_N) can be expressed as:

g1(x,Q2)=g1τ2+TMCs(x,Q2)+g1(x,Q2)HT,subscript𝑔1𝑥superscript𝑄2superscriptsubscript𝑔1𝜏2TMCs𝑥superscript𝑄2subscript𝑔1subscript𝑥superscript𝑄2HTg_{1}(x,Q^{2})=g_{1}^{\rm\tau 2+TMCs}(x,Q^{2})+g_{1}(x,Q^{2})_{\mathrm{HT}}~{},italic_g start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) = italic_g start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_τ 2 + roman_TMCs end_POSTSUPERSCRIPT ( italic_x , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) + italic_g start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUBSCRIPT roman_HT end_POSTSUBSCRIPT , (1)

In this equation, g1τ2+TMCs(x,Q2)superscriptsubscript𝑔1𝜏2TMCs𝑥superscript𝑄2g_{1}^{\rm\tau 2+TMCs}(x,Q^{2})italic_g start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_τ 2 + roman_TMCs end_POSTSUPERSCRIPT ( italic_x , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) represents the leading twist (τ=2𝜏2\tau=2italic_τ = 2) contribution to the g1subscript𝑔1g_{1}italic_g start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT structure function, including the effects of target mass corrections (TMCs). The term g1(x,Q2)HTsubscript𝑔1subscript𝑥superscript𝑄2HTg_{1}(x,Q^{2})_{\mathrm{HT}}italic_g start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUBSCRIPT roman_HT end_POSTSUBSCRIPT accounts for the contributions from higher-twist QCD operators as well.

The g1(x,Q2)τ2subscript𝑔1superscript𝑥superscript𝑄2𝜏2g_{1}(x,Q^{2})^{\tau 2}italic_g start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT italic_τ 2 end_POSTSUPERSCRIPT is the well-known (logarithmic in Q2superscript𝑄2Q^{2}italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT) NLO perturbative QCD contribution, and can be written as,

g1(x,Q2)τ2subscript𝑔1superscript𝑥superscript𝑄2𝜏2\displaystyle g_{1}(x,Q^{2})^{\tau 2}italic_g start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT italic_τ 2 end_POSTSUPERSCRIPT =\displaystyle== 12q,q¯nfeq2{[1+αs2πδCq]δq(x,Q2)\displaystyle\frac{1}{2}\sum_{q,\bar{q}}^{n_{f}}e_{q}^{2}\left\{\left[1+\frac{% \alpha_{s}}{2\pi}\delta{C}_{q}\right]\right.\otimes\delta q(x,Q^{2})divide start_ARG 1 end_ARG start_ARG 2 end_ARG ∑ start_POSTSUBSCRIPT italic_q , over¯ start_ARG italic_q end_ARG end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_n start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT end_POSTSUPERSCRIPT italic_e start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT { [ 1 + divide start_ARG italic_α start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT end_ARG start_ARG 2 italic_π end_ARG italic_δ italic_C start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT ] ⊗ italic_δ italic_q ( italic_x , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) (2)
+αs2π 2δCgδg(x,Q2)}.\displaystyle+\left.\frac{\alpha_{s}}{2\pi}\>2\delta C_{g}\otimes\delta g(x,Q^% {2})\right\}\,.+ divide start_ARG italic_α start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT end_ARG start_ARG 2 italic_π end_ARG 2 italic_δ italic_C start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT ⊗ italic_δ italic_g ( italic_x , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) } .

In equation above, αssubscript𝛼𝑠\alpha_{s}italic_α start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT represents the strong coupling constant at NLO accuracy, and δq𝛿𝑞\delta qitalic_δ italic_q, δq¯𝛿¯𝑞\delta\bar{q}italic_δ over¯ start_ARG italic_q end_ARG, and δg𝛿𝑔\delta gitalic_δ italic_g correspond to the polarized quark, antiquark, and gluon density functions, respectively. In Eq. 2, the symbol tensor-product\otimes represents the convolution integral. The terms δCq𝛿subscript𝐶𝑞\delta C_{q}italic_δ italic_C start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT and δCg𝛿subscript𝐶𝑔\delta C_{g}italic_δ italic_C start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT denote the corresponding hard scattering Wilson coefficients, which are calculable in perturbative QCD and can be found, for example, in Refs. Lampe:1998eu ; Arbabifar:2013tma . The explicit twist-2 expression for the g1subscript𝑔1g_{1}italic_g start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT in the presence of TMCs can be written as Dong:2006jm ; Dong:2008zg ; Dong:2007iv ; Dong:2007zzc ; Sidorov:2006fi :

g1τ2+TMCs(x,Q2)superscriptsubscript𝑔1𝜏2TMCs𝑥superscript𝑄2\displaystyle g_{1}^{\tau 2+\mathrm{TMCs}}(x,Q^{2})italic_g start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_τ 2 + roman_TMCs end_POSTSUPERSCRIPT ( italic_x , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) (3)
=\displaystyle== xg1τ2(ξ,Q2;M=0)ξ(1+4M2x2/Q2)3/2+4M2x2Q2(x+ξ)ξ(1+4M2x2/Q2)2ξ1dξξg1τ2(ξ,Q2;M=0)𝑥superscriptsubscript𝑔1𝜏2𝜉superscript𝑄2M0𝜉superscript14superscriptM2superscript𝑥2superscript𝑄2324superscriptM2superscript𝑥2superscript𝑄2𝑥𝜉𝜉superscript14superscriptM2superscript𝑥2superscript𝑄22superscriptsubscript𝜉1𝑑superscript𝜉superscript𝜉superscriptsubscript𝑔1𝜏2superscript𝜉superscript𝑄2M0\displaystyle\frac{xg_{1}^{\tau 2}(\xi,Q^{2};{\mathrm{M}}=0)}{\xi(1+4{\mathrm{% M}}^{2}x^{2}/Q^{2})^{3/2}}+\frac{4{\mathrm{M}}^{2}x^{2}}{Q^{2}}\frac{(x+\xi)}{% \xi(1+4{\mathrm{M}}^{2}x^{2}/Q^{2})^{2}}\int_{\xi}^{1}\frac{d\xi^{\prime}}{\xi% ^{\prime}}g_{1}^{\tau 2}(\xi^{\prime},Q^{2};{\mathrm{M}}=0)divide start_ARG italic_x italic_g start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_τ 2 end_POSTSUPERSCRIPT ( italic_ξ , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ; roman_M = 0 ) end_ARG start_ARG italic_ξ ( 1 + 4 roman_M start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT 3 / 2 end_POSTSUPERSCRIPT end_ARG + divide start_ARG 4 roman_M start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG divide start_ARG ( italic_x + italic_ξ ) end_ARG start_ARG italic_ξ ( 1 + 4 roman_M start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ∫ start_POSTSUBSCRIPT italic_ξ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT divide start_ARG italic_d italic_ξ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_ARG start_ARG italic_ξ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_ARG italic_g start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_τ 2 end_POSTSUPERSCRIPT ( italic_ξ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ; roman_M = 0 )
\displaystyle-- 4M2x2Q2(24M2x2/Q2)2(1+4M2x2/Q2)5/2×ξ1dξξξ1dξ′′ξ′′g1τ2(ξ′′,Q2;M=0),4superscriptM2superscript𝑥2superscript𝑄224superscriptM2superscript𝑥2superscript𝑄22superscript14superscriptM2superscript𝑥2superscript𝑄252superscriptsubscript𝜉1𝑑superscript𝜉superscript𝜉superscriptsubscriptsuperscript𝜉1𝑑superscript𝜉′′superscript𝜉′′superscriptsubscript𝑔1𝜏2superscript𝜉′′superscript𝑄2M0\displaystyle\frac{4{\mathrm{M}}^{2}x^{2}}{Q^{2}}\frac{(2-4{\mathrm{M}}^{2}x^{% 2}/Q^{2})}{2(1+4{\mathrm{M}}^{2}x^{2}/Q^{2})^{5/2}}\times\int_{\xi}^{1}\frac{d% \xi^{\prime}}{\xi^{\prime}}\int_{\xi^{\prime}}^{1}\frac{d\xi^{\prime\prime}}{% \xi^{\prime\prime}}g_{1}^{\tau 2}(\xi^{\prime\prime},Q^{2};{\mathrm{M}}=0)\,,divide start_ARG 4 roman_M start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG divide start_ARG ( 2 - 4 roman_M start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG start_ARG 2 ( 1 + 4 roman_M start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT 5 / 2 end_POSTSUPERSCRIPT end_ARG × ∫ start_POSTSUBSCRIPT italic_ξ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT divide start_ARG italic_d italic_ξ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_ARG start_ARG italic_ξ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_ARG ∫ start_POSTSUBSCRIPT italic_ξ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT divide start_ARG italic_d italic_ξ start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT end_ARG start_ARG italic_ξ start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT end_ARG italic_g start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_τ 2 end_POSTSUPERSCRIPT ( italic_ξ start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ; roman_M = 0 ) ,

in which, MM{\mathrm{M}}roman_M denotes the nucleon mass. The parameter ξ𝜉\xiitalic_ξ in the above equation refers to the Nachtmann variable Nachtmann:1973mr , which is defined as:

ξ=2x1+1+4M2x2/Q2.𝜉2𝑥114superscriptM2superscript𝑥2superscript𝑄2\xi=\frac{2x}{1+\sqrt{1+4{\mathrm{M}}^{2}x^{2}/Q^{2}}}~{}.italic_ξ = divide start_ARG 2 italic_x end_ARG start_ARG 1 + square-root start_ARG 1 + 4 roman_M start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG . (4)

The maximum kinematic value of the Nachtmann variable ξ𝜉\xiitalic_ξ presented above is less than unity, implying that both the polarized and unpolarized target mass-corrected leading-twist structure functions do not vanish at x=1𝑥1x=1italic_x = 1. Finally, the dynamical higher-twist effects are given by:

g1(x,Q2)HT=h(x,Q2)Q2+𝒪(Λ4Q4).subscript𝑔1subscript𝑥superscript𝑄2HT𝑥superscript𝑄2superscript𝑄2𝒪superscriptΛ4superscript𝑄4g_{1}(x,Q^{2})_{\mathrm{HT}}=\frac{h(x,Q^{2})}{Q^{2}}+{\cal O}\left(\frac{% \Lambda^{4}}{Q^{4}}\right)~{}\,.italic_g start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUBSCRIPT roman_HT end_POSTSUBSCRIPT = divide start_ARG italic_h ( italic_x , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG start_ARG italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG + caligraphic_O ( divide start_ARG roman_Λ start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG start_ARG italic_Q start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG ) . (5)

This contribution is fitted to the data along with the fit parameters.

Given that AAK24 is limited to the NLO accuracy, consistent with most global polarized PDF analyses available in literature, we utilize the hard scattering coefficients computed to the same NLO order. In addition, nfsubscript𝑛𝑓n_{f}italic_n start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT is the number of active flavors (nf=3subscript𝑛𝑓3n_{f}=3italic_n start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT = 3 in our analysis).

In the leading-twist approximation and at the NLO accuracy, the g2(x,Q2)subscript𝑔2𝑥superscript𝑄2g_{2}(x,Q^{2})italic_g start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_x , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) spin-dependent structure-function which contains the transverse spin carried by the nucleons is given in terms of the twist-2 component of g1(x,Q2)subscript𝑔1𝑥superscript𝑄2g_{1}(x,Q^{2})italic_g start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) via the Wandzura-Wilczek relation, as detailed in Refs. Piccione:1997zh . The g2(x,Q2)subscript𝑔2𝑥superscript𝑄2g_{2}(x,Q^{2})italic_g start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_x , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) can be written as,

g2(x,Q2)=g1(x,Q2)+x1dyyg1(y,Q2).subscript𝑔2𝑥superscript𝑄2subscript𝑔1𝑥superscript𝑄2superscriptsubscript𝑥1𝑑𝑦𝑦subscript𝑔1𝑦superscript𝑄2g_{2}({x},Q^{2})=-g_{1}({x},Q^{2})+\int_{x}^{1}\frac{dy}{y}\,g_{1}({y},Q^{2})\,.italic_g start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_x , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) = - italic_g start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) + ∫ start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT divide start_ARG italic_d italic_y end_ARG start_ARG italic_y end_ARG italic_g start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_y , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) . (6)

It is important to note that these leading-twist approximation results are, strictly speaking, valid in the Bjorken limit Q2superscript𝑄2Q^{2}\to\inftyitalic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT → ∞ for a given fixed value of x𝑥xitalic_x. At the finite values of photon virtuaity Q2superscript𝑄2Q^{2}italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT the power-suppressed corrections to the spin-dependent structure function can make important contributions, particularly in specific kinematic regions. The simplest of these are the target mass corrections (TMCs) in which studied in details in our prior analysis Khanpour:2017cha , and it has also been examined in other studies found in the literature, such as Sato:2016tuz ; Courtoy:2022kca ; Mirjalili:2022cal . In the following analysis, we consider the TMC and HT corrections as presented in Eqs. 3 and 5.

The unpolarized structure function F1(x,Q2)subscript𝐹1𝑥superscript𝑄2F_{1}(x,Q^{2})italic_F start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ), which is essential for extracting the polarized PDFs from the measured asymmetries, can be expressed as follows:

F1(x,Q2)=(1+γ2)2x(1+R(x,Q2))F2(x,Q2),subscript𝐹1𝑥superscript𝑄21superscript𝛾22𝑥1𝑅𝑥superscript𝑄2subscript𝐹2𝑥superscript𝑄2F_{1}({x},Q^{2})=\frac{({1}+\gamma^{2})}{2x(1+R({x},Q^{2}))}F_{2}({x},Q^{2})\,,italic_F start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) = divide start_ARG ( 1 + italic_γ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG start_ARG 2 italic_x ( 1 + italic_R ( italic_x , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) ) end_ARG italic_F start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_x , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) , (7)

in which is determined per kinetic variables x𝑥xitalic_x, R𝑅Ritalic_R and the nucleon unpolarized structure function F2(x,Q2)subscript𝐹2𝑥superscript𝑄2F_{2}(x,Q^{2})italic_F start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_x , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ). In Eq. 7, the factor γ2superscript𝛾2\gamma^{2}italic_γ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT defined as 4M2x2Q24superscript𝑀2superscript𝑥2superscript𝑄2\frac{4M^{2}x^{2}}{Q^{2}}divide start_ARG 4 italic_M start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG, and R𝑅Ritalic_R is the ratio of the longitudinal σLsubscript𝜎𝐿\sigma_{L}italic_σ start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT to the transverse σTsubscript𝜎𝑇\sigma_{T}italic_σ start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT virtual photoproduction cross sections; R=σLσT𝑅subscript𝜎𝐿subscript𝜎𝑇R=\frac{\sigma_{L}}{\sigma_{T}}italic_R = divide start_ARG italic_σ start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT end_ARG start_ARG italic_σ start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT end_ARG Whitlow:1990gk , which in our calculations, we use the experimental values measurement by NMC collaboration NewMuon:1995zeb .

The virtual photoproduction longitudinal and transverse asymmetries in the Bjorken limit can be expressed as ratios of the spin-dependent (g1(x,Q2)subscript𝑔1𝑥superscript𝑄2g_{1}({x},Q^{2})italic_g start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) and g2(x,Q2)subscript𝑔2𝑥superscript𝑄2g_{2}({x},Q^{2})italic_g start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_x , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT )) and unpolarized (F1(x,Q2)subscript𝐹1𝑥superscript𝑄2F_{1}({x},Q^{2})italic_F start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) and F2(x,Q2)subscript𝐹2𝑥superscript𝑄2F_{2}({x},Q^{2})italic_F start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_x , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT )) structure functions.

Finally, having at hand the spin-dependent and unpolarized structure functions, one can define the A1subscript𝐴1A_{1}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and A2subscript𝐴2A_{2}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT asymmetries for polarized DIS processes as follows:

A1(x,Q2)=g1(x,Q2)F1(x,Q2)(1+γ2),subscript𝐴1𝑥superscript𝑄2subscript𝑔1𝑥superscript𝑄2subscript𝐹1𝑥superscript𝑄21superscript𝛾2A_{1}({x},Q^{2})=\frac{g_{1}({x},Q^{2})}{F_{1}({x},Q^{2})}({1}+\gamma^{2})\,,italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) = divide start_ARG italic_g start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG start_ARG italic_F start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG ( 1 + italic_γ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) , (8)

for A1subscript𝐴1A_{1}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT, and

A2(x,Q2)=γg1(x,Q2)+g2(x,Q2)F1(x,Q2),subscript𝐴2𝑥superscript𝑄2𝛾subscript𝑔1𝑥superscript𝑄2subscript𝑔2𝑥superscript𝑄2subscript𝐹1𝑥superscript𝑄2A_{2}({x},Q^{2})=\gamma\frac{g_{1}({x},Q^{2})+g_{2}(x,Q^{2})}{F_{1}({x},Q^{2})% }\,,italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_x , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) = italic_γ divide start_ARG italic_g start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) + italic_g start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_x , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG start_ARG italic_F start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG , (9)

for the A2subscript𝐴2A_{2}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT asymmetries, respectively.

We should stress here that, a very large data sample of the A2subscript𝐴2A_{2}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT asymmetry has been incorporated into the AAK24 QCD analysis. We will provide a comprehensive discussion of this dataset in the forthcoming section. The effect of the A2subscript𝐴2A_{2}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT asymmetry data on the extracted polarized PDFs are studied and clearly discussed in Sec. 6 as well.

As previously emphasized, one of the primary objectives of this study is to incorporate the SIDIS data alongside the polarized DIS data, thereby enabling the examination of light sea-quark decomposition and to provide well-constraint gluon PDFs. In order to increase the statistics of the SIDIS observables, we utilize the most recent and precise experimental data measured by JLAB17, COMPASS16, and COMPASS17 experiments.

At the leading-twist approximation and at NLO accuracy, the spin-dependent structure function g1Nh(x,z,Q2)superscriptsubscript𝑔1𝑁𝑥𝑧superscript𝑄2g_{1N}^{h}({x},{z},Q^{2})italic_g start_POSTSUBSCRIPT 1 italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( italic_x , italic_z , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) for the SIDIS processes can be calculated in terms of the polarized PDFs and the Fragmentation Function (FFs). It reads,

g1Nh,twist2(x,z,Q2)superscriptsubscript𝑔1𝑁𝑡𝑤𝑖𝑠𝑡2𝑥𝑧superscript𝑄2\displaystyle g_{1N}^{h,{twist-2}}({x},{z},Q^{2})italic_g start_POSTSUBSCRIPT 1 italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h , italic_t italic_w italic_i italic_s italic_t - 2 end_POSTSUPERSCRIPT ( italic_x , italic_z , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) =\displaystyle== 12q,q¯nfeq2{[δq(1+αs(Q2)2πδCqq)Dqh\displaystyle\frac{1}{2}\sum_{q,\bar{q}}^{n_{f}}e_{q}^{2}\left\{\phantom{\int% \limits_{a}^{b}}\left[\delta q\left({1}+\otimes\frac{\alpha_{s}(Q^{2})}{2\pi}% \delta C_{qq}\otimes\right){D}_{q}^{h}\right.\right.divide start_ARG 1 end_ARG start_ARG 2 end_ARG ∑ start_POSTSUBSCRIPT italic_q , over¯ start_ARG italic_q end_ARG end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_n start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT end_POSTSUPERSCRIPT italic_e start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT { [ italic_δ italic_q ( 1 + ⊗ divide start_ARG italic_α start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG start_ARG 2 italic_π end_ARG italic_δ italic_C start_POSTSUBSCRIPT italic_q italic_q end_POSTSUBSCRIPT ⊗ ) italic_D start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT
+δqαs(Q2)2πδCgq(1)Dghtensor-producttensor-product𝛿𝑞subscript𝛼𝑠superscript𝑄22𝜋𝛿superscriptsubscript𝐶𝑔𝑞1superscriptsubscript𝐷𝑔\displaystyle+\delta q\otimes\frac{\alpha_{s}(Q^{2})}{{2}\pi}\delta{C}_{gq}^{(% 1)}\otimes{D}_{g}^{h}+ italic_δ italic_q ⊗ divide start_ARG italic_α start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG start_ARG 2 italic_π end_ARG italic_δ italic_C start_POSTSUBSCRIPT italic_g italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT ⊗ italic_D start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT
+δgαs(Q2)2πδCqg(1)Dqh](x,z,Q2)},\displaystyle+\left.\left.\phantom{\int}\delta g\otimes\frac{\alpha_{s}(Q^{2})% }{{2}\pi}\delta C_{qg}^{(1)}\otimes{D}_{q}^{h}\right]({x},{z},Q^{2})\phantom{% \int\limits_{a}^{b}}\right\}~{},+ italic_δ italic_g ⊗ divide start_ARG italic_α start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG start_ARG 2 italic_π end_ARG italic_δ italic_C start_POSTSUBSCRIPT italic_q italic_g end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT ⊗ italic_D start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ] ( italic_x , italic_z , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) } ,
,absent\displaystyle\,,, (10)

in which the δq𝛿𝑞\delta qitalic_δ italic_q, δq¯𝛿¯𝑞\delta\bar{q}italic_δ over¯ start_ARG italic_q end_ARG and δg𝛿𝑔\delta gitalic_δ italic_g are the polarized quark, antiquark and gluon density functions, the Cqqsubscript𝐶𝑞𝑞C_{qq}italic_C start_POSTSUBSCRIPT italic_q italic_q end_POSTSUBSCRIPT, Cqgsubscript𝐶𝑞𝑔C_{qg}italic_C start_POSTSUBSCRIPT italic_q italic_g end_POSTSUBSCRIPT and Cggsubscript𝐶𝑔𝑔C_{gg}italic_C start_POSTSUBSCRIPT italic_g italic_g end_POSTSUBSCRIPT denote the corresponding hard scattering Wilson coefficients. The Dqhsubscriptsuperscript𝐷𝑞D^{h}_{q}italic_D start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT and Dghsubscriptsuperscript𝐷𝑔D^{h}_{g}italic_D start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT are the corresponding FFs of a given hadron.

The unpolarized structure function of F1Nh(x,z,Q2)superscriptsubscript𝐹1𝑁𝑥𝑧superscript𝑄2F_{1N}^{h}({x},{z},Q^{2})italic_F start_POSTSUBSCRIPT 1 italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( italic_x , italic_z , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) for a SIDIS process is given by,

F1Nh,twist2(x,z,Q2)superscriptsubscript𝐹1𝑁𝑡𝑤𝑖𝑠𝑡2𝑥𝑧superscript𝑄2\displaystyle F_{1N}^{h,{twist-2}}({x},{z},Q^{2})italic_F start_POSTSUBSCRIPT 1 italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h , italic_t italic_w italic_i italic_s italic_t - 2 end_POSTSUPERSCRIPT ( italic_x , italic_z , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) =\displaystyle== 12q,q¯nfeq2{[q(1+αs(Q2)2πCqq)Dqh\displaystyle\frac{1}{2}\sum_{q,\bar{q}}^{n_{f}}e_{q}^{2}\left\{\phantom{\int% \limits_{a}^{b}}\left[q\left(1+\otimes\frac{\alpha_{s}(Q^{2})}{{2}\pi}C_{qq}% \otimes\right)D_{q}^{h}\right.\right.divide start_ARG 1 end_ARG start_ARG 2 end_ARG ∑ start_POSTSUBSCRIPT italic_q , over¯ start_ARG italic_q end_ARG end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_n start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT end_POSTSUPERSCRIPT italic_e start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT { [ italic_q ( 1 + ⊗ divide start_ARG italic_α start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG start_ARG 2 italic_π end_ARG italic_C start_POSTSUBSCRIPT italic_q italic_q end_POSTSUBSCRIPT ⊗ ) italic_D start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT
+qαs(Q2)2πCgq(1)Dghtensor-producttensor-product𝑞subscript𝛼𝑠superscript𝑄22𝜋superscriptsubscript𝐶𝑔𝑞1superscriptsubscript𝐷𝑔\displaystyle+q\otimes\frac{\alpha_{s}(Q^{2})}{2\pi}{C}_{gq}^{(1)}\otimes D_{g% }^{h}+ italic_q ⊗ divide start_ARG italic_α start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG start_ARG 2 italic_π end_ARG italic_C start_POSTSUBSCRIPT italic_g italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT ⊗ italic_D start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT
+gαs(Q2)2πCqg(1)Dqh](x,z,Q2)}\displaystyle+\left.\left.\phantom{\int}{g}\otimes\frac{\alpha_{s}(Q^{2})}{2% \pi}{C}_{qg}^{(1)}\otimes D_{q}^{h}\right]({x},{z},Q^{2})\phantom{\int\limits_% {a}^{b}}\right\}~{}\,+ italic_g ⊗ divide start_ARG italic_α start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG start_ARG 2 italic_π end_ARG italic_C start_POSTSUBSCRIPT italic_q italic_g end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT ⊗ italic_D start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ] ( italic_x , italic_z , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) }

where the q𝑞qitalic_q and g𝑔gitalic_g are the unpolarized PDFs, the Dqhsubscriptsuperscript𝐷𝑞D^{h}_{q}italic_D start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT and Dghsubscriptsuperscript𝐷𝑔D^{h}_{g}italic_D start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT are the quark and gluon FFs of hadron. Finlay, the Cq(g)g(q)subscript𝐶𝑞𝑔𝑔𝑞C_{q(g)g(q)}italic_C start_POSTSUBSCRIPT italic_q ( italic_g ) italic_g ( italic_q ) end_POSTSUBSCRIPT are the hard scattering Wilson coefficients available in literature deFlorian:1997zj ; Arbabifar:2013tma . In the following analysis, we use the publicly available unpolarized PDFs of CT18 group Hou:2019efy and for the quark and gluon FFs of hadron we employ the DSS DSS results.

The observable for the SIDIS process that we used in AAK24 QCD analysis is the A1Nhsubscriptsuperscript𝐴1𝑁A^{h}_{1N}italic_A start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 1 italic_N end_POSTSUBSCRIPT asymmetry which can be calculated as ratio of spin-dependent structure function g1Nh(x,z,Q2)superscriptsubscript𝑔1𝑁𝑥𝑧superscript𝑄2g_{1N}^{h}({x},{z},Q^{2})italic_g start_POSTSUBSCRIPT 1 italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( italic_x , italic_z , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) and g2Nh(x,z,Q2)superscriptsubscript𝑔2𝑁𝑥𝑧superscript𝑄2g_{2N}^{h}({x},{z},Q^{2})italic_g start_POSTSUBSCRIPT 2 italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( italic_x , italic_z , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) to the unpolarized structure function F1Nh(x,z,Q2)superscriptsubscript𝐹1𝑁𝑥𝑧superscript𝑄2F_{1N}^{h}({x},{z},Q^{2})italic_F start_POSTSUBSCRIPT 1 italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( italic_x , italic_z , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ). It is given by Leader:2010rb ,

A1Nh(x,z,Q2)=g1Nh(x,z,Q2)γ2g2Nh(x,z,Q2)F1Nh(x,z,Q2),superscriptsubscript𝐴1𝑁𝑥𝑧superscript𝑄2superscriptsubscript𝑔1𝑁𝑥𝑧superscript𝑄2superscript𝛾2superscriptsubscript𝑔2𝑁𝑥𝑧superscript𝑄2superscriptsubscript𝐹1𝑁𝑥𝑧superscript𝑄2A_{1N}^{h}(x,z,Q^{2})=\frac{g_{1N}^{h}(x,z,Q^{2})-\gamma^{2}g_{2N}^{h}(x,z,Q^{% 2})}{F_{1N}^{h}(x,z,Q^{2})},italic_A start_POSTSUBSCRIPT 1 italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( italic_x , italic_z , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) = divide start_ARG italic_g start_POSTSUBSCRIPT 1 italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( italic_x , italic_z , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) - italic_γ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT 2 italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( italic_x , italic_z , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG start_ARG italic_F start_POSTSUBSCRIPT 1 italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( italic_x , italic_z , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG , (12)

As in the case of polarized DIS, the factor γ2superscript𝛾2\gamma^{2}italic_γ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT is defined as 4M2x2Q24superscript𝑀2superscript𝑥2superscript𝑄2\frac{4M^{2}x^{2}}{Q^{2}}divide start_ARG 4 italic_M start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG. In this work, we do not include HT and TMC corrections for the case of SIDIS observables. The HT and TMC corrections are expected to be less significant due to the kinematic region and the precision of the currently available SIDIS data sets Leader:2010rb .

Having defined the polarized observables A1(x,Q2)subscript𝐴1𝑥superscript𝑄2A_{1}({x},Q^{2})italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_x , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ), A2(x,Q2)subscript𝐴2𝑥superscript𝑄2A_{2}({x},Q^{2})italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_x , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ), A1Nh(x,z,Q2)superscriptsubscript𝐴1𝑁𝑥𝑧superscript𝑄2A_{1N}^{h}({x},{z},Q^{2})italic_A start_POSTSUBSCRIPT 1 italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( italic_x , italic_z , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) and the structure functions necessary for the AAK24 QCD-based analysis, our next step in the upcoming section is to present an overview of the polarized inclusive DIS and SIDIS datasets that are subject to analysis in this study.

3 Polarized inclusive DIS and SIDIS datasets

As we discussed earlier, in AAK24 QCD analysis, we use two types of datasets from polarized DIS and SIDIS experiments on proton, neutron and deuteron targets which mostly come from relevant experiments that were done at DESY, SLAC, JLAB and CERN. Our objective is to comprehensively incorporate all available data from polarized DIS and SIDIS experiments, ensuring coverage across a wide range of kinematic variables in terms of both x𝑥xitalic_x and Q2superscript𝑄2Q^{2}italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT. In the following, we will discuss in details the datasets analyzed in our work.

From a statistical perspective, the wealth of experimental datasets used in our analysis are for the inclusive polarized DIS experiments for A1subscript𝐴1A_{1}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT asymmetry for lepton scattering from stationary targets with various combinations of targets and lepton spin.

The published experimental datasets for the A1subscript𝐴1A_{1}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and g1F1subscript𝑔1subscript𝐹1\frac{g_{1}}{F_{1}}divide start_ARG italic_g start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_F start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG observables are comprehensively presented in Table. 1. This table provides details on the measured observables, their kinematic coverage in terms of x𝑥xitalic_x and Q2superscript𝑄2Q^{2}italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT, as well as references to the publications where these datasets can be found.

These datasets have been measured and published by various experimental collaborations, including COMPASS COMPASS:2015mhb ; COMPASS:2016jwv ; COMP1 , JLAB CLAS:2017qga ; JLabn , CLAS CLA1pd , EMC EMCp , SMC SMCpd , SLAC E142n ; E143pd ; E154n ; E155p ; E155d , and HERMES HERMn ; HERMpd .

It is worth noting that, in AAK24 QCD analysis, we have analyzed the latest A1subscript𝐴1A_{1}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT asymmetry data from the COMPASS experiment at CERN, encompassing measurements on proton, neutron, and deuteron targets, COMPASS16 COMPASS:2015mhb and COMPASS17 COMPASS:2016jwv . For the case of COMPASS16 in which cover the range from 1 to 190 GeV2 in the photon virtuality and from 0.00250.00250.00250.0025 to 0.70.70.70.7 in the Bjorken scaling variable x𝑥xitalic_x, they improved the statistical precision of g1psuperscriptsubscript𝑔1𝑝g_{1}^{p}italic_g start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT by about a factor of two in the newly explored low-x𝑥xitalic_x region x<0.02𝑥0.02x<0.02italic_x < 0.02.

The COMPASS17 data were taken at 160 GeV2 beam energy and the results are presented for the kinematic range 1 to 100 GeV2 in photon virtuality, and 4×103<x<0.74superscript103𝑥0.74\times 10^{-3}<x<0.74 × 10 start_POSTSUPERSCRIPT - 3 end_POSTSUPERSCRIPT < italic_x < 0.7 in the Bjorken scaling variable, and W>4𝑊4W>4italic_W > 4 GeV2 in the mass of the hadronic final state.

As indicated in Table. 1, the recently incorporated data from the COMPASS experiment are particularly valuable as they offer measurements at the lowest accessible values of x𝑥xitalic_x, (approximately x4×103similar-to𝑥4superscript103x\sim 4\times 10^{-3}italic_x ∼ 4 × 10 start_POSTSUPERSCRIPT - 3 end_POSTSUPERSCRIPT). Consequently, these data play a crucial role in accurately determining the polarization information of gluons and sea quarks through inclusive polarized DIS data analysis.

In addition to the COMPASS data, we have incorporated the recent A1subscript𝐴1A_{1}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT asymmetry data from the JLAB experiment, focusing on measurements with the proton target, as detailed in Ref. CLAS:2017qga . This inclusion further enriches our dataset for analysis.

The CLAS experiment CLA1pd , contributes significantly to our dataset, providing nearly 60 percent of the total data points. In total, we include 996 data points for the A1subscript𝐴1A_{1}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and g1F1subscript𝑔1subscript𝐹1\frac{g_{1}}{F_{1}}divide start_ARG italic_g start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_F start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG asymmetry measurements. The E155 E155p ; E155d , HERMES HERMn and JLAB JLabn experiment present special cases of g1/F1subscript𝑔1subscript𝐹1g_{1}/F_{1}italic_g start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT / italic_F start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT.

Experiment Process Ndatasubscript𝑁dataN_{\rm data}italic_N start_POSTSUBSCRIPT roman_data end_POSTSUBSCRIPT xminsubscript𝑥minx_{\rm min}italic_x start_POSTSUBSCRIPT roman_min end_POSTSUBSCRIPT xmaxsubscript𝑥maxx_{\rm max}italic_x start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT Qmin2superscriptsubscript𝑄min2Q_{\rm min}^{{2}}italic_Q start_POSTSUBSCRIPT roman_min end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT [GeV2] Qmax2superscriptsubscript𝑄max2Q_{\rm max}^{{2}}italic_Q start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT [GeV2] Observable χ2superscript𝜒2\chi^{2}italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT χ2superscript𝜒2\chi^{2}italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT (Exc. A2subscript𝐴2A_{2}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT)
COMPASS16 COMPASS:2015mhb DIS(p)  44 0.0045 0.58 1.02 82.04 A1psuperscriptsubscript𝐴1𝑝A_{1}^{p}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT 40.188 39.361
COMPASS17 COMPASS:2016jwv DIS(d)  43 0.0045 0.569 1.03 74.1 A1dsuperscriptsubscript𝐴1𝑑A_{1}^{d}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT 34.269 33.737
JLAB17  CLAS:2017qga DIS(p)  10 0.15733 0.58355 1.0522 4.1293 A1psuperscriptsubscript𝐴1𝑝A_{1}^{p}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT 36.312 36.478
EMC EMCp DIS(p)  10 0.015 0.466 3.5 29.5 A1psuperscriptsubscript𝐴1𝑝A_{1}^{p}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT 4.100 4.686
SMC SMCpd DIS(p)  12 0.005 0.48 1.3 58 A1psuperscriptsubscript𝐴1𝑝A_{1}^{p}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT 4.315 4.526
SMC SMCpd DIS(d)  12 0.005 0.479 1.3 54.8 A1dsuperscriptsubscript𝐴1𝑑A_{1}^{d}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT 21.428 21.261
COMPASS COMP1 DIS(p)  15 0.0046 0.568 1.1 62.1 A1psuperscriptsubscript𝐴1𝑝A_{1}^{p}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT 11.601 11.186
SLAC/E142 E142n DIS(n)   8 0.035 0.466 1.1 5.5 A1nsuperscriptsubscript𝐴1𝑛A_{1}^{n}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT 7.404 4.445
SLAC/E143 E143pd DIS(p)  28 0.031 0.749 1.27 9.52 A1psuperscriptsubscript𝐴1𝑝A_{1}^{p}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT 21.111 20.757
SLAC/E143 E143pd DIS(d)  28 0.031 0.749 1.27 9.52 A1dsuperscriptsubscript𝐴1𝑑A_{1}^{d}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT 52.322 50.999
SLAC/E154 E154n DIS(n)  11 0.017 0.564 1.2 15 A1nsuperscriptsubscript𝐴1𝑛A_{1}^{n}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT 3.805 3.803
SLAC/E155 E155p DIS(p)  24 0.015 0.75 1.22 34.72 g1pF1psuperscriptsubscript𝑔1𝑝superscriptsubscript𝐹1𝑝\frac{g_{1}^{p}}{F_{1}^{p}}divide start_ARG italic_g start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT end_ARG start_ARG italic_F start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT end_ARG 29.226 34.812
SLAC/E155 E155d DIS(d)  24 0.015 0.75 1.22 34.72 g1dF1dsuperscriptsubscript𝑔1𝑑superscriptsubscript𝐹1𝑑\frac{g_{1}^{d}}{F_{1}^{d}}divide start_ARG italic_g start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT end_ARG start_ARG italic_F start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT end_ARG 21.652 21.450
HERMES HERMn DIS(n)   9 0.033 0.464 1.22 5.25 A1nsuperscriptsubscript𝐴1𝑛A_{1}^{n}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT 3.939 3.460
HERMES HERMn DIS(p)  19 0.028 0.66 1.01 7.36 g1pF1psuperscriptsubscript𝑔1𝑝superscriptsubscript𝐹1𝑝\frac{g_{1}^{p}}{F_{1}^{p}}divide start_ARG italic_g start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT end_ARG start_ARG italic_F start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT end_ARG 16.324 16.092
HERMES HERMpd DIS(p)  15 0.0264 0.7248 1.12 12.21 A1psuperscriptsubscript𝐴1𝑝A_{1}^{p}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT 10.284 9.947
HERMES HERMpd DIS(d)  15 0.0264 0.7248 1.12 12.21 A1dsuperscriptsubscript𝐴1𝑑A_{1}^{d}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT 19.408 19.569
JLab-Hall A JLabn DIS(n)   3 0.33 0.6 2.71 4.38 g1nF1nsuperscriptsubscript𝑔1𝑛superscriptsubscript𝐹1𝑛\frac{g_{1}^{n}}{F_{1}^{n}}divide start_ARG italic_g start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT end_ARG start_ARG italic_F start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT end_ARG 2.396 1.840
CLAS CLA1pd DIS(p) 151 0.1088 0.5916 1.01 4.96 A1psuperscriptsubscript𝐴1𝑝A_{1}^{p}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT 140.949 145.795
CLAS CLA1pd DIS(d) 482 0.1366 0.57 1.01 4.16 A1dsuperscriptsubscript𝐴1𝑑A_{1}^{d}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT 434.995 441.439
 TOTAL: 951 916.038 925.652
Table 1: The published experimental data sets for A1Nsuperscriptsubscript𝐴1𝑁A_{1}^{N}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_N end_POSTSUPERSCRIPT and g1NF1Nsuperscriptsubscript𝑔1𝑁superscriptsubscript𝐹1𝑁\frac{g_{1}^{N}}{F_{1}^{N}}divide start_ARG italic_g start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_N end_POSTSUPERSCRIPT end_ARG start_ARG italic_F start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_N end_POSTSUPERSCRIPT end_ARG asymmetry measurements from different experiments for different polarized DIS target. The measured observables, their kinematic range in x𝑥xitalic_x and Q2superscript𝑄2Q^{2}italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT, and the published references are also shown. The individual χ2superscript𝜒2\chi^{2}italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT values per data point for both the analysis with and without A2subscript𝐴2A_{2}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT asymmetry data are also presented. The details of the kinematical cuts applied on the data and the method on χ2superscript𝜒2\chi^{2}italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT minimization are explained in the text.

In Fig. 1, the kinematic coverage of the A1Nsuperscriptsubscript𝐴1𝑁A_{1}^{N}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_N end_POSTSUPERSCRIPT asymmetry measurement dataset is represented in the x𝑥xitalic_x and Q2superscript𝑄2Q^{2}italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT plane. The data points are categorized by different experiments, as illustrated in the presentation of the dataset. As one can see, the newly added data from COMPASS experiment offer measurements at the lowest accessible values of Bjorken variable x𝑥xitalic_x, (x4×103similar-to𝑥4superscript103x\sim 4\times 10^{-3}italic_x ∼ 4 × 10 start_POSTSUPERSCRIPT - 3 end_POSTSUPERSCRIPT).

Refer to caption
Figure 1: The kinematic coverage in the x𝑥xitalic_x and Q2superscript𝑄2Q^{2}italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT plane of the A1Nsuperscriptsubscript𝐴1𝑁A_{1}^{N}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_N end_POSTSUPERSCRIPT asymmetry measurement data set. The data points are classified by different experiments.

In the AAK24 global polarized PDF analysis, we incorporate, for the first time, all available worldwide data on the inclusive DIS of leptons on proton, neutron, and deuteron targets for the measurement of the A2subscript𝐴2A_{2}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT asymmetry measurement. This comprehensive inclusion expect to enhance the coverage and precision of the analysis.

This inclusion encompasses all datasets from various experiments, including SLAC/E143 E143:1995xmc ; E143pd ; E155:2002iec ; E155:1999eug ; E154:1997eyc ; E154:1997eyc , SMC SpinMuonSMC:1997ixm , JLab-Hall A JLabn , as well as recent high-precision asymmetry measurements from HERMES HERMES:2011xgd . The datasets for the A2subscript𝐴2A_{2}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT asymmetry measurement are presented in detail in Table. 2. This table provides the kinematic coverage in terms of x𝑥xitalic_x and Q2superscript𝑄2Q^{2}italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT, as well as references to the corresponding publications. The A2Nsuperscriptsubscript𝐴2𝑁A_{2}^{N}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_N end_POSTSUPERSCRIPT asymmetry measurements include 140 data points. To assess the impact of these data sets, two distinct analyses were conducted in this work, with and without the inclusion of A2Nsuperscriptsubscript𝐴2𝑁A_{2}^{N}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_N end_POSTSUPERSCRIPT asymmetry measurements. A detailed discussion on this matter is provided in Sec. 6.

Experiment Process Ndatasubscript𝑁dataN_{\rm data}italic_N start_POSTSUBSCRIPT roman_data end_POSTSUBSCRIPT xminsubscript𝑥minx_{\rm min}italic_x start_POSTSUBSCRIPT roman_min end_POSTSUBSCRIPT xmaxsubscript𝑥maxx_{\rm max}italic_x start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT Qmin2superscriptsubscript𝑄min2Q_{\rm min}^{2}italic_Q start_POSTSUBSCRIPT roman_min end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT [GeV2] Qmax2superscriptsubscript𝑄max2Q_{\rm max}^{2}italic_Q start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT [GeV2] Observable χ2superscript𝜒2\chi^{2}italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT
SLAC/E143 E143:1995xmc DIS(p)  12 0.038 0.595 1.49 8.85 A2psuperscriptsubscript𝐴2𝑝A_{2}^{p}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT 7.260
SLAC/E143 E143:1995xmc DIS(d)  12 0.038 0.594 1.49 8.84 A2dsuperscriptsubscript𝐴2𝑑A_{2}^{d}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT 11.441
SLAC/E143 E143pd DIS(p)  12 0.038 0.595 1.49 8.85 A2psuperscriptsubscript𝐴2𝑝A_{2}^{p}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT 7.414
SLAC/E143 E143pd DIS(d)  12 0.038 0.595 1.49 8.86 A2dsuperscriptsubscript𝐴2𝑑A_{2}^{d}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT 11.870
SLAC/E143 E143pd DIS(n)  12 0.038 0.595 1.49 8.86 A2nsuperscriptsubscript𝐴2𝑛A_{2}^{n}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT 7.414
SLAC/E155 E155:2002iec DIS(p)  8 0.038 0.78 1.1 8.4 A2psuperscriptsubscript𝐴2𝑝A_{2}^{p}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT 9.130
SLAC/E155 E155:2002iec DIS(d)  8 0.038 0.78 1.1 8.4 A2dsuperscriptsubscript𝐴2𝑑A_{2}^{d}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT 14.759
SLAC/E155 E155:1999eug DIS(p)  10 0.022 0.839 1.15 27.18 A2psuperscriptsubscript𝐴2𝑝A_{2}^{p}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT 46.832
SLAC/E155 E155:1999eug DIS(d)  10 0.022 0.839 1.15 27.18 A2dsuperscriptsubscript𝐴2𝑑A_{2}^{d}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT 21.102
SLAC/E154 E154:1997eyc DIS(n)  17 0.017 0.564 1.2 15 A2psuperscriptsubscript𝐴2𝑝A_{2}^{p}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT 0.736
SMC SpinMuonSMC:1997ixm DIS(d)  4 0.0108 0.228 2.6 18 A2dsuperscriptsubscript𝐴2𝑑A_{2}^{d}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT 6.224
HERMES HERMES:2011xgd DIS(p)  20 0.039 0.678 1.09 10.35 A2psuperscriptsubscript𝐴2𝑝A_{2}^{p}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT 20.394
JLab-Hall A JLabn DIS(n)   3 0.33 0.6 2.71 4.38 A2nsuperscriptsubscript𝐴2𝑛A_{2}^{n}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT 3.369
 TOTAL: 140 113.703
Table 2: Same as Table. 1 but this time for the A2Nsuperscriptsubscript𝐴2𝑁A_{2}^{N}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_N end_POSTSUPERSCRIPT asymmetry measurement.

In Fig. 2, the kinematic coverage of the A2subscript𝐴2A_{2}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT asymmetry measurement are presented in the x𝑥xitalic_x and Q2superscript𝑄2Q^{2}italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT plane. The data points are categorized by different experiments as well.

Refer to caption
Figure 2: The kinematic coverage in the x𝑥xitalic_x and Q2superscript𝑄2Q^{2}italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT plane for the A2Nsuperscriptsubscript𝐴2𝑁A_{2}^{N}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_N end_POSTSUPERSCRIPT asymmetry measurement in which, for the first time, analyzed in AAK24 QCD analysis. The data points are classified by different experiments.

In our analysis, we have also incorporated all available data for the SIDIS process, covering various detected hadrons, including pions (π±superscript𝜋plus-or-minus\pi^{\pm}italic_π start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT), kaons (k±superscript𝑘plus-or-minusk^{\pm}italic_k start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT), and light-charged hadrons (h±superscriptplus-or-minush^{\pm}italic_h start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT). This comprehensive addition further enhances the richness and diversity of our dataset for analysis. The published data points are related to the A1N,hsuperscriptsubscript𝐴1𝑁A_{1}^{N,h}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_N , italic_h end_POSTSUPERSCRIPT asymmetry, as described in Eq. 12, originating from the SIDIS processes. The SIDIS observables are instrumental in enabling the consideration of flavor decomposition within the sea-quark density. These observables provide valuable information about the different flavors of sea quarks within the nucleon.

The measurements of the A1N,hsuperscriptsubscript𝐴1𝑁A_{1}^{N,h}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_N , italic_h end_POSTSUPERSCRIPT asymmetry, encompassing data from various experiments, incorporate datasets from SMC SMC98 , COMPASS08 COMP1_h , and COMPASS09 COMP1_pK . Furthermore, the recently obtained and highly precise proton data from the COMPASS10 experiment compass10 , have been included in our QCD analysis.

The analyzed SIDIS data have been comprehensively summarized in Table. 3. This table provides a classification of the data based on each experiment and different detected hadrons, which include pions (π±superscript𝜋plus-or-minus\pi^{\pm}italic_π start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT) and kaons (k±superscript𝑘plus-or-minusk^{\pm}italic_k start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT), and light-charged hadrons (h±superscriptplus-or-minush^{\pm}italic_h start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT). In the AAK24 QCD fit, a total of 160 data points from the SIDIS process have been utilized.

Experiment Process Ndatasubscript𝑁dataN_{\rm data}italic_N start_POSTSUBSCRIPT roman_data end_POSTSUBSCRIPT xminsubscript𝑥minx_{\rm min}italic_x start_POSTSUBSCRIPT roman_min end_POSTSUBSCRIPT xmaxsubscript𝑥maxx_{\rm max}italic_x start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT Qmin2superscriptsubscript𝑄min2Q_{\rm min}^{{2}}italic_Q start_POSTSUBSCRIPT roman_min end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT [GeV2] Qmax2superscriptsubscript𝑄max2Q_{\rm max}^{2}italic_Q start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT [GeV2] Observable χ2superscript𝜒2\chi^{2}italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT χ2superscript𝜒2\chi^{2}italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT (Exc. A2subscript𝐴2A_{2}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT)
SMC SMC98 SIDIS(p,h+superscripth^{+}italic_h start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT)  12 0.005 0.48 10 10 A1p,h+superscriptsubscript𝐴1𝑝superscriptA_{1}^{p,h^{+}}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p , italic_h start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT 24.973 26.246
SMC SMC98 SIDIS(p,hsuperscripth^{-}italic_h start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT)  12 0.005 0.48 10 10 A1p,hsuperscriptsubscript𝐴1𝑝superscriptA_{1}^{p,h^{-}}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p , italic_h start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT 12.096 12.497
SMC SMC98 SIDIS(d,h+superscripth^{+}italic_h start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT)  12 0.005 0.48 10 10 A1d,h+superscriptsubscript𝐴1𝑑superscriptA_{1}^{d,h^{+}}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_d , italic_h start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT 6.131 6.514
SMC SMC98 SIDIS(d,hsuperscripth^{-}italic_h start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT)  12 0.005 0.48 10 10 A1d,hsuperscriptsubscript𝐴1𝑑superscriptA_{1}^{d,h^{-}}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_d , italic_h start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT 16.697 17.004
COMPASS08 COMP1_h SIDIS(d,h+superscripth^{+}italic_h start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT)  12 0.0052 0.482 1.17 60.2 A1d,h+superscriptsubscript𝐴1𝑑superscriptA_{1}^{d,h^{+}}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_d , italic_h start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT 16.108 17.906
COMPASS08 COMP1_h SIDIS(d,hsuperscripth^{-}italic_h start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT)  12 0.0052 0.482 1.17 60.2 A1d,hsuperscriptsubscript𝐴1𝑑superscriptA_{1}^{d,h^{-}}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_d , italic_h start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT 18.210 18.178
COMPASS09 COMP1_pK SIDIS(d,π+superscript𝜋\pi^{+}italic_π start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT)  10 0.0052 0.24 1.16 32.8 A1d,π+superscriptsubscript𝐴1𝑑superscript𝜋A_{1}^{d,\pi^{+}}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_d , italic_π start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT 11.598 13.578
COMPASS09 COMP1_pK SIDIS(d,πsuperscript𝜋\pi^{-}italic_π start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT)  10 0.0052 0.24 1.16 32.8 A1d,πsuperscriptsubscript𝐴1𝑑superscript𝜋A_{1}^{d,\pi^{-}}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_d , italic_π start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT 13.476 12.296
COMPASS09 COMP1_pK SIDIS(d,k+superscript𝑘k^{+}italic_k start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT)  10 0.0052 0.24 1.16 32.8 A1d,k+superscriptsubscript𝐴1𝑑superscript𝑘A_{1}^{d,k^{+}}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_d , italic_k start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT 26.377 29.621
COMPASS09 COMP1_pK SIDIS(d,ksuperscript𝑘k^{-}italic_k start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT)  10 0.0052 0.24 1.16 32.8 A1d,ksuperscriptsubscript𝐴1𝑑superscript𝑘A_{1}^{d,k^{-}}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_d , italic_k start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT 26.354 29.791
COMPASS10 compass10 SIDIS(p,π+superscript𝜋\pi^{+}italic_π start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT)  12 0.0052 0.48 1.16 55.6 A1p,π+superscriptsubscript𝐴1𝑝superscript𝜋A_{1}^{p,\pi^{+}}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p , italic_π start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT 15.294 15.853
COMPASS10 compass10 SIDIS(p,πsuperscript𝜋\pi^{-}italic_π start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT)  12 0.0052 0.48 1.16 55.6 A1p,πsuperscriptsubscript𝐴1𝑝superscript𝜋A_{1}^{p,\pi^{-}}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p , italic_π start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT 13.491 15.303
COMPASS10 compass10 SIDIS(p,k+superscript𝑘k^{+}italic_k start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT)  12 0.0052 0.48 1.16 55.6 A1p,k+superscriptsubscript𝐴1𝑝superscript𝑘A_{1}^{p,k^{+}}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p , italic_k start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT 20.450 22.077
COMPASS10 compass10 SIDIS(p,ksuperscript𝑘k^{-}italic_k start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT)  12 0.0052 0.48 1.16 55.6 A1p,ksuperscriptsubscript𝐴1𝑝superscript𝑘A_{1}^{p,k^{-}}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p , italic_k start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT 6.673 7.648
 TOTAL: 160 227.934 244.520
Table 3: The published data points for the A1N,hsuperscriptsubscript𝐴1𝑁A_{1}^{{N},{h}}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_N , italic_h end_POSTSUPERSCRIPT asymmetry measured by SMC and COMPASS collaborations on different targets. The measured observables, the number of data points, their kinematic range in x𝑥xitalic_x and Q2superscript𝑄2Q^{2}italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT, and the published references are also shown. The individual χ2superscript𝜒2\chi^{2}italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT values per data point for both the analysis with and without A2subscript𝐴2A_{2}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT asymmetry data are also presented.

The kinematic coverage of the A1N,hsuperscriptsubscript𝐴1𝑁A_{1}^{N,h}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_N , italic_h end_POSTSUPERSCRIPT asymmetry measurement are presented in Fig. 3 in the x𝑥xitalic_x and Q2superscript𝑄2Q^{2}italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT plane. The data points are categorized by different experiments as well.

Refer to caption
Figure 3: The kinematic coverage in the x𝑥xitalic_x and Q2superscript𝑄2Q^{2}italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT plane of the A1N,hsuperscriptsubscript𝐴1𝑁A_{1}^{N,h}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_N , italic_h end_POSTSUPERSCRIPT asymmetry measurement data set. The data points are classified by different experiments.

All the datasets analyzed in AAK24 polarized PDFs are subject to the cut on the Q21superscript𝑄21Q^{2}\geq 1italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ≥ 1 GeV2, in which below Q2=1superscript𝑄21Q^{2}=1italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = 1 GeV2 the perturbative QCD is not reliable. For all polarized DIS experiments where they are available, we fit directly the asymmetries of as A1subscript𝐴1A_{1}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and A2subscript𝐴2A_{2}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT.

4 AAK24 parametrization of polarized PDFs

In this section, we will present and discuss the AAK24 parametrization of polarized PDFs and the assumptions that we have considered. For the generic parametrization of polarized PDFs, assigned to each parton at a specified input scale of Q02=1GeV2superscriptsubscript𝑄021𝐺𝑒superscript𝑉2Q_{0}^{2}=1{GeV}^{2}italic_Q start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = 1 italic_G italic_e italic_V start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT, we choose the following standard functional form,

xδq=𝒩qηqxaq(1x)bq(1+cqx0.5+dqx),𝑥𝛿𝑞subscript𝒩𝑞subscript𝜂𝑞superscript𝑥subscript𝑎𝑞superscript1𝑥subscript𝑏𝑞1subscript𝑐𝑞superscript𝑥0.5subscript𝑑𝑞𝑥\displaystyle x\>\delta q={\cal N}_{q}\eta_{q}x^{a_{q}}({1}-x)^{b_{q}}(1+c_{q}% x^{0.5}+d_{q}x)\,,italic_x italic_δ italic_q = caligraphic_N start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT italic_η start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT italic_x start_POSTSUPERSCRIPT italic_a start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT end_POSTSUPERSCRIPT ( 1 - italic_x ) start_POSTSUPERSCRIPT italic_b start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT end_POSTSUPERSCRIPT ( 1 + italic_c start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT italic_x start_POSTSUPERSCRIPT 0.5 end_POSTSUPERSCRIPT + italic_d start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT italic_x ) , (13)

which is defined in terms of five shape fit parameters ηqsubscript𝜂𝑞\eta_{q}italic_η start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT, aqsubscript𝑎𝑞a_{q}italic_a start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT, bqsubscript𝑏𝑞b_{q}italic_b start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT, cqsubscript𝑐𝑞c_{q}italic_c start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT and dqsubscript𝑑𝑞d_{q}italic_d start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT, and the normalization factor 𝒩qsubscript𝒩𝑞{\cal N}_{q}caligraphic_N start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT as well. The normalization factor 𝒩𝒩{\cal N}caligraphic_N in Eq. 13 is calculated as follows:

1𝒩q1subscript𝒩𝑞\displaystyle\frac{1}{{\cal N}_{q}}divide start_ARG 1 end_ARG start_ARG caligraphic_N start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT end_ARG =\displaystyle== (1+dqaqaq+bq+1)B(aq,bq+1)1subscript𝑑𝑞subscript𝑎𝑞subscript𝑎𝑞subscript𝑏𝑞1𝐵subscript𝑎𝑞subscript𝑏𝑞1\displaystyle\left(1+d_{q}\frac{a_{q}}{a_{q}+b_{q}+1}\right)\,B\left(a_{q},b_{% q}+1\right)( 1 + italic_d start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT divide start_ARG italic_a start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT end_ARG start_ARG italic_a start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT + italic_b start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT + 1 end_ARG ) italic_B ( italic_a start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT , italic_b start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT + 1 ) (14)
+cqB(aq+12,bq+1),subscript𝑐𝑞𝐵subscript𝑎𝑞12subscript𝑏𝑞1\displaystyle+c_{q}B\left(a_{q}+\frac{1}{2},b_{q}+1\right)~{},+ italic_c start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT italic_B ( italic_a start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT + divide start_ARG 1 end_ARG start_ARG 2 end_ARG , italic_b start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT + 1 ) ,

which is chosen such that the ηqsubscript𝜂𝑞\eta_{q}italic_η start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT represents the first moments of δq(x,Q02)𝛿𝑞𝑥superscriptsubscript𝑄02\delta q(x,Q_{0}^{2})italic_δ italic_q ( italic_x , italic_Q start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ). The B(a,b)𝐵𝑎𝑏B(a,b)italic_B ( italic_a , italic_b ) indicates to the Euler beta function.

Since the AAK24 QCD analysis considers inclusive DIS data and the SIDIS observables as well, we attempt to consider the flavor asymmetric polarized PDFs and try to fit the PDFs δq=δu+δu¯,δd+δd¯,δu¯,δd¯,δs¯𝛿𝑞𝛿𝑢𝛿¯𝑢𝛿𝑑𝛿¯𝑑𝛿¯𝑢𝛿¯𝑑𝛿¯𝑠\delta q=\delta u+\delta\bar{u},\delta d+\delta\bar{d},\delta\bar{u},\delta% \bar{d},\delta\bar{s}italic_δ italic_q = italic_δ italic_u + italic_δ over¯ start_ARG italic_u end_ARG , italic_δ italic_d + italic_δ over¯ start_ARG italic_d end_ARG , italic_δ over¯ start_ARG italic_u end_ARG , italic_δ over¯ start_ARG italic_d end_ARG , italic_δ over¯ start_ARG italic_s end_ARG and δg𝛿𝑔\delta gitalic_δ italic_g at the input scale.

The observables obtained from the SIDIS processes can, in principle, provide the capability to distinguish between the strange quark density (δs𝛿𝑠\delta sitalic_δ italic_s) and its antiquark counterpart (δs¯𝛿¯𝑠\delta\bar{s}italic_δ over¯ start_ARG italic_s end_ARG) through separate production of K+superscript𝐾K^{+}italic_K start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT and Ksuperscript𝐾K^{-}italic_K start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT mesons. However, in practice, our current fit lacks the sensitivity to achieve a clear separation between these distributions. Therefore, in line with the polarized PDF analyses conducted in Refs. Sato:2016tuz ; Arbabifar:2013tma , we make the assumption of a flavor-symmetric strange quark distribution, meaning that δs(x,Q2)=δs¯(x,Q2)𝛿𝑠𝑥superscript𝑄2𝛿¯𝑠𝑥superscript𝑄2\delta{s}({x},Q^{2})=\delta{\bar{s}}({x},Q^{2})italic_δ italic_s ( italic_x , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) = italic_δ over¯ start_ARG italic_s end_ARG ( italic_x , italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ). To control the behavior of δu¯𝛿¯𝑢\delta\bar{u}italic_δ over¯ start_ARG italic_u end_ARG and δd¯𝛿¯𝑑\delta\bar{d}italic_δ over¯ start_ARG italic_d end_ARG distributions over the small x𝑥xitalic_x region, additional constraints on the shape parameters are applied by assuming aδu+δu¯=aδu¯subscript𝑎𝛿𝑢𝛿¯𝑢subscript𝑎𝛿¯𝑢a_{\delta u+\delta\bar{u}}=a_{\delta\bar{u}}italic_a start_POSTSUBSCRIPT italic_δ italic_u + italic_δ over¯ start_ARG italic_u end_ARG end_POSTSUBSCRIPT = italic_a start_POSTSUBSCRIPT italic_δ over¯ start_ARG italic_u end_ARG end_POSTSUBSCRIPT, and aδd+δd¯=aδd¯=aδssubscript𝑎𝛿𝑑𝛿¯𝑑subscript𝑎𝛿¯𝑑subscript𝑎𝛿𝑠a_{\delta d+\delta\bar{d}}=a_{\delta\bar{d}}=a_{\delta s}italic_a start_POSTSUBSCRIPT italic_δ italic_d + italic_δ over¯ start_ARG italic_d end_ARG end_POSTSUBSCRIPT = italic_a start_POSTSUBSCRIPT italic_δ over¯ start_ARG italic_d end_ARG end_POSTSUBSCRIPT = italic_a start_POSTSUBSCRIPT italic_δ italic_s end_POSTSUBSCRIPT as well.

In general, polarized PDF analyses utilize two well-known sum rules that relate the first moments of polarized PDFs to the F𝐹Fitalic_F and D𝐷Ditalic_D quantities, which are evaluated in neutron and hyperon β𝛽\betaitalic_β-decays PDG under the assumption of SU(2) and SU(3) flavor symmetries:

a3subscript𝑎3\displaystyle a_{3}italic_a start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT =\displaystyle== ΔΣuΔΣd=F+D,ΔsubscriptΣ𝑢ΔsubscriptΣ𝑑𝐹𝐷\displaystyle\Delta\Sigma_{u}-\Delta\Sigma_{d}=F+D,roman_Δ roman_Σ start_POSTSUBSCRIPT italic_u end_POSTSUBSCRIPT - roman_Δ roman_Σ start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT = italic_F + italic_D , (15)
a8subscript𝑎8\displaystyle a_{8}italic_a start_POSTSUBSCRIPT 8 end_POSTSUBSCRIPT =\displaystyle== ΔΣu+ΔΣd2ΔΣs=3FD.ΔsubscriptΣ𝑢ΔsubscriptΣ𝑑2ΔsubscriptΣ𝑠3𝐹𝐷\displaystyle\Delta\Sigma_{u}+\Delta\Sigma_{d}-2\Delta\Sigma_{s}=3F-D.roman_Δ roman_Σ start_POSTSUBSCRIPT italic_u end_POSTSUBSCRIPT + roman_Δ roman_Σ start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT - 2 roman_Δ roman_Σ start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT = 3 italic_F - italic_D . (16)

Here, a3subscript𝑎3a_{3}italic_a start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT and a8subscript𝑎8a_{8}italic_a start_POSTSUBSCRIPT 8 end_POSTSUBSCRIPT represent non-singlet combinations of the first moments of the polarized parton distributions, corresponding to the non-singlet q3subscript𝑞3q_{3}italic_q start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT and q8subscript𝑞8q_{8}italic_q start_POSTSUBSCRIPT 8 end_POSTSUBSCRIPT distributions:

q3subscript𝑞3\displaystyle q_{3}italic_q start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT =\displaystyle== (δu+δu¯)(δd+δd¯),𝛿𝑢𝛿¯𝑢𝛿𝑑𝛿¯𝑑\displaystyle(\delta u+\delta\overline{u})-(\delta d+\delta\overline{d}),( italic_δ italic_u + italic_δ over¯ start_ARG italic_u end_ARG ) - ( italic_δ italic_d + italic_δ over¯ start_ARG italic_d end_ARG ) , (17)
q8subscript𝑞8\displaystyle q_{8}italic_q start_POSTSUBSCRIPT 8 end_POSTSUBSCRIPT =\displaystyle== (δu+δu¯)+(δd+δd¯)2(δs+δs¯).𝛿𝑢𝛿¯𝑢𝛿𝑑𝛿¯𝑑2𝛿𝑠𝛿¯𝑠\displaystyle(\delta u+\delta\overline{u})+(\delta d+\delta\overline{d})-2(% \delta s+\delta\overline{s}).( italic_δ italic_u + italic_δ over¯ start_ARG italic_u end_ARG ) + ( italic_δ italic_d + italic_δ over¯ start_ARG italic_d end_ARG ) - 2 ( italic_δ italic_s + italic_δ over¯ start_ARG italic_s end_ARG ) . (18)

A recent reanalysis of the F𝐹Fitalic_F and D𝐷Ditalic_D parameters, with updated β𝛽\betaitalic_β-decay constants PDG , yielded F=0.464±0.008𝐹plus-or-minus0.4640.008F=0.464\pm 0.008italic_F = 0.464 ± 0.008 and D=0.806±0.008𝐷plus-or-minus0.8060.008D=0.806\pm 0.008italic_D = 0.806 ± 0.008. We have utilized these updated values in our current analysis. However, since our focus is not on flavor symmetry and we have δu¯δd¯δs𝛿¯𝑢𝛿¯𝑑𝛿𝑠\delta\bar{u}\neq\delta\bar{d}\neq\delta sitalic_δ over¯ start_ARG italic_u end_ARG ≠ italic_δ over¯ start_ARG italic_d end_ARG ≠ italic_δ italic_s, we can use the combination of Eqs. 15 and 16 as follows:

Δu+Δu¯Δ𝑢Δ¯𝑢\displaystyle\Delta u+\Delta\bar{u}roman_Δ italic_u + roman_Δ over¯ start_ARG italic_u end_ARG =\displaystyle== 0.9275+Δs+Δs¯,0.9275Δ𝑠Δ¯𝑠\displaystyle 0.9275+\Delta s+\Delta\bar{s},0.9275 + roman_Δ italic_s + roman_Δ over¯ start_ARG italic_s end_ARG ,
Δd+Δd¯Δ𝑑Δ¯𝑑\displaystyle\Delta d+\Delta\bar{d}roman_Δ italic_d + roman_Δ over¯ start_ARG italic_d end_ARG =\displaystyle== 0.3415+Δs+Δs¯.0.3415Δ𝑠Δ¯𝑠\displaystyle-0.3415+\Delta s+\Delta\bar{s}.- 0.3415 + roman_Δ italic_s + roman_Δ over¯ start_ARG italic_s end_ARG . (19)

We apply these relations in our analysis, thus excluding the parameters that define the first moments of (δu+δu¯)𝛿𝑢𝛿¯𝑢(\delta u+\delta\bar{u})( italic_δ italic_u + italic_δ over¯ start_ARG italic_u end_ARG ) and (δd+δd¯)𝛿𝑑𝛿¯𝑑(\delta d+\delta\bar{d})( italic_δ italic_d + italic_δ over¯ start_ARG italic_d end_ARG ) (i.e., ηu+u¯subscript𝜂𝑢¯𝑢\eta_{u+\bar{u}}italic_η start_POSTSUBSCRIPT italic_u + over¯ start_ARG italic_u end_ARG end_POSTSUBSCRIPT and ηd+d¯subscript𝜂𝑑¯𝑑\eta_{d+\bar{d}}italic_η start_POSTSUBSCRIPT italic_d + over¯ start_ARG italic_d end_ARG end_POSTSUBSCRIPT) from the analysis and obtaining them using Eq. 19. The impact of symmetry breaking on the first moments of polarized PDFs has been discussed in detail in the literature su3 ; deFlorian:2009vb .

It is important to mention that the calculations take into account Isospin symmetry for both the proton and neutron. This symmetry is defined as follows:

δup=δdn;δu¯p=δd¯n;δdp=δunformulae-sequence𝛿superscript𝑢𝑝𝛿superscript𝑑𝑛formulae-sequence𝛿superscript¯𝑢𝑝𝛿superscript¯𝑑𝑛𝛿superscript𝑑𝑝𝛿superscript𝑢𝑛\displaystyle\delta u^{p}=\delta d^{n}\,\,\,;\,\delta\bar{u}^{p}=\delta\bar{d}% ^{n}\,\,\,;\,\delta d^{p}=\delta u^{n}italic_δ italic_u start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT = italic_δ italic_d start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT ; italic_δ over¯ start_ARG italic_u end_ARG start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT = italic_δ over¯ start_ARG italic_d end_ARG start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT ; italic_δ italic_d start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT = italic_δ italic_u start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT
δd¯p=δu¯n;δsp=δsn;δs¯p=δs¯n.formulae-sequence𝛿superscript¯𝑑𝑝𝛿superscript¯𝑢𝑛formulae-sequence𝛿superscript𝑠𝑝𝛿superscript𝑠𝑛𝛿superscript¯𝑠𝑝𝛿superscript¯𝑠𝑛\displaystyle\delta\bar{d}^{p}=\delta\bar{u}^{n}\,\,\,;\,\delta s^{p}=\delta s% ^{n}\,\,\,;\,\delta\bar{s}^{p}=\delta\bar{s}^{n}\,.italic_δ over¯ start_ARG italic_d end_ARG start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT = italic_δ over¯ start_ARG italic_u end_ARG start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT ; italic_δ italic_s start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT = italic_δ italic_s start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT ; italic_δ over¯ start_ARG italic_s end_ARG start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT = italic_δ over¯ start_ARG italic_s end_ARG start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT . (20)

Finally, we should highlight here that for some shape parameters as presented in Eq. 13 may lack strong constraints from available data. In such cases, additional assumptions become necessary. We will provide a detailed discussion of these assumptions, as considered in the AAK24 analysis, in Sec. 6.1.

5 Minimization strategy

In this section, we will introduce the AAK24 methodology employed for extracting polarized PDFs and for estimating their uncertainties, specifically through the standard Hessian method.

The AAK24 best-fit parameters are determined through the minimization of the χ2superscript𝜒2\chi^{2}italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT function, which is defined as follows:

χglobal2superscriptsubscript𝜒global2\displaystyle\chi_{\rm global}^{2}italic_χ start_POSTSUBSCRIPT roman_global end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT =\displaystyle== nwnχn2,subscript𝑛subscript𝑤𝑛superscriptsubscript𝜒𝑛2\displaystyle\sum_{n}{w}_{n}\chi_{n}^{2},∑ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT italic_w start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT italic_χ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ,
χn2superscriptsubscript𝜒𝑛2\displaystyle\chi_{n}^{2}italic_χ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT =\displaystyle== (1NnΔNn)2+i(NnA1,2,idataA1,2,itheoryNnΔA1,2,idata)2,superscript1subscript𝑁𝑛Δsubscript𝑁𝑛2subscript𝑖superscriptsubscript𝑁𝑛superscriptsubscript𝐴12𝑖datasuperscriptsubscript𝐴12𝑖theorysubscript𝑁𝑛Δsuperscriptsubscript𝐴12𝑖data2\displaystyle\left(\frac{{1}-{N}_{n}}{\Delta N_{n}}\right)^{2}+\sum_{i}\left(% \frac{{N}_{n}A_{1,2,i}^{\rm data}-A_{1,2,i}^{\rm theory}}{N_{n}\Delta A_{1,2,i% }^{\rm data}}\right)^{2},( divide start_ARG 1 - italic_N start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT end_ARG start_ARG roman_Δ italic_N start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + ∑ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( divide start_ARG italic_N start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT italic_A start_POSTSUBSCRIPT 1 , 2 , italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT roman_data end_POSTSUPERSCRIPT - italic_A start_POSTSUBSCRIPT 1 , 2 , italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT roman_theory end_POSTSUPERSCRIPT end_ARG start_ARG italic_N start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT roman_Δ italic_A start_POSTSUBSCRIPT 1 , 2 , italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT roman_data end_POSTSUPERSCRIPT end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ,

In the equation above, the quantity A1,2,idatasuperscriptsubscript𝐴12𝑖dataA_{1,2,i}^{\rm data}italic_A start_POSTSUBSCRIPT 1 , 2 , italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT roman_data end_POSTSUPERSCRIPT represents the measured value of the observable for the data point i𝑖iitalic_i from the experimental data set n𝑛nitalic_n. A1,2,itheorysuperscriptsubscript𝐴12𝑖theoryA_{1,2,i}^{\rm theory}italic_A start_POSTSUBSCRIPT 1 , 2 , italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT roman_theory end_POSTSUPERSCRIPT refers to the corresponding theoretical value for the ithsuperscript𝑖thi^{\rm th}italic_i start_POSTSUPERSCRIPT roman_th end_POSTSUPERSCRIPT experimental data point, and ΔA1,2,idataΔsuperscriptsubscript𝐴12𝑖data\Delta A_{1,2,i}^{\rm data}roman_Δ italic_A start_POSTSUBSCRIPT 1 , 2 , italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT roman_data end_POSTSUPERSCRIPT denotes the experimental uncertainty, which is typically the combination of statistical and systematic uncertainties computed in quadrature. Nnsubscript𝑁𝑛{N}_{n}italic_N start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT is an overall normalization factor for the data of experimental data set n𝑛nitalic_n, and the ΔNnΔsubscript𝑁𝑛{\Delta{N}_{n}}roman_Δ italic_N start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT is the experimental normalization uncertainty. We allow for a relative normalization shift Nnsubscript𝑁𝑛{N}_{n}italic_N start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT between different experimental data sets within uncertainties of ΔNnΔsubscript𝑁𝑛{\Delta{N}_{n}}roman_Δ italic_N start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT quoted by the experiments.

For conducting the fitting process and deriving the AAK24 best-fit parameters, we utilize the publicly available MINUIT package, which has been developed by CERN James:1975dr . A very good total χ2/d.o.f=1.0528formulae-sequencesuperscript𝜒2𝑑𝑜𝑓1.0528\chi^{2}/d.o.f=1.0528italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / italic_d . italic_o . italic_f = 1.0528 value is achieved in these calculations. Table. 4 presents the numeric values of the AAK24 fit parameters obtained through the QCD analysis. The χ2superscript𝜒2\chi^{2}italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT values for the total and individual data sets are presented in Tables. 1, 2 and 3 respectively. The detailed discussions on the obtained best fit parameters, and the fit quality will be presented in Sec. 6.1.

At this point, we are ready to introduce the Hessian method for the estimation of polarized PDF uncertainty. The standard Hessian method employed in the AAK24 QCD analysis was originally developed in Refs. Martin:2002aw ; Pumplin:2001ct . This method essentially aims to approximate the posterior distribution as a Gaussian distribution within the parameter space. To do this end, one first needs to find the set of parameters 𝐚𝟎subscript𝐚0\bf{a}_{0}bold_a start_POSTSUBSCRIPT bold_0 end_POSTSUBSCRIPT corresponding to the maximum a posteriori estimate. These the set of parameters could be obtained by minimizing a certain the χ2superscript𝜒2\chi^{2}italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT function in Eq. 5. In a neighborhood of 𝐚𝟎subscript𝐚0\bf{a}_{0}bold_a start_POSTSUBSCRIPT bold_0 end_POSTSUBSCRIPT parameters, the χ2(𝐚)superscript𝜒2𝐚\chi^{2}(\bf{a})italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( bold_a ) function can be expanded up to the quadratic term in its parameters and the posterior approximated by a multi-dimensional Gaussian function,

p(a|m)exp(12χ2(a,m))exp(12ΔaTHΔa).proportional-to𝑝conditionalam12superscript𝜒2a𝑚proportional-to12Δsuperscripta𝑇𝐻Δa\displaystyle p(\textbf{{a}}|\textbf{{m}})\propto\exp\Big{(}-\frac{{1}}{{2}}% \chi^{{2}}(\textbf{a},{m})\Big{)}\propto\exp\Big{(}-\frac{1}{2}\Delta\textbf{a% }^{T}\,H\,{\Delta}\textbf{a}\Big{)}\,.italic_p ( a | m ) ∝ roman_exp ( - divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( a , italic_m ) ) ∝ roman_exp ( - divide start_ARG 1 end_ARG start_ARG 2 end_ARG roman_Δ a start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT italic_H roman_Δ a ) .

In the last term, Δa=aa0Δaasubscripta0\Delta\textbf{a}=\textbf{a}-\textbf{a}_{0}roman_Δ a = a - a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, the Hessian matrix elements are given by

Hij=122χ2(a)aiaj|a=a0,i,j=1,npar.formulae-sequencesubscript𝐻𝑖𝑗evaluated-at12superscript2superscript𝜒2asuperscript𝑎𝑖superscript𝑎𝑗asubscripta0𝑖𝑗1subscript𝑛parH_{ij}=\frac{{1}}{{2}}\left.\frac{\partial^{{2}}\chi^{{2}}(\textbf{a})}{% \partial{a}^{{i}}\partial{a}^{j}}\right|_{\textbf{a}=\textbf{a}_{0}}\,,\quad i% ,j=1,\dots{{n}}_{\rm par}\,.italic_H start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG 2 end_ARG divide start_ARG ∂ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( a ) end_ARG start_ARG ∂ italic_a start_POSTSUPERSCRIPT italic_i end_POSTSUPERSCRIPT ∂ italic_a start_POSTSUPERSCRIPT italic_j end_POSTSUPERSCRIPT end_ARG | start_POSTSUBSCRIPT a = a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT , italic_i , italic_j = 1 , … italic_n start_POSTSUBSCRIPT roman_par end_POSTSUBSCRIPT . (23)

and exp(12χ2(a0))12superscript𝜒2subscripta0\exp(-\frac{1}{2}{\chi}^{2}(\textbf{a}_{0}))roman_exp ( - divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) ) has been absorbed into the normalization of the posterior.

Additionally, the tolerances can be fine-tuned by setting a criterion that the χ2superscript𝜒2{\chi}^{2}italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT must increase above its minimum by a specified Δχ2Δsuperscript𝜒2\Delta{\chi}^{2}roman_Δ italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT value, as discussed in Ref. Accardi:2021ysh . In the AAK24 QCD analysis, fits will be conducted on a statistically consistent set of polarized DIS data, and thus, there is no need to employ a tolerance criterion. Nevertheless, it is important to calculate the tolerance explicitly in order to account for any actual deviations from Gaussian behavior and to obtain precise estimates of the variance. In the results that will be presented in the following Sec. 6, we take the tolerance criterion into account by setting Δχ2=T=1Δsuperscript𝜒2𝑇1\Delta{\chi}^{2}=T=1roman_Δ italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = italic_T = 1.

6 AAK24 fit results

This section includes the main results and findings of AAK24 polarized PDFs obtained from QCD analysis of polarized DIS and SIDIS data. We first present the AAK24 polarized PDFs, the fit parameters, and the results for different Parton species. Then we compare the AAK24 NLO theoretical predictions with the analyzed data. We also compare our results with other polarized PDFs available in the literature by discussing in detail the similarity and differences between these different polarized PDF sets.

6.1 AAK24 polarized PDFs

This section presents the AAK24 polarized PDFs at the input scale, the extracted best-fit parameters, and assumptions considered for the shape fit parameter as well.

The best fit parameters for AAK24 are presented in Table. 4. Upon examining the numbers presented in this table, several remarks are pertinent. In general, a QCD analysis of polarized DIS experiments should incorporate additional constraints on the shape parameters. This necessity arises from the affects of available data to fully constrain all parameters and the limited kinematic coverage of the data. Given that the existing experimental data sets lack the power to comprehensively determine the behavior of polarized PDFs, certain assumptions must be taken into consideration.

Firstly, to control the small x𝑥xitalic_x behavior of the xδu¯𝑥𝛿¯𝑢x\delta\bar{u}italic_x italic_δ over¯ start_ARG italic_u end_ARG, xδd¯𝑥𝛿¯𝑑x\delta\bar{d}italic_x italic_δ over¯ start_ARG italic_d end_ARG and xδs=xδs¯𝑥𝛿𝑠𝑥𝛿¯𝑠x\delta s=x\delta\bar{s}italic_x italic_δ italic_s = italic_x italic_δ over¯ start_ARG italic_s end_ARG, we assume aδu+δu¯=aδu¯subscript𝑎𝛿𝑢𝛿¯𝑢subscript𝑎𝛿¯𝑢a_{\delta u+\delta\bar{u}}=a_{\delta\bar{u}}italic_a start_POSTSUBSCRIPT italic_δ italic_u + italic_δ over¯ start_ARG italic_u end_ARG end_POSTSUBSCRIPT = italic_a start_POSTSUBSCRIPT italic_δ over¯ start_ARG italic_u end_ARG end_POSTSUBSCRIPT and aδd+δd¯=aδd¯=aδs¯subscript𝑎𝛿𝑑𝛿¯𝑑subscript𝑎𝛿¯𝑑subscript𝑎𝛿¯𝑠a_{\delta d+\delta\bar{d}}=a_{\delta\bar{d}}=a_{\delta\bar{s}}italic_a start_POSTSUBSCRIPT italic_δ italic_d + italic_δ over¯ start_ARG italic_d end_ARG end_POSTSUBSCRIPT = italic_a start_POSTSUBSCRIPT italic_δ over¯ start_ARG italic_d end_ARG end_POSTSUBSCRIPT = italic_a start_POSTSUBSCRIPT italic_δ over¯ start_ARG italic_s end_ARG end_POSTSUBSCRIPT. Since we have determined that the best-fitted values of bδu¯subscript𝑏𝛿¯𝑢b_{\delta\bar{u}}italic_b start_POSTSUBSCRIPT italic_δ over¯ start_ARG italic_u end_ARG end_POSTSUBSCRIPT, bδd¯subscript𝑏𝛿¯𝑑b_{\delta\bar{d}}italic_b start_POSTSUBSCRIPT italic_δ over¯ start_ARG italic_d end_ARG end_POSTSUBSCRIPT, and bδs=δs¯subscript𝑏𝛿𝑠𝛿¯𝑠b_{\delta s=\delta\bar{s}}italic_b start_POSTSUBSCRIPT italic_δ italic_s = italic_δ over¯ start_ARG italic_s end_ARG end_POSTSUBSCRIPT are not strongly constrained by the QCD fit and are very similar, all approximately equal to 10, we have chosen to maintain them as fixed values at 10.

Furthermore, we have observed that the parameter cqsubscript𝑐𝑞c_{q}italic_c start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT for certain parton species, such as xδu+xδu¯𝑥𝛿𝑢𝑥𝛿¯𝑢x\delta u+x\delta\bar{u}italic_x italic_δ italic_u + italic_x italic_δ over¯ start_ARG italic_u end_ARG, xδd+xδd¯𝑥𝛿𝑑𝑥𝛿¯𝑑x\delta d+x\delta\bar{d}italic_x italic_δ italic_d + italic_x italic_δ over¯ start_ARG italic_d end_ARG, xδs¯𝑥𝛿¯𝑠x\delta\bar{s}italic_x italic_δ over¯ start_ARG italic_s end_ARG, and xδg𝑥𝛿𝑔x\delta gitalic_x italic_δ italic_g, is extremely close to zero. As a result, we have opted to set these parameters to a fixed value of 0. Additionally, it is worth noting that some of the c𝑐citalic_c and d𝑑ditalic_d shape parameters for the xδd¯𝑥𝛿¯𝑑x\delta\bar{d}italic_x italic_δ over¯ start_ARG italic_d end_ARG density have been determined with relatively large errors. Incorporating additional high-energy polarized data sets from various types of experiments has the potential to alleviate some of the extra assumptions we have made regarding the shape parameters of the polarized distribution. To be more specific, the inclusion of proton-proton collision data from RHIC can have a significant impact on our understanding of the distributions within the nucleon, particularly in relation to the gluon PDFs.

Flavor ηqsubscript𝜂𝑞\eta_{q}italic_η start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT aqsubscript𝑎𝑞a_{q}italic_a start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT bqsubscript𝑏𝑞b_{q}italic_b start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT cqsubscript𝑐𝑞c_{q}italic_c start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT dqsubscript𝑑𝑞d_{q}italic_d start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT
δu+δu¯𝛿𝑢𝛿¯𝑢\delta u+\delta\bar{u}italic_δ italic_u + italic_δ over¯ start_ARG italic_u end_ARG 0.863 0.4126±plus-or-minus\pm± 0.0079 2.4636±plus-or-minus\pm± 0.0174 0.0superscript0.00.0^{*}0.0 start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT 40.9115±+1.3840.9115\pm+1.3840.9115 ± + 1.38
δd+δd¯𝛿𝑑𝛿¯𝑑\delta d+\delta\bar{d}italic_δ italic_d + italic_δ over¯ start_ARG italic_d end_ARG -0.4053 0.2487±plus-or-minus\pm± 0.0059 4.5203±plus-or-minus\pm± 0.00272 0.0superscript0.00.0^{*}0.0 start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT 100.9866±10.92plus-or-minus100.986610.92100.9866\pm 10.92100.9866 ± 10.92
δu¯𝛿¯𝑢\delta\bar{u}italic_δ over¯ start_ARG italic_u end_ARG 0.0639±plus-or-minus\pm±0.175 0.4126±plus-or-minus\pm± 0.0079 10.0 12.5012±plus-or-minus\pm± 0.923 39.5083±8.09plus-or-minus39.50838.09-39.5083\pm 8.09- 39.5083 ± 8.09
δd¯𝛿¯𝑑\delta\bar{d}italic_δ over¯ start_ARG italic_d end_ARG -0.0288±plus-or-minus\pm± 0.355 0.2487±plus-or-minus\pm± 0.0059 10.0 102.436±plus-or-minus\pm±52.07 800.569±284.799plus-or-minus800.569284.799800.569\pm 284.799800.569 ± 284.799
δs¯𝛿¯𝑠\delta\bar{s}italic_δ over¯ start_ARG italic_s end_ARG -0.0433±plus-or-minus\pm± 0.0017 0.2487±plus-or-minus\pm± 0.0059 10.0 0.0superscript0.00.0^{*}0.0 start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT 28.8485±0.0027plus-or-minus28.84850.0027-28.8485\pm 0.0027- 28.8485 ± 0.0027
δg𝛿𝑔\delta gitalic_δ italic_g -0.3788±plus-or-minus\pm±0.0026 2.9019 ±plus-or-minus\pm± 0.0908 10.0 0.0superscript0.00.0^{*}0.0 start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT 3.271±0.0875plus-or-minus3.2710.0875-3.271\pm 0.0875- 3.271 ± 0.0875
Table 4: The best-fitted values and their associated statistical errors are provided at the input scale Q02=1superscriptsubscript𝑄021Q_{0}^{2}=1italic_Q start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = 1 GeV2. It is important to note that parameters marked with an asterisk (*) have been held fixed.

As mentioned earlier, we also take into account the TMC and HT corrections, which play important roles in certain kinematic regions of x𝑥xitalic_x and Q2. The values of the HT corrections extracted from the data in our analysis are presented in Table. 5.

Table 5: The values of HT corrections extracted from the data in a model independent way.
xisubscript𝑥𝑖x_{i}italic_x start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT     hp(xi)[GeV2]superscript𝑝subscript𝑥𝑖delimited-[]superscriptGeV2h^{p}(x_{i})~{}[\mathrm{GeV}^{2}]italic_h start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT ( italic_x start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) [ roman_GeV start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ]        xisubscript𝑥𝑖x_{i}italic_x start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT     hn(xi)[GeV2]superscript𝑛subscript𝑥𝑖delimited-[]superscriptGeV2h^{n}(x_{i})~{}[\mathrm{GeV}^{2}]italic_h start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT ( italic_x start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) [ roman_GeV start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ]
0.028     0.03092 ±plus-or-minus\pm± 0.00817       0.028     0.0418 ±plus-or-minus\pm± 0.0170
0.100     -0.08545 ±plus-or-minus\pm± 0.00386       0.100     0.0556 ±plus-or-minus\pm± 0.0140
0.200     -0.04048 ±plus-or-minus\pm± 0.00270       0.300     -0.01539 ±plus-or-minus\pm± 0.0108
0.350     -0.03238 ±plus-or-minus\pm± 0.00274       0.300     0.008924 ±plus-or-minus\pm± 0.00581
0.600     -0.07021 ±plus-or-minus\pm± 0.00443       0.500     -0.008407 ±plus-or-minus\pm± 0.00374

In Fig. 4, we display the AAK24 polarized PDFs at the input scale for all parton species including xδu+xδu¯𝑥𝛿𝑢𝑥𝛿¯𝑢x\delta u+x\delta\bar{u}italic_x italic_δ italic_u + italic_x italic_δ over¯ start_ARG italic_u end_ARG, xδd+xδd¯𝑥𝛿𝑑𝑥𝛿¯𝑑x\delta d+x\delta\bar{d}italic_x italic_δ italic_d + italic_x italic_δ over¯ start_ARG italic_d end_ARG, xδu¯𝑥𝛿¯𝑢x\delta\bar{u}italic_x italic_δ over¯ start_ARG italic_u end_ARG, xδd¯𝑥𝛿¯𝑑x\delta\bar{d}italic_x italic_δ over¯ start_ARG italic_d end_ARG, xδs¯𝑥𝛿¯𝑠x\delta\bar{s}italic_x italic_δ over¯ start_ARG italic_s end_ARG and xδg𝑥𝛿𝑔x\delta gitalic_x italic_δ italic_g. The results for the case where we have excluded the measurements of A2psuperscriptsubscript𝐴2𝑝A_{2}^{p}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT, A2dsuperscriptsubscript𝐴2𝑑A_{2}^{d}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT, and A2nsuperscriptsubscript𝐴2𝑛A_{2}^{n}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT asymmetries are also presented for the purpose of comparison. The ability of determining the polarized sea-quark distributions are due to the impact of SIDIS data in which added to our data sample. The error bands, obtained through the Hessian method, are also shown as well. Upon examining the distributions in this figure, it becomes evident that the sea quark densities xδu¯𝑥𝛿¯𝑢x\delta\bar{u}italic_x italic_δ over¯ start_ARG italic_u end_ARG and xδd¯𝑥𝛿¯𝑑x\delta\bar{d}italic_x italic_δ over¯ start_ARG italic_d end_ARG exhibit relatively large uncertainties. Additionally, the input strange sea-quark distribution xδs=xδs¯𝑥𝛿𝑠𝑥𝛿¯𝑠x\delta s=x\delta{\bar{s}}italic_x italic_δ italic_s = italic_x italic_δ over¯ start_ARG italic_s end_ARG displays the typical sign change.

Refer to caption
Refer to caption
Refer to caption
Refer to caption
Refer to caption
Refer to caption
Figure 4: The AAK24 quark and gluon helicity distributions at the input scale Q02=1subscriptsuperscript𝑄201Q^{2}_{0}=1italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 1 GeV2 with and without the A2subscript𝐴2A_{2}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT asymmetry measurements. The uncertainty bands, determined using the Hessian method, are also shown as well.

Let us now discuss the results obtained by excluding the A2subscript𝐴2A_{2}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT asymmetry measurements to examine the effect of these data sets on both the central values and extracted uncertainties. As observed in Fig. 4, the inclusion of the A2subscript𝐴2A_{2}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT asymmetry data not only affects the central values of certain parton species but also impacts the error bands. These effects are most pronounced in the case of xδu¯𝑥𝛿¯𝑢x\delta\bar{u}italic_x italic_δ over¯ start_ARG italic_u end_ARG, xδd¯𝑥𝛿¯𝑑x\delta\bar{d}italic_x italic_δ over¯ start_ARG italic_d end_ARG, xδs𝑥𝛿𝑠x\delta sitalic_x italic_δ italic_s, and gluon densities. The most significant changes are noticeable in the case of the gluon PDF in terms of central value, and in xδd¯𝑥𝛿¯𝑑x\delta\bar{d}italic_x italic_δ over¯ start_ARG italic_d end_ARG and xδs𝑥𝛿𝑠x\delta sitalic_x italic_δ italic_s in terms of error bands. In terms of individual parton flavors, it can be observed that for xδu+xδu¯𝑥𝛿𝑢𝑥𝛿¯𝑢x\delta u+x\delta\bar{u}italic_x italic_δ italic_u + italic_x italic_δ over¯ start_ARG italic_u end_ARG and xδd+xδd¯𝑥𝛿𝑑𝑥𝛿¯𝑑x\delta d+x\delta\bar{d}italic_x italic_δ italic_d + italic_x italic_δ over¯ start_ARG italic_d end_ARG, the two results exhibit a similar shape, with a slight reduction of error bands for xδd+xδd¯𝑥𝛿𝑑𝑥𝛿¯𝑑x\delta d+x\delta\bar{d}italic_x italic_δ italic_d + italic_x italic_δ over¯ start_ARG italic_d end_ARG. In the case of xδu¯𝑥𝛿¯𝑢x\delta\bar{u}italic_x italic_δ over¯ start_ARG italic_u end_ARG, when compared to the analysis excluding A2subscript𝐴2A_{2}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT, it is evident that the distribution is slightly increased at large values of x𝑥xitalic_x and slightly decreased over the medium to small range of x𝑥xitalic_x, with the same size of error bands spanning the entire range of x𝑥xitalic_x. Regarding the xs𝑥𝑠xsitalic_x italic_s density, in general, the shape remains nearly unchanged while the error band is increased. In the case of the gluon PDF, the inclusion of the A2subscript𝐴2A_{2}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT asymmetry data significantly affects the central value, while the error band remains unchanged. Finally, for xδd¯𝑥𝛿¯𝑑x\delta\bar{d}italic_x italic_δ over¯ start_ARG italic_d end_ARG, the central values are the same, but the analysis including the A2subscript𝐴2A_{2}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT asymmetry measurements leads to smaller error bands.

6.2 Comparison with other polrized PDFs

In this section, we present the AAK24 polarized PDFs and provide detailed comparisons with other results available in the literature, namely the results from DSSV14 deFlorian:2014yva , AKS14 Arbabifar:2013tma , LSS10 Leader:2010rb , and JAM17 Sato:2016tuz all evaluated at the input scale Q02=1subscriptsuperscript𝑄201Q^{2}_{0}=1italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 1 GeV2 and at higher value of 10 GeV2. The comparisons are displayed in Figs. 5 and 6 at both the input scale of Q2=1superscript𝑄21Q^{2}=1italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = 1 GeV2 and a higher value of 10 GeV2.

Examining the xδu+xδu¯𝑥𝛿𝑢𝑥𝛿¯𝑢x\delta u+x\delta\bar{u}italic_x italic_δ italic_u + italic_x italic_δ over¯ start_ARG italic_u end_ARG distribution, generally speaking, one can observe that all of the fits are in very good agreement with each other. In the case of the δd+δd¯𝛿𝑑𝛿¯𝑑\delta d+\delta\bar{d}italic_δ italic_d + italic_δ over¯ start_ARG italic_d end_ARG distribution, all the groups agree very well, except for the DSSV14 result, which is smaller than those of the others at the input scale. However, the differences are within the uncertainty. As for Q2=10superscript𝑄210Q^{2}=10italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = 10 GeV2, the AAK24 polarized PDFs, JAM17 and LSS10 agree over the entire x𝑥xitalic_x region. However, for all other parton species, clear differences can be observed among the results of various groups.

A significant difference can be observed in the case of the sea-quark densities xδu¯𝑥𝛿¯𝑢x\delta\bar{u}italic_x italic_δ over¯ start_ARG italic_u end_ARG and xδd¯𝑥𝛿¯𝑑x\delta\bar{d}italic_x italic_δ over¯ start_ARG italic_d end_ARG. For the case of xδu¯𝑥𝛿¯𝑢x\delta\bar{u}italic_x italic_δ over¯ start_ARG italic_u end_ARG, the AAK24 results are larger than those of other groups for x<0.1𝑥0.1x<0.1italic_x < 0.1 at the input scale. However, the difference becomes less pronounced at the higher scale. Despite the similar results between AAK24 and JAM17, for the xδd¯𝑥𝛿¯𝑑x\delta\bar{d}italic_x italic_δ over¯ start_ARG italic_d end_ARG density at both input scale and higher values of Q2 = 10 GeV2, the other results exhibit different behaviors at

Examining the strange distribution, we can observe that we obtained very similar results compared to those of JAM17 at the input scale. For the case of the gluon density xδg𝑥𝛿𝑔x\delta gitalic_x italic_δ italic_g, the difference between all groups are more pronounced in this case. Like for the case of LSS10, our gluon density is smaller than those of others for medium value of x𝑥xitalic_x; (x0.1similar-to𝑥0.1x\sim 0.1italic_x ∼ 0.1)

Focusing on the JAM17 Sato:2016tuz polarized PDFs, one can observe significant differences in their gluon density xδg𝑥𝛿𝑔x\delta gitalic_x italic_δ italic_g, xδu¯𝑥𝛿¯𝑢x\delta\bar{u}italic_x italic_δ over¯ start_ARG italic_u end_ARG, xδd¯𝑥𝛿¯𝑑x\delta\bar{d}italic_x italic_δ over¯ start_ARG italic_d end_ARG, and strange distribution compared to all other distributions presented in Fig. 5. Distinctions are also more apparent for xδu¯𝑥𝛿¯𝑢x\delta\bar{u}italic_x italic_δ over¯ start_ARG italic_u end_ARG and xδd¯𝑥𝛿¯𝑑x\delta\bar{d}italic_x italic_δ over¯ start_ARG italic_d end_ARG, and xδg𝑥𝛿𝑔x\delta gitalic_x italic_δ italic_g. It is essential to highlight that the AAK24 and JAM17 are very similar in the case of x(δu+δu¯)𝑥𝛿𝑢𝛿¯𝑢x(\delta u+\delta\bar{u})italic_x ( italic_δ italic_u + italic_δ over¯ start_ARG italic_u end_ARG ) and x(δd+δd¯)𝑥𝛿𝑑𝛿¯𝑑x(\delta d+\delta\bar{d})italic_x ( italic_δ italic_d + italic_δ over¯ start_ARG italic_d end_ARG ) and strange distribution.

These findings hold at a higher scale of 10 GeV2 as well, as depicted in Fig. 6. Despite the similarity in x(δu+δu¯)𝑥𝛿𝑢𝛿¯𝑢x(\delta u+\delta\bar{u})italic_x ( italic_δ italic_u + italic_δ over¯ start_ARG italic_u end_ARG ) and x(δd+δd¯)𝑥𝛿𝑑𝛿¯𝑑x(\delta d+\delta\bar{d})italic_x ( italic_δ italic_d + italic_δ over¯ start_ARG italic_d end_ARG ), differences become more pronounced for the gluon density xδg𝑥𝛿𝑔x\delta gitalic_x italic_δ italic_g, sea quarks and strange distribution xδs𝑥𝛿𝑠x\delta sitalic_x italic_δ italic_s. It is crucial to emphasize that the JAM17 analysis employs the iterative Monte Carlo fitting technique, generating stable fits to polarized PDFs with statistically rigorous uncertainties.

Refer to caption
Refer to caption
Refer to caption
Refer to caption
Refer to caption
Refer to caption
Figure 5: The AAK24 quark and gluon helicity distributions at the input scale Q02=1subscriptsuperscript𝑄201Q^{2}_{0}=1italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 1 GeV2, with uncertainty bands now clearly corresponding to the different PDF sets. The comparisons include DSSV14 deFlorian:2014yva (black dashed), LSS10 Leader:2010rb (blue dashed dotted), JAM17 Sato:2016tuz (green dashed dot dotted), and AKS14 Arbabifar:2013tma (red short dashed dotted). The uncertainty bands for each set are shown with the same color as their corresponding line in the legend.
Refer to caption
Refer to caption
Refer to caption
Refer to caption
Refer to caption
Refer to caption
Figure 6: Same as Fig. 5 but this time at the higher value of 10 GeV2.

6.3 The AAK24 fit quality and the data/theory comparison

We will now shift our focus to discussing the fit quality of AAK24 and conducting a detailed comparison of NLO theory predictions with the data analyzed in this study.

The quality of the AAK24 global QCD analysis is summarized in Tables. 1, 2 and 3. For both the A1subscript𝐴1A_{1}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and g1F1subscript𝑔1subscript𝐹1\frac{g_{1}}{F_{1}}divide start_ARG italic_g start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_F start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG observables, we have presented both the total and individual χ2superscript𝜒2\chi^{2}italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT values in Table 1. It is noteworthy that in the total, we have considered 951 data points, resulting in a total χ2superscript𝜒2\chi^{2}italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT of 1170.172, which indicates a relatively good fit within the framework of QCD. When examining individual χ2superscript𝜒2\chi^{2}italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT values, it becomes evident that in most cases, there is a strong agreement between the data and theory. This is particularly evident in the case of the CLAS data, which contributes significantly to the total dataset with 136 data points for A1psuperscriptsubscript𝐴1𝑝A_{1}^{p}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT and 452 for A1dsuperscriptsubscript𝐴1𝑑A_{1}^{d}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT. For the newly added COMPASS and JLAB data, our observations indicate that the fit quality is generally good for COMPASS16 and COMPASS17. However, for the JLAB17 dataset, a less favorable χ2superscript𝜒2\chi^{2}italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT value is obtained, suggesting a relatively worse agreement between the data and theoretical predictions in this case.

The corresponding χ2superscript𝜒2\chi^{2}italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT for the A2subscript𝐴2A_{2}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT asymmetry observables, with a total of Ndat=140subscript𝑁𝑑𝑎𝑡140N_{dat}=140italic_N start_POSTSUBSCRIPT italic_d italic_a italic_t end_POSTSUBSCRIPT = 140 data points, are available in Table. 2. It is worth noting that not only do we achieve a very favorable total χ2superscript𝜒2\chi^{2}italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT value, but we also observe good individual χ2superscript𝜒2\chi^{2}italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT values for datasets originating from various experiments analyzed in this work. This suggests a strong agreement between the data and theoretical predictions across different experimental sources.

Finally, Table. 3 provides the χ2superscript𝜒2\chi^{2}italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT values for the A1Nhsuperscriptsubscript𝐴1𝑁A_{1N}^{h}italic_A start_POSTSUBSCRIPT 1 italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT SIDIS observables that have been examined in this study as well, encompassing a total of Ndat=160subscript𝑁𝑑𝑎𝑡160N_{dat}=160italic_N start_POSTSUBSCRIPT italic_d italic_a italic_t end_POSTSUBSCRIPT = 160 data points. As is evident from this table, the individual χ2superscript𝜒2\chi^{2}italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT values for each experiment exhibit a consistent and stable pattern, and the overall total values indicate a good agreement between the data and the theoretical predictions.

As discussed earlier, an additional analysis was performed where the A2subscript𝐴2A_{2}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT asymmetry observables (data in Table. 2) were excluded from the QCD fit. The corresponding χ2superscript𝜒2\chi^{2}italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT values per data point are presented in the last column of Tables. 1 and 3. The numbers in these tables suggest that the exclusion of A2subscript𝐴2A_{2}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT asymmetry data could potentially increase the total values of χ2superscript𝜒2\chi^{2}italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT. Table. 1 shows that the inclusion of A2subscript𝐴2A_{2}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT data increases the total value of χ2superscript𝜒2\chi^{2}italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT from 1170.172 to 1304.509. Similarly, Table. 3 indicates an increase in the total value of χ2superscript𝜒2\chi^{2}italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT from 227.934 to 244.520.

In summary, the global χ2/Ndatsuperscript𝜒2subscript𝑁𝑑𝑎𝑡\chi^{2}/N_{dat}italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / italic_N start_POSTSUBSCRIPT italic_d italic_a italic_t end_POSTSUBSCRIPT value for the entire dataset of 1266 data points in the AAK24 global QCD analysis is calculated to be 1.0528. This value indicates a stable and acceptable fit quality, demonstrating that the theoretical model aligns well with the experimental data. However, when the A2subscript𝐴2A_{2}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT asymmetry observables are excluded, the χ2/Ndatsuperscript𝜒2subscript𝑁𝑑𝑎𝑡\chi^{2}/N_{dat}italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / italic_N start_POSTSUBSCRIPT italic_d italic_a italic_t end_POSTSUBSCRIPT value increases to 1.0647. This result underscores the significance of including the A2subscript𝐴2A_{2}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT asymmetry data in QCD analyses of polarized PDFs. It suggests that the inclusion of A2subscript𝐴2A_{2}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT data leads to a better-constrained set of polarized PDFs with improved agreement between the data and theory, emphasizing the importance of these measurements in refining our understanding of polarized PDFs.

We can now proceed to compare the NLO theory predictions obtained for the asymmetry observables used in this study with the experimental datasets, particularly the newly added polarized data from COMPASS16 COMPASS:2015mhb and COMPASS17 COMPASS:2016jwv experiments, focusing on the A1psuperscriptsubscript𝐴1𝑝A_{1}^{p}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT and A1dsuperscriptsubscript𝐴1𝑑A_{1}^{d}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT asymmetries. Such comparisons are displayed in Figs. 7 and  8 for the A1psuperscriptsubscript𝐴1𝑝A_{1}^{p}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT and A1dsuperscriptsubscript𝐴1𝑑A_{1}^{d}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT asymmetry data from the COMPASS16 and COMPASS17 experiments, respectively. The comparisons are made as functions of x𝑥xitalic_x for two different selections of Q2superscript𝑄2Q^{2}italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT, specifically 15 and 13 GeV2. As these plots illustrate, the predictions from the AAK24 polarized PDFs are in excellent agreement with both the COMPASS16 and COMPASS17 datasets across the entire range of x𝑥xitalic_x. These results align with the individual χ2superscript𝜒2\chi^{2}italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT values presented in Table. 1, further confirming the quality of the fit and the compatibility between the theoretical model and the experimental data.

In Figs. 9 and 10, we present the polarized proton structure functions xg2p𝑥subscriptsuperscript𝑔𝑝2xg^{p}_{2}italic_x italic_g start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT calculated using the Wandzura-Wilczek relation in Eq. 6. The results from NMA23 (dashed dashed dotted) Nematollahi:2023dvj and KTA17 (dashed dotted) Khanpour:2017cha are also shown. The experimental data from the E155 collaboration E155:2002iec are included for comparison. As can be seen, the xg2p𝑥subscriptsuperscript𝑔𝑝2xg^{p}_{2}italic_x italic_g start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT structure functions calculated using the Wandzura-Wilczek approximation provide good agreement with both the data and other results available in the literature.

To provide additional support for our discussion regarding the comparison between data and theory, we present detailed comparisons of our theoretical predictions with a different selection of A1psuperscriptsubscript𝐴1𝑝A_{1}^{p}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT asymmetry data as outlined in Table. 1. As one can see from Fig. 11, the agreement between the theoretical model and various selections of A1psuperscriptsubscript𝐴1𝑝A_{1}^{p}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT asymmetry data is good. Such agreement between the theoretical predictions and a variety of experimental datasets for A1psuperscriptsubscript𝐴1𝑝A_{1}^{p}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT asymmetry provides evidence of the validity and reliability of the AAK24 analysis in describing the polarized PDFs in nucleons.

Refer to caption
Figure 7: Comparison of the AAK24 theory prediction with the A1psuperscriptsubscript𝐴1𝑝A_{1}^{p}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT asymmetry data from COMPASS16 experiment COMPASS:2015mhb .
Refer to caption
Figure 8: Same as Fig. 7 but this time for the A1dsuperscriptsubscript𝐴1𝑑A_{1}^{d}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT asymmetry data from COMPASS17 experiment COMPASS:2016jwv .
Refer to caption
Figure 9: The polarized proton structure functions xg2p𝑥subscriptsuperscript𝑔𝑝2xg^{p}_{2}italic_x italic_g start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT as a function of x𝑥xitalic_x at Q2=5superscript𝑄25Q^{2}=5italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = 5 GeV2. The results from NMA23 (dashed dashed dotted) Nematollahi:2023dvj and KTA17 (dashed dotted) Khanpour:2017cha are also shown. The experimental data from the E155 collaboration E155:2002iec are presented for comparison.
Refer to caption
Figure 10: Same as Fig. 9, but this time the comparisons are shown for Q2superscript𝑄2Q^{2}italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT=3 GeV2.
Refer to caption
Figure 11: Comparison of the AAK24 theory prediction with some selection of A1psuperscriptsubscript𝐴1𝑝A_{1}^{p}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT, A1nsuperscriptsubscript𝐴1𝑛A_{1}^{n}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT, and A1dsuperscriptsubscript𝐴1𝑑A_{1}^{d}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT asymmetry data presented in Table. 1.

We are now in a position to compare the theoretical predictions from AAK24 with the corresponding experimental results for the SIDIS asymmetries. This comparison is illustrated in Fig. 12 for selected data from SMC and COMPASS collaborations. Our calculations include the theory predictions for proton, neutron, and deuteron, taking into account the light charged hadron, kaon, and pion FFs. As shown in the figures, the agreement between the data and theory is excellent across the entire range of x𝑥xitalic_x and for all hadron species. This alignment corroborates the favorable χ2superscript𝜒2\chi^{2}italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT values presented in Table. 3. It further confirms the consistency between the theoretical framework and the experimental data, affirming the model ability to accurately describe the observed phenomena in the realm of polarized SIDIS asymmetries.

Refer to caption
Figure 12: Comparison of AAK24 QCD results for the SIDIS asymmetries with the corresponding SIDIS observables. The calculations include the theory predictions for proton, neutron, and deuteron, taking into account the light charged hadron, kaon, and pion FFs.

To provide more detailed discussions of the results, we provide a comparison in Fig. 13 between the AAK24 theoretical predictions for the A2subscript𝐴2A_{2}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT asymmetry for proton, neutron, and deuteron targets and the corresponding data analyzed in this study. As mentioned earlier, this marks the first inclusion of A2subscript𝐴2A_{2}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT asymmetry data in a global QCD analysis of polarized PDFs. Upon inspection, the data/theory agreement is generally good for most experiments. However, some deviations between the data and theory are evident. For example, discrepancies can be observed for the A2dsuperscriptsubscript𝐴2𝑑A_{2}^{d}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT asymmetry from the SMC experiment and the A2psuperscriptsubscript𝐴2𝑝A_{2}^{p}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT asymmetry from HERMES.

Refer to caption
Figure 13: Comparison of AAK24 QCD results for the A2subscript𝐴2A_{2}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT asymmetries with a selected set of corresponding data. The calculations include the theory predictions for proton, neutron, and deuteron

In summary, the AAK24 results presented in this section, which integrated recent experimental data from JLAB17, COMPASS16, and COMPASS17, including the novel inclusion of A2subscript𝐴2A_{2}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT asymmetry measurements, exhibited favorable overall agreement with experimental data. These findings also align well with other polarized PDF determinations, notably JAM17, DSSV14, LSS10, and AKS14, all of which also accounted for SU(2) and SU(3) symmetry breaking.

7 Conclusions

In this article, we have presented AAK24 polarized PDFs, a next-to-leading order (NLO) QCD analysis of polarized data from both polarize DIS and SIDIS experiments on the nucleon. Throughout the analysis, we considered the scenario of SU(2) and SU(3) symmetry breaking, specifically δu¯δd¯δs¯𝛿¯𝑢𝛿¯𝑑𝛿¯𝑠\delta\bar{u}\neq\delta\bar{d}\neq\delta\bar{s}italic_δ over¯ start_ARG italic_u end_ARG ≠ italic_δ over¯ start_ARG italic_d end_ARG ≠ italic_δ over¯ start_ARG italic_s end_ARG, while treating δs¯𝛿¯𝑠\delta\bar{s}italic_δ over¯ start_ARG italic_s end_ARG and δs𝛿𝑠\delta sitalic_δ italic_s as equal since the available experimental data are not enough to distinguish them. Additionally, we have included the nonperturbative target mass corrections (TMC) as well as higher-twist terms (HT) which are particularly important for the polarized DIS. Our study addressed the crucial role of semi-inclusive data in determining the polarized sea quark distributions. The recent experimental data from JLAB17, COMPASS16, and COMPASS17 experiments have been incorporated into our data samples. Additionally, this analysis marks simultaneous incorporation of SIDIS and A2subscript𝐴2A_{2}italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT asymmetry measurements as well.

We have thoroughly examined the impact of these new data sets on both the central values and the extracted uncertainties, as well as their effect on the overall fit quality. Our findings are presented and discussed to provide a comprehensive understanding of their implications. To quantify the uncertainties in the polarized PDFs, we employed the standard Hessian method. The main results of AAK24 QCD analysis exhibit an overall good agreement with the analyzed experimental data, aligning well with other polarized PDFs determinations which considered the SU(2) and SU(3) symmetry breaking as well.

A Fortran package containing AAK24 polarized PDFs at NLO accuracy can be obtained via e-mail from the authors.

Data Availability Statement

This manuscript has associated data in a data in a repository. [Authors comment: This data sets generated during the current study are available from the corresponding author on responsible request.]

Acknowledgments

The authors appreciate the financial support from the Iran National Science Foundation (INSF) under grant number 4013570. S. A. T. and H. K. are also grateful to the School of Particles and Accelerators, Institute for Research in Fundamental Sciences (IPM). F. A. acknowledges the Farhangian University for the provided support to conduct this research. H. Khanpour appreciates the financial support from NAWA under grant number BPN/ULM/2023/1/00160 and from the IDUB programme at the AGH University. We also thank Dimiter Stamenov for providing the necessary LSS10 grids, and Rodolfo Sassot for supplying the required DSSV14 grids data.

References

  • (1) A. Deur, S. J. Brodsky and G. F. De Téramond, Rept.Prog.Phys. 82 (2019) 076201
  • (2) C. A. Aidala, S. D. Bass, D. Hasch and G. K. Mallot, Rev. Mod. Phys. 85, 655-691 (2013),
  • (3) J. J. Ethier and E. R. Nocera, Ann. Rev. Nucl. Part. Sci. 70, 43-76 (2020)
  • (4) R. Abdul Khalek, A. Accardi, J. Adam, D. Adamiak, W. Akers, M. Albaladejo, A. Al-bataineh, M. G. Alexeev, F. Ameli and P. Antonioli, et al. Nucl. Phys. A 1026, 122447 (2022),
  • (5) D. P. Anderle, T. J. Hou, H. Xing, M. Yan, C. P. Yuan and Y. Zhao, JHEP 08, 034 (2021),
  • (6) C. Adolph et al. [COMPASS], Phys. Lett. B 753, 18-28 (2016)
  • (7) D. De Florian, G. A. Lucero, R. Sassot, M. Stratmann and W. Vogelsang, Phys. Rev. D 100, no.11, 114027 (2019)
  • (8) T. Khan, T. Liu and R. S. Sufian, Phys. Rev. D 108, no.7, 074502 (2023)
  • (9) D. Adamiak et al. [Jefferson Lab Angular Momentum], Phys. Rev. D 104, no.3, L031501 (2021)
  • (10) F. Arbabifar, S. Atashbar Tehrani and H. Khanpour, Phys. Rev. C 108, no.3, 035203 (2023)
  • (11) F. Taghavi-Shahri, H. Khanpour, S. Atashbar Tehrani and Z. Alizadeh Yazdi, Phys. Rev. D 93, no.11, 114024 (2016)
  • (12) M. Hirai et al. [Asymmetry Analysis], Nucl. Phys. B 813, 106-122 (2009)
  • (13) E. Leader, A. V. Sidorov and D. B. Stamenov, Phys. Rev. D 91, no.5, 054017 (2015)
  • (14) H. Khanpour, S. T. Monfared and S. Atashbar Tehrani, Phys. Rev. D 95, no.7, 074006 (2017),
  • (15) N. Sato et al. [Jefferson Lab Angular Momentum], Phys. Rev. D 93, no.7, 074005 (2016),
  • (16) J. Blumlein and H. Bottcher, Nucl. Phys.  B 841 (2010) 205
  • (17) D. de Florian, R. Sassot, M. Stratmann and W. Vogelsang, Phys. Rev.  D 80 (2009) 034030
  • (18) D. Adamiak, N. Baldonado, Y. V. Kovchegov, W. Melnitchouk, D. Pitonyak, N. Sato, M. D. Sievert, A. Tarasov and Y. Tawabutr, [arXiv:2308.07461 [hep-ph]].
  • (19) F. Arbabifar, A. N. Khorramian and M. Soleymaninia, Phys. Rev. D 89, no.3, 034006 (2014)
  • (20) H. Khanpour, S. T. Monfared and S. Atashbar Tehrani, Phys. Rev. D 95, no.7, 074006 (2017)
  • (21) H. Nematollahi, A. Mirjalili and S. Atashbar Tehrani, Phys. Rev. D 107, no.5, 054033 (2023)
  • (22) D. de Florian, R. Sassot, M. Stratmann and W. Vogelsang, Phys. Rev. Lett. 113, no.1, 012001 (2014)
  • (23) E. Leader, A. V. Sidorov and D. B. Stamenov, Phys. Rev. D 82, 114018 (2010)
  • (24) N. Sato et al. [Jefferson Lab Angular Momentum], Phys. Rev. D 93, no.7, 074005 (2016)
  • (25) C. Amsler et al. (Particle Data Group), Phys. Lett. B 667 (2008) 1.
  • (26) H. J. Lipkin, Phys. Lett.  B 214, 429 (1988); Phys. Lett.  B 230, 135 (1989); F. E. Close and R. G. Roberts, Phys. Rev. Lett.  60, 1471 (1988); M. Roos, Phys. Lett.  B 246, 179 (1990); Z. Dziembowski and J. Franklin, J. Phys. G 17, 213 (1991). P. G. Ratcliffe, Phys. Lett.  B 242, 271 (1990); Phys. Lett.  B 365, 383 (1996); arXiv:hep-ph/0012133; S. L. Zhu, G. Sacco, and M. J. Ramsey-Musolf, Phys. Rev.  D 66, 034021 (2002); E. Leader and D. B. Stamenov, Phys. Rev.  D 67, 037503 (2003).
  • (27) C. Adolph et al. [COMPASS], Phys. Lett. B 769, 34-41 (2017)
  • (28) R. Fersch et al. [CLAS], Phys. Rev. C 96, no.6, 065208 (2017)
  • (29) B. Lampe and E. Reya, Phys. Rept. 332, 1-163 (2000),
  • (30) Y. B. Dong, Phys. Lett. B 641, 272 (2006).
  • (31) Y. B. Dong, Phys. Rev. C 78, 028201 (2008)
  • (32) Y. B. Dong, Phys. Lett. B 653, 18 (2007)
  • (33) Y. B. Dong and D. Y. Chen, Nucl. Phys. A 791, 342 (2007).
  • (34) A. V. Sidorov and D. B. Stamenov, Mod. Phys. Lett. A 21, 1991 (2006)
  • (35) O. Nachtmann, Nucl. Phys. B 63, 237 (1973).
  • (36) A. Piccione and G. Ridolfi, Nucl. Phys. B 513, 301-316 (1998)
  • (37) A. Courtoy, A. S. Miramontes, H. Avakian, M. Mirazita and S. Pisano, Phys. Rev. D 106, no.1, 014027 (2022),
  • (38) A. Mirjalili and S. Atashbar Tehrani , Phys. Rev. D 105, no.7, 074023 (2022)
  • (39) L. W. Whitlow, S. Rock, A. Bodek, E. M. Riordan and S. Dasu, Phys. Lett. B 250, 193-198 (1990)
  • (40) M. Arneodo et al. [New Muon], Phys. Lett. B 364, 107-115 (1995)
  • (41) D. de Florian, M. Stratmann and W. Vogelsang, Phys. Rev. D 57, 5811-5824 (1998),
  • (42) A. D. Martin, R. G. Roberts, W. J. Stirling and R. S. Thorne, Eur. Phys. J. C 28, 455-473 (2003)
  • (43) T. J. Hou, J. Gao, T. J. Hobbs, K. Xie, S. Dulat, M. Guzzi, J. Huston, P. Nadolsky, J. Pumplin and C. Schmidt, et al. Phys. Rev. D 103, no.1, 014013 (2021)
  • (44) D. de Florian, R. Sassot, and M. Stratmann, Phys. Rev. D 75, 114010 (2007) Phys. Rev. D 76, 074033 (2007)
  • (45) J. Ashman et al. [European Muon Collaboration], J. Ashman et al. [European Muon Collaboration], Nucl. Phys. B 328 (1989) 1.
  • (46) B. Adeva et al. [Spin Muon Collaboration], Phys. Rev. D 58 (1998) 112001.
  • (47) M. G. Alekseev et al. [COMPASS Collaboration], Phys. Lett. B 690, 466 (2010)V. Y. Alexakhin et al. [COMPASS Collaboration], Phys. Lett. B 647 (2007) 8
  • (48) P. L. Anthony et al. [E142 Collaboration], Phys. Rev. D 54 (1996) 6620
  • (49) K. Abe et al. [E143 collaboration], Phys. Rev. D 58 (1998) 112003
  • (50) K. Abe et al. [E154 Collaboration], Phys. Rev. Lett. 79 (1997) 26
  • (51) P. L. Anthony et al. [E155 Collaboration], Phys. Lett. B 493 (2000) 19
  • (52) P. L. Anthony et al. [E155 Collaboration], Phys. Lett. B 463 (1999) 339
  • (53) SMC Collaboration, B. Adeva et al., Phys. Lett. B 420 (1998) 180.
  • (54) K. Ackerstaff et al. [HERMES Collaboration], Phys. Lett. B 404 (1997) 383 A. Airapetian et al. [HERMES Collaboration], Phys. Lett. B 442 (1998) 484
  • (55) A. Airapetian et al. [HERMES Collaboration], Phys. Rev. D 75 (2007) 012007
  • (56) K. V. Dharmawardane et al. [CLAS Collaboration], Phys. Lett. B 641 (2006) 11
  • (57) X. Zheng et al. (JLab/Hall A Collaboration), Phys. Rev. Lett.92, 012004 (2004);
  • (58) K. Abe et al. [E143], Phys. Rev. Lett. 76, 587-591 (1996)
  • (59) P. L. Anthony et al. [E155], Phys. Lett. B 553, 18-24 (2003)
  • (60) P. L. Anthony et al. [E155], Phys. Lett. B 458, 529-535 (1999)
  • (61) K. Abe et al. [E154], Phys. Lett. B 404, 377-382 (1997)
  • (62) D. Adams et al. [Spin Muon (SMC)], Phys. Lett. B 396, 338-348 (1997)
  • (63) A. Airapetian et al. [HERMES], Eur. Phys. J. C 72, 1921 (2012)
  • (64) M.G. Alekseev et al. (COMPASS Collaboration), Phys. Lett. B660, 458 (2008).
  • (65) M.G. Alekseev et al. (COMPASS Collaboration), Phys. Lett. B680, 217 (2009).
  • (66) M. G. Alekseev et al. [COMPASS Collaboration], Phys. Lett.  B 693 (2010) 227
  • (67) F. James and M. Roos, Comput. Phys. Commun.  10, 343 (1975); F. James, “MINUIT Function Minimization and Error Analysis: Reference Manual Version 94.1,” CERN-D-506, CERN-D506.
  • (68) J. Pumplin, D. Stump, R. Brock, D. Casey, J. Huston, J. Kalk, H. L. Lai and W. K. Tung, Phys. Rev. D 65, 014013 (2001),
  • (69) A. Accardi, T. J. Hobbs, X. Jing and P. M. Nadolsky, Eur. Phys. J. C 81, no.7, 603 (2021),