The Sivers asymmetry in charged Kaon and $\Lambda$ hyperon produced SIDIS process at electron ion colliders

We investigate the single transverse-spin asymmetry with a $\sin (\phi_h-\phi_S)$ modulation in the charged Kaon produced and in $\Lambda$ hyperon produced SIDIS process within the theoretical framework of transverse momentum dependent (TMD) factorization at the next-to-leading-logarithmic order. The asymmetry is contributed by the convolution of Sivers function and the unpolarized fragmentation function $D_1$ for the produced hadron. The parametrization for the proton Qiu-Sterman function, which is closely related to the Sivers function, is adopted to numerically estimate the Sivers asymmetry at the kinematical region of Electron Ion Collider (EIC) and Electron Ion Collider in China (EicC). The TMD evolution of the TMD parton distribution functions are considered by employing the nonperturbative Sudakov form factor. It is found that the predicted Sivers asymmetries $A_{UT}^{\sin (\phi_h-\phi_S)}$ as functions of $x$, $z$ and $P_{hT}$ are sizable at the kinematical configurations of both EIC and EicC. The strange constituent of the produced charged Kaon and $\Lambda$ hyperon in the final state can be a promising probe of the sea quark Sivers function as well as the flavor dependence in the proton target. Therefore, it is important to utilize the future EIC facilities to constrain the sea quark distribution functions as well as the validity of the generalized universality of the Sivers function.


I. INTRODUCTION
Since the measurement by the European Muon Collaboration [1,2] showed that the spin fraction carried by the internal quarks is much smaller than the spin of the proton, which is contradict to the conventional theoretical prediction that the constituent quark spin contributes the total proton spin, numerous studies have been carried out to explore the nucleon spin structure from both theoretical and experimental aspects.Among the spin-related observables, the transverse single spin asymmetries (TSSAs) can be the key access to the information of transverse momentum structure of nucleon, which is encoded in the transverse momentum dependent parton distribution functions (TMD PDFs).In leading twist there are eight TMD PDFs, each one describes a distribution of three-dimensional motion of partons with specified polarization inside the nucleon.Particularly, the time-reversal odd (T-odd) Sivers function f ⊥ 1T (x, p T ) [3,4] denotes the asymmetric distribution of unpolarized quarks inside a transversely polarized nucleon, which arises from the correlation between the quark transverse momentum and the nucleon transverse spin.Due to its T-odd property, Sivers function as well as its chiral-odd partner the Boer-Mulders function has been assumed to be forbidden by the naive time-reversal invariance of QCD [5], the very existence of the two T-odd distribution functions was not so obvious.However, the situation has changed since the calculations in Refs.[6][7][8], which showed that the T-odd distributions can actually survive using spectator model calculations incorporating gluon exchange between the struck quark and the spectator.In Ref. [9], the time-reversal-invariant argument was reexamined and showed the gauge-link in the operator definition of the correlators guarantee the T-odd distribution functions to be nonzero.More importantly, the presence of the gauge-link indicates that Sivers function or the Boer-Mulders function has opposite sign between semi-inclusive deeply inelastic scattering (SIDIS) and Drell-Yan processes [6,7,9] which is a significant prediction by QCD.The verification of the sign change is one of the most fundamental tests of QCD prediction, and it is also the main pursue of the existing and future Drell-Yan facilities.
The transverse single spin asymmetry can be utilized to extract the information of the Sivers function, and has been intensively investigated in the past two decades from both experimental and theoretical aspects.The first non-zero Sivers asymmetry was measured by the HERMES Collaboration at DESY in electroproduction of charge pions off the transversely polarized hydrogen target [10].Updated measurements on the Sivers asymmetry in pion produced as well as those in Kaon and p/p produced in three-dimensional kinematic bin and enlarged phase space were reported in Refs.[11,12].The COMPASS Collaboration at CERN also measured the Sivers asymmetries in charged hadrons produced SIDIS process through muon beam scattering off the transversely polarized proton and deuteron targets [13][14][15][16][17][18].In addition, the data on the weighted Sivers asymmetry are also released in Ref. [19], allowing for the extractions of the Sivers function and its first transverse moment.The Hall A Collaboration at Jefferson Lab presented the measurement of TSSA in charged pion produced SIDIS process with a transversely polarized 3 He target [20,21].Besides the SIDIS process, COMPASS also measured the Sivers asymmetry in Drell-Yan process via π N collision [22].Measurement of TSSAs sensitive to Sivers function in the W ± boson produced in proton-proton collisions has also been performed by the STAR experiments at RHIC [23].The data from these measurements have been applied to extract the Sivers function using parametrizations and phenomenological approaches [24][25][26][27][28][29][30][31][32][33][34][35].
From the theoretical aspect, intensive studies on the quark Sivers function were performed using the QCD-inspired models, such as the spectator model [6,8,36], the light-cone quark model [37,38], the light-front quark-diquark model [39,40], the non-relativistic constituent quark model [41], the MIT bag model [42,43], and the Holographic QCD [44].The sea quark Sivers function has been estimated from the light-cone wave function in Refs.[45][46][47].However, so far only the valence quark Sivers functions are constrained in the valence region with relatively large uncertainties in the transverse momentum space.One of the reasons is that it is difficult to describe the corresponding physical observables since there are complicated TMD effects.The TMD evolution effects are encoded in the Sudakovlike form factor which also includes details of the non-perturbative QCD dynamics.Therefore, this part of the Sudakovlike form factor can not be calculated from the perturbative QCD and is mostly unknown.Another reason is that the knowledge of the scale dependence of the TSSAs is very limited since the measurements are mostly performed in the fixed-target experiments with similar hard scales.With the expected high energy and high precision of the Electron Ion Collider (EIC) [48,49] and the Electron Ion Collider in China (EicC) [50], the precise knowledge on the TMD distribution functions may be gained, not only for the valence quarks, but also for sea quarks and gluons.Concerning the sea quark TMDs, the charged Kaon produced or the Lambda hyperon produced in SIDIS can be recognized as an ideal probe to the sea quark distribution of the proton due to the strange constituent quark inside the kaon and the Lambda hyperon.Therefore, through the Sivers asymmetry in K ± produced and in Lambda produced off transversely polarized nucleon at EIC and EicC, there might be an opportunity to obtain the information of the Sivers distribution function of the sea quark as well as its flavor dependence.
The purpose of this work is to evaluate the Sivers asymmetry in ep ↑ → eK ± X and in ep ↑ → eΛX at the kinematical region of EIC and EicC.The theoretical tool adopted in this study is the TMD factorization formalism [51][52][53][54], which has been widely applied to various high energy processes, such as SIDIS [51,53,[55][56][57][58], e + e − annihilation [53,59,60], Drell-Yan [53,61], and W/Z produced in hadron collision [52,53,62].In this framework, the differential cross section can be written as the convolution of the well-defined TMD PDFs and/or fragmentation functions (FFs) at the small transverse momentum region P hT ≪ Q as an approximation.The energy dependence of the TMD PDFs and FFs is encoded in the TMD evolution equations, their solutions are usually given in b space, which is conjugate to the transverse momentum space [52,53] through Fourier transformation.After solving the TMD evolution equations, the scale dependence of the TMDs may be included in the exponential form of the so-called Sudakov-like form factor [52,53,56,63].The Sudakov-like form factor can be further separated into perturbatively calculable part and the nonperturbative part, the latter one can not be calculated through perturbative theory and may be obtained by fitting experimental data.We will consider the TMD evolution effects of the corresponding TMDs to obtain the numerical results for the Sivers asymmetry in charged Kaon K ± and Λ hyperon produced SIDIS process.
The rest of the paper is organized as follows.In Sec.II, we provide the theoretical framework of Sivers asymmetry A sin(φ h −φS ) UT in the charged Kaon produced and Λ hyperon produced in SIDIS process within the TMD factorization formalism.In Sec.III, we perform the numerical estimation of the Sivers asymmetry at the kinematical region of EIC and EicC.In Sec.IV, we summarize the work and discuss the results.

II. FORMALISM OF THE SIVERS ASYMMETRY IN SIDIS PROCESS
In this section, we will set up the detailed formalism of the Sivers asymmetry with a modulation of sin(φ h − φ S ) in the charged Kaon or Lambda produced SIDIS process with an unpolarized electron beam scattered off a transversely polarized proton target e(ℓ) + p(P where ℓ and ℓ ′ represent the four-momenta of the incoming and outgoing electrons, P is the four-momentum of the target proton, the up-arrow represents the proton is transversely polarized, P h stands for the four-momentum of the final-state hadron, which can be charged Kaon K ± or Lambda hyperon.The four-momentum of the exchanged virtual photon is q = ℓ − ℓ ′ and the usual defined energy scale is Q 2 = −q 2 .We denote the masses of the proton target and the final-state hadron by M and M h .To express the differential cross section as well as the physical observables, we adopt the following Lorentz invariants where S is the squared center of mass energy, x represents the Bjorken variable, y represents the lepton (quark) energy momentum transferring fraction, and z represents the longitudinal momentum fraction of the final fragmented hadron to the parent quark.
The reference frame applied in our study is shown in Fig. 1.According to the Trento convention [64], the z-axis is defined by the direction of the virtual photon momentum.The azimuthal angle φ h of the outgoing hadron (charged Kaon or Lambda) is defined by with ℓ µ T = g µν ⊥ ℓ ν and P µ hT = g µν ⊥ P hν being the transverse components of ℓ and P h respect to z-axis.The tensor g µν ⊥ is with γ = 2Mx Q .The azimuthal angle φ S of the proton spin vector S is defined by replacing P h by S in Eq. ( 4), and the transverse component of S is S µ T = g µν ⊥ S ν similar to the definition of P hT .Assuming one photon exchange, the model-independent differential cross section can be written as a set of structure functions with the general form as [65] where F UU,T stands for the unpolarized structure function, F is the transverse spin-dependent structure function contributed by the Sivers function, ǫ is the ratio of the longitudinal flux and the transverse flux of the photon which has the definition ǫ = 4 γ 2 y 2 and the ellipsis denotes other structure functions, which will not be considered in this work.The three subscripts in the structure functions F XY,Z stand for the polarization of the lepton beam (X), the target proton (Y ) and the virtual photon (Z) with U being unpolarized, T being transversely polarized.The Sivers asymmetry is defined as the ratio of the difference between the spin-dependent differential cross sections and the unpolarized differential cross section where The structure functions in Eq. ( 7) can be expressed as the convolution of the corresponding TMD PDFs and FFs as [65] Here, f 1 (x, p T ) is the unpolarized TMD PDF, and is the unpolarized TMD FF, which depends on the longitudinal momentum fraction z and the transverse momentum k T of the finalstate quark.ĥ = P hT |P hT | is the unit vector along P hT .The notation C represents the convolution among the transverse momenta Substituting Eq. ( 10) into Eq.( 8), we can expand the unpolarized structure function F UU as hereafter, the term with a tilde denotes it is in the b space.Similarly, substituting Eq. ( 10) into Eq.( 9), we can obtain expansions for the spin-dependent structure function F The Sivers function in b space can also be obtained from momentum space f ⊥ 1 (x, p T ) to b space by Fourier Transformation as where T q,F (x, x) is the Qiu-Sterman (QS) function.Therefore, the Sivers asymmetry can be rewritten as One should note that the energy dependence of the TMD structure functions were not encoded in the above formalism, which will be studied in details in the following subsections.

A. TMD evolution effects
First we set up the basic formalism of the TMD evolution effects for TMD PDFs and FFs, which mainly serves to solve the energy dependence of the TMD PDFs f 1 (x, p T ), f ⊥ 1T (x, p T ) and the TMD FF D 1 (z, k T ).Since the complicated convolution among the transverse momenta can be transformed into a simple product after performing the Fourier Transformation, it is convenient to solve the energy dependence in the b space.
Particularly, there are two different energy dependencies while the µ dependence is given by the renormalization group equation with α s being the running strong coupling at the energy scale µ, K being the CS evolution kernel, and γ K , γ F and γ D being the anomalous dimensions.Hereafter, we will assume Solving these TMD evolution equations, one can obtain the solution of the energy dependence for TMD parton distribution functions and the fragmentation functions, of which the solution has the general form as where F and D is the factor related to the hard scattering and depend on the factorization schemes, S is the Sudakovlike form factor. Eq. ( 21) and Eq.(22) show that the energy evolution of TMD PDFs and TMD fragmentation functions from an initial energy µ B to another energy Q is encoded in the Sudakov-like form factor S by the exponential form exp(−S).By performing the reverse Fourier Transformation of the TMDs in b space, the TMDs in momentum space can be obtained, thus it is of great importance to study the b space behavior of the TMDs.In the small b region (b ≪ 1/Λ QCD ), the b dependence of TMDs is perturbative and can be calculated by perturbative QCD.However, the dependence in large b region turns to be nonperturbative, since the operators are separated by a large distance.
To include the evolution effect in this region, a nonperturbative Sudakov-like form factor S NP is introduced and is usually given in a parameterized form.The parameters of S NP can be determined by analyzing experimental data, given the lack of non-perturbative calculations.In order to combine the information from both the small b region and the large b region, a matching procedure is applied with a parameter b max serving as the boundary between the two regions.
as given in the original CSS prescription [52].The prescription also allows for a smooth transition from perturbative to nonperturbative regions and avoids the Landau pole singularity in α s (µ B ).The typical value of b max is chosen around 1.5 GeV −1 to guarantee that b * is always in the perturbative region.With the constraint of b * , we can calculate TMDs within a small b region.
In the small b region, the TMDs can be expressed as the convolutions of the perturbatively calculable hard coefficients and the corresponding collinear counterparts at fixed energy µ B , which could be the collinear PDFs/FFs or the multiparton correlation functions [51,67] where µ B = c 0 /b * and c 0 = 2e −γE and the Euler constant γ E ≈ 0.577 [51], the ⊗ stands for the convolution in the momentum fraction x Ĉj←q where C coefficients in the formula has different values in different processes, and its specific value will be given in the subsequent calculation.In addition, the sum i runs over all parton flavors.Now we can combine all the above information to get the expression for TMD distribution function and the fragmentation function in b space as The Sudakov-like form factor S can be separated into the perturbatively calculable part S pert (Q; b * ) and the nonper- The perturbative part S pert (Q; b * ) has a general form and can be expanded as the series of ( αs π ) [31,58,[68][69][70].
In Eq. ( 31) the coefficients A and B can be expanded as following In our calculation we take A (n) up to A (2) and B (n) up to B (1) [52, 56, 58, 68, 71, 72], For the non-perturbative Sudakov-like form factor S NP (Q; b), it cannot be obtained from perturbation calculation, and it is usually extracted from the experimental data.Inspired by Refs.[71,73], a widely used parametrization (the EIKV parametrization) of S NP for TMD PDFs or fragmentation functions was proposed [31,67,71,[73][74][75]] where the initial energy is Q 2 0 = 2.4 GeV 2 , and the factor 1/2 in front of g 2 comes from the fact that only one hadron is involved for the parametrization of S pdf/ff NP .The parameter g pdf/ff 1 in Eq. ( 37) depends on the type of TMDs, which can be regarded as the width of the intrinsic transverse momentum for the relevant TMDs at the initial energy scale Q 0 [29,56,72].Assuming a Gaussian form for the dependence for the transverse momentum, one can obtain where p 2 ⊥ Q0 and k 2 ⊥ Q0 represent the averaged intrinsic transverse momenta squared for TMD PDFs and FFs at the initial scale Q 0 , respectively.g 2 = 0.16 [31] has the same form in all forms of TMDs.
As the information on the Sudakov-like form factor for the Kaon fragmentation functions can not be determined, we assume that the evolution of the TMD distribution function and the fragmentation function for producing the K meson from the initial energy scale µ to another energy scale Q follows the Gaussian form of g 2 (b) in Eq. ( 37), so we can obtain the nonperturbative Sudakov-like form factor for the PDF and FF in the produced of the K mesons as Combining all the steps mentioned above together, the scale-dependent TMD PDFs and FFs in b space as functions of x (or z), b, and Q can be rewritten as B. The solution of the unpolarized structure function In the following we solve the denominator of the Sivers asymmetry in details, which is the unpolarized structure function F UU in Eq. ( 8).Since we have expanded F UU in Eq. ( 11), we will directly follow Eq. ( 11) to give the complete expression for F UU .We can divide Eq. ( 11) into three parts to solve the problem: the weighted Bessel function, the distribution function f q/p 1 (x, b) and the fragmentation function Dh/q 1 (z, b).The weighted Bessel function can be obtained as According to Eq. ( 41) and Eq. ( 42), f q 1 (x, b) and Dh/q 1 (z, b) can be expanded as where the Sudakov-like form factor is expressed in Eq. ( 31), Eq. ( 39) and Eq. ( 40).The hard scattering coefficients F (α s (Q)) and D(α s (Q)) can be set equal to 1.In addition, the C coefficient in Eq. ( 44) and Eq. ( 45) can be expressed as [76] C (SIDIS) where α s represents the strong coupling, and the expansion at the next-to-leading order can be written as follows In Eq. ( 50), Q 2 is the running energy scale and n f = 5, Λ QCD = 0.225 GeV.The splitting functions P q←q and P g←q in Eq. ( 48) and Eq. ( 49) have the general form where C F = 4/3, T R = 1/2, and the subscript symbol "+" denotes the following prescription Combining the Bessel function part, the PDF part f q/p 1 (x, b) and the FF part Dq 1 (z, b), we can expand the denominator of Sivers asymmetry as where the nonperturbative Sudakov-like form factor S SIDIS NP (Q; b) receives the contribution from the unpolarized TMD PDF and FF: Here, g f1 1 and g D1 1 are obtained from Eq. ( 38), so C. The solution of the transverse spin-dependent structure function Similar to the previous subsection, we further obtain the expression for the transverse spin-dependent structure function F sin(φ h −φs) UT .Eq. ( 14) can be separated to three parts to solve it, namely, the weighted Bessel function, the Qiu-Sterman (QS) function T q,F (x ′ , x ′′ ) and the fragmentation function Dh/q 1 (z, b).The solution for the weighted Bessel function is similar to Eq. ( 43) The QS function part T q,F (x ′ , x ′′ ) in the b space has a similar solution form after solving the TMD evolution equations In the small b region, the QS function can be calculated from the convolution of the hard scattering coefficients and the collinear counterpart utilizing the perturbative QCD where ⊗ denotes convolution and can be expended as thus the QS function T q,F (x ′ , x ′′ ) and the fragmentation function Dh/q 1 (z, b) can be expanded as where the Sudakov-like form factor is expressed in Eq. ( 31), Eq. ( 39) and Eq.(40).And the hard scattering coefficient 48), ( 49), ( 51) and (52).The expression of ∆C in Eq. ( 60) is [77 where N c = 3. Combining all the ingredients together, the numerator part of the Sivers asymmetry can be rewritten as where the nonperturbative Sudakov-like form factor S SIDIS NP Sivers (Q; b) is the combination of and the one for the Sivers function the one for the unpolarized TMD FF where g Sivers

III. NUMERICAL ESTIMATE
In this section, we present the numerical estimate for the Sivers asymmetry in charged Kaon and Λ hyperon produced SIDIS process at the kinematical configurations of EIC and EicC.
In order to obtain the results for the Sivers asymmetry, we shall have the collinear distribution functions and fragmentation functions as input.For the unpolarized proton collinear distribution function f 1 (x, µ B ), we apply the parametrization from MSTW2008 [78].For the collinear unpolarized fragmentation function D 1 (z) of charged Kaon, the DSS parameterization results at NLO accuracy is used [79].For the collinear unpolarized fragmentation function D Λ 1 (z) of Λ hyperon, we adopt the model results from the diquark spectator model [80] The values of the free parameters in Eq. ( 65) are taken from Ref. [80].Since the model result is obtained at the initial energy of 0.23 GeV 2 , to make it applicable to a more general energy range, we use the QCDNUM evolution package [81] to evolve the unpolarized fragmentation function D 1 (z) from the initial energy of 0.23 GeV 2 to another energy scale.
In our calculation we adopt the Qiu-Sterman (QS) function T q,F (x, x; µ B ) from the parameterization in Ref. [31] T q,F (x, x, µ B ) = N q (α q + β q ) (αq+βq) The parameters obtained by fitting to the data in SIDIS processes from HERMES, COMPASS and JLab [11,15,17,20], are listed in  We should notice that in Ref. [31], the QS function is assumed to be proportional to the collinear unpolarized distribution function f 1 (x, µ B ) in which the DGLAP evolution effect is not included.To investigate the impact of the QS function DGLAP evolution effect on the Sivers asymmetry, we take Q 2 0 = 2.4 GeV 2 as the initial energy and evolve the QS function in Eq. ( 66) into another energy.We adopt two different approaches to evolve the QS function: one is to assume the QS function follows the same evolution effect as that for unpolarized distribution function, the other one is to change the evolution kernel in the QCDNUM evolution package to include the QS evolution kernel by considering the homogenous terms [the terms containing T q,F (x, x, µ B )] in the evolution kernel as an approximation [82]: In Eq. ( 37) and Eq. ( 38), the free parameter g 1 and the universal parameter g 2 contain information about the evolution of TMDs and are the key parameters that determine the evolution of TMDs from one initial energy µ to another Q.Here we adopt the results given in Ref. [31] for the mean transverse momentum squared For the universal parameter g 2 in the nonperturbative Sudakov-like form factor, the specific value g 2 = 0.16 is also given in Ref. [31].
The kinematical region that is available at EIC is chosen as follows [48] 0.001 < x < 0.4, 0.07 < y < 0.9, 0.2 < z < 0.8, As for the EicC, the following kinematical cuts are adopted 0.005 < x < 0.5, 0.07 <y < 0.9, 0.2 < z < 0.7, where W 2 = (P + q) 2 ≈ 1−x x Q 2 is invariant mass of the virtual photon-nucleon system.Since TMD factorization is proved to be valid to describe the physical observables in the region P hT ≪ Q, P hT < 0.5 GeV is chosen to guarantee the validity of TMD factorization.Combining Eq. ( 7), Eq. ( 54) and Eq. ( 63) and the kinematical regions of EIC and EicC, we can calculate the single-spin dependent Sivers asymmetry of charged Kaon and Λ hyperon produced SIDIS process within the EIC and EicC kinematical region.
The results are shown in Fig. 2. The six rows in the figure depict the results of Sivers asymmetry of K + production, K − production, and Λ hyperon production in the EIC and EicC kinematical region, respectively.The left, central and right panels denote the Sivers asymmetry as the function of x, z and P hT , respectively.In the figure, the black solid line represents the Sivers asymmetry that is obtained when the QS evolution kernel is included to evolve Eq. ( 66).The red dashed line represents the Sivers asymmetry obtained by using the evolution kernel for the unpolarized distribution function from QCDNUM evolution package to evolve Eq. ( 66).
As can be seen from Fig. 2, the Sivers asymmetries in the charged Kaon produced and in Λ hyperon produced SIDIS process are sizable in the EIC and EicC kinematical region.Thus, the measurements of the Sivers asymmetry in these facilities provide an ideal tool to obtain the sea quark Sivers function as well as the flavor dependence of the Sivers function.Also, it seems that the magnitude of the asymmetry in EicC is larger than that in EIC.In addition, Fig. 2 shows that the value of the solid line is larger than the dashed line, which indicates that the Sivers asymmetry obtained using the evolution kernel of QS function is greater than that obtained using the unpolarized distribution function evolution kernel.One should note that the DGLAP evolution in the TMD effects may play some role in the future phenomenological analysis of the Sivers asymmetry.
We can also find from Fig. 2 that both the z-dependent and P hT -dependent Sivers asymmetries in K + produced SIDIS process at the kinematical region of EIC and EicC are negative, and the magnitude of Sivers asymmetry increases with increasing z or increasing P hT .While for x dependent Sivers asymmetry of K + production in the EIC and EicC kinematical region, there is a node around x = 0.15.On the other hand, for the Sivers asymmetry of K − production, it is always positive in all cases.The Sivers asymmetry increases with x at small x region and decreases along x at large x region, and gradually increases with z and P hT .In addition, by comparing the magnitudes of the Sivers asymmetry of K + produced with that of K − produced, we can see that the former one is significantly larger than the latter one, since for K + meson the constituent quarks are u and s, and for K − meson are ū and s quarks.Therefore there is relatively large contribution of valence u quark Sivers function for K + , while the contributions for K − are both from sea quark Sivers function.It is known from Ref. [31] that the QS function of the valence quark is significantly larger than the that of the sea quarks, so the Sivers asymmetry in K + produced process is greater than that in K − produced process.Finally, the Sivers asymmetry of Λ hyperon production in both the EIC and EicC kinematics is positive, and the tendency of x, z and P hT dependent asymmetry is generally consistent with that of K − meson.In addition, the Sivers asymmetry in Λ hyperon produced process is larger than that in charged Kaon K ± produced process.
Therefore, the future higher precision EICs can provide a unique opportunity to extract the proton Sivers function of valence quark and sea quark from the charged Kaon K + production, to analyze the flavor dependence of the Sivers distribution function among sea quarks from K − produced SIDIS process, to investigate the flavor dependence from Λ hyperon produced SIDIS process.Combining the experimental data from Drell-Yan process, it also provides an ideal tool to study the sign change of the Sivers function and the generalized universality of the T-odd distribution functions.
where K T represents the transverse momentum of the final state hadron with respect to the fragmentation quark, which has the relation K T = −zk T with k T being the final-state quark transverse momentum respect to z axis.The δ-function Fourier Transformation was performed in the fourth line.The TMD distribution function f1 (x, b) and TMD fragmentation function D1 (z, b) in the b space can be obtained by performing a Fourier Transformation from momentum space to b space µ and ζ F (ζ D ) of the TMD PDF F (x, b) and the TMD FF D(z, b) in b space according to TMD factorization.µ is the renormalization scale related to the corresponding collinear PDFs/FFs, and ζ F (ζ D ) is the energy scale serving as a cutoff to regularize the light-cone singularity in the operator definition of the TMDs.The µ and ζ F (ζ D ) dependencies are encoded in different TMD evolution equations.The energy evolution for the ζ F (ζ D ) dependence is encoded in the Collins-Soper (CS) equation[52,53,66] then the TMD PDFs and FFs can be written as F (x, b; Q) and D(z, b; Q) for simplicity.
Furthermore, one can define a b-dependent function b * , which has the property b * ≈ b in small b region and b * ≈ b max in large b region b

FIG. 2 :
FIG. 2:The Sivers asymmetry in semi-inclusive charged Kaon and Λ hyperon produced SIDIS process at the kinematics of EIC and EicC as functions of x (left panels), z (middle panels), and P hT (right panels). Table.I.

TABLE I :
[31]es of the parameters in the parameterization of the Qiu-Sterman function in Ref.[31]