integration - Methods to attack integrals that include $(1+x)^{a}ln^{b}(1+x)$ in the integrand

I am looking for systematic methods to attack the following class of integrals involving logarithmic functions

$$\begin{aligned} I_{0} &= \int_{0}^{1}(1+x)^{a}\ln^{m}(1+x)\,\mathrm{d}x \\ I_{1} &= \int_{0}^{1} x^{a}(1+x)^{b}\ln^{m}x\ln^{n}(1+x)\,\mathrm{d}x \\ I_{2} &= \int_{0}^{1} (1-x)^{a}(1+x)^{b}\ln^{m}(1-x)\ln^{n}(1+x)\,\mathrm{d}x \\ I_{3} &= \int_{0}^{1}x^{a}(1-x)^{b}(1+x)^{c}\ln^{\ell}x\ln^{m}(1-x)\ln^{n}(1+x)\,\mathrm{d}x\end{aligned}$$

where $a, b, c$ are real numbers (usually integers and half-integers) and $\ell, m, n$ are whole numbers such that the integrals converge. Of course, the first three are special cases of $I_{3}$, but I thought it would be best to tackle the easier cases first before moving to the fully general case in case different techniques are needed.

Motivation. The beta function can be written in the form

$$\frac{\Gamma(a)\Gamma(b)}{\Gamma(a+b)} = \int_{0}^{1}t^{a-1}(1-t)^{b-1}\,\mathrm{d}t.$$

From here, it is straightforward to compute integrals of the type

$$\int_{0}^{1}x^{a}(1-x)^{b}\ln^{m}x\ln^{n}(1-x)\,\mathrm{d}x$$

through various series expansions or derivatives of the beta function. Thus one is not afraid of any powers on the logarithms, as one merely has to keep enough terms in the expansion, or compute enough derivatives.

Attempt at $I_{0}$. The presence of $1+x$ as a power and in the logarithm means that one cannot use the beta function directly. I am able to get away with it with the integral $I_{0}$ for now by shifting the boundaries and proceeding as

$$\begin{aligned} \int_{0}^{1}(1+x)^{a+\epsilon}\,\mathrm{d}x &= \sum_{m=0}^{\infty}\frac{\epsilon^{m}}{m!}\int_{0}^{1}(1+x)^{a}\ln^{m}(1+x)\,\mathrm{d}x \\ &= \int_{1}^{2}x^{a+\epsilon}\,\mathrm{d}x = \frac{2^{a+1+\epsilon} - 1}{a+1+\epsilon} = \frac{2^{a+1}e^{\epsilon\ln 2} - 1}{a+1+\epsilon} \end{aligned}$$

and keeping the appropriate order coefficient by expanding the exponential in its Taylor series and the denominator in a power series. For example,

$$\begin{aligned} \int_{0}^{1}(1+x)^{1+\epsilon}\,\mathrm{d}x &= \sum_{m=0}^{\infty}\frac{\epsilon^{m}}{m!}\int_{0}^{1}(1+x)\ln^{m}(1+x)\,\mathrm{d}x = \frac{4e^{\epsilon\ln 2} - 1}{2+\epsilon} \\ &\approx \left[4\left(1 + \epsilon\ln 2 + \frac{1}{2}\ln^{2}2\,\epsilon^{2}\right) - 1\right]\frac{1}{2}\left(1 - \frac{\epsilon}{2} + \frac{\epsilon^{2}}{4}\right) \\ &\approx \frac{1}{2}\left(3 + 4\ln 2\,\epsilon + 2\ln^{2}2\,\epsilon^{2}\right)\left(1 - \frac{\epsilon}{2} + \frac{\epsilon^{2}}{4}\right) \end{aligned}$$

so if one requires the $\epsilon^{2}$ coefficient, he merely keep all the terms, multiply by $2! = 2$ to account for the factorial in the original expansion, and arrive at

$$\int_{0}^{1}(1+x)\ln^{2}(1+x)\,\mathrm{d}x = \frac{3}{4} - 2\ln 2 + 2\ln^{2}2$$

with the integrals

$$\int_{0}^{1}(1+x)\ln(1+x)\,\mathrm{d}x = -\frac{3}{4} + 2\ln 2 \quad\text{and}\quad \int_{0}^{1}(1+x)\,\mathrm{d}x = \frac{3}{2}$$

found for "free" if we were not lazy with the multiplication. So we have a method for doing $I_{0}$ for any powers desired.

Attempts at $I_{1}$. The above technique does not seem to generalize to the other integrals. For example, I am currently trying to derive the result

$$\int_{0}^{1}x(1+x)\ln^{2}x\ln(1+x)\,\mathrm{d}x = -\frac{1}{4}\zeta(3) + \frac{4}{27}\ln 2 - \frac{5}{36}\zeta(2) + \frac{37}{72} \tag{1}$$

and many other integrals of a similar form.

(Edit: a method for $(1)$ has been found, see $(2)$ below)

$$\int_{0}^{1}(1+x)\ln x\ln^{2}(1+x)\,\mathrm{d}x = -\frac{1}{8}\zeta(3) + 5\ln 2 - 2\ln^{2}2 + \frac{1}{4}\zeta(2) - \frac{23}{8} \tag{2}$$

The class of integrals represented by $I_{1}$ seem difficult, as the zeta functions are familiar from the beta function, the $\ln 2$ is familiar from the treatment of $I_{0}$, but how an answer like this combines these constants is elusive. Hypergeometric functions seem out of reach, as I seem to lose information about the answer from the Pochhammer symbols. To demonstrate this, one may derive from the Euler integral representation of ${}_{2}F_{1}$

$$\begin{aligned} \int_{0}^{1}x^{\epsilon}(1+x)^{\delta}\,\mathrm{d}x &= \frac{1}{1+\epsilon}\,{}_{2}F_{1}(-\delta, 1+\epsilon; 2+\epsilon; -1) \\ \int_{0}^{1}(1-x)^{\epsilon}(1+x)^{\delta}\,\mathrm{d}x &= \frac{1}{1+\epsilon}\,{}_{2}F_{1}(-\delta, 1; 2+\epsilon; -1) \\ \int_{0}^{1}x^{\epsilon}(1-x)^{\delta}(1+x)^{\lambda}\,\mathrm{d}x &= \frac{\Gamma(1+\epsilon)\Gamma(1+\delta)}{\Gamma(2+\delta+\epsilon)}\,{}_{2}F_{1}(-\lambda, 1+\epsilon; 2+\delta+\epsilon; -1). \end{aligned}$$

To take the first one as an example, we can write this expression as

$$\frac{1}{1+\epsilon}\,{}_{2}F_{1}(-\delta, 1+\epsilon; 2+\epsilon; -1) = \sum_{k=0}^{\infty}\frac{(-1)^{k}}{k!}\frac{(-\delta)_{k}}{1+k+\epsilon} = \sum_{k=0}^{\infty}{\delta\choose k}\frac{1}{1+k+\epsilon}.$$

Thus one can use power series again to find the desired coefficient of $\epsilon$, but the $(-\delta)_{k}$ makes things difficult. If we wanted to find the coefficient of $\epsilon\delta^{2}$, for example as needed in $(2)$, we would need to sum up an infinite series, since every term $k \geq 2$ will contribute to this coefficient. This particular series does numerically converge to the correct answer, but I cannot find a way to relate it to well-known constants such as $\zeta(3)$ and $\ln 2$.

I have also searched up some other integrals of similar forms on this site, and have found that many of the answers proceed with pages of horrendous algebra involving double series, harmonic numbers, etc. which I would like to avoid if possible as I am not very familiar with manipulating harmonic numbers as they appear in series. Of course, I will accept manipulations of harmonic numbers, etc. if it is the only straightforward way to proceed.

Lastly, the references below may be of some use. The second reference may well be the key to solving this problem but I would like some clarification as dealing with the $\ln(1+x)$ still seems problematic...

1. Devoto, A. and Duke, D. Table of integrals and formulae for Feynman diagram calculations. La Rivista del Nuovo Cimento, 1984.




2. Vermaseren, J. A. M. Harmonic sums, Mellin transforms and integrals. International Journal of Modern Physics A, 2012.




3. Sofo, A. Integrals of logarithmic and hypergeometric functions. Communications in Mathematics, 2016.