In a (great) paper "A theorem for the closed-form evaluation of the first generalized Stieltjes constant at rational arguments and some related summations" by Iaroslav V. Blagouchine, the following integral representation of the first Stieltjes constant $\gamma_1$ is given (on page 539): $$\gamma_1=-\left[\gamma-\frac{\ln2}2\right]\ln2+i\int_0^\infty\frac{dx}{e^{\pi x}+1}\left\{\frac{\ln(1-ix)}{1-ix}-\frac{\ln(1+ix)}{1+ix}\right\}.\tag1$$ It's possible to get rid of imaginary numbers in this formula, and rewrite it in terms of only real-valued functions: $$2\int_0^\infty\frac{\arctan x}{1+x^2}\frac{dx}{e^{\pi x}+1}-\int_0^\infty\frac{x\ln(1+x^2)}{1+x^2}\frac{dx}{e^{\pi x}+1}=\gamma_1+\left[\gamma-\frac{\ln2}2\right]\ln2.\tag2$$
Question:
Is it possible to find closed forms separately for each integral on the left-hand side of $(2)$?
Just a partial answer for now, but a promising one. Let: $$\begin{eqnarray*} I(z)&\triangleq&\int_{0}^{+\infty}\frac{\arctan\frac{t}{z}}{(t^2+1)(e^{\pi t}+1)}\,dt\\ &=&\int_{0}^{+\infty}\frac{\arctan\frac{t}{z}}{(t^2+1)(e^{\pi t}-1)}\,dt+2\int_{0}^{+\infty}\frac{\arctan\frac{t}{z}}{(t^2+1)(e^{2\pi t}-1)}\,dt\\&=&2\int_{0}^{+\infty}\left(\frac{\arctan\frac{2t}{z}}{4t^2+1}+\frac{\arctan\frac{t}{z}}{t^2+1}\right)\frac{dt}{e^{2\pi t}-1}.\end{eqnarray*}$$ We have: $$ I'(z) = -2\int_{0}^{+\infty}\left(\frac{2}{(4t^2+1)(4t^2+z^2)}+\frac{1}{(t^2+1)(t^2+z^2)}\right)\frac{t\,dt}{e^{2\pi t}-1}$$ that is treatable through the Abel-Plana's formula with $$ f(u) = \frac{1}{(u+1)(u+z)}.$$ We have: $$ \sum_{n\geq 0}f(n) = \frac{\psi(z)-\psi(1)}{z-1},\qquad \int_{0}^{+\infty}f(u)\,du = \frac{\log z}{z-1},$$ and since for any $z>0$ we have: $$I(z)\leq \frac{1}{z}\int_{0}^{+\infty}\frac{t\,dt}{e^{\pi t}+1}=\frac{1}{12 z}$$ $I(1)$ just depends on: $$\int_{1}^{+\infty}\frac{\log z}{z^2-1}\,dz = \frac{\pi^2}{8}$$ and: $$ J = \int_{1}^{+\infty}\frac{\psi(z)-\psi(1)}{z^2-1}\,dz = \sum_{n\geq 0}\frac{1}{n+1}\int_{1}^{+\infty}\frac{dz}{(z+1)(z+n)}=\sum_{n\geq 0}\frac{\log\frac{n+1}{2}}{n^2-1}.$$ Since: $$ J = \log 2 +\frac{1}{4}-\sum_{n\geq 2}\frac{\log 2}{n^2-1}+\sum_{n\geq 2}\frac{\log(n+1)}{n^2-1}=\frac{1+\log 2}{4}+\sum_{n\geq 2}\frac{\log(n+1)}{n^2-1}$$ it is sufficient to compute: $$ \sum_{n\geq 2}\frac{\log(n+1)}{n^2-1}=\frac{d}{d\alpha}\left.\sum_{n\geq 2}\frac{(n+1)^{\alpha-1}}{(n-1)}\right|_{\alpha=0}.$$