In equation (9) of this paper, it is claimed that the limiting behaviour $$ \int_0^\infty \frac{1-\cos(kx_0)}{s+Dk^\alpha}dk \sim \frac{\Gamma(2-\alpha)\sin(\pi(2-\alpha)/2)x_0^{\alpha-1}}{(\alpha-1)D} \tag{1}\label{papereq} $$
holds as $s\to0$, when $1<\alpha<2$. However, by using the series expansion of $\cos(kx_0)$ $$ \cos(kx_0) = \sum_{n=0}^\infty\frac{(-1)^n(kx_0)^{2n}}{(2n)!} ~, $$ along with the identity (Appendix p152, The H-function with Applications in Statistics and Other Disciplines) $$ \frac{z^\beta}{1+az^\alpha}=a^{-\beta/\alpha}H_{1,1}^{1,1}\left[az^\alpha\Bigg|\begin{matrix}(\beta/\alpha,1)\\(\beta/\alpha,1)\end{matrix}\right] ~, $$ we get $$ \int_0^\infty \frac{1-\cos(kx_0)}{s+Dk^\alpha}dk\\ = \sum_{n=1}^\infty\frac{(-1)^{n+1}x_0^{2n}(D/s)^{-2n/\alpha}}{(2n)!\alpha s} \int_0^\infty dk_1k_1^{1/\alpha-1}H_{1,1}^{1,1}\left[\frac{D}{s}k_1\Bigg|\begin{matrix}(2n/\alpha,1)\\(2n/\alpha,1)\end{matrix}\right] ~, $$ and using the Mellin transform of the Fox $H$-function (equation (2.8), The H-Function: Theory and Applications) gives $$ \int_0^\infty \frac{1-\cos(kx_0)}{s+Dk^\alpha}dk\\ = \sum_{n=1}^\infty\frac{(-1)^{n+1}x_0^{2n}(D/s)^{-2n/\alpha}}{(2n)!\alpha s} (D/s)^{-1/\alpha} \Gamma(2n/\alpha+1/\alpha)\Gamma(1-2n/\alpha-1/\alpha) ~. $$ Taking the first term of the series gives us $$ \int_0^\infty \frac{1-\cos(kx_0)}{s+Dk^\alpha}dk \sim \frac{x_0^{2}(D/s)^{-3/\alpha}}{(2)!\alpha s}\Gamma(3/\alpha)\Gamma(1-3/\alpha) ~. \tag{2} $$
However, this expression cannot be reconciled with $\eqref{papereq}$, especially considering that the powers of $D$ and $x_0$ have differing dependences on $\alpha$ between the two equations. This leads one to believe that this process of obtaining asymptotic approximations doesn't work with these expressions.
What might be the correct way to obtain asymptotic approximations for this expression, as $s\to0$?