The value of logsine integrals $$\int_0^\frac{\pi}4 \log^k(\sin(\theta))d\theta, \ \ \int_0^\frac{\pi}4 \log^k(\cos(\theta))d\theta$$ Where $k=1,2$ are well known. Moreover, this answer by nospoon provides (essentially through the nice structure of the Fourier series of $\log\left(2\left|\sin\frac{x}{2}\right|\right)$) two closed forms for quite non-elementary integrals:
\begin{eqnarray*} \int_{0}^{\pi/4} \log^3(\sin\theta)\,d\theta &=& 3\,\text{Im}\,\text{Li}_4(1-i)-\frac{25\pi^3}{256}\log(2)+\frac{3}{2}\log(2)\,\text{Im}\,\text{Li}_3\left(\tfrac{1+i}{2}\right)\\ && -\frac{3K}{8}\log^2(2)-\frac{17\pi}{64}\log^3(2)-\frac{3\pi}{8}\zeta(3)+\frac{3}{4}\beta(4)\tag{1} \end{eqnarray*} and by differentiating the Euler Beta function $\int_0^\frac{\pi}2 \sin^a(\theta)d\theta$ three times, the logsine integral $\int_0^\frac{\pi}2 \log^3(\sin(\theta))d\theta$ is trivially evaluated, from which and reflection $\theta\to \frac{\pi}2-\theta$ the result below follows: \begin{eqnarray}\label{intlog3}\notag \int_{0}^{\pi/4} \log^3(\cos\theta)\,d\theta &=& -\frac{7\pi^3}{256}\log(2)-\frac{15\pi}{64}\log^3(2)+\frac{3K}{8}\log^2(2)-\frac{3\pi}{8}\zeta(3)\\ && -\frac{3}{4}\beta(4)-\frac{3}{2}\log(2)\,\text{Im}\,\text{Li}_3\left(\tfrac{1+i}{2}\right)-3\,\text{Im}\,\text{Li}_4(1-i).\tag{2} \end{eqnarray}
Q: I am interested in a closed form evaluation (in terms of Euler sums) for the integrals $$\color{blue}{ \int_{0}^{\pi/4}\log^4(\sin\theta)\,d\theta,\qquad \int_{0}^{\pi/4}\log^4(\cos\theta)\,d\theta }$$ whose sum is clearly given by $\frac{19\pi^5}{480}+\frac{\pi^3}{4}\log^2(2)+\frac{\pi}{2}\log^4(2)+3\pi\zeta(3)\log(2)$ thanks to Euler's Beta function again. Are their value already known in the literature? If so, does the evaluation procedure exploit the convolution identity $$ \log^2\left(2\sin\frac{x}{2}\right)\stackrel{L^2(0,\pi)}{=}\frac{\pi^2}{12}+\sum_{n\geq 1}\cos(nx)\frac{H_n+H_{n-1}}{n}\quad?\tag{3} $$
Based on my limited knowledge on logsine integrals, value of these two are not given in literature until this July when "paper $1$" (see link below) was published, and based on whose associated algorithm an equivalent $_pF_q$ representation of these integrals is given in "paper $2$". The solution does not depend on the mentioned convolutional identity but MZV theory instead. Indeed, let $x\to \tan^{-1}(u)$, these integrals are converted to $4$-admissible logarithmic integrals (for its definition, see "paper $1$"). Using the algorithm given by Paper $1$ one immediately get the result: $$\small \int_0^{\frac{\pi }{4}} \log ^4(\sin (x)) \, dx=\frac{1}{4} C \log ^3(2)-3 \beta (4) \log (2)-12 \Im\left(\text{Li}_5\left(\frac{1}{2}+\frac{i}{2}\right)\right)-\frac{3}{2} \log ^2(2) \Im\left(\text{Li}_3\left(\frac{1}{2}+\frac{i}{2}\right)\right)+6 \log (2) \Im\left(\text{Li}_4\left(\frac{1}{2}+\frac{i}{2}\right)\right)-3 S+\frac{3}{2} \pi \zeta (3) \log (2)+\frac{2093 \pi ^5}{30720}+\frac{43}{128} \pi \log ^4(2)+\frac{23}{64} \pi ^3 \log ^2(2)$$ $$\small\int_0^{\frac{\pi }{4}} \log ^4(\cos (x)) \, dx=-\frac{1}{4} C \log ^3(2)+3 \beta (4) \log (2)+12 \Im\left(\text{Li}_5\left(\frac{1}{2}+\frac{i}{2}\right)\right)+\frac{3}{2} \log ^2(2) \Im\left(\text{Li}_3\left(\frac{1}{2}+\frac{i}{2}\right)\right)-6 \log (2) \Im\left(\text{Li}_4\left(\frac{1}{2}+\frac{i}{2}\right)\right)+3 S+\frac{3}{2} \pi \zeta (3) \log (2)-\frac{877 \pi ^5}{30720}+\frac{21}{128} \pi \log ^4(2)-\frac{7}{64} \pi ^3 \log ^2(2)$$ Here $S=\Im L(4,1|i,1)=\sum _{k=1}^{\infty } \frac{(-1)^{k-1} H_{2 k-2}}{(2 k-1)^4}$ is an irreducible colored MZV, which is also expressible via hypergeometric functions by binomial expansion, due to the following identity given by paper $2$: $$\small \sqrt{2} \, _6F_5\left(\frac{1}{2},\frac{1}{2},\frac{1}{2},\frac{1}{2},\frac{1}{2},\frac{1}{2};\frac{3}{2},\frac{3}{2},\frac{3}{2},\frac{3}{2},\frac{3}{2};\frac{1}{2}\right)=-\frac{1}{4} S-\frac{1}{8} \beta (4) \log (2)-\Im\left(\text{Li}_5\left(\frac{1}{2}+\frac{i}{2}\right)\right)+\frac{1}{16} \pi \zeta (3) \log (2)+\frac{2093 \pi ^5}{368640}+\frac{1}{512} \pi \log ^4(2)+\frac{23 \pi ^3 \log ^2(2)}{3072}$$ For verification of $2$ integrals above, download the Mathematica package associated to "paper $1$" here and make use of the function MZIntegrate. Note that one may derive an Euler-sum identity using the convolutional identity mentioned and orthogonal relations with help of values of quadratic logsines above. Nevertheless, they themselves are easily transformed into MZVs directly.