I believe a closed form for this integral exists since WFA is able to compute the antiderivative. It does not come out with a closed form when I give it bounds on the integral for some reason though. This integral is tricky. I tried Feynman's Trick by paramterising with $$I(a)=\int_0^1 \ln(ax)\operatorname{Li}_2(x)dx$$
however, $a=0$ does not give a value so this cannot be used. It would be very useful though since after differentiating we get $\frac{\operatorname{Li}_2(x)}{ax}$ as the integrand which is just $\operatorname{Li}_3(x) a^{-1}$ but of course this is not possible because of the $a=0$ problem. Are there any other ways to approach this? I think integration by parts may work but I'm not sure how to deal with the integrals that result from it.