La primera integral no es realmente acerca de la adición de un arbitrario regulador en la forma de exponencial de amortiguación $e^{-\epsilon x}$, como podemos ver en algo como
$$\sum_n n =\lim_{\epsilon\rightarrow 0} \sum_n n e^{-\epsilon n}=\lim_{\epsilon\rightarrow 0}\frac{1}{4 \cosh^2{\epsilon/2}}=\lim_{\epsilon\rightarrow 0}\left(\frac{1}{\epsilon^2}-\frac{1}{12}+O(\epsilon^2)\right)\stackrel{magic}{\sim} -\frac{1}{12}.$$
Por el contrario, el punto es no elegir algo arbitrario que tiene el límite correcto, pero la resolución de un problema más general dejando $p \to \tilde p = p + i \epsilon \in\mathbb C $. El resultado es entonces perfectamente válido para $\text{Im}(\tilde p)=\epsilon>0$, y el límite es a menudo escrita como
$$I_0 = \frac{i}{p+i 0},$$
cual es la forma que podemos ver a menudo en QFT. El "$+i0$" tiene gran importancia, como veremos un poco más adelante.
Métodos similares a lo OP hizo con $I(\alpha)$ funcionaría, si se producen constantes las respuestas. Pero, OP mostraron que aunque el $\epsilon\rightarrow 0$ límite es el mismo, la respuesta es no. Con $I_0$ sin embargo, no tenemos que introducir nada nuevo, acabamos de promover la $p$ a una variable compleja, y resolver la misma integral que originalmente tenía.
Espero que esto explique por qué este particular integral "obras". No estoy diciendo que no hay muchas otras dudosa integrales en la física, aunque!
Ahora un poco de una larga digresión acerca de por qué queremos tener $i0$ términos en el denominador y por qué tiene sentido físicamente, lo que puede saltar libremente si usted está satisfecho con mi respuesta a la pregunta hasta ahora.
Vamos a usar la función de Green para resolver el oscilador armónico forzado ecuación. Sólo para recordar, si tenemos un lineal de educación a distancia (o PDE) en el formulario $$\mathcal L_x[y]=f(x).$$ Then we can try writing the solution as a convolution with the forcing: $y(x)=\int \mathrm d x'G(x,x')f(x')$, and inserting this into the differential equation we get $$f(x)=\mathcal L_x[y]=\mathcal L[\int \mathrm d x' G(x,x')f(x')]=\int \mathrm d x'\mathcal L_x[ G(x,x')]f(x'),$$ which means that $L_x[ G(x,x')] = \delta(x-x').$ In most cases $G(x,x')=G(x-x')$, y para el oscilador armónico, tenemos
$$\partial_{t}^2G(t-t') + \omega_0^2 G(t-t')=\delta(t-t'),$$
que podemos resolver a través de la transformada de Fourier. Tenemos
$$G(t-t')=-\frac{1}{2\pi}\int_{\mathbb R} d\omega \frac{e^{i\omega(t-t')}}{\omega^2-\omega_0^2}.$$
Esto no convergen debido a los polos en el eje real. La solución es la adición de $+i0$ a los polacos y a cambio de ellos sobre el eje real. Ahora, cuando $t>t'$ debemos cerrar el contorno sobre el eje real, donde tenemos los 2 y los residuos de todo es grande. Pero cuando $t<t'$, para cerrar en la mitad inferior del plano, y el resultado es cero ya que el integrando es analítica de allí. Y, en general, la adición de $+i0$ restringe la función de Green para $G(t,t')\to G(t,t')\theta(t-t')$, el cual es una declaración acerca de la causalidad, y desde entonces se han
$$y(t)=\int_{-\infty}^{\infty}\mathrm d t G(t,t')\theta(t-t') f(t')=\int_{-\infty}^{t}\mathrm d t G(t,t') f(t').$$
El mismo truco se utiliza en la electrodinámica clásica para solucionar $(\nabla^2 -1/c^2 \partial_t^2)\phi = -\frac{1}{\epsilon_0}\rho$ y garantiza la causal de propagación a la velocidad de la luz. Casi todos los propagador (es decir, la función de Green) en QFT ha $\pm i0$ - a pesar de que estos términos pueden ser omitidos en la escritura, ellos están ahí. En el oscilador armónico caso, el resultado es básicamente el mismo como la adición de amortiguación para el problema original, que es físicamente mejor, ya que si tenemos $f(t)=f_0 e^{i\Omega t}$, $y(t)$ va a explotar si tenemos resonancia $\omega_0=\Omega$ - esto es exactamente lo que sucede en la integral de Fourier, donde estamos, básicamente, la integración a través de una descomposición espectral de $f(t)$, lo cual también contiene la 'resonancia'. El punto principal sobre el $i0$ es la relación de causalidad, por supuesto.