Dimensional regularization and expansion of gamma function

Question

In my calculations, I used dimensional regularization, i.e. replace $d\rightarrow d-\epsilon$ and calculated the divergent integral. Then, I would like to expand the answer into seriers by $\epsilon$ around $\epsilon=0$. But I obtained strange result. I start from the following integral (where I denote $d=3-\epsilon$): $$\int_{0}^{\infty}dp\frac{p^{3-1}}{p^2+m^2},$$ which is divergent. Then, I have calculated the integral $$I(\epsilon)=\int_{0}^{\infty}dp\frac{p^{d-1}}{p^2+m^2}=\frac{m^{d-2}}{2}\Gamma(d/2)\Gamma(1-d/2),$$ which is convergent for $d<2$. Also, I also have integral over angles, which is in $d$-dimensional case can be written as $$\frac{2\pi^{d/2}}{\Gamma(d/2)}.$$ So, my final answer is $$I(\epsilon)\propto\Gamma(\epsilon/2-1/2).$$ Using Wolfram Mathematia, I find the expansion around $\epsilon=0$. My expectation was that the divergence of my integral will be appear like a pole, $1/\epsilon$. But from the expansion I see no one singular term.

Answer 1

The key problem is that dimensional regularization is insensitive to higher divergences, i.e. if the integral diverges linear (my case) the divergence is not a pole $1/\epsilon$.

The question is near to close. Now the main problem is convergence of considered integral. Indeed, one can calculate the integral in $d$-dimensional space and expressed the answer in term of Euler $\beta$-function. Then, if one set $d=3$ the answer is finite while the initial integral is linearly divergent.