En primer lugar observamos que la integral puede escribirse como
$$\frac12 \int_{-\infty}^{\infty} dx \frac{\cos{a x}}{1+x^2} - \frac12 \frac{d}{da} \int_{-\infty}^{\infty} dx \frac{\cos{a x}}{1+x^2} $$
Así que considere la integral de contorno
$$\oint_{C(a)} dz \frac{e^{i a z}}{1+z^2} $$
donde $C(a)$ es un semicírculo de radio $R$ en la mitad superior del plano al $a\gt 0$ y en la mitad inferior del plano al $a \lt 0$. Por ejemplo, cuando se $a \gt 0$, la integral es
$$\int_{-R}^R dx \frac{e^{i a x}}{1+x^2} + i R \int_0^{\pi} d\phi \, e^{i \phi} \frac{e^{i a R \cos{\phi}} e^{-a R \sin{\phi}}}{1+R^2 e^{i 2 \phi}} $$
Como $R \to \infty$, la magnitud de la segunda integral es delimitada por
$$\frac{2 R}{R^2-1} \int_0^{\pi/2} d\phi \, e^{-2 a R \phi/\pi} = \frac{\pi}{a (R^2-1)} \left ( 1-e^{-a R} \right )$$
en el que claramente se desvanece en este límite. Por otro lado, el contorno de la integral es también igual a $i 2 \pi$ veces el residuo en la pole $z=i$, por lo que
$$\int_{-\infty}^{\infty} dx \frac{\cos{a x}}{1+x^2} = i 2 \pi \frac{e^{-a}}{2 i} = \pi \, e^{-a} \quad (a\gt 0)$$
(Aquí se usa la simetría de el integrando para obtener el original de la integral de la parte de atrás.) Para $a \lt 0$, se cierra en la mitad inferior del plano (o simplemente el uso de la uniformidad del coseno) y encontrar que
$$\int_{-\infty}^{\infty} dx \frac{\cos{a x}}{1+x^2} = \pi \, e^{-|a|} $$
A continuación, tomamos la derivada de esta integral (por separado para los casos de $a \gt 0$$a \lt 0$) encuentra que el original de la integral es
$$\int_0^{\infty} dx \frac{\cos{a x} + x \sin{a x}}{1+x^2} = \frac{\pi}{2} \left (1+\operatorname{sgn}{a} \right )e^{-|a|} = \pi \theta(a) e^{-|a|}$$
donde $\theta$ es la función escalón unitario, que es igual a $1/2$ al $a=0$.