Integral $\int_{-1}^1\frac1x\sqrt{\frac{1+x}{1-x}}\ln\left(\frac{2\,x^2+2\,x+1}{2\,x^2-2\,x+1}\right) \mathrm dx$

I will transform the integral via a substitution, break it up into two pieces and recombine, perform an integration by parts, and perform another substitution to get an integral to which I know a closed form exists. From there, I use a method I know to attack the integral, but in an unusual way because of the 8th degree polynomial in the denominator of the integrand.

First sub $t=(1-x)/(1+x)$, $dt=-2/(1+x)^2 dx$ to get

$$2 \int_0^{\infty} dt \frac{t^{-1/2}}{1-t^2} \log{\left (\frac{5-2 t+t^2}{1-2 t +5 t^2} \right )} $$

Now use the symmetry from the map $t \mapsto 1/t$. Break the integral up into two as follows:

\begin{align} & 2 \int_0^{1} dt \frac{t^{-1/2}}{1-t^2} \log{\left (\frac{5-2 t+t^2}{1-2 t +5 t^2} \right )} + 2 \int_1^{\infty} dt \frac{t^{-1/2}}{1-t^2} \log{\left (\frac{5-2 t+t^2}{1-2 t +5 t^2} \right )} \\ &= 2 \int_0^{1} dt \frac{t^{-1/2}}{1-t^2} \log{\left (\frac{5-2 t+t^2}{1-2 t +5 t^2} \right )} + 2 \int_0^{1} dt \frac{t^{1/2}}{1-t^2} \log{\left (\frac{5-2 t+t^2}{1-2 t +5 t^2} \right )} \\ &= 2 \int_0^{1} dt \frac{t^{-1/2}}{1-t} \log{\left (\frac{5-2 t+t^2}{1-2 t +5 t^2} \right )} \end{align}

Sub $t=u^2$ to get

$$4 \int_0^{1} \frac{du}{1-u^2} \log{\left (\frac{5-2 u^2+u^4}{1-2 u^2 +5 u^4} \right )}$$

Integrate by parts:

$$\left [2 \log{\left (\frac{1+u}{1-u} \right )} \log{\left (\frac{5-2 u^2+u^4}{1-2 u^2 +5 u^4} \right )}\right ]_0^1 \\- 32 \int_0^1 du \frac{\left(u^5-6 u^3+u\right)}{\left(u^4-2 u^2+5\right) \left(5 u^4-2 u^2+1\right)} \log{\left (\frac{1+u}{1-u} \right )}$$

One last sub: $u=(v-1)/(v+1)$ $du=2/(v+1)^2 dv$, and finally get

$$8 \int_0^{\infty} dv \frac{(v^2-1)(v^4-6 v^2+1)}{v^8+4 v^6+70v^4+4 v^2+1} \log{v}$$

With this form, we may finally conclude that a closed form exists and apply the residue theorem to obtain it. To wit, consider the following contour integral:

$$\oint_C dz \frac{8 (z^2-1)(z^4-6 z^2+1)}{z^8+4 z^6+70z^4+4 z^2+1} \log^2{z}$$

where $C$ is a keyhole contour about the positive real axis. This contour integral is equal to (I omit the steps where I show the integral vanishes about the circular arcs)

$$-i 4 \pi \int_0^{\infty} dv \frac{8 (v^2-1)(v^4-6 v^2+1)}{v^8+4 v^6+70v^4+4 v^2+1} \log{v} + 4 \pi^2 \int_0^{\infty} dv \frac{8 (v^2-1)(v^4-6 v^2+1)}{v^8+4 v^6+70v^4+4 v^2+1}$$

It should be noted that the second integral vanishes; this may be easily seen by exploiting the symmetry about $v \mapsto 1/v$.

On the other hand, the contour integral is $i 2 \pi$ times the sum of the residues about the poles of the integrand. In general, this requires us to find the zeroes of the eight degree polynomial, which may not be possible analytically. Here, on the other hand, we have many symmetries to exploit, e.g., if $a$ is a root, then $1/a$ is a root, $-a$ is a root, and $\bar{a}$ is a root. For example, we may deduce that

$$z^8+4 z^6+70z^4+4 z^2+1 = (z^4+4 z^3+10 z^2+4 z+1) (z^4-4 z^3+10 z^2-4 z+1)$$

which exploits the $a \mapsto -a$ symmetry. Now write

$$z^4+4 z^3+10 z^2+4 z+1 = (z-a)(z-\bar{a})\left (z-\frac{1}{a}\right )\left (z-\frac{1}{\bar{a}}\right )$$

Write $a=r e^{i \theta}$ and get the following equations:

$$\left ( r+\frac{1}{r}\right ) \cos{\theta}=-2$$ $$\left (r^2+\frac{1}{r^2}\right) + 4 \cos^2{\theta}=10$$

From these equations, one may deduce that a solution is $r=\phi+\sqrt{\phi}$ and $\cos{\theta}=1/\phi$, where $\phi=(1+\sqrt{5})/2$ is the golden ratio. Thus the poles take the form

$$z_k = \pm \left (\phi\pm\sqrt{\phi}\right) e^{\pm i \arctan{\sqrt{\phi}}}$$

Now we have to find the residues of the integrand at these 8 poles. We can break this task up by computing:

$$\sum_{k=1}^8 \operatorname*{Res}_{z=z_k} \left [\frac{8 (z^2-1)(z^4-6 z^2+1) \log^2{z}}{z^8+4 z^6+70z^4+4 z^2+1}\right ]=\sum_{k=1}^8 \operatorname*{Res}_{z=z_k} \left [\frac{8 (z^2-1)(z^4-6 z^2+1)}{z^8+4 z^6+70z^4+4 z^2+1}\right ] \log^2{z_k}$$

Here things got very messy, but the result is rather unbelievably simple:

$$\operatorname*{Res}_{z=z_k} \left [\frac{8 (z^2-1)(z^4-6 z^2+1)}{z^8+4 z^6+70z^4+4 z^2+1}\right ] = \text{sgn}[\cos{(\arg{z_k})}]$$

EDIT

Actually, this is a very simple computation. Inspired by @sos440, one may express the rational function of $z$ in a very simple form:

$$\frac{8 (z^2-1)(z^4-6 z^2+1)}{z^8+4 z^6+70z^4+4 z^2+1} = -\left [\frac{p'(z)}{p(z)} + \frac{p'(-z)}{p(-z)} \right ]$$

where

$$p(z)=z^4+4 z^3+10 z^2+4 z+1$$

The residue of this function at the poles are then easily seen to be $\pm 1$ according to whether the pole is a zero of $p(z)$ or $p(-z)$.

END EDIT

That is, if the pole has a positive real part, the residue of the fraction is $+1$; if it has a negative real part, the residue is $-1$.

Now consider the log piece. Expanding the square, we get 3 terms:

$$\log^2{|z_k|} - (\arg{z_k})^2 + i 2 \log{|z_k|} \arg{z_k}$$

Summing over the residues, we find that because of the $\pm1$ contributions above, that the first and third terms sum to zero. This leaves the second term. For this, it is crucial that we get the arguments right, as $\arg{z_k} \in [0,2 \pi)$. Thus, we have

$$\begin{align}I= \int_0^{\infty} dv \frac{8 (v^2-1)(v^4-6 v^2+1)}{v^8+4 v^6+70v^4+4 v^2+1} \log{v} &= \frac12 \sum_{k=1}^8 \text{sgn}[\cos{(\arg{z_k})}] (\arg{z_k})^2 \\ &= \frac12 [2 (\arctan{\sqrt{\phi}})^2 + 2 (2 \pi - \arctan{\sqrt{\phi}})^2 \\ &- 2 (\pi - \arctan{\sqrt{\phi}})^2 - 2 (\pi + \arctan{\sqrt{\phi}})^2]\\ &= 2 \pi^2 -4 \pi \arctan{\sqrt{\phi}} \\ &= 4 \pi \, \text{arccot}{\sqrt{\phi}}\\\end{align}$$


NEW ANSWER. I found yet another way of solving this problem. My new solution does not use contour integration, and is based on the following observation: for $|z| \leq 1$,

$$ - \int_{-1}^{1} \frac{1}{x} \sqrt{\frac{1+x}{1-x}} \log(1 - zx) \, dz= \pi \sin^{-1} z - \pi \log \left( \tfrac{1}{2}+\tfrac{1}{2}\sqrt{1-z^{2}} \right) . $$

As I want to keep both the old answer and the new answer, I posted my new solution to other page. You can check it here.


OLD ANSWER. Okay here is another solution. It is also related to my generalization.

We claim the following proposition:

Proposition. If $0 < r < 1$ and $r < s$, then $$ I(r, s) := \int_{-1}^{1} \frac{1}{x} \sqrt{\frac{1+x}{1-x}} \log \left( \frac{1 + 2rsx + (r^{2} + s^{2} - 1)x^{2}}{1 - 2rsx + (r^{2} + s^{2} - 1)x^{2}} \right) \, dx = 4\pi \arcsin r. \tag{1} $$

Assuming this proposition, all that we have to do is to solve the non-linear system of equations

$$ 2rs = 2 \quad \text{and} \quad r^{2} + s^{2} - 1 = 2. $$

The unique solution satisfying the condition of the proposition is $r = \phi - 1$ and $s = \phi$. So by $\text{(1)}$ we have

\begin{align*} \int_{-1}^{1} \frac{1}{x} \sqrt{\frac{1+x}{1-x}} \log \left( \frac{1 + 2x + 2x^{2}}{1 - 2x + 2x^{2}} \right) \, dx & = I(\phi-1, \phi) \\ &= 4\pi \arcsin (\phi - 1) = 4\pi \operatorname{arccot} \sqrt{\phi}. \end{align*}

Thus it remains to prove the proposition.


Proof of Proposition. We divide the proof into several steps.

Step 1. (Case reduction by analytic continuation) We first remark that given $r$ and $s$, we always have

$$ \min_{-1 \leq x \leq 1} \{ 1 \pm 2rsx + (r^{2} + s^{2} - 1)x^{2} \} > 0. \tag{2} $$

Indeed, it is not hard to check if we utilize the following equality

$$ 1 \pm 2rsx + (r^{2} + s^{2} - 1)x^{2} = (1 \pm rsx)^{2} - (1 - r^{2})(1 - s^{2}) x^{2}. $$

Then $\text{(2)}$ shows that the integrand of $I(r, s)$ remains holomoprhic under small perturbation of $s$ in $\Bbb{C}$. So it allows us to extend $s \mapsto I(r, s)$ as a holomorphic function on some open set containing the line segment $(r, \infty) \subset \Bbb{C}$. Then by the principle of analytic continuation, it is sufficient to prove that $\text{(1)}$ holds for $r < s < 1$.

Step 2. (Integral representation of $I$) Put $r = \sin \alpha$ and $s = \sin \beta$, where $ 0 < \alpha < \beta < \frac{\pi}{2}$. Then

\begin{align*} I(r, s) &= \int_{-1}^{1} \frac{1+x}{x\sqrt{1-x^{2}}} \log \left( \frac{1 + 2rsx + (r^{2} + s^{2} - 1)x^{2}}{1 - 2rsx + (r^{2} + s^{2} - 1)x^{2}} \right) \, dx \\ &= \int_{0}^{1} \frac{2}{x\sqrt{1-x^{2}}} \log \left( \frac{1 + 2rsx + (r^{2} + s^{2} - 1)x^{2}}{1 - 2rsx + (r^{2} + s^{2} - 1)x^{2}} \right) \, dx \qquad (\because \text{ parity}) \\ &= \int_{1}^{\infty} \frac{2}{\sqrt{x^{2}-1}} \log \left( \frac{x^{2} + 2rsx + (r^{2} + s^{2} - 1)}{x^{2} - 2rsx + (r^{2} + s^{2} - 1)} \right) \, dx \qquad (x \mapsto x^{-1}) \\ &= \int_{0}^{1} \frac{2}{t} \log \left( \frac{\left(t+t^{-1}\right)^{2} + 4rs\left(t+t^{-1}\right) + 4(r^{2} + s^{2} - 1)}{\left(t+t^{-1}\right)^{2} - 4rs\left(t+t^{-1}\right) + 4(r^{2} + s^{2} - 1)} \right) \, dt, \end{align*}

where in the last line we utilized the substitution $x = \frac{1}{2}(t + t^{-1})$. If we introduce the quartic polynomial \begin{align*} p(t) = t^{4} + 4rst^{3} + (4r^{2}+4s^{2}-2)t^{2} + 4rst + 1, \end{align*}

then by the property $p(1/t) = t^{-4}p(t)$, we can simplify

\begin{align*} I(r, s) &= 2 \int_{0}^{1} \frac{\log p(t) - \log p(-t)}{t} \, dt = \int_{0}^{\infty} \frac{\log p(t) - \log p(-t)}{t} \, dt \\ &= - \int_{0}^{\infty} \left( \frac{p'(t)}{p(t)} + \frac{p'(-t)}{p(-t)} \right) \log t \, dt = - \frac{1}{2} \Re \int_{-\infty}^{\infty} \left( \frac{p'(z)}{p(z)} + \frac{p'(-z)}{p(-z)} \right) \log z \, dz, \end{align*}

where we choose the branch cut of $\log$ in such a way that it avoids the upper-half plane

$$\Bbb{H} = \{ z \in \Bbb{C} : \Im z > 0 \}.$$

Step 3. (Residue calculation) Since

$$ f(z) := \left( \frac{p'(z)}{p(z)} + \frac{p'(-z)}{p(-z)} \right) \log z = O\left(\frac{\log z}{z^{2}} \right) \quad \text{as } z \to \infty, $$

by replacing the contour of integration by a semicircle of sufficiently large radius, it follows that

\begin{align*} I(r, s) = - \frac{1}{2} \Re \left\{ 2 \pi i \sum_{z_{0} \in \Bbb{H}} \operatorname{Res}_{z = z_{0}} f(z) \right\} = \pi \Im \sum_{z_{0} \in \Bbb{H}} \operatorname{Res}_{z = z_{0}} f(z). \end{align*}

(It turns out that $f(z)$ has only logarithmic singularity at the origin. So it does not account for the value of $I(r, s)$.) But by a simple calculation, together with the condition $ 0 < \alpha < \beta < \frac{\pi}{2}$, we easily notice that the zeros of $p(z)$ are exactly

$$ e^{\pm i(\alpha + \beta)} \quad \text{and} \quad -e^{\pm i(\alpha - \beta)}. $$

Now let $Z_{+}$ be the set of zeros of $p(z)$ in $\Bbb{H}$ and $Z_{-}$ be the set of zeros of $p(z)$ in $-\Bbb{H}$. Then

$$ Z_{+} = \{ e^{i(\beta+\alpha)}, -e^{-i(\beta - \alpha)} \} \quad \text{and} \quad Z_{-} = \{ e^{-i(\beta+\alpha)}, -e^{i(\beta- \alpha)} \}. $$

This in particular shows that

$$ \frac{p'(z)}{p(z)}\log z = \sum_{z_{0} \in Z_{+}} \frac{\log z}{z - z_{0}} + \text{holomorphic function on } \Bbb{H} $$

and

$$ \frac{p'(-z)}{p(-z)}\log z = -\sum_{z_{0} \in -Z_{-}} \frac{\log z}{z - z_{0}} + \text{holomorphic function on } \Bbb{H}. $$

So it follows that

\begin{align*} I(r, s) &= \pi \Im \left\{ \sum_{z_{0} \in Z_{+}} \log z_{0} - \sum_{z_{0} \in -Z_{-}} \log z_{0} \right\} \\ &= \pi \Im \left\{ \log e^{i(\beta+\alpha)} + \log e^{i(\pi-\beta+\alpha)} - \log e^{i(\pi-\beta-\alpha)} - \log e^{i(\beta-\alpha)} \right\} \\ &= \pi \Im \left\{ i(\beta+\alpha) + i(\pi-\beta+\alpha) - i(\pi-\beta-\alpha) - i(\beta-\alpha) \right\} \\ &= 4\pi \alpha = 4\pi \arcsin r. \end{align*}

This completes the proof.