We apply the formalism developed earlier [1, 2] for studying transverse momentum dependent parton distribution functions (TMDs) at small Bjorken x to construct the small-x asymptotics of the quark Sivers function. First, we explicitly construct the complete fundamental “polarized Wilson line” operator to sub-sub-eikonal order: this object can be used to study a variety of quark TMDs at small x. We then express the quark Sivers function in terms of dipole scattering amplitudes containing various components of the “polarized Wilson line” and show that the dominant (eikonal) term which contributes to the quark Sivers function at small x is the spin-dependent odderon, confirming the re- cent results of Dong, Zheng and Zhou [3]. Our conclusion is also similar to the case of the gluon Sivers function derived by Boer, Echevarria, Mulders and Zhou [4] (see also [5]). We also analyze the sub-eikonal corrections to the quark Sivers function using the constructed “polarized Wilson line” operator. We derive new small-x evolution equations re-summing double-logarithmic powers of αs ln2(1/x) with αs the strong coupling constant. We solve the corresponding novel evolution equations in the large-Nc limit, obtaining a sub-eikonal correction to the spin-dependent odderon contribution. We conclude that the quark Sivers function at small x receives contributions from two terms and is given by$$ {f}_{1T}^{\perp q}\left(x,{k}_T^2\right)={C}_O\left(x,{k}_T^2\right)\frac{1}{x}+{C}_1\left({k}_T^2\right){\left(\frac{1}{x}\right)}^0+\cdots $$
f
1
T
⊥
q
x
k
T
2
=
C
O
x
k
T
2
1
x
+
C
1
k
T
2
1
x
0
+
⋯
with the function CO(x,$$ {k}_T^2 $$
k
T
2
) varying slowly with x and the ellipsis denoting the subasymptotic and sub-sub-eikonal (order-x) corrections.