Here is my trial, which is partially successful but still not fully answering to your question.
Using some coordinate change, I derived that
$$ I_{n} := \int_{0}^{\infty} \mathrm{erfc}^{n}(x) \, dx = \frac{1}{\sqrt{n}} \left( \frac{2}{\sqrt{\pi}} \right)^{n-1} \int_{T^{n-1}} \int_{0}^{\infty} s^{n-1}e^{-(|x|^{2}-1)s^{2}} \mathrm{erfc}(s) \, ds d\sigma_{x}, \tag{*} $$
where $T^{n-1} = \{ x \in [0, \infty)^{n} : x_{1} + \cdots + x_{n} = \sqrt{n} \}$ is an $(n-1)$-simplex and $d\sigma_{x}$ is the surface measure on $T^{n-1}$.
Case $n = 1$. Since $T^{0} = \{1\}$ and $d\sigma_{x} = \delta(x-1)$ is a unit mass located at $x = 1$, the equation $\text{(*)}$ gives no new information.
Case $n = 2.$ It is easy to find that
$$ \int_{0}^{\infty} s e^{-(|x|^{2}-1)s^{2}} \mathrm{erfc}(s) \, ds = \frac{1}{2|x|(|x|+1)}$$
and by suitable parametrization of the line segment $T^{1}$, we get
\begin{align*}
I_{2}
&= \sqrt{\frac{2}{\pi}} \int_{T^{1}} \frac{1}{2|x|(|x|+1)} \, d\sigma_{x} \\
&= \sqrt{\frac{2}{\pi}} \int_{0}^{\sqrt{2}} \frac{\sqrt{2} \, dt}{2\sqrt{t^{2} + (\sqrt{2}-t)^{2}} ( \sqrt{t^{2} + (\sqrt{2}-t)^{2}} + 1)} \\
&= \sqrt{\frac{2}{\pi}} \int_{0}^{\frac{\pi}{4}} \frac{d\theta}{1+\cos\theta}
= \sqrt{\frac{2}{\pi}} \tan\left(\frac{\pi}{8}\right)
= \frac{2 - \sqrt{2}}{\sqrt{\pi}}.
\end{align*}
We can also write
$$ I_{2} = \boxed{ \displaystyle \frac{2}{\sqrt{\pi}} - \frac{2\sqrt{2}}{\pi^{3/2}} \arctan(\infty) }. $$
Case $n = 3$. We have
$$ \int_{0}^{\infty} s^{2}e^{-(|x|^{2}-1)s^{2}} \mathrm{erfc}(s) \, ds = \frac{1}{2\sqrt{\pi}(|x|^{2}-1)} \left( \frac{\arctan\sqrt{|x|^{2}-1}}{\sqrt{|x|^{2}-1}} - \frac{1}{|x|^{2}} \right). $$
Now we know that $T^{2}$ is a regular triangle with side length $\sqrt{6}$. Thus by introducing polar coordinate change, we get
\begin{align*}
I_{3}
&= \frac{2}{\pi^{3/2}\sqrt{3}} \int_{T^{2}} \frac{1}{|x|^{2}-1} \left( \frac{\arctan\sqrt{|x|^{2}-1}}{\sqrt{|x|^{2}-1}} - \frac{1}{|x|^{2}} \right) \, d\sigma_{x} \\
&= \frac{4\sqrt{3}}{\pi^{3/2}} \int_{0}^{\frac{\pi}{3}} \int_{0}^{\frac{\sec\theta}{\sqrt{2}}} \frac{1}{r} \left( \frac{\arctan r}{r} - \frac{1}{r^{2} + 1} \right) \, drd\theta \\
&= \frac{4\sqrt{3}}{\pi^{3/2}} \int_{0}^{\frac{\pi}{3}} \left( 1- \frac{\arctan\left(\sec\theta / \sqrt{2}\right)}{\sec\theta / \sqrt{2}} \right) \, d\theta \\
&= \frac{4\sqrt{3}}{\pi^{3/2}} \int_{\frac{1}{\sqrt{2}}}^{\sqrt{2}} \left( 1 - \frac{\arctan u}{u} \right) \frac{du}{u\sqrt{2u^{2} - 1}}.
\end{align*}
I haven't tried the last integral, but Mathematica suggests that
$$ \int_{\frac{1}{\sqrt{2}}}^{\sqrt{2}} \left( 1 - \frac{\arctan u}{u} \right) \frac{du}{u\sqrt{2u^{2} - 1}} = \frac{\sqrt{3}\pi}{4} -\sqrt{\frac{3}{2}} \arctan\left( \sqrt{2} \right).$$
This gives
$$ I_{3} = \boxed{ \displaystyle \frac{3}{\sqrt{\pi}} - \frac{6\sqrt{2}}{\pi^{3/2}} \arctan\left( \sqrt{2} \right) }, $$
which can be numerically checked.
I am trying to generalize this calculation for higher dimensions, but it seems somewhat daunting. For example, if $n = 4$ we have to evaluate
$$ \int_{0}^{\infty} \mathrm{erfc}^{4}(x) \, dx = \frac{1}{\pi^{3/2}} \int_{T^{3}} \frac{2|x|+1}{(|x|+1)^{2}|x|^{3}} \, d\sigma_{x}. $$
Nevertheless, using the formula
$$ \int_{T^{n-1}} f(|x|) \, d\sigma_{x} = n! \int_{0}^{\sqrt{\frac{1}{n-1}}} \int_{0}^{\sqrt{\frac{n}{n-2}}s_{1}} \cdots \int_{0}^{\sqrt{\frac{3}{1}}s_{n-2}} f\left( \sqrt{1+s_{1}^{2} + \cdots + s_{n-1}^{2}} \right) \, ds_{n-1} \cdots ds_{1} $$
together with aid of Mathematica, I was able to evaluate $I_{4}$ and it was
$$ I_{4} = \boxed{ \displaystyle \frac{4}{\sqrt{\pi}} - \frac{24\sqrt{2}}{\pi^{3/2}} \arctan \left( \frac{1}{\sqrt{8}} \right) }. $$
These series of observations lead us to believe that $I_{n}$ is of the form
$$ I_{n} = \frac{n}{\sqrt{\pi}} - \frac{n!\sqrt{2}}{\pi^{3/2}} \arctan \alpha_{n} $$
for some reasonably nice $\alpha_{n}$ (with $\alpha_{1} = 0$, $\alpha_{2} = \infty$, $\alpha_{3} = \sqrt{2}$ and $\alpha_{4} = \frac{1}{\sqrt{8}}$), but inverse symbolic calculations show that this seems no longer the case for $n \geq 5$.
Indeed, for $n = 5$ we have
$$ I_{5} = \frac{5}{\sqrt{\pi}} - \frac{80\sqrt{2}}{\pi^{3/2}} \arctan \sqrt{\frac{2}{3}} + \frac{240\sqrt{2}}{\pi^{5/2}} A\left( \sqrt{\frac{5}{3}}, \sqrt{\frac{3}{2}}, \sqrt{\frac{4}{5}} \right), $$
where $A(p, q, r)$ is the generalized Ahmed's integral. I have no idea whether it will reduce to a familiar closed form or not.
$$I:=\int_{0}^{\infty}\frac{\ln{(x)}}{\sqrt{x}\,\sqrt{x+1}\,\sqrt{2x+1}}\mathrm{d}x.$$
After first multiplying and dividing the integrand by 2, substitute $x=\frac{t}{2}$:
$$I=\int_{0}^{\infty}\frac{2\ln{(x)}}{\sqrt{2x}\,\sqrt{2x+2}\,\sqrt{2x+1}}\mathrm{d}x=\int_{0}^{\infty}\frac{\ln{\left(\frac{t}{2}\right)}}{\sqrt{t}\,\sqrt{t+2}\,\sqrt{t+1}}\mathrm{d}t.$$
Next, substituting $t=\frac{1}{u}$ yields:
$$\begin{align}
I
&=-\int_{0}^{\infty}\frac{\ln{(2u)}}{\sqrt{u}\sqrt{u+1}\sqrt{2u+1}}\mathrm{d}u\\
&=-\int_{0}^{\infty}\frac{\ln{(2)}}{\sqrt{u}\sqrt{u+1}\sqrt{2u+1}}\mathrm{d}u-\int_{0}^{\infty}\frac{\ln{(u)}}{\sqrt{u}\sqrt{u+1}\sqrt{2u+1}}\mathrm{d}u\\
&=-\int_{0}^{\infty}\frac{\ln{(2)}}{\sqrt{u}\sqrt{u+1}\sqrt{2u+1}}\mathrm{d}u-I\\
\implies I&=-\frac{\ln{(2)}}{2}\int_{0}^{\infty}\frac{\mathrm{d}x}{\sqrt{x}\sqrt{x+1}\sqrt{2x+1}}.
\end{align}$$
Making the sequence of substitutions $x=\frac{u-1}{2}$, then $u=\frac{1}{t}$, and finally $t=\sqrt{w}$, puts this integral into the form of a beta function:
$$\begin{align}
\int_{0}^{\infty}\frac{\mathrm{d}x}{\sqrt{x}\sqrt{x+1}\sqrt{2x+1}}
&=\int_{1}^{\infty}\frac{\mathrm{d}u}{\sqrt{u-1}\sqrt{u+1}\sqrt{u}}\\
&=\int_{1}^{\infty}\frac{\mathrm{d}u}{\sqrt{u^2-1}\sqrt{u}}\\
&=\int_{1}^{0}\frac{t^{3/2}}{\sqrt{1-t^2}}\frac{(-1)}{t^2}\mathrm{d}t\\
&=\int_{0}^{1}\frac{\mathrm{d}t}{\sqrt{t}\,\sqrt{1-t^2}}\\
&=\frac12\int_{0}^{1}\frac{\mathrm{d}w}{w^{3/4}\,\sqrt{1-w}}\\
&=\frac12\operatorname{B}{\left(\frac14,\frac12\right)}\\
&=\frac12\frac{\Gamma{\left(\frac12\right)}\Gamma{\left(\frac14\right)}}{\Gamma{\left(\frac34\right)}}\\
&=\frac{\pi^{3/2}}{2^{1/2}\Gamma^2{\left(\frac34\right)}}
\end{align}$$
Hence,
$$I=-\frac{\ln{(2)}}{2}\frac{\pi^{3/2}}{2^{1/2}\Gamma^2{\left(\frac34\right)}}=-\frac{\pi^{3/2}\,\ln{(2)}}{2^{3/2}\,\Gamma^2{\left(\frac34\right)}}.~~~\blacksquare$$
Possible Alternative: You could also derive the answer from the complete elliptic integral of the first kind instead of from the beta function by making the substitution $t=z^2$ instead of $t=\sqrt{w}$.
$$\begin{align}
\int_{0}^{\infty}\frac{\mathrm{d}x}{\sqrt{x}\sqrt{x+1}\sqrt{2x+1}}
&=\int_{1}^{\infty}\frac{\mathrm{d}u}{\sqrt{u-1}\sqrt{u+1}\sqrt{u}}\\
&=\int_{1}^{\infty}\frac{\mathrm{d}u}{\sqrt{u^2-1}\sqrt{u}}\\
&=\int_{1}^{0}\frac{t^{3/2}}{\sqrt{1-t^2}}\frac{(-1)}{t^2}\mathrm{d}t\\
&=\int_{0}^{1}\frac{\mathrm{d}t}{\sqrt{t}\,\sqrt{1-t^2}}\\
&=2\int_{0}^{1}\frac{\mathrm{d}z}{\sqrt{1-z^4}}\\
&=2\,K{(-1)}\\
&=\frac{\Gamma^2{\left(\frac14\right)}}{2\sqrt{2\pi}}\\
&=\frac{\pi^{3/2}}{2^{1/2}\Gamma^2{\left(\frac34\right)}}.
\end{align}$$
Best Answer
Start with integration by parts (IBP) by setting $u=\ln^3(1+x)$ and $dv=\dfrac{\ln x}{x}\ dx$ yields \begin{align} I&=-\frac32\int_0^1\frac{\ln^2(1+x)\ln^2 x}{1+x}\ dx\\ &=-\frac32\int_1^2\frac{\ln^2x\ln^2 (x-1)}{x}\ dx\quad\Rightarrow\quad\color{red}{x\mapsto1+x}\\ &=-\frac32\int_{\large\frac12}^1\left[\frac{\ln^2x\ln^2 (1-x)}{x}-\frac{2\ln^3x\ln(1-x)}{x}+\frac{\ln^4x}{x}\right]\ dx\quad\Rightarrow\quad\color{red}{x\mapsto\frac1x}\\ &=-\frac32\int_{\large\frac12}^1\frac{\ln^2x\ln^2 (1-x)}{x}\ dx+3\int_{\large\frac12}^1\frac{\ln^3x\ln(1-x)}{x}\ dx-\left.\frac3{10}\ln^5x\right|_{\large\frac12}^1\\ &=-\frac32\color{red}{\int_{\large\frac12}^1\frac{\ln^2x\ln^2 (1-x)}{x}\ dx}+3\int_{\large\frac12}^1\frac{\ln^3x\ln(1-x)}{x}\ dx-\frac3{10}\ln^52. \end{align} Applying IBP again to evaluate the red integral by setting $u=\ln^2(1-x)$ and $dv=\dfrac{\ln^2 x}{x}\ dx$ yields \begin{align} \color{red}{\int_{\large\frac12}^1\frac{\ln^2x\ln^2 (1-x)}{x}\ dx}&=\frac13\ln^52+\frac23\color{blue}{\int_{\large\frac12}^1\frac{\ln^3x\ln (1-x)}{1-x}\ dx}. \end{align}
Let us integrating the indefinite form of the blue integral. \begin{align} \color{blue}{\int\frac{\ln^3x\ln (1-x)}{1-x}\ dx}=&-\int\sum_{n=1}^\infty H_nx^n\ln^3x\ dx\\ =&-\sum_{n=1}^\infty H_n\int x^n\ln^3x\ dx\\ =&-\sum_{n=1}^\infty H_n\frac{\partial^3}{\partial n^3}\left[\int x^n\ dx\right]\\ =&-\sum_{n=1}^\infty H_n\frac{\partial^3}{\partial n^3}\left[\frac{x^{n+1}}{n+1}\right]\\ =&-\sum_{n=1}^\infty H_n\left[\frac{x^{n+1}\ln^3x}{n+1}-\frac{3x^{n+1}\ln^2x}{(n+1)^2}+\frac{6x^{n+1}\ln x}{(n+1)^3}-\frac{6x^{n+1}}{(n+1)^4}\right]\\ =&-\ln^3x\sum_{n=1}^\infty \frac{H_{n+1}x^{n+1}}{n+1}+\ln^3x\sum_{n=1}^\infty \frac{x^{n+1}}{(n+1)^2}+3\ln^2x\sum_{n=1}^\infty \frac{H_{n+1}x^{n+1}}{(n+1)^2}\\&-3\ln^2x\sum_{n=1}^\infty \frac{x^{n+1}}{(n+1)^3}-6\ln x\sum_{n=1}^\infty \frac{H_{n+1}x^{n+1}}{(n+1)^3}+6\ln x\sum_{n=1}^\infty \frac{x^{n+1}}{(n+1)^4}\\&+6\sum_{n=1}^\infty \frac{H_{n+1}x^{n+1}}{(n+1)^4}-6\sum_{n=1}^\infty \frac{x^{n+1}}{(n+1)^5}\\ =&\ -\sum_{n=1}^\infty\left[\frac{H_nx^{n}\ln^3x}{n}-\frac{x^{n}\ln^3x}{n^2}-\frac{3H_nx^{n}\ln^2x}{n^2}+\frac{3x^{n}\ln^2x}{n^3}\right.\\& \left.\ +\frac{6H_nx^{n}\ln x}{n^3}-\frac{6x^{n}\ln x}{n^4}-\frac{6H_nx^{n}}{n^4}+\frac{6x^{n}}{n^5}\right]\\ =&\ -\color{blue}{\mathbf{H}_{1}(x)}\ln^3x+\operatorname{Li}_2(x)\ln^3x+3\color{blue}{\mathbf{H}_{2}(x)}\ln^2x-3\operatorname{Li}_3(x)\ln^2x\\&\ -6\color{blue}{\mathbf{H}_{3}(x)}\ln x+6\operatorname{Li}_4(x)\ln x+6\color{blue}{\mathbf{H}_{4}(x)}-6\operatorname{Li}_5(x). \end{align} Therefore \begin{align} \color{blue}{\int_{\Large\frac12}^1\frac{\ln^3x\ln (1-x)}{1-x}\ dx} =&\ 6\color{blue}{\mathbf{H}_{4}(1)}-6\operatorname{Li}_5(1)-\left[\color{blue}{\mathbf{H}_{1}\left(\frac12\right)}\ln^32-\operatorname{Li}_2\left(\frac12\right)\ln^32\right.\\&\left.\ +3\color{blue}{\mathbf{H}_{2}\left(\frac12\right)}\ln^22-3\operatorname{Li}_3\left(\frac12\right)\ln^22+6\color{blue}{\mathbf{H}_{3}\left(\frac12\right)}\ln 2\right.\\&\ -6\operatorname{Li}_4(x)\ln 2+6\color{blue}{\mathbf{H}_{4}(x)}-6\operatorname{Li}_5(x)\bigg]\\ =&\ 12\zeta(5)-\pi^2\zeta(3)+\frac{3}8\zeta(3)\ln^22-\frac{\pi^4}{120}\ln2-\frac{1} {4}\ln^52\\&\ -6\color{blue}{\mathbf{H}_{4}\left(\frac12\right)}+6\operatorname{Li}_4\left(\frac12\right)\ln 2+6\operatorname{Li}_5\left(\frac12\right). \end{align} Using the similar approach as calculating the blue integral, then \begin{align} \int\frac{\ln^3x\ln (1-x)}{x}\ dx&=-\int\sum_{n=1}^\infty \frac{x^{n-1}}{n}\ln^3x\ dx\\ &=-\sum_{n=1}^\infty \frac{1}{n}\int x^{n-1}\ln^3x\ dx\\ &=-\sum_{n=1}^\infty \frac{1}{n}\frac{\partial^3}{\partial n^3}\left[\int x^{n-1}\ dx\right]\\ &=-\sum_{n=1}^\infty \frac{1}{n}\frac{\partial^3}{\partial n^3}\left[\frac{x^{n}}{n}\right]\\ &=-\sum_{n=1}^\infty \frac{1}{n}\left[\frac{x^{n}\ln^3x}{n}-\frac{3x^{n}\ln^2x}{n^2}+\frac{6x^{n}\ln x}{n^3}-\frac{6x^{n}}{n^4}\right]\\ &=\sum_{n=1}^\infty \left[-\frac{x^{n}\ln^3x}{n^2}+\frac{3x^{n}\ln^2x}{n^3}-\frac{6x^{n}\ln x}{n^4}+\frac{6x^{n}}{n^5}\right]\\ &=6\operatorname{Li}_5(x)-6\operatorname{Li}_4(x)\ln x+3\operatorname{Li}_3(x)\ln^2x-\operatorname{Li}_2(x)\ln^3x. \end{align} Hence $$ \int_{\large\frac{1}{2}}^1\frac{\ln^3x\ln (1-x)}{x}\ dx=\frac{\pi^2}{6}\ln^32-\frac{21}{8}\zeta(3)\ln^22-6\operatorname{Li}_4\left(\frac{1}{2}\right)\ln2-6\operatorname{Li}_5\left(\frac{1}{2}\right)+6\zeta(5). $$ Combining altogether, we have
Continuing my answer in: A sum containing harmonic numbers $\displaystyle\sum_{n=1}^\infty\frac{H_n}{n^3\,2^n}$, we have \begin{align} \color{blue}{\mathbf{H}_{3}\left(x\right)}=&\frac12\zeta(3)\ln x-\frac18\ln^2x\ln^2(1-x)+\frac12\ln x\left[\color{blue}{\mathbf{H}_{2}\left(x\right)}-\operatorname{Li}_3(x)\right]\\&+\operatorname{Li}_4(x)-\frac{\pi^2}{12}\operatorname{Li}_2(x)-\frac12\operatorname{Li}_3(1-x)\ln x+\frac{\pi^4}{60}.\tag1 \end{align} Dividing $(1)$ by $x$ and then integrating yields $$\small\begin{align} \color{blue}{\mathbf{H}_{4}\left(x\right)}=&\frac14\zeta(3)\ln^2 x-\frac18\int\frac{\ln^2x\ln^2(1-x)}x\ dx+\frac12\int\frac{\ln x}x\bigg[\color{blue}{\mathbf{H}_{2}\left(x\right)}-\operatorname{Li}_3(x)\bigg]\ dx\\&+\operatorname{Li}_5(x)-\frac{\pi^2}{12}\operatorname{Li}_3(x)-\frac12\int\frac{\operatorname{Li}_3(1-x)\ln x}x\ dx+\frac{\pi^4}{60}\ln x\\ =&\frac14\zeta(3)\ln^2 x+\frac{\pi^4}{60}\ln x+\operatorname{Li}_5(x)-\frac{\pi^2}{12}\operatorname{Li}_3(x)-\frac18\color{red}{\int\frac{\ln^2x\ln^2(1-x)}x\ dx}\\&+\frac12\left[\color{purple}{\sum_{n=1}^\infty\frac{H_{n}}{n^2}\int x^{n-1}\ln x\ dx}-\color{green}{\int\frac{\operatorname{Li}_3(x)\ln x}x\ dx}-\color{orange}{\int\frac{\operatorname{Li}_3(1-x)\ln x}x\ dx}\right].\tag2 \end{align}$$ Evaluating the red integral using the same technique as the previous one yields \begin{align} \color{red}{\int\frac{\ln^2x\ln^2(1-x)}x\ dx}&=\frac13\ln^3x\ln^2(1-x)-\frac23\color{blue}{\int\frac{\ln(1-x)\ln^3 x}{1-x}\ dx}. \end{align} Evaluating the purple integral yields \begin{align} \color{purple}{\sum_{n=1}^\infty\frac{H_{n}}{n^2}\int x^{n-1}\ln x\ dx}&=\sum_{n=1}^\infty\frac{H_{n}}{n^2}\frac{\partial}{\partial n}\left[\int x^{n-1}\ dx\right]\\ &=\sum_{n=1}^\infty\frac{H_{n}}{n^2}\left[\frac{x^n\ln x}{n}-\frac{x^n}{n^2}\right]\\ &=\color{blue}{\mathbf{H}_{3}(x)}\ln x-\color{blue}{\mathbf{H}_{4}(x)}. \end{align} Evaluating the green integral using IBP by setting $u=\ln x$ and $dv=\dfrac{\operatorname{Li}_3(x)}{x}\ dx$ yields \begin{align} \color{green}{\int\frac{\operatorname{Li}_3(x)\ln x}x\ dx}&=\operatorname{Li}_4(x)\ln x-\int\frac{\operatorname{Li}_4(x)}x\ dx\\ &=\operatorname{Li}_4(x)\ln x-\operatorname{Li}_5(x). \end{align} Evaluating the orange integral using IBP by setting $u=\operatorname{Li}_3(1-x)$ and $dv=\dfrac{\ln x}{x}\ dx$ yields \begin{align} \color{orange}{\int\frac{\operatorname{Li}_3(1-x)\ln x}x\ dx}&=\frac12\operatorname{Li}_3(1-x)\ln^2 x+\frac12\color{maroon}{\int\frac{\operatorname{Li}_2(1-x)\ln^2 x}{1-x}\ dx}. \end{align} Applying IBP again to evaluate the maroon integral by setting $u=\operatorname{Li}_2(1-x)$ and $$ dv=\dfrac{\ln^2 x}{1-x}\ dx\quad\Rightarrow\quad v=2\operatorname{Li}_3(x)-2\operatorname{Li}_2(x)\ln x-\ln(1-x)\ln^2x, $$ we have $$\small{\begin{align} \color{maroon}{\int\frac{\operatorname{Li}_2(1-x)\ln^2 x}{1-x}\ dx}=&\left[2\operatorname{Li}_3(x)-2\operatorname{Li}_2(x)\ln x-\ln(1-x)\ln^2x\right]\operatorname{Li}_2(1-x)\\ &-2\int\frac{\operatorname{Li}_3(x)\ln x}{1-x}\ dx+2\int\frac{\operatorname{Li}_2(x)\ln x}{1-x}\ dx+\color{blue}{\int\frac{\ln(1-x)\ln^3 x}{1-x}\ dx}. \end{align}}$$
We use the generating function for the generalized harmonic numbers evaluate the above integrals involving polylogarithm.
Dividing generating function of $\color{blue}{\mathbf{H}^{(k)}(x)}$ by $x$ and then integrating yields
Repeating the process above yields
where it is easy to show by using IBP that
and
Now, all unknown terms have been obtained. Putting altogether to $(2)$, we have $$\small{\begin{align} \color{blue}{\mathbf{H}_{4}(x)} =&\ \frac1{10}\zeta(3)\ln^2 x+\frac{\pi^4}{150}\ln x-\frac{\pi^2}{30}\operatorname{Li}_3(x)-\frac1{60}\ln^3x\ln^2(1-x)+\frac65\operatorname{Li}_5(x)\\&-\frac15\left[\operatorname{Li}_3(x)-\operatorname{Li}_2(x)\ln x-\frac12\ln(1-x)\ln^2x\right]\operatorname{Li}_2(1-x)-\frac15\operatorname{Li}_4(x)\\&-\frac35\operatorname{Li}_4(x)\ln x+\frac15\operatorname{Li}_3(x)\ln x+\frac15\operatorname{Li}_3(x)\ln^2x-\frac1{10}\operatorname{Li}_3(1-x)\ln^2 x\\&-\frac1{15}\operatorname{Li}_2(x)\ln^3x-\frac15\color{blue}{\mathbf{H}_{2}^{(3)}(x)}+\frac15\color{blue}{\mathbf{H}_{2}^{(2)}(x)} +\frac15\color{blue}{\mathbf{H}_{1}^{(3)}(x)}\ln x\\&-\frac15\color{blue}{\mathbf{H}_{1}^{(2)}(x)}\ln x+\frac25\color{blue}{\mathbf{H}_{3}(x)}\ln x-\frac15\color{blue}{\mathbf{H}_{2}(x)}\ln^2x+\frac1{15}\color{blue}{\mathbf{H}_{1}(x)}\ln^3x+C.\tag3 \end{align}}$$ The next step is finding the constant of integration. Setting $x=1$ to $(3)$ yields $$\small{\begin{align} \color{blue}{\mathbf{H}_{4}(1)} &=-\frac{\pi^2}{30}\operatorname{Li}_3(1)+\frac65\operatorname{Li}_5(1)-\frac15\operatorname{Li}_4(1)-\frac15\color{blue}{\mathbf{H}_{2}^{(3)}(1)}+\frac15\color{blue}{\mathbf{H}_{2}^{(2)}(1)}+C\\ 3\zeta(5)+\zeta(2)\zeta(3)&=-\frac{\pi^2}{30}\operatorname{Li}_3(1)+\frac{19}{30}\operatorname{Li}_5(1)+\frac{3}{5}\operatorname{Li}_3(1)+C\\ C&=\frac{\pi^4}{450}+\frac{\pi^2}{5}\zeta(3)-\frac35\zeta(3)+3\zeta(5). \end{align}}$$ Thus $$\small{\begin{align} \color{blue}{\mathbf{H}_{4}(x)} =&\ \frac1{10}\zeta(3)\ln^2 x+\frac{\pi^4}{150}\ln x-\frac{\pi^2}{30}\operatorname{Li}_3(x)-\frac1{60}\ln^3x\ln^2(1-x)+\frac65\operatorname{Li}_5(x)\\&-\frac15\left[\operatorname{Li}_3(x)-\operatorname{Li}_2(x)\ln x-\frac12\ln(1-x)\ln^2x\right]\operatorname{Li}_2(1-x)-\frac15\operatorname{Li}_4(x)\\&-\frac35\operatorname{Li}_4(x)\ln x+\frac15\operatorname{Li}_3(x)\ln x+\frac15\operatorname{Li}_3(x)\ln^2x-\frac1{10}\operatorname{Li}_3(1-x)\ln^2 x\\&-\frac1{15}\operatorname{Li}_2(x)\ln^3x-\frac15\color{blue}{\mathbf{H}_{2}^{(3)}(x)}+\frac15\color{blue}{\mathbf{H}_{2}^{(2)}(x)} +\frac15\color{blue}{\mathbf{H}_{1}^{(3)}(x)}\ln x\\&-\frac15\color{blue}{\mathbf{H}_{1}^{(2)}(x)}\ln x+\frac25\color{blue}{\mathbf{H}_{3}(x)}\ln x-\frac15\color{blue}{\mathbf{H}_{2}(x)}\ln^2x+\frac1{15}\color{blue}{\mathbf{H}_{1}(x)}\ln^3x\\&+\frac{\pi^4}{450}+\frac{\pi^2}{5}\zeta(3)-\frac35\zeta(3)+3\zeta(5)\tag4 \end{align}}$$ and setting $x=\frac12$ to $(4)$ yields \begin{align} \color{blue}{\mathbf{H}_{4}\left(\frac12\right)}=&\ \frac{\ln^52}{40}-\frac{\pi^2}{36}\ln^32+\frac{\zeta(3)}{2}\ln^22-\frac{\pi^2}{12}\zeta(3)\\&+\frac{\zeta(5)}{32}-\frac{\pi^4}{720}\ln2+\operatorname{Li}_4\left(\frac12\right)\ln2+2\operatorname{Li}_5\left(\frac12\right).\tag5 \end{align}
Finally, we obtain
which again matches @Cleo's answer.
References :
$[1]\ $ Harmonic number
$[2]\ $ Polylogarithm