An integral over three Legendre polynomials

[alfredska]

Homework Statement

I encountered the following integral in my research, and I've yet to find an analytic solution:
$$I(n_1,n_2,n_3) = \int_{-1}^{1} d(\cos\theta_1) \int_{-1}^{1} d(\cos\theta_2) P_{n_1}(\cos\theta_1) P_{n_2}[\cos(\theta_1-\theta_2)] P_{n_3}(\cos\theta_2)$$
where $P_n(x)$ is the n^th Legendre polynomial.

Homework Equations

Orthogonal polynomial relation:
$$\int_{-1}^{1} d(\cos\theta) P_{n_1}(\cos\theta) P_{n_2}(\cos\theta) = \int_{0}^{\pi} d\theta (\sin\theta) P_{n_1}(\cos\theta) P_{n_2}(\cos\theta) = \frac{2}{2n_1+1}\delta_{n_1,n_2}$$

The Attempt at a Solution

Mathematica does not have much trouble in evaluating the integral exactly (Integrate, not NIntegrate) for a given set of $\{n_1, n_2, n_3\}$.

k[n1_,n2_,n3_]:=Integrate[Sin[x]Sin[y]LegendreP[n1,Cos[x]]LegendreP[n2,Cos[x-y]]LegendreP[n3,Cos[y]],{x,0,Pi},{y,0,Pi}]

I can see already that $I(n_1, n_2, n_3)$ is not diagonal.

The diagonal terms contain the factor $$\frac{2^n}{(2n+1)^2}$$
but that is not all. They also contain something that looks like a truncated double factorial: $(2n-1)!$, where the truncation point is not clear yet. And still there's more to the series that I haven't figured out:
$$I(n,n,n) = 2^n \left[\frac{2}{2n+1}\right]^2 \times \left\{1,\frac{1}{2},\frac{1}{3},\frac{1}{5},\frac{4}{35},\frac{4}{63},\frac{8}{231},\frac{8}{429}, \frac{64}{6435}, ... \right\}$$
for $n=0,1,2,...$

The off-diagonal terms appear to all contain a $\pi^2$ term, and some power of 2.

I suspect that the solution will need to be written in the form
$$I(n_1,n_2,n_3) = \delta_{n_1,n_2}\delta_{n_2,n_3}f(n_1) + (1-\delta_{n_1,n_2}\delta_{n_2,n_3})\pi^2 g(n_1,n_2,n_3)$$

 

[lippyka]

but I haven't been able to find the exact form of f and g. Any help or hints would be greatly appreciated.