The transverse momentum dependent (TMD) parton distribution functions (PDFs), also called as unintegrated PDFs, describe the probability of finding a quark with certain light-cone momentum fraction \(x\)  and transverse momentum \(k_T\)  in a hadron [1]}.
Among them the Sivers function [2]} is of particular interests.
Due to the time-reversal-odd (T-odd) spin-orbit correlation of partons inside the nucleon,
the Sivers function describes the distortion in the distribution of unpolarized partons inside a transversely polarized nucleon.
Thereby it can give rise to the single-spin asymmetries (SSAs) in leading-twist in various semi-inclusive high energy processes.
During the last two decades, the SSAs related to Sivers function in semi-inclusive deep inelastic scattering (SIDIS) have been measured by the HERMES [3]}, [4]}, COMPASS [5]}, [6]}, [7]}, [8]}, and JLab Hall A [9]}, [10]} Collaborations.
Recently, the measurement on the Sivers effects in Drell-Yan and W/Z production have also been performed by COMPASS [11]} and STAR [12]}.
The data from these experiments were utilized by
different groups [13]}, [14]}, [15]}, [16]}, [17]}, [18]}, [19]}, [20]}, [21]}, [22]}, [23]}, [24]}
to extract the quark Sivers functions of the proton without or with TMD evolution.
On the other hand, various model calculations have been applied to calculate the quark Sivers functions, such as the the spectator model [25]}, [26]}, [27]}, the light-cone quark model [28]}, [29]}, the non-relativistic constituent quark model [30]}, and the MIT bag model [31]}, [32]}.
