Conjecture $\int_0^1\frac{\ln^2\left(1+x+x^2\right)}x dx\stackrel?=\frac{2\pi}{9\sqrt3}\psi^{(1)}(\tfrac13)-\frac{4\pi^3}{27\sqrt3}-\frac23\zeta(3)$ I'm interested in the following definite integral:
$$I=\int_0^1\frac{\ln^2\!\left(1+x+x^2\right)}x\,dx.\tag1$$
The corresponding antiderivative can be evaluated with Mathematica, but even after simplification is quite clumsy. It matches results of numerical integration, and its correctness can potentially be verified by hand using differentiation. So, we are assured that a closed form exists for $I$, albeit complicated one.
My program for a numerical search for closed forms found a much simpler candidate:

$$I\stackrel{\color{gray}?}=\frac{2\pi}{9\sqrt3}\psi^{\small(1)}\!\left(\tfrac13\right)-\frac{4\pi^3}{27\sqrt3}-\frac23\zeta(3).\tag2$$

Note that the trigamma value here can be expressed in terms of the dilogarithm of complex argument (see formula $(5)$ here) or of the $2^{nd}$ order harmonic number of fractional argument:
$$\begin{align}\psi^{\small(1)}\!\left(\tfrac13\right)&=\frac{2\pi^2}3+2\sqrt3\,\Im\,\operatorname{Li}_2\!\left[(-1)^{\small1/3}\right],\tag3\\\psi^{\small(1)}\!\left(\tfrac13\right)&=\frac{\pi^2}6+9-H^{\small(2)}_{\small1/3}.\tag4\end{align}$$
Can we prove $(2)$, preferably not going through the huge intermediate antiderivative?

One possible direction that I thought of is to factor the polynomial under the logarithm:
$$I=\int_0^1\Big[\ln\!\left(x+(-1)^{\small1/3}\right)+\ln\!\left(x-(-1)^{\small2/3}\right)\Big]^2x^{-1}\,dx.\tag5$$
After expanding the square brackets, Mathematica can find a simpler antiderivative for it. Can we reach $(2)$ following this direction manually?
 A: By some simple manipulations the integral gets reduced to calculating 
$$\frac{1}{2}\int_1^3 \frac{\log ^2(x)}{x-1} \ dx+\frac{1}{2}\int_1^3\frac{\log ^2(x)}{ (x-1) \sqrt{4 x-3}} \ dx,$$ 
and both residual integrals are pretty easy to finish at this point. 
A: *

*Replace $\ln(1+x+x^2)$ by $\ln(1-x^3)-\ln(1-x)$.

*Two of the resulting integrals are easy to compute:
$$\int_0^1\frac{\ln^2\left(1-x^3\right)dx}{x}=\frac23\zeta(3),\quad \int_0^1\frac{\ln^2\left(1-x\right)dx}{x}=2\zeta(3).$$

*Mathematica computes and fullsimplifies the remaining nontrivial integral $\int_0^1\frac{\ln\left(1-x\right)\ln\left(1-x^3\right)dx}{x}$ to a one-line expression containing a sum of two trilogarithms $\operatorname{Li}_3(z_1)+\operatorname{Li}_3(z_2)$.

*It so happens that $z_1+z_2=1$, hence thanks to Landen's identity the above sum is equal to $$-\operatorname{Li}_3\left(\frac{z_1}{z_1-1}\right)+\zeta\left(3\right)+\text{elementary}.$$ 

*Finally, it so happens that $\frac{z_1}{z_1-1}=e^{2\pi i/3}$ and the corresponding trilogarithmic value is known in terms of $\zeta(3)$. 
Altogether this should lead to your answer.
