Following is the proof that,
\begin{equation*}
\displaystyle \int_0^{\pi/4} \ln(\tan x)\ln(\cos x-\sin x)dx=\dfrac{G\ln 2}{2}
\end{equation*}
Knowing that,
\begin{equation*}
\displaystyle \beta(3)=\sum_{n=1}^{\infty} \dfrac{(-1)^n}{(2n+1)^3}
\end{equation*}
then,
\begin{equation*}
\displaystyle\int_0^1 \dfrac{\arctan x\ln x}{x}dx=-\beta(3)
\end{equation*}
(series expansion of $\dfrac{\arctan x}{x}$)
Let,
\begin{align*}
\displaystyle A&=\int_0^1 \dfrac{x\arctan x\ln x}{1+x^2}dx\\
\displaystyle B&=\int_0^1 \dfrac{\ln x \ln(1+x^2)}{1+x^2}dx\\
\displaystyle C&=\int_0^1 \dfrac{\arctan x\ln x}{1+x}dx
\end{align*}
Let $R$ the function defined on $[0;1]$ such that for all $x$ in $[0;1]$,
\begin{equation*}
\displaystyle R(x)=\int_0^x \dfrac{\ln t}{1+t^2}dt=\int_0^1 \dfrac{x\ln(tx)}{1+t^2x^2}dt
\end{equation*}
Let $\epsilon$ such that $0<\epsilon<1$.
\begin{align*}
\varphi(\epsilon)&=\int_0^{1-\epsilon}\dfrac{\ln x\ln(1-x)}{1+x^2}dx\\
&=\Big[\left(R(x)-R(1)\right)\ln(1-x)\Big]_0^{1-\epsilon}+\int_0^{1-\epsilon}\dfrac{\left(R(x)-R(1)\right)}{1-x}dx\\
&=\displaystyle \left(R(1-\epsilon)-R(1)\right)\ln(\epsilon)+R(1)\ln(\epsilon)+\int_0^{1-\epsilon}\dfrac{R(x)}{1-x}dx\\
&=\displaystyle \left(R(1-\epsilon)-R(1)\right)\ln(\epsilon)+R(1)\ln(\epsilon)+\int_0^{1-\epsilon} \left(\dfrac{x\ln(tx)}{(1+t^2x^2)(1-x)}dt\right)dx\\
&=\displaystyle \left(R(1-\epsilon)-R(1)\right)\ln(\epsilon)+R(1)\ln(\epsilon)+\int_0^{1-\epsilon} \left(\int_0^1 \dfrac{x\ln x}{(1+t^2x^2)(1-x)}dt\right)dx+\int_0^1 \left(\int_0^{1-\epsilon} \dfrac{x\ln t}{(1+t^2x^2)(1-x)}dx\right)dt\\
&=\displaystyle \left(R(1-\epsilon)-R(1)\right)\ln(\epsilon)+R(1)\ln(\epsilon)+\int_0^{1-\epsilon} \left[\dfrac{\ln x\arctan(tx)}{1-x}\right]_{t=0}^{t=1}dx+\\
&\int_0^1 \left[\dfrac{\ln t\ln(1+t^2x^2)}{2(1+t^2)}+\dfrac{t\ln t\arctan(tx)}{1+t^2}-\dfrac{\ln t\arctan(tx)}{t}-\dfrac{\ln t\ln(1-x)}{1+t^2}\right]_{x=0}^{x=1-\epsilon} dt\\
&=\displaystyle \left(R(1-\epsilon)-R(1)\right)\ln(\epsilon)+R(1)\ln(\epsilon)+\int_0^{1-\epsilon} \dfrac{\ln x\arctan(x)}{1-x}dx+\int_0^1 \dfrac{\ln t\ln\left(1+t^2(1-\epsilon)^2\right)}{2(1+t^2)}dt+\\
&\int_0^{1-\epsilon}\dfrac{t\ln t\arctan\left(t(1-\epsilon)\right)}{1+t^2}dt-\int_0^1\dfrac{\ln t\arctan\left(t(1-\epsilon)\right)}{t}dt-\int_0^1\dfrac{\ln t\ln \epsilon}{1+t^2}dt\\
&=\displaystyle \left(R(1-\epsilon)-R(1)\right)\ln(\epsilon)+\int_0^{1-\epsilon} \dfrac{\ln x\arctan(x)}{1-x}dx+\int_0^1 \dfrac{\ln t\ln\left(1+t^2(1-\epsilon)^2\right)}{2(1+t^2)}dt+\\
&\int_0^{1-\epsilon}\dfrac{t\ln t\arctan\left(t(1-\epsilon)\right)}{1+t^2}dt-\int_0^1\dfrac{\ln t\arctan\left(t(1-\epsilon)\right)}{t}dt
\end{align*}
\begin{equation*}
\displaystyle\lim_{\epsilon\rightarrow 0}\varphi(\epsilon)=\int_0^{1} \dfrac{\ln x\arctan(x)}{1-x}dx+\int_0^1 \dfrac{\ln t\ln\left(1+t^2\right)}{2(1+t^2)}dt+\int_0^{1}\dfrac{t\ln t\arctan t}{1+t^2}dt-\int_0^1\dfrac{\ln t\arctan t}{t}dt
\end{equation*}
Thus,
\begin{equation}
(1)\boxed{\displaystyle\int_0^{1}\dfrac{\ln x\ln(1-x)}{1+x^2}dx=A+\dfrac{1}{2}B+\beta(3)+\int_0^1\dfrac{\ln x\arctan x }{1-x}dx }
\end{equation}
In the following integral apply the change of variable $y=\dfrac{1-x}{1+x}$,
\begin{align*}
\displaystyle\int_0^1\dfrac{\ln x\arctan(x)}{1-x}dx&=\int_0^1 \dfrac{\left(\ln(1+x)-\ln(1-x)\right)\arctan\left(\dfrac{x-1}{x+1}\right)}{x}dx+\int_0^1 \dfrac{\left(\ln(1-x)-\ln(1+x)\right)\arctan\left(\dfrac{x-1}{x+1}\right)}{1+x}dx\\
&=\displaystyle \int_0^1 \dfrac{\ln(1+x)\arctan x}{x}dx-\int_0^1 \dfrac{\ln(1-x)\arctan x}{x}dx-\dfrac{\pi}{4}\int_0^1\dfrac{\ln\left(\tfrac{1+x}{1-x}\right)}{x}dx+\\
&\int_0^1\dfrac{\ln\left(\tfrac{1-x}{1+x}\right)\arctan\left(\tfrac{x-1}{x+1}\right)}{1+x}dx
\end{align*}
and,
\begin{align*}
\displaystyle \int_0^1 \dfrac{\ln(1+x)\arctan x}{x}dx&=\Big[\ln x\ln(1+x)\arctan x\Big]_0^1-\int_0^1\ln x\left(\dfrac{\ln(1+x)}{1+x^2}+\dfrac{\arctan x}{1+x}\right)dx\\
&=\displaystyle -\int_0^1 \dfrac{\ln x\ln(1+x) }{1+x^2}dx-C
\end{align*}
\begin{align*}
\displaystyle \int_0^1 \dfrac{\ln(1-x)\arctan x}{x}dx&=\Big[\ln x\ln(1-x)\arctan x\Big]_0^1-\int_0^1\ln x\left(\dfrac{\ln(1-x)}{1+x^2}-\dfrac{\arctan x}{1-x}\right)dx\\
&=\displaystyle \int_0^1 \dfrac{\ln x\arctan x }{1-x}dx-\int_0^1 \dfrac{\ln x\ln(1-x) }{1+x^2}dx
\end{align*}
In the following integral apply the change of variable $y=\dfrac{1-x}{1+x}$,
\begin{align*}
\displaystyle \int_0^1\dfrac{\ln\left(\tfrac{1+x}{1-x}\right)}{x}dx&=-2\int_0^1\dfrac{\ln x}{1-x^2}dx\\
&=\dfrac{\pi^2}{4}
\end{align*}
In the following integral apply the change of variable $y=\dfrac{1-x}{1+x}$,
\begin{equation*}
\displaystyle \int_0^1\dfrac{\ln\left(\tfrac{1-x}{1+x}\right)\arctan\left(\tfrac{x-1}{x+1}\right)}{1+x}dx=-C
\end{equation*}
Thus,
\begin{equation*}
(2)\boxed{\displaystyle \int_0^1\dfrac{\ln x\arctan(x)}{1-x}dx=\dfrac{1}{2}\int_0^1 \dfrac{\ln x\ln(1-x)}{1+x^2}dx-\dfrac{1}{2}\int_0^1 \dfrac{\ln x\ln(1+x)}{1+x^2}dx-C-\dfrac{\pi^3}{32}}
\end{equation*}
Let $S$ the function defined on $[0;1]$ such that for all $x$ in $[0;1]$,,
\begin{align*}
\displaystyle S(x)&=\int_0^x\dfrac{\ln x}{1+x^2}dt\\
&=\int_0^1\dfrac{x\ln(tx)}{1+t^2x^2}dt
\end{align*}
Note that $S(1)=-G$, $G$ being the Catalan constant,
\begin{align*}
\int_0^1\dfrac{\ln x\ln(1+x)}{1+x^2}dx&=\Big[S(x)\ln(1+x)\Big]_0^1-\int_0^1 \dfrac{S(x)}{1+x}dx\\
&=-G\ln 2-\int_0^1 \int_0^1 \dfrac{x\ln(tx)}{(1+t^2x^2)(1+x)}dtdx\\
&=-G\ln 2-\int_0^1 \int_0^1 \dfrac{x\ln x}{(1+t^2x^2)(1+x)}dtdx-\!\!\int_0^1 \int_0^1 \dfrac{x\ln t}{(1+t^2x^2)(1+x)}dtdx\\
&=-G\ln 2-\!\!\!\int_0^1\left[\dfrac{\ln x\arctan(tx)}{1+x}\right]_{t=0}^{t=1}dx-\\
&\int_0^1\!\!\left[\dfrac{\ln t\ln(1+x^2t^2)}{2(1+t^2)}\!-\!\!\dfrac{\ln t\ln(1+x)}{1+t^2}\!-\!\!\dfrac{t\ln t\arctan(tx)}{1+t^2}\!+\!\!\dfrac{\ln t\arctan(tx)}{t}\right]_{x=0}^{x=1}\!dt\\
&=-G\ln 2-\int_0^1\dfrac{\ln x\arctan x}{1+x}dx-\dfrac{1}{2}\int_0^1\dfrac{\ln t\ln(1+t^2)}{1+t^2}dt+\\
&\ln 2\int_0^1\dfrac{\ln t}{1+t^2}dt+\int_0^1\dfrac{t\ln t\arctan t}{1+t^2}dt-\int_0^1\dfrac{\ln t\arctan t}{t}dt\\
\end{align*}
Thus,
\begin{equation*}
(3)\boxed{\displaystyle \int_0^1\dfrac{\ln x\ln(1+x)}{1+x^2}dx=A-\dfrac{1}{2}B-C-2G\ln 2+\beta(3)}
\end{equation*}
Let $T$ the function defined on $[0;1]$ such that for all $x$ in $[0;1]$,
\begin{align*}
\displaystyle T(y)&=\int_0^y \dfrac{t\ln(t)}{1+t^2}dt\\
&=\displaystyle\int_0^1 \dfrac{ty^2\ln(ty)}{1+t^2y^2}dt
\end{align*}
Note that $T(1)=-\dfrac{\pi^2}{48}$
\begin{align*}
\displaystyle A&=\Big[T(y)\arctan y\Big]_0^1-\int_0^1 \dfrac{S(y)}{1+y^2}dy\\
&=\displaystyle-\dfrac{\pi^3}{192}-\int_0^1\int_0^1\dfrac{ty^2\ln(ty)}{(1+t^2y^2)(1+y^2)}dtdy\\
\displaystyle &=-\dfrac{\pi^3}{192}-\int_0^1\int_0^1\dfrac{ty^2\ln y }{(1+t^2y^2)(1+y^2)}dtdy-\int_0^1\int_0^1\dfrac{ty^2\ln t }{(1+t^2y^2)(1+y^2)}dtdy\\
\displaystyle &=-\dfrac{\pi^3}{192}-\dfrac{1}{2}\int_0^1\left[\dfrac{\ln y \ln(1+t^2y^2)}{1+y^2}\right]_{t=0}^{t=1}dy-\\
&\dfrac{1}{2}\int_0^1 \left[\dfrac{\ln t \arctan y+\ln t\arctan(ty)}{1+t}-\dfrac{\ln t \arctan(ty)-\ln t\arctan y}{t-1}\right]_{y=0}^{y=1}dt
\end{align*}
Thus,
\begin{align*}
(4)\boxed{\displaystyle A=-\dfrac{1}{64}\pi^3-\dfrac{1}{2}B-\dfrac{1}{2}C-\dfrac{1}{2}\int_0^1\dfrac{\ln x\arctan x}{1-x}dx}
\end{align*}
Plug $\displaystyle \int_0^{1}\dfrac{\ln x\ln(1-x)}{1+x^2}dx$ from (1) and
$\displaystyle \int_0^{1}\dfrac{\ln x\ln(1+x)}{1+x^2}dx$ from (4) into (2),
Thus,
\begin{align*}
(5)\boxed{\displaystyle \int_0^1\dfrac{\ln x\arctan x}{1-x}dx=B-C+2G\ln 2-\dfrac{1}{16}\pi^3}
\end{align*}
Plug $\displaystyle \int_0^1\dfrac{\ln x\arctan x}{1-x}dx$ from (5) into (4)
Thus,
\begin{align*}
(6)\boxed{A=\dfrac{1}{64}\pi^3-B-G\ln 2}
\end{align*}
In the following integral apply the change of variable $y=\arctan x$,
\begin{align*}
(7) \boxed{\displaystyle \int_0^{\pi/4} \ln(\tan x)\ln(\cos x-\sin x)dx=\int_0^1\dfrac{\ln x\ln(1-x)}{1+x^2}dx-\dfrac{B}{2}}
\end{align*}
and, from (1), it follows that,
\begin{align*}
(8) \boxed{\displaystyle \int_0^{\pi/4} \ln(\tan x)\ln(\cos x-\sin x)dx=A+\beta(3)-\int_0^1\dfrac{\ln x\arctan x}{1-x}dx}
\end{align*}
Plug $A$ from (6) and $\displaystyle \int_0^1\dfrac{\ln \arctan x}{1-x}dx$ from (5) into (7),
\begin{align*}
(8) \boxed{\displaystyle \int_0^{\pi/4} \ln(\tan x)\ln(\cos x-\sin x)dx=-\dfrac{3}{64}\pi^3+G\ln 2+\beta(3)-C}
\end{align*}
$C$ have been already evaluated (see Evaluating $\int_0^1 \frac{\arctan x \log x}{1+x}dx$ )
\begin{equation}
\boxed{\displaystyle C=\dfrac{G\ln 2}{2}-\dfrac{\pi^3}{64}}
\end{equation}
and, knowing that,
\begin{equation}
\beta(3)=\dfrac{\pi^3}{32}
\end{equation}
Therefore,
\begin{equation}
\boxed{\displaystyle \int_0^{\pi/4} \ln(\tan x)\ln(\cos x-\sin x)dx=\dfrac{G\ln 2}{2}}
\end{equation}