26

If $\displaystyle a_0=\frac12$ and $\displaystyle a_{n+1}=\frac{1-\sqrt{1-a_n}}{1+\sqrt{1-2a_n}}$, show that

$$\sum_{n=0}^\infty2^na_{n+1}\sqrt{a_n}=\frac{\Gamma^2(\frac14)-4\cdot\Gamma^2(\frac34)}{8\sqrt{2\pi}}$$

The closed form for this series involves the square of gamma function, hence I try to connect it with integrals, which requires us first to find the explicit form from the recursion equation. But this recusion equation is highly non-linear. I try to multiply $1-\sqrt{1-2a_n}$ to rationalize the denominator but seems no help. I also tried some non-linear sub, such as

$$\tan(x-y)=\frac{\tan x-\tan y}{1+\tan x\tan y}$$

where $x=\frac\pi4, \tan y=\sqrt{1-a_n}$, but the $2a_n$ term inside the square root kills this attempt. If let $a_n=\sin^2\theta_n$, then we get

$$1-\cos^2\theta_{n+1}=\sin^2\theta_{n+1}=\frac{1-\cos\theta_n}{1+\sqrt{\cos2\theta_n}}$$

Is there any hint? Thank you!

MathFail
  • 21,529
  • 2
    Can you tell us where the problem comes from? – Martin R Jun 19 '23 at 18:47
  • Someone shares the problem he collected in a group. – MathFail Jun 19 '23 at 18:50
  • This constant in the RHS has appeared somewhere else before it seems, though I don't know if that means there's actually a connection or if it's a coincidence: https://math.stackexchange.com/q/3070035/1104384 – Bruno B Jun 19 '23 at 19:27
  • Mathematica evaluates RHS as 0.2008398... whereas the first term in the sum is $(1-1/\sqrt{2})/\sqrt{2}\approx 0.207$, and the terms are all non-negative, so if Mathematica is correct about the numerical value, then the equality cannot be true. – Kevin.S Jun 19 '23 at 19:37
  • 2
    Found a better post to link I think: https://math.stackexchange.com/a/229265/1104384 – Bruno B Jun 19 '23 at 19:37
  • @Kevin.S For me Wolfram is saying that the RHS is over $0.3$. – Bruno B Jun 19 '23 at 19:48
  • @BrunoB I think you accidentally put $\pi$ inside the square root sign in the denominator, which is different from the conjectured closed form. – Kevin.S Jun 19 '23 at 19:50
  • @Kevin.S Ah my bad! MathFail are you sure that $\pi$ is outside the square root? Just to make sure. – Bruno B Jun 19 '23 at 19:53
  • Sorry for the typo, the $\pi$ is inside the square root, nice catch! @Kevin.S – MathFail Jun 19 '23 at 19:57
  • @BrunoB Sorry for the typo, I have fixed it. I check it numerically for $\sum_{n=0}^8 2^n a_{n+1}\sqrt{a_n}=0.35597646$ and the RHS is $0.35597932$. It should be good now. – MathFail Jun 19 '23 at 20:00
  • There is a daily comment vote limit?? Anyway, hopefully there is a link between the $a_n$s and the integrand from the post I linked, though the $2^n$ is puzzling. – Bruno B Jun 19 '23 at 20:01
  • 2
    I add that if someone manages to prove that the sum is equivalent to $$\frac{\sqrt{\pi}}{4}\sum_{n=0}^\infty(-1)^n \frac{\Gamma\left(\frac{2n+3}{4}\right)}{\Gamma\left(\frac{2n+5}{4}\right)}$$ then the proof is complete, just by using this result: https://math.stackexchange.com/q/4717394/1174256 – Zima Jun 19 '23 at 20:28
  • 1
    This is false, since all terms of $a_n$ in my OP are postive, but your referenced series is alternating which includes negative terms. @Zima – MathFail Jun 19 '23 at 21:32
  • @MathFail what do you mean by that? The two series have the same result, hence they are equal. – Zima Jun 19 '23 at 21:33
  • 1
    Equality means nothing, because you can always make two convergent series equal to each other, for example, suppose two series converge to $a$ and $b$, then $a=\frac ab b$, but it doesn't mean they have the same structure. Here the OP has all positive terms, your referenced one are alternating, they have different structure, and this difference is intrinsic @Zima – MathFail Jun 19 '23 at 21:35
  • 1
    @MathFail I honestly don't understand what you are saying. I just pointed out that maybe someone could transform your series into my series, so the result would be proved. I know my series is alternating, but series can be manipulated, so why are you against it? – Zima Jun 19 '23 at 21:37
  • 1
    $\sum (\frac12)^n=1=\frac1e\sum\frac1{n!}$ they are equal, but they have different structures and no way to manipulate. @Zima – MathFail Jun 19 '23 at 21:43
  • How is the sn function defined? The link uses sl @ParamanandSingh – MathFail Jun 21 '23 at 03:06
  • Typo $a_n=sn^2(u_n, 2)$ and $u_0=\frac{\Gamma^2(1/4)}{4\sqrt{2\pi}}$ @ParamanandSingh – MathFail Jun 21 '23 at 03:43
  • I use this https://en.wikipedia.org/wiki/Jacobi_elliptic_functions @ParamanandSingh – MathFail Jun 21 '23 at 03:58

2 Answers2

15

I got the basic idea for this post from Find the limit of $4^n a_n$, for the recurrent sequence $a_{n+1}=\frac{1-\sqrt{1-a_n}}{1+\sqrt{1+a_n}}$


The recurrence relation for $a_n$ can be more properly understood in terms of Jacobian elliptic function. I use the notation $\text{sn} (u, k) $ with $k$ as modulus. The Wikipedia entry uses $\text{sn} (u, m) $ with $m=k^2$ as parameter. We have by definition $$u=\int_0^{\text{sn}(u,k)}\frac{dx}{\sqrt{(1-x^2)(1-k^2x^2)}}\tag{1}$$ The function $\text{sn} (u, k) $ satisfies the half argument formula $$\text{sn} ^2(u/2,k)=\frac{1-\sqrt{1-\text{sn}^2(u,k)}}{1+\sqrt {1-k^2\text{sn}^2(u,k)} }\tag{2}$$ Putting $k^2=2$ and setting $a_n=\text{sn} ^2(u_n,\sqrt{2})$ we see that recursion for $a_n$ leads to $u_{n+1}=u_n/2$ via $(2)$ and hence $u_n=u_0/2^n$.

Since $a_0=1/2$ we have $\text{sn} (u_0,\sqrt{2})=1/\sqrt{2}$ and hence $$u_0=\int_0^{1/\sqrt{2}}\frac{dx}{\sqrt{(1-x^2)(1-2x^2)}}=\int_0^{\pi/4}(1-2\sin^2x)^{-1/2}\,dx$$ The integral above equals $$\frac{1}{2}\int_0^{\pi/2}(\cos x) ^{-1/2}\,dx=\frac{1}{4}B(1/4,1/2)=\frac{\Gamma^2(1/4)}{4\sqrt{2\pi}}$$ The series in question can now be rewritten as $$\sum_{n\geq 0}2^n\text{sn}(u_0/2^n)\text{sn}^2(u_0/2^{n+1})$$ and the expected sum equals $(u_0/2)-\pi/8u_0$.

The term $\pi/4u_0$ is actually the value of another integral $$\int_0^{\pi/4}(1-2\sin^2x)^{1/2}\,dx$$ which equals $B(3/4,1/2)/4=\Gamma^2(3/4)/\sqrt{2\pi}$ and hence we can observe that the expected sum of series is $$\frac{1}{2}\left(\int_0^{\pi/4}\frac{dx}{\sqrt{1-2\sin^2x}}-\int_{0}^{\pi/4}\sqrt{1-2\sin^2x}\,dx\right)=\int_0^{\pi/4}\frac{\sin^2x}{\sqrt{1-2\sin^2x}}\,dx$$


User K B Dave mentions in comments about the Jacobi Epsilon function $$\mathcal{E} (u, k) =\int_0^{\text{sn}(u,k)}\sqrt {\frac{1-k^2t^2}{1-t^2}}\,dt=\int_0^u\text{dn}^2(t,k)\,dt$$ which satisfies the relationship $$\mathcal{E} (u+v, k)=\mathcal{E}(u, k)+\mathcal{E}(v,k)-k^2\text {sn}(u, k) \text{sn} (v, k) \text{sn} (u+v, k) $$ Replacing $u, v$ both by $u/2$ and $k^2$ by $2$ we get $$2\operatorname{sn}^2(u/2)\operatorname{sn}u= 2\mathcal{E}(u/2)-\mathcal{E} (u) $$ Putting $u=u_0/2^n$ and multiplying the equation by $2^{n-1}$ we get $$2^n\operatorname{sn}^2(u_0/2^{n+1})\operatorname{sn}(u_0/2^n)=2^n\mathcal{E}(u_0/2^{n+1})-2^{n-1}\mathcal{E}(u_0/2^n)$$ and hence our series evaluates to $$\lim_{n\to\infty} 2^n\mathcal{E}(u_0/2^{n+1})-\frac{1}{2}\mathcal{E}(u_0)$$ The second term above is $$\frac{1}{2}\int_0^{\pi/4}\sqrt{1-2\sin^2t}\,dt$$ and it remains to prove that $2^n\mathcal{E}(u_0/2^{n+1})\to u_0/2$ which equivalently requires us to prove that $\mathcal{E} (u) /u\to 1$ as $u\to 0$. This is luckily an easy consequence of fundamental theorem of calculus.

I should have remembered this Jacobian Epsilon function because I had asked about a proof of its addition formula some years ago. Damn!!

  • It is not a telescoping series. Or maybe there are other ways to telescope it. – MathFail Jun 21 '23 at 06:00
  • 1
    $2 \mathcal{E}(t/2) - \mathcal{E}(t)=k^2,\mathrm{sn}^2\tfrac{t}{2},\mathrm{sn},t$, with $\mathcal{E}$ the Jacobi epsilon function, so... – K B Dave Jul 08 '23 at 22:48
  • Thanks a lot @KBDave. Sometimes the solution is right there in front of our eyes but is not visible unless one gets some guidance. – Paramanand Singh Jul 09 '23 at 02:06
  • 1
    A remark: this question doesn't rely on the extra symmetry of lemniscatic $k^2=2$, save for the reductions of the complete integrals to gamma products. An example that does: if $x_0=1$, $x_{2n+1}=2(x_{2n}+\sqrt{x_{2n}^2-1})$, $x_{2n+2}=\tfrac{1}{2}(x_{2n+1}+\sqrt{x_{2n+1}^2+4})$, then $\sum_{n=0}^{\infty}\tfrac{(-1)^n}{4^{\lfloor n/2\rfloor}x_n}=\tfrac{32\pi}{\Gamma(\tfrac{1}{4})^4}$. – K B Dave Jul 09 '23 at 03:31
  • 1
    @KBDave: that calls for an interesting separate question. If time permits, maybe you can ask it as a self answered question with all relevant context. – Paramanand Singh Jul 09 '23 at 04:33
8

Check @ParamanandSingh's initial conjecture, and show that it is NOT a telescoping series. Use this Reference Convension (RC) for Jacobi elliptic function. Let $a_n=\text{sn}^2(u_n,k=\sqrt2)$, for short, we suppress the index $k$ and denote $\text{sn}(u_n)=\text{sn}(u_n,k=\sqrt2)$, the recursion equation becomes:

$$\text{sn}(u_{n+1})=\frac{1-\sqrt{1-\text{sn}^2(u_n)}}{1+\sqrt{1-2\text{sn}^2(u_n)}}$$

Use the property (23), (24) in RC, we get

$$\text{sn}^2(u_{n+1})=\frac{1-\text{cn}(u_n)}{1+\text{dn}(u_n)}$$

Use (69) in RC, we get

$$\text{sn}^2(u_{n+1})=\text{sn}^2(\frac{u_{n}}2)\Longrightarrow u_n=\frac{1}{2}u_{n-1}\Rightarrow \boxed{a_n=\text{sn}^2(\frac{u_{0}}{2^n})}$$

where $$u_0=\int_0^{\pi/4}\frac1{\sqrt{1-2\sin^2x}}dx=\int_0^{\pi/4}\cos^{-\frac12}(2x)dx=\frac14B\left(\frac12,\frac14\right)=\frac1{4\sqrt{2\pi}}\Gamma^2(\frac14)$$

the series becomes

$$S=\sum_{n=0} 2^n \text{sn}^2(\frac{u_n}2)\text{sn}(u_{n})$$

Use (69)

$$S=\sum_{n=0} 2^n\frac{1-\text{cn}(u_n)}{1+\text{dn}(u_n)} \text{sn}(u_{n})=\sum_{n=0} \frac{2^n\text{sn}(u_{n})}{1+\text{dn}(u_n)}-\frac{2^{n-1}\cdot2\text{sn}(u_n)\text{cn}(u_{n})}{1+\text{dn}(u_n)}\tag{*}$$

Use (63)

$$2\text{sn}(u_n)\text{cn}(u_{n})=\frac{\text{sn}(u_{n-1})(1-2\text{sn}^4(u_n))}{\text{dn}(u_n)}$$

and (65)

$$1+\text{dn}(u_{n-1})=\frac{2-4\text{sn}^2(u_n)}{1-2\text{sn}^4(u_n)}$$

Plug into (*), we get

$$S=\sum_{n=0} \frac{2^n\text{sn}(u_{n})}{1+\text{dn}(u_n)}-\frac{2^{n-1}\cdot\text{sn}(u_{n-1})}{1+\text{dn}(u_{n-1})}\cdot\frac{2(1-2\text{sn}^2(u_{n}))}{\text{dn}(u_n)(1+\text{dn}(u_n))}$$

Use (24)

$$S=\sum_{n=0} \frac{2^n\text{sn}(u_{n})}{1+\text{dn}(u_n)}-\frac{2^{n-1}\cdot\text{sn}(u_{n-1})}{1+\text{dn}(u_{n-1})}\cdot\color{red}{\frac{2\text{dn}(u_{n})}{1+\text{dn}(u_n)}}\tag{**}$$

Numerically checked, (*)$=$(**).

For $n=1, n=2, n=3$, the first half of (**) are

$$~~~~~~~~~~~~~~~~~~~~~~~~~~~~~\sum_{n=1}^3 \frac{2^n\text{sn}(u_{n})}{1+\text{dn}(u_n)}=0.658552 + 0.655704 + 0.655526$$

the second half of (**) are

$$\sum_{n=1}^3 \frac{2^{n-1}\cdot\text{sn}(u_{n-1})}{1+\text{dn}(u_{n-1})}\cdot\color{red}{\frac{2\text{dn}(u_{n})}{1+\text{dn}(u_n)}}=0.553774 + 0.623159 + 0.646898$$

Unfortunately, it is NOT a telescoping series, due to the annoying red term.

MathFail
  • 21,529
  • The expected sum is not the first term of the series, whence there is no reason to expect a telescoping property. – Gary Jun 21 '23 at 06:18
  • If you read ParamanandSingh's comment (it is deleted now) and he first conjectured as a telescoping series, and I verified it is not. – MathFail Jun 21 '23 at 06:21
  • 1
    I think there should be a recursion of the form $b_n=a_{n+1}\sqrt{a_n}+2b_{n+1}$ where $2^nb_n\to 0$. This gives the sum of series as $b_0$. My guess is that $b_n=\int_0^{\sqrt{a_n}}\frac{x^2}{\sqrt{(1-x^2)(1-x^2/a_n)}},dx$ but I don't have any way to verify this guess. It may be wrong but some sort of similar integral should work. – Paramanand Singh Jun 21 '23 at 07:49
  • Yes, it means to find another way to telescope it, $2^na_{n+1}\sqrt{a_n}=2^nb_n-2^{n+1}b_{n+1}$. If this works, then the series $\sum_{n=0}2^na_{n+1}\sqrt{a_n}=b_0$ @ParamanandSingh – MathFail Jun 21 '23 at 16:13
  • This kind of recursion is already used to evaluate certain elliptic integrals and hence I am hopeful here. See this blog post of mine https://paramanands.blogspot.com/2009/08/pi-and-the-agm-evaluating-elliptic-integrals-contd.html?m=0 (search for $L(a, b) $ there). – Paramanand Singh Jun 21 '23 at 16:23
  • @ParamanandSingh This idea doesn't work, I check numerically, $b_0=0.355979, b_1=0.141370$ and $a_0=0.5, a_1=0.292893$. hence $$b_0-2b_1\neq a_1\sqrt{a_0}$$ – MathFail Jun 21 '23 at 16:46
  • Thanks for checking it. I don't have further ideas at the moment. – Paramanand Singh Jun 21 '23 at 17:30
  • I have added a bounty to draw more attention and hopefully get a nice solution. I also think that maybe you can remove gamma-function tag and add elliptic-functions (so that experts from that area can look at it). But I leave that decision to you. – Paramanand Singh Jun 26 '23 at 08:33
  • Yes, I have edited the tag. Thank you very much! @ParamanandSingh – MathFail Jun 26 '23 at 12:25
  • The problem is finally solved thanks to commenter K B Dave. See my updated answer. – Paramanand Singh Jul 09 '23 at 01:47