Siguiendo los comentarios, pongo aquí algunos avances en la 2ª identidad. Los principales resultados son (1) y (2), y queda por demostrar (3). Probablemente terminaré el cálculo más tarde o pondré una recompensa si nadie se le ocurre una idea más sencilla.
Segunda identidad : Utilizando la misma notación anterior, podemos reescribir la integral como $$ A=\int_0^{\pi}\left(\int_0^{\pi}\ln I(x,y)\;dy\right)dx=\left(\int_0^{\frac{2\pi}{3}}+\int_{\frac{2\pi}{3}}^{\pi}\right)\left(\int_0^{\pi}\ln I(x,y)\;dy\right)dx.$$
- Para $x\in(0,2\pi/3)$ tenemos $|2\cos\frac{x}{2}|>1$ y podemos escribir (usando (2) de lo anterior) \begin{align} F_+(x)=\int_0^{\pi}\ln I(x,y)\;dy=\\=\int_0^{\pi}\left[\ln\left(1+\frac{1}{2\cos\frac{x}{2}}e^{-i(y-x/2)}\right)+\ln\left(1+\frac{1}{2\cos\frac{x}{2}}e^{i(y-x/2)}\right)+\ln 4\cos^2\frac{x}{2}\right]dy=\\ =i\left[\mathrm{Li}_2\left(\frac{e^{-ix/2}}{2\cos\frac{x}{2}}\right)-\mathrm{Li}_2\left(-\frac{e^{-ix/2}}{2\cos\frac{x}{2}}\right)+\mathrm{Li}_2\left(-\frac{e^{ix/2}}{2\cos\frac{x}{2}}\right)-\mathrm{Li}_2\left(\frac{e^{ix/2}}{2\cos\frac{x}{2}}\right)\right]+\\+\pi \ln 4\cos^2\frac{x}{2}. \end{align} En el último paso hemos utilizado las evaluaciones estándar de $\mathrm{Li}_2(z)$ : \begin{align} \int_0^{\pi}\ln\left(1+re^{\pm iy}\right)dy=\pm i\Bigl[\mathrm{Li}_2(r)-\mathrm{Li}_2(-r)\Bigr],\qquad |r|<1. \end{align}
- Del mismo modo, para $x\in(2\pi/3,\pi)$ tenemos $|2\cos\frac{x}{2}|<1$ para que \begin{align} F_-(x)=\int_0^{\pi}\ln I(x,y)\;dy=\qquad\\=\int_0^{\pi}\left[\ln\left(1+2\cos\frac{x}{2}e^{-i(y-x/2)}\right)+\ln\left(1+2\cos\frac{x}{2}e^{i(y-x/2)}\right)\right]dy=\qquad\\ =i\left[\mathrm{Li}_2\left(2\cos\frac{x}{2}e^{-ix/2}\right)-\mathrm{Li}_2\left(-2\cos\frac{x}{2}e^{-ix/2}\right)+\mathrm{Li}_2\left(-2\cos\frac{x}{2}e^{ix/2}\right)-\mathrm{Li}_2\left(2\cos\frac{x}{2}e^{ix/2}\right)\right]. \end{align} $$=i\left[\mathrm{Li}_2\left(1+e^{-ix}\right)-\mathrm{Li}_2\left(1+e^{ix}\right)-\mathrm{Li}_2\left(-1-e^{-ix}\right)+\mathrm{Li}_2\left(-1-e^{ix}\right)\right].$$
Las funciones $F_+(x)$ y $F_-(x)$ pueden continuarse analíticamente fuera de sus respectivos ámbitos de definición. Resulta que estas continuaciones coinciden en un neobarrio abierto del intervalo real $x\in(0,\pi)$ . Después de todo, esto no es muy sorprendente, ya que ambos $F_{\pm}(x)$ representan la misma integral $\int_0^{\pi}\ln I(x,y)\;dy$ . Para demostrar que, efectivamente, $F_+(x)=F_-(x)$ para $x\in(0,\pi)$ se deben utilizar las identidades de tipo $\mathrm{Li}_2(z)+\mathrm{Li}_2(z^{-1})=-\frac{\pi^2}{6}-\frac12 \ln^2(-z)$ . [O simplemente verificar que las derivadas (dadas por funciones elementales) de $F_{\pm}(x)$ coinciden y calculan ambas funciones en algún punto concreto].
Ahora bien, como $F_+(x)=F_-(x)$ en $x\in(0,\pi)$ podemos escribir \begin{align} A=\int_0^{\frac{2\pi}{3}}F_+(x)dx+\int_{\frac{2\pi}{3}}^{\pi}F_-(x)dx=\int_0^{\pi}F_-(x)dx=\\= i\int_0^{\pi}\left[\mathrm{Li}_2\left(1+e^{-ix}\right)-\mathrm{Li}_2\left(1+e^{ix}\right)-\mathrm{Li}_2\left(-1-e^{-ix}\right)+\mathrm{Li}_2\left(-1-e^{ix}\right)\right]\tag{1}dx. \end{align} N.B.: Dilogaritmos $\mathrm{Li}_2(z)$ se definen en su hoja principal: el plano de corte $z\in\mathbb{C}\backslash[1,\infty)$ .
Vamos a mostrar cómo calcular las dos primeras integrales de (1). Esto se puede hacer de la siguiente manera. En primer lugar, se utiliza el identidad $\mathrm{Li}_2(z)+\mathrm{Li}_2(1-z)=-\ln z\ln(1-z)+\frac{\pi^2}{6}$ para reescribir estas integrales como \begin{align} A_1=i\int_0^{\pi}\left[\mathrm{Li}_2\left(1+e^{-ix}\right)-\mathrm{Li}_2\left(1+e^{ix}\right)\right]dx=\\= i\int_0^{\pi}\left[\mathrm{Li}_2\left(-e^{ix}\right)-\mathrm{Li}_2\left(-e^{-ix}\right)-i(\pi-x)\ln4\cos^2\frac{x}{2}\right]dx. \end{align} Para los dos primeros trimestres, introduzca $z=e^{ix}$ y deformar los semicírculos de integración a los diámetros. Esto da como resultado $$i\int_0^{\pi}\left[\mathrm{Li}_2\left(-e^{ix}\right)-\mathrm{Li}_2\left(-e^{-ix}\right)\right]dx=2\int_{-1}^{1}\frac{\mathrm{Li}_2(z)}{z}dz= \frac{7}{2}\zeta(3),$$ donde la última igualdad se obtiene expandiendo $\mathrm{Li}_2(z)$ en series, calculando las integrales y utilizando esa $\zeta(3)=\sum_{k=0}^{\infty}\frac{1}{(2k+1)^3}+\frac{1}{8}\zeta(3)$ . A continuación, en $\int_0^{\pi}(\pi-x)\ln4\cos^2\frac{x}{2}\,dx$ , integrando una vez por partes para matar el término lineal, obtenemos la mitad de la integral inicial $A_1$ . Por lo tanto, $$A_1=i\int_0^{\pi}\left[\mathrm{Li}_2\left(1+e^{-ix}\right)-\mathrm{Li}_2\left(1+e^{ix}\right)\right]dx=7\zeta(3)\tag{2}$$ y ahora queda por demostrar que $$A_2=i\int_0^{\pi}\left[\mathrm{Li}_2\left(-1-e^{ix}\right)-\mathrm{Li}_2\left(-1-e^{-ix}\right)\right]dx=\frac{7}{3}\zeta(3).\tag{3} $$