Este es un reto integral con inesperados giros y vueltas en su evaluación. En última instancia, sin embargo, todo funciona.
Comenzamos por la integración por partes, algo que normalmente no tienen que hacer, pero los infinitos cancelar:
$$\int_0^{\infty} dx \frac{x}{x^2+a^2}\log{\left(\frac{x+1}{x-1}\right)} = \\ \lim_{N \rightarrow \infty}\left[\frac{1}{2} \log{(x^2+a^2)}\log{\left(\frac{x+1}{x-1}\right)} \right]_0^N + \frac{1}{2} \int_0^{\infty} dx \log{(x^2+a^2)}\left(\frac{1}{x-1}-\frac{1}{x+1} \right)$$
Ahora, el límite en el primer término en el lado derecho de la ecuación anterior converge a $-i \pi \, \log{a}$. El segundo término, la integral, tiene una singularidad de la que se pueden deformar el intervalo de integración; la contribución de esta deformación es
$$\lim_{\epsilon \rightarrow 0} i \frac{1}{2} \epsilon \int_0^1 d\phi e^{i \phi} \frac{\log{(1+a^2)}}{\epsilon e^{i \phi}} = i \frac{\pi}{2} \log{(1+a^2)}$$
Por lo tanto la integral es ahora igual a
$$PV \int_0^{\infty} dx \frac{\log{(x^2+a^2)}}{x^2-1} + i \frac{\pi}{2} \log{\left ( 1+ \frac{1}{a^2} \right )}$$
Para evaluar esto, vamos a $I(a)$ ser la parte real de la expresión anterior (la integral), y se diferencian $I(a)$ con respecto al $a$ (que supongo que es válido el paso, no voy a demostrar aquí) y obtener
$$\frac{\partial I}{\partial a} = 2 a \, PV \int_0^{\infty} \frac{dx}{x^2+a^2} \frac{1}{x^2-1} = a PV \int_{-\infty}^{\infty} \frac{dx}{x^2+a^2} \frac{1}{x^2-1}$$
Tenga en cuenta que esta integral es igual a
$$a \oint_C \frac{dz}{z^2+a^2} \frac{1}{z^2-1}$$
donde $C$ es un arco semicircular en la mitad superior del plano, justo encima del eje real (por lo que podemos ignorar los polos $z= \pm 1$, ya que la integral es un valor principal de Cauchy). Esta integral es igual a $i 2 \pi$ veces la suma de los residuos en los polos dentro de $C$. El único polo dentro de$C$$z=ia$, por lo que la integral es
$$a PV \int_{-\infty}^{\infty} \frac{dx}{x^2+a^2} \frac{1}{x^2-1} = a \, i 2 \pi \frac{1}{i 2 a} \frac{1}{(-a^2-1)} = -\frac{\pi}{a^2+1}$$
Integrando con respecto a $a$, nos encontramos con que
$$I(a) = K - \pi \arctan{a}$$
donde $K$ es una constante de integración. Determinamos $K$ igual $I(0)$. Así que el problema ahora se reduce a la evaluación de
$$I(0) =2 \int_0^{\infty} dx \frac{\log{x}}{x^2-1}$$
Tenga en cuenta que la singularidad en $x=1$ es desmontable en esta integral y por lo tanto no es necesario utilizar un valor principal de Cauchy. Podemos evaluar esta integral por una vez más atractivo para el teorema de los residuos, pero esta vez, consideramos que
$$\oint_{C'} dz \frac{\log^2{z}}{z^2-1}$$
donde $C'$ es un ojo de la cerradura de contorno con respecto al eje real positivo. Mediante la integración de todo este contorno, y observando que el integrando se desvanece suficientemente rápido como el radio de la sección circular de $C'$ aumenta sin límite, obtenemos
$$\oint_{C'} dz \frac{\log^2{z}}{z^2-1} = -i 4 \pi \int_0^{\infty} dx \frac{\log{x}}{x^2-1} + 4 \pi^2 \int_0^{\infty} dx \frac{1}{x^2-1}$$
Esto es igual a, por el teorema de los residuos, $i 2 \pi$ veces la suma de los residuos de los polos de el integrando del complejo integral dentro de $C'$. Como el único polo es al $z=-1$, vemos que
$$\begin{align}\oint_{C'} dz \frac{\log^2{z}}{z^2-1} &= i 2 \pi \frac{\log^2{(-1)}}{2 (-1)} \\ &= i 2 \pi \frac{\pi^2}{2}\end{align}$$
Ahora, la parte real de la integral anterior se divide en un valor principal de Cauchy y un pedazo de sangría acerca de la singularidad en $x=1$. El valor principal de Cauchy es cero:
$$\begin{align}PV \int_0^{\infty} dx \frac{1}{x^2-1} &= \lim_{\epsilon \rightarrow 0} \left [\int_0^{1-\epsilon} dx \frac{1}{x^2-1} + \int_{1+\epsilon}^{\infty} dx \frac{1}{x^2-1}\right]\\ &= \lim_{\epsilon \rightarrow 0} \left [\int_0^{1-\epsilon} dx \frac{1}{x^2-1} + \int_0^{1/(1+\epsilon)} \left (-\frac{dx}{x^2} \right ) \frac{1}{(1/x^2)-1} \right ]\\ &= \lim_{\epsilon \rightarrow 0} \left [\int_0^{1-\epsilon} dx \frac{1}{x^2-1} - \int_0^{1-\epsilon} \frac{dx}{x^2-1} \right ] \\ &= 0\end{align}$$
La sangría en el contorno, sin embargo, produce una contribución; deje $x=1+\epsilon e^{i \phi}$$\phi \in [\pi,0]$:
$$4 \pi^2 i \epsilon \int_{-\pi}^0 d\phi \frac{e^{i \phi}}{2 \epsilon e^{i \phi}} = i \frac{\pi}{2} 4 \pi^2$$
así que
$$-i 4 \pi \int_0^{\infty} dx \frac{\log{x}}{x^2-1} = i 2 \pi \frac{\pi^2}{2} - i \frac{\pi}{2} 4 \pi^2 = -i 2 \pi \frac{\pi^2}{2}$$
Por lo tanto
$$K = 2 \int_0^{\infty} dx \frac{\log{x}}{x^2-1} = \frac{\pi^2}{2}$$
y
$$I(a) = \frac{\pi^2}{2} - \pi \arctan{a} = \pi \left ( \frac{\pi}{2} - \arctan{a}\right ) = \pi \arctan{\frac{1}{a}}$$
y
$$\int_0^{\infty} dx \frac{x}{x^2+a^2}\log{\left(\frac{x+1}{x-1}\right)} = \pi \arctan{\frac{1}{a}}+ i \frac{\pi}{2} \log{\left ( 1+ \frac{1}{a^2} \right )}$$
Nota la parte imaginaria, la cual difiere de la declaración original del problema. Mediante el uso de los valores absolutos, esta parte imaginaria se va y se la declaró resultado es true.