Integral ${\large\int}_0^1\frac{\ln^2\ln\left(\frac1x\right)}{1+x+x^2}dx$

Solution 1:

You may write $$\begin{align*} {\large\int}_0^1\frac{\ln^2\ln\left(\frac1x\right)}{1+x+x^2}dx&={\large\int}_0^\infty \left(1-e^{-t}\right)\frac{\ln^2t}{1-e^{-3t}}e^{-t}dt\\ &=\sum_{n=0}^\infty{\large\int}_0^\infty \left(e^{-t}-e^{-2t}\right)e^{-3nt}\ln^2t\:dt\\ &=\sum_{n=0}^\infty \left({\large\int}_0^\infty e^{-(3n+1)t}\ln^2t\:dt-{\large\int}_0^\infty e^{-(3n+2)t}\ln^2t\:dt\right)\\ &=\sum_{n=0}^\infty \left. \partial_{s}^2 \left({\large\int}_0^\infty t^se^{-(3n+1)t}\:dt-{\large\int}_0^\infty t^s e^{-(3n+2)t}\:dt\right)\right|_{s=0}\\ &=\left. \partial_{s}^2 \left(\Gamma(s+1) \left(\sum_{n=0}^\infty \frac{1}{(3n+1)^{s+1}}-\sum_{n=0}^\infty \frac{1}{(3n+2)^{s+1}}\right)\right)\right|_{s=0}\\ &=\left. \partial_{s}^2 \left(\frac{\Gamma(s+1)}{3^{s+1}} \left(\zeta\left(s+1,\frac13\right)-\zeta\left(s+1,\frac23\right)\right)\right)\right|_{s=0}\\ \end{align*}$$ Then, using the Laurent series expansion of the Hurwitz zeta function near $1$, $$ \zeta(s+1,a) = \frac1s+\sum_{k=0}^{\infty}\frac{(-1)^{k}}{k!} \gamma_{k}(a) s^{k}, \quad a>0, \, s\to 0, $$ with $\displaystyle \gamma_{0}(a)= - \psi(a)=-\Gamma'(a)/\Gamma(a)$, we get

$$ \color{#0052A3}{\frac{\pi^3\sqrt{3}}{54}+\frac{\sqrt{3}\pi}{9} (\gamma+\ln 3)^2+\frac{2}{3} (\gamma+\ln 3)\! \left(\gamma_1\!\!\left(\frac{1}{3}\right)-\gamma_1\!\!\left(\frac{2}{3}\right)\right)+\frac{1}{3}\!\!\left(\gamma_2\!\!\left(\frac{1}{3}\right)-\gamma_2\!\!\left(\frac{2}{3}\right)\right)}. $$

Observe that, as pointed out by Vladimir Reshetnikov, you can make the following substitution

$$ \color{#3366FF}{\gamma_1\left(\tfrac13\right)-\gamma_1\left(\tfrac23\right)=\frac\pi{\sqrt 3}\left(6\ln\Gamma\left(\tfrac13\right)-\gamma+\frac{\ln 3}2-4\ln(2\pi)\right)}. $$