For what $a,b$ the integral
$$\int^1_0 x^a(1-x)^b\operatorname{Li}_2 (x)\, \mathrm dx$$
has a closed form solution? I tried to solve it by expanding dilogarithm function, or by reducing it to linear combination of similar integrals, but I was unable to sum up series containing product of harmonic numbers. I wonder whether a closed formula is indeed available.
Reduction. As @Raffaele mentioned, by expanding dilog, integrating term by term (Beta integral) and summing up, one reduces the integral to a generalized hypergeometric function, i.e. $$\int^1_0 x^a (1-x)^b\operatorname{Li}_2 (x)\, \mathrm dx=B(a+2,b+1) \, _4F_3(1,1,1,a+2;2,2,a+b+3;1)$$
Claim. This hypergeometric function can be expressed by Multiple Zeta Values of level $1,2,4$ (up to Gamma factors) whenever:
And can be expressed via polygamma function whenever:
Proof. For proof of the MZV part (mainly by Fourier Legendre expansion), see Theorem $1, 2$ of this paper (there seems to be no general formula though). For case $1$ of polygamma part, recall reflection formula of dilog one have a general formula that
For case $2$ of polygamma part, see Theorem $3$ of paper linked above. For case $3$, due to second entry of this page the integral is evaluable whenever $b+2a=-2$. As @pisco mentioned, by repeated IBP and modulo polygamma integrals $\int_0^1 x^* (1-x)^* \log^*(1-x) \mathrm dx$ one may evaluate all $I(a+m,b+n)$ once $I=I(a,b)$ is known; this completes the proof of case $b+2a\in\mathbb Z$. The case $a+2b\in \mathbb Z$ is direct by reflection again (see case $1$).
Examples. Case $1$ of MZV part:
Case $2$ of MZV part:
Case $3$ of MZV part:
Polygamma part:
Notes.