An integral of a cosine multiplied by a log of a cosine My integration is rather rusty these days so I am keen to get any available help. I would like to find a general expression for the following integral (preferably in closed form, but a series representation would also be helpful)
$$I_{l,m}=\int_0^{2π}{\rm d}x\left(\int_0^{2π}{\rm d}y\ln(1-\cos(x-y))\cos\left(\frac{lx+my}{2}\right)\right)$$
where $l$ and $m$ are both integers and $\ln$ is the natural logarithm. Mathematica will evaluate the integral for individual pairs of integers $(l,m)$ but not the general integral. 
A few observations:


*

*When $l$ and $m$ are both odd, $I_{l,m}$ appears to take the form $Q+Q\ln2$, where $Q$ is a rational number, but I don’t have expressions for these in terms of $l$ and $m$.

*When $l+m=0$ with $l$ even and non-zero, $I_{l,m}$ appears to take the form $(-8π^2)/|l|$, although I haven’t proven this.
When $l=m=0$ 

*Otherwise $I_{l,m}$ appears to vanish.


Therefore it’s really the case when $l$ and $m$ are both odd that I’m interested in. The weakness of the log singularity when $x=y$ means that the integral is well-defined (I believe), but I don’t think you can use the Taylor expansion for $\ln(1+x)$ because the resulting summation doesn’t commute with the integrals (I don’t think the dominated convergence theorem is valid here). Thanks in advance for any help.
 A: Using the identity $1-\cos a=2\sin^2a/2$,
\begin{align}
 I_{l,m}&=\int_0^{2\pi}{\rm d}x\left(\int_0^{2\pi}{\rm d}y\ln(1-\cos(x-y))\cos\left(\frac{lx+my}{2}\right)\right)\\
 &=\ln 2\int_0^{2\pi}\int_0^{2\pi}\cos\left(\frac{lx+my}{2}\right)\,dxdy+\int_0^{2\pi}\int_0^{2\pi}\ln\left(\sin^2\left( \frac{x-y}{2} \right)\right)\cos\left(\frac{lx+my}{2}\right)\,dxdy\\
 &=4\ln 2\int_0^{\pi}\int_0^{\pi}\cos\left(lu+mv\right)\,dudv+4\int_0^{\pi}\int_0^{\pi}\ln\left(\sin^2(u-v)\right)\cos\left(lu+mv\right)\,dudv
\end{align}
The first integral is evaluated using the identity $\cos\left(lu+mv\right)=\cos lu\cos mv-\sin lu\sin mv$:
\begin{align}
J_{l,m}&= \int_0^{\pi}\int_0^{\pi}\cos\left(lu+mv\right)\,dudv\\
&=-\frac{4}{l m}
\end{align}
if $l$ and $m$ are both odd integers.
To evaluate the second integral we remark that
\begin{align}
 K_{l,m}&=\int_0^{\pi}\int_0^{\pi}\ln\left(\sin^2(u-v)\right)\cos\left(lu+mv\right)\,dudv\\
 &=2\int_0^{\pi}\int_0^{\pi}\ln\left(\sin\left|u-v\right|\right)\cos\left(lu+mv\right)\,dudv
\end{align}
From the Fourier expansion of $\ln\sin x$,
\begin{equation}
 \log(\sin x)=-\ln 2-\sum_{n\ge 1}\frac{\cos(2nx)}{n}
\end{equation} 
valid for $0<x<\pi$, we have
\begin{align}
 K_{l,m}=-2\ln 2 &\int_0^{\pi}\int_0^{\pi}\cos\left(lu+mv\right)\,dudv\\& -2\sum_{n\ge1}\frac{1}{n}\int_0^{\pi}\int_0^{\pi}\cos2n(u-v)\cos\left(lu+mv\right)\,dudv\\
 &=-2\ln 2 J_{l,m}-8\sum_{n\ge1}\frac{1}{n}\frac{4n^2-lm}{\left( 4n^2-l^2 \right)\left( 4n^2-m^2 \right)}
\end{align}
Here again, $l$ and $m$ are both supposed to be odd integers. Then,
\begin{equation}
 I_{l,m}=16\frac{\ln 2}{lm}-32\sum_{n\ge1}\frac{1}{n}\frac{4n^2-lm}{\left( 4n^2-l^2 \right)\left( 4n^2-m^2 \right)}
\end{equation} 
Decomposing the terms of the summation as
\begin{equation}
\frac{1}{n}\frac{4n^2-lm}{\left( 4n^2-l^2 \right)\left( 4n^2-m^2 \right)}=\frac{1}{l+m}\left[\frac{1}{n\left( 4n^2-l^2 \right)}+\frac{1}{n\left( 4n^2-m^2 \right)}\right]
\end{equation}
and, introducing the Psi function, we find
\begin{align}
\sum_{n\ge1}\frac{1}{n\left( 4n^2-l^2 \right)}&=\sum_{n\ge1}\frac{1}{2l}\left[\frac{1}{n-l/2}-\frac{1}{n}+\frac{1}{n+l/2}-\frac{1}{n}\right]\\
&=-\frac{1}{2l}\left[\psi\left( \frac{l}{2}+1 \right)+\psi\left( -\frac{l}{2}+1 \right)-2\gamma\right]
\end{align}
and a similar expression for $m$. Now, the recurrence identities and some simplifications finally give
\begin{equation}
I_{l,m}=16\frac{\ln 2}{lm}+\frac{32}{l+m}\left[\frac{1}{l}\psi\left( \frac{l}{2} \right)+\frac{1}{m}\psi\left( \frac{m}{2} \right)+\frac{\gamma l+1}{l^2}+\frac{\gamma m+1}{m^2}\right]
\end{equation} 
