How to evaluate improper integral $\int_{0}^{\infty}\frac{\tan^{-1}{x}}{e^{ax}-1}dx$? I'm trying to evaluate the improper integral,
$$I(a)=\int_{0}^{\infty}\frac{\tan^{-1}{x}}{e^{ax}-1}\mathrm{d}x,~~~\text{where }a\in\mathbb{R}^+.$$
Does this integral have a simple closed form expression? And if so, how best to obtain it?

My attempt
My first idea was to integrate by parts using $f=\tan^{-1}{x}$ and $dg=\frac{\mathrm{d}x}{e^{ax}-1}$, with the hope that the resulting integral would be amenable to solution by differentiating under the integral sign. So I found the indefinite integral,
$$\int\frac{\tan^{-1}{x}}{e^{ax}-1}\mathrm{d}x=\tan^{-1}{x}\left(\frac{1}{a}\log{(1-e^{ax})}-x\right)-\int\left(\frac{\log{(1-e^{ax})}}{a(x^2+1)}-\frac{x}{x^2+1}\right)\mathrm{d}x,$$
but then I realized that this would result in an imaginary boundary term for the corresponding definite integral over $[0,\infty)$ since,
$$\begin{cases}\lim_{x\to\infty}\tan^{-1}{x}\left(\frac{1}{a}\log{(1-e^{ax})}-x\right)=i\frac{\pi^2}{2a},\\
\lim_{x\to0}\tan^{-1}{x}\left(\frac{1}{a}\log{(1-e^{ax})}-x\right)=0.\end{cases}$$
I want to avoid complex variables if at all possible, so I don't know if I want to continue down this route.
Can anyone offer any hints or suggestions?

CORRECTION: As Vladimir pointed out, the correct anti-derivative is actually,
$$\int\frac{\tan^{-1}{x}}{e^{ax}-1}\mathrm{d}x=\tan^{-1}{x}\left(\frac{1}{a}\log{(e^{ax}-1)}-x\right)-\int\left(\frac{\log{(e^{ax}-1)}}{a(x^2+1)}-\frac{x}{x^2+1}\right)\mathrm{d}x.$$
Then, since,
$$\begin{cases}\lim_{x\to\infty}\tan^{-1}{x}\left(\frac{1}{a}\log{(e^{ax}-1)}-x\right)=\frac{\pi}{2}\cdot0=0,\\
\lim_{x\to0}\tan^{-1}{x}\left(\frac{1}{a}\log{-e^{ax}-1)}-x\right)=0,\end{cases}$$
we have,
$$I(a)=-\int_{0}^{\infty}\left(\frac{\log{(e^{ax}-1)}}{a(x^2+1)}-\frac{x}{x^2+1}\right)\mathrm{d}x.$$

I hope it is OK that I chime in, but I had played around with this one sometime back using the Abel-Plana formula. A formula which can come in handy.  It isn't as nice as Achille and Felix solutions, but it does show how cool the AP formula can be. 
In case you're unfamiliar, the Abel-Plana formula states:
$$\sum_{n=0}^{\infty}f(n)=1/2f(0)+\int_{0}^{\infty}f(x)dx+i\int_{0}^{\infty}\frac{f(ix)-f(-ix)}{e^{2\pi x}-1}dx$$
Now, let $\displaystyle f(n)=\frac{1}{(n+t)^{a}}$
Note the sum $\sum_{n=0}^{\infty}\frac{1}{(n+t)^{a}}=\zeta(a,t)$
and we have:  $$\sum_{n=0}^{\infty}\frac{1}{(n+t)^{a}}=\frac{1}{2t^{a}}+\frac{1}{(s-1)t^{s-1}}+i\int_{0}^{\infty}\frac{(t+ix)^{-a}-(t-ix)^{-a}}{e^{2\pi x}-1}dx \tag{1}$$.
Note that $$\frac{1}{(t+ix)^{a}}-\frac{1}{(t-ix)^{a}}=\frac{2}{i(t^{2}+x^{2})^{a/2}}\sin(a\cdot \tan^{-1}(x/t))$$
diff (1) w.r.t 'a' and get:
$$\zeta'(a,t)=\frac{-\log(t)}{2t^{a}}-\frac{t^{1-a}[1+(a-1)\log(t)]}{(a-1)^{2}}+i\int_{0}^{\infty}\frac{(t-ix)^{-a}\log(t-ix)-(t+ix)^{-a}\log(t+ix)}{e^{2\pi x}-1}dx \tag{2}$$
The integral on the right can be written as:
$$2\Im \int_{0}^{\infty}\frac{(t+ix)^{-a}\log(t+ix)}{e^{2\pi x}-1}dx$$
Now, let $a=0$ in (2):
$$\zeta'(0,t)=(t-1/2)\log(t)-t+i\int_{0}^{\infty}\frac{\log(t-ix)-\log(t+ix)}{e^{2\pi x}-1}dx$$
Notice what the right side is beginning to look like?. Remember the complex arctan identity.
so, we have $$\zeta'(0,t)=(t-1/2)\log(t)-t+2\int_{0}^{\infty}\frac{\tan^{-1}(x/t)}{e^{2\pi x}-1}dx$$
We can also use an identity that states $\displaystyle \zeta'(0,t)=\log\Gamma(t)-1/2\log(2\pi)$
and finally arrive at:
$$\log\Gamma(t)-1/2\log(2\pi)=(t-1/2)\log(t)-t+2\int_{0}^{\infty}\frac{\tan^{-1}(x/t)}{e^{2\pi x}-1}dx$$
But, to get this last arctan integral into your form, I used the same sub Felix used and let $t\to \frac{a}{2\pi}$
$$\frac{2\pi}{a}\int_{0}^{\infty}\frac{\tan^{-1}(\frac{2\pi u}{a})}{e^{2\pi u}-1}du=\frac{\pi}{a}\left[\log\Gamma(\frac{a}{2\pi})-\frac{a}{2\pi}\log(\frac{a}{2\pi})+1/2\log(\frac{a}{2\pi})+\frac{a}{2\pi}-1/2\log(2\pi)\right]$$
i.e. Let $a=2\pi$ and get the famous:
$\displaystyle \int_{0}^{\infty}\frac{\tan^{-1}(x)}{e^{2\pi x}-1}dx=1/2-1/4\log(2\pi)$
