Ramanujan's take on Chudnovsky series for 1/π(PI): Part 1

Be the first to leave a comment!
We have discussed a proof of Chudnovsky series for $1/\pi$ in this post based on Ramanujan's ideas as presented in this post. However a serious look at one of the pages from his lost notebook suggests that Ramanujan used a slightly different approach to obtain Chudnovsky type series and he also performed all the desired calculations needed to get the series in explicit form. This is what we intend to discuss in the current post.

A page from Ramanujan's Lost Notebook

On page 211 of his Lost Notebook (original manuscript) Ramanujan provides a table of relations between his functions $Q,R$ and values of a cryptic expression alongside certain positive integers
Page 211 from the Lost Notebook
Page 211 from the Lost Notebook
Bruce C. Berndt and his collaborators tried to make some guesses and discovered correctly the functions $Q,R$ as well the cryptic expression involved and its relation to the integers listed at the beginning of each row in the table. Based on their research let us demystify the above page and to that end we assume $n$ to be a positive integer and define \begin{align} P(q)&= 1-24\sum_{j=1}^{\infty}\frac{jq^j}{1-q^j}\tag{1}\\ Q(q)&= 1+240\sum_{j=1}^{\infty}\frac{j^3q^j}{1-q^j}\tag{2}\\ R(q)&= 1-504\sum_{j=1}^{\infty}\frac{j^5q^j}{1-q^j}\tag{3} \end{align} and $$P_n=P(-e^{-\pi\sqrt{n}}),Q_n=Q(-e^{-\pi\sqrt{n}}),R_n=R(-e^{-\pi\sqrt{n}}).\tag{4}$$ The table by Ramanujan above then consists of relations between $Q_n,R_n$ for different values of $n$ listed in first column. The cryptic expression in last column is $(\sqrt{n}P_n -(6/\pi))/\sqrt{Q_n}$. Ramanujan barely writes this expression with just first term $1$ for $P_n$ and two terms for $Q_n$. Let us then rewrite the table by Ramanujan as \begin{array}{|c|c|c|} \hline n& \text{Relation between } Q_n\text{ and } R_n&\frac{\sqrt{n} P_n-(6/\pi)}{\sqrt{Q_n}}\\ \hline 11& (8^3+27)Q_n^3-8^3R_n^2=0&\sqrt{2}\\ \hline 19&(8^3+1)Q_n^3-8^3R_n^2=0&\sqrt{6}\\ \hline 27&(40^3+9)Q_n^3-40^3R_n^2=0&3\sqrt{\frac{6}{5}}\\ \hline 43&(80^3+1)Q_n^3-80^3R_n^2=0&6\sqrt{\frac{3}{5}}\\ \hline 67&(440^3+1)Q_n^3-440^3R_n^2=0&19\sqrt{\frac{6}{55}}\\ \hline 163&(53360^3+1)Q_n^3-53360^3R_n^2=0&362\sqrt{\frac{3}{3335}}\\ \hline 35&((60+28\sqrt{5})^3+27)Q_n^3-(60+28\sqrt{5})^3R_n^2=0&(2+\sqrt{5})\sqrt{\frac{2}{\sqrt{5}}} \\ \hline \end{array} The table thus deals with values of Ramanujan's functions $P(-q),Q(-q),R(-q)$ with nome $q=e^{-\pi\sqrt{n}}$ for certain integer values of $n$. These functions can be evaluated in terms of elliptic moduli and integrals associated with nome $q$ as follows \begin{align} P(-q)&=\left(\frac{2K(k)}{\pi}\right)^2\left(\frac{6E(k)}{K(k)}+4k^2-5\right)\tag{5a}\\ Q(-q)&=\left(\frac{2K(k)}{\pi}\right)^4\left(1-4G^{-24}\right)\tag{5b}\\ R(-q)&=\left(\frac{2K(k)}{\pi}\right)^6\left(1-2k^2\right)\left(1+8G^{-24}\right)\tag{5c}\\ \end{align} where $G=(2kk')^{-1/12}$ is one of Ramanujan's class invariants, $k\in(0,1)$ is the elliptic modulus, $k'=\sqrt{1-k^2}$ is the complementary modulus and elliptic integrals $K,E$ are given by $$K(k)=\int_0^{\pi/2}\frac{dx}{\sqrt{1-k^2\sin^2x}},E(k)=\int_0^{\pi/2}\sqrt{1-k^2\sin^2x}\,dx\tag{6}$$ and nome $q$ is related to elliptic integrals via $$q=\exp\left(-\pi\frac{K(k')}{K(k)}\right).$$ The function $P(q)$ is crucial in obtaining series for $1/\pi$ as decribed by Ramanujan in his famous paper "Modular equations and approximations to $\pi$" and I have discussed his approach in detail in my blog posts (see posts starting from here). The approach described in these posts deals with values of $P(q^2)$ and its connection with Dedekind's eta function and elliptic integral $K=K(k)$ \begin{align} \eta(q)&=q^{1/24}\prod_{j=1}^{\infty}(1-q^j)\tag{7a}\\ P(q^2)&=12q\frac{d}{dq}\log\eta(q^2)\tag{7b}\\ P(q^2)&= \left(\frac{2K}{\pi}\right)^2(1-2k^2)+\frac{3kk'^2}{2}\frac{d}{dk} \left(\frac{2K}{\pi}\right)^2\tag{7c} \end{align} We also have the relation $$2P(q^2)-P(-q)=\left(\frac{2K}{\pi}\right)^2(1-2k^2)\tag{8}$$ and thus the values of $P_n=P(-q)$ given in lost notebook can be used to calculate the values of $P(q^2)$ and one can follow the procedure mentioned in Ramanujan's paper to obtain series for $1/\pi$ by choosing a suitable hypergeometric series for $(2K/\pi)^2$. However it is better to deal with the function $P(-q)$ directly and while doing so we will also see the usefulness of the relationship between $Q_n,R_n$ given in table above.

Evaluation of $P_n=P(-q) $

Let us observe that $$P(-q) =24q\frac {d} {dq} \log\eta(-q) =q\frac{d} {dq} \log\eta^{24}(-q) =q\frac{d} {dq} \log(-\eta^{24}(-q))$$ and we have \begin{align} -\eta^{24}(-q)&=q\prod_{j\geq 1}(1-(-q)^j)^{24}\notag\\ &=q\prod_{j\geq 1}(1+q^{2j-1})^{24}(1-q^{2j})^{24}\notag\\ &=\left(q^{-1}\prod_{j\geq 1}(1+q^{2j-1})^{24}\right) \left(q^2\prod_{j\geq 1}(1-q^{2j}) ^{24}\right) \notag\\ &=2^6(2kk')^{-2}\eta^{24}(q^2)\notag\\ &=2^4(kk')^{-2}\cdot 2^{-8}(2K/\pi)^{12}(kk')^{4}\notag\\ &=2^{-4}\left(\frac{2K}{\pi}\right)^{12} (kk') ^2\notag \end{align} Using logarithmic differentiation and noting that $$\frac{dq} {dk} =\frac{\pi^2q}{2kk'^2K^2}$$ we arrive at $$P(-q) =\left(\frac {2K}{\pi}\right)^2\left(\frac{6E}{K}+4k^2-5\right)$$ mentioned earlier. A similar procedure leads us to $$P(q^2)= \left(\frac {2K}{\pi}\right)^2\left(\frac{3E}{K}+k^2-2\right)$$ and then we get $$2P(q^2)-P(-q)= \left(\frac {2K}{\pi}\right)^2\left(1-2k^2\right)$$ which is also mentioned earlier.

Let us now take another elliptic modulus $l$ which is of degree $n$ over $k$ so that its nome is $q^n$ and the corresponding elliptic integral is denoted by $L$ and $L'/L=nK'/K$. Then we have $$\eta^{24}(-q^n)=-2^{-4}(2L/\pi)^{12}(ll')^{2}$$ Next we have $$nP(-q^n) =q\frac{d} {dq} \log\eta^{24}(-q^n)$$ so that \begin{align} nP(-q^n) - P(-q) &=q\frac{d} {dq} \log\frac{\eta^{24}(-q^n)}{\eta^{24}(-q)} \notag\\ &=q\frac{dk} {dq} \frac{d} {dk} \log\left\{\left(\frac{L} {K} \right) ^{12}\left(\frac{ll'} {kk'} \right) ^2\right\} \notag\\ &=\left(\frac{2K}{\pi}\right) ^2kk'^2\frac{d} {dk} \log\left\{\left(\frac{L} {K} \right) ^{6}\frac{ll'} {kk'} \right\} \notag\\ &=\frac{4KL}{\pi^2}mkk'^2\frac{d} {dk} \log\frac{ll'} {m^6kk'}\notag \end{align} where $m=K/L$ is the multiplier. Since $m$ as well as $l$ are algebraic functions of $k$ it follows that $$nP(-q^n) - P(-q) =\frac{4KL} {\pi^2}B_n(l, k)\tag{9} $$ where $B_n(l, k) $ is an algebraic function of $l, k$. If we put $l=k'$ so that $L=K'$ and $K'/K=1/\sqrt{n}$ the $q=e^{-\pi/\sqrt {n}} $ and we get $$nP(-e^{-\pi\sqrt{n}} ) - P(-e^{-\pi/\sqrt{n}} ) =\left(\frac{2L}{\pi}\right)^2B_n\tag{10}$$ where $L$ corresponds to nome $e^{-\pi\sqrt{n}} $ and $B_n$ is some algebraic number dependent on $n$.

Next note that $$\frac{\eta^{24}(-e^{-\pi\sqrt{n}})} {\eta^{24}(-e^{-\pi/\sqrt{n}}) } =\left(\frac{K'} {K} \right) ^{12}=n^{-6}$$ and logarithmic differentiation with respect to $n$ gives us $$nP(-e^{-\pi\sqrt{n}}) +P(-e^{-\pi/\sqrt{n}}) =\frac{12\sqrt{n}}{\pi}\tag{11}$$ Using equations $(10),(11)$ we see that $$P(-e^{-\pi\sqrt{n}}) =\frac{6}{\pi\sqrt{n}}+\left(\frac{2K}{\pi} \right) ^2\frac{B_n}{2n}$$ where elliptic integral $K$ corresponds to nome $e^{-\pi\sqrt{n}} $. Ramanujan perhaps calculated certain values of this expression for different values of $n$ and figured out that multiplying the above expression by $\sqrt{n} $ and further dividing the result by $\sqrt{Q(-q)} $ to get rid of factor $(2K/\pi) ^2 $ leads to very simple algebraic numbers. And this probably led to the table of values of the expression $(\sqrt{n} P_n-(6/\pi)) /\sqrt{Q_n} $.

Ramanujan gave the value $$\frac{1}{\sqrt{Q_n}}\left(\sqrt{n}P_n-\frac{6}{\pi}\right)=\sqrt{2}$$ for $n=11$ and we show the calculations needed to obtain this value. We need two ingredients here with the first one being the modular equation of degree $11$ given by Ramanujan $$\sqrt{kl} +\sqrt{k'l'} +2(4klk'l')^{1/6}=1\tag{12}$$ The equation is symmetric in $k, l$ and hence works when modulus $l$ is of degree $11$ over modulus $k$ as well as when $k$ is of degree $11$ over $l$. Let us assume the latter here and let $k, l$ correspond to nomes $e^{-\pi\sqrt{11}},e^{-\pi/\sqrt{11}} $ respectively so that $l=k'$ and then the equation $(12)$ reduces to $$2\sqrt{kk'}+2(2kk')^{1/3}=1$$ Let us further observe that $G=G_{11}=(2kk')^{-1/12}$ and hence the above equation can be written in terms of $G$ as $$2(G^{-12}/2)^{1/2}+2G^{-4}=1$$ or $$\sqrt{2}+2G^{2}=G^6$$ so that $a=G^2=G_{11}^2$ is the root of $$f(x) = x^3-2x-\sqrt{2}\tag{13}$$ It can be checked that this equation has no roots in $\mathbb{Q} (\sqrt{2})$ and hence $a$ is of degree $3$ over $\mathbb{Q} (\sqrt{2})$ and hence of degree $6$ over $\mathbb{Q} $.

The minimal polynomial for $a$ over rationals is then $$(x^3-2x-\sqrt{2})(x^3-2x+\sqrt{2})$$ or $$x^6-4x^4+4x^2-2$$ so that $b=a^2=G_{11}^4$ is a root of $$g(x) =x^3-4x^2+4x-2\tag{14}$$ Next ingredient we need is the expression for $nP(q^{2n})-P(q^2)$ for $n=11$ provided by Ramanujan as $$nP(q^{2n})-P(q^2)=\frac{8KL}{\pi^2}\left\{2(1+kl+k'l')+\sqrt{kl} +\sqrt{k'l'} - \sqrt{klk'l'} \right\}\tag{15}$$ and here we can again assume $k$ being of degree $11$ over $l$ so that $l$ corresponds to $q$ and $k$ to $q^{n} $ and $K'/K=nL'/L$. Putting $q=e^{-\pi/\sqrt{11}}$ so that $l=k'$ and $L/K=K'/K=\sqrt {11}$ we get $$nP(q^{2n})-P(q^2)=\sqrt{11}\left(\frac{2K}{\pi}\right)^2(4(1+2kk')+4\sqrt{kk'}-2kk')$$ or $$nP(q^{2n})-P(q^2)=\sqrt{11}\left(\frac{2K}{\pi}\right)^2(4(1+G^{-12})+2\sqrt{2}G^{-6}-G^{-12})$$ or $$nP(q^{2n})-P(q^2)=\sqrt{11}\left(\frac{2K}{\pi}\right)^2(4+3b^{-3}+2\sqrt{2}a^{-3})$$ Using $\sqrt{2}=a^3-2a$ we get $$nP(q^{2n})-P(q^2)=\sqrt{11}\left(\frac{2K}{\pi}\right)^2\frac{6b^3-4b^2+3}{b^3}$$ We need another identity $$nP(q^{2n})+P(q^2)=\frac{6\sqrt{n}}{\pi}$$ which holds for $q=e^{-\pi/\sqrt{n}} $. Using this identity together with previous equation we get $$2P(e^{-2\pi\sqrt {11}})=\frac{6}{\pi\sqrt{11}}+\frac{1}{\sqrt{11}} \left(\frac{2K}{\pi}\right)^2\frac{6b^3-4b^2+3}{b^3}\tag{16}$$ Next we use the identity $$2P(q^2)-P(-q)=\left(\frac{2K}{\pi}\right)^2(1-2k^2)$$ with $q=e^{-\pi\sqrt{11}}$ to get $$P(-q)=\frac{6}{\pi\sqrt{11}}+\left(\frac {2K}{\pi}\right)^2\left(\frac{6b^3-4b^2+3}{b^3\sqrt{11}}-(1-2k^2)\right)$$ or $$\sqrt{n} P_n-\frac{6}{\pi}=\left(\frac {2K}{\pi}\right)^2\left(\frac{6b^3-4b^2+3}{b^3}-\frac{\sqrt{11(b^6-1)}}{b^3}\right)\tag{17}$$ as $$(1-2k^2)^2=1-G^{-24}=1-b^{-6}$$ To evaluate $Q_n$ we note that $$Q(-q) =\left(\frac{2K}{\pi}\right) ^4(1-4G^{-24})=\left(\frac{2K}{\pi}\right) ^4\frac{b^6-4}{b^6}$$ and hence $$\frac{1}{\sqrt{Q_n}}\left(\sqrt{n}P_n-\frac{6}{\pi}\right)=\frac{6b^3-4b^2+3-\sqrt{11(b^6-1)}}{\sqrt {b^6-4}}.$$ The above expression can be simplified using equation $(14)$ which says that $$b^3-4b^2+4b-2=0$$ or $$b(b^2+4)=4b^2+2.$$ On squaring the above equation we get $$b^2(b^2+4)^2=(4b^2+2)^2$$ so that $b^2=G_{11}^8$ is a root of $$x(x+4)^2=(4x+2)^2$$ or $$x^3-8x^2-4=0.$$ This means that $$\sqrt {b^6-4}=\sqrt{8b^4}=2\sqrt {2}b^2.$$ We can thus prove that $$\frac{6b^3-4b^2+3-\sqrt{11(b^6-1)}}{\sqrt{b^6-4}}=\sqrt {2}$$ if $$\sqrt{11(b^6-1)}=6b^3-8b^2+3$$ or $$11(b^6-1)=(6(b^3-4b^2+4b-2)+16b^2-24b+15)^2$$ or $$11(8b^4+3)=(16b^2-24b+15)^2$$ or $$88b^4+33=256b^4+576b^2+480b^2+225-768b^3-720b$$ or $$168b^4-768b^3+1056b^2-720b+192=0$$ or $$7b^4-32b^3+44b^2-30b+8=0$$ or $$(7b-4)(b^3-4b^2+4b-2)=0$$ and this is true via equation $(14)$.

The case for $n=27$ can be handled with more ease using modular equations of degree $3$. Let the moduli $k, l$ and elliptic integrals $K, L$ correspond to nomes $q^3,q$ respectively so that $k$ is of degree $3$ over $l$. Ramanujan established that $$3P(q^6)-P(q^2)=\frac{4KL}{\pi^2}\left\{1+kl+k'l'\right\}\tag{18}$$ Using Landen transformation we can replace $q^2$ by $q$, $K$ by $(1+k)K$, $L$ by $(1+l)L$, $k$ by $2\sqrt{k}/(1+k)$, $l$ by $2\sqrt{l}/(1+l)$ in above formula to get $$3P(q^3)-P(q)=\frac{4(1+k)(1+l)KL}{\pi^2}\left\{1+\frac{4\sqrt {kl}} {(1+k)(1+l)} +\frac{(1-k)(1-l)}{(1+k)(1+l)}\right\}$$ and after some simplification we get $$3P(q^3)-P(q)=\frac{8KL}{\pi^2}\left\{1+\sqrt{kl}\right\}^2$$ Next we change $q$ to $-q$. This changes $K$ to $k'K$ and $L$ to $l'L$. To see how this transforms the expression $\sqrt{kl} $ let us observe that $$\sqrt{kl} =4q\frac{\psi(q^2)\psi(q^6)}{\phi(q)\phi(q^3)}$$ where $$\phi(q) =\sum_{n\in\mathbb{Z}} q^{n^2}$$ and $$\psi(q) =\sum_{n\geq 0}q^{n(n+1)/2}$$ are Ramanujan theta functions. The function $\phi(q) $ is same as $\vartheta_3(q)$ of Jacobi and one can now observe that changing $q$ to $-q$ changes $\sqrt {kl} $ to $-\sqrt {kl} /\sqrt{k'l'} $ and hence we get $$3P(-q^3)-P(-q)=\frac{8KL}{\pi^2}\left\{\sqrt{k'l'}-\sqrt{kl}\right\}^2\tag{19}$$ We are going to use this equation repeatedly in what follows. We will also need the following formulas related to multiplier $m=L/K$ established by Ramanujan: \begin{align} m&=\left\{1+4\left(\frac{(kk')^3}{ll'}\right)^{1/4}\right\} ^{1/2}\tag{20a}\\ k^{2} &= \frac{(m - 1)^{3}(m + 3)}{16m}\tag{20b} \\ l^{2} &= \frac{(m - 1)(m + 3)^{3}}{16m^{3}}\tag{20c} \\ k'^{2} &= \frac{(m + 1)^{3}(3 - m)}{16m}\tag{20d} \\ l'^{2} &= \frac{(m + 1)(3 - m)^{3}}{16m^{3}}\tag{20e} \end{align} Next let us now use moduli related to $n=27$ and let $k, l, r, s$ be elliptic moduli and $K, L, R, S$ be elliptic integrals corresponding to nomes $q^{27},q^9,q^3,q=\exp(-\pi/3\sqrt{3})$ respectively so that $s=k', r=l'$ and further $k, l, r$ are of degree $3$ over $l, r, s$ respectively. Also we have $$l=\frac{\sqrt{3}-1}{2\sqrt{2}},l'=\frac{\sqrt{3}+1}{2\sqrt{2}}$$ so that $ll'=1/4$. Using $(19)$ repeatedly we have $$3P(-q^3)-P(-q)=\frac{8RS}{\pi^2}(\sqrt{r's'}-\sqrt{rs})^2=\frac{72KL}{\pi^2}(\sqrt{k'l'}-\sqrt{kl} )^2\tag{21}$$ and $$3P(-q^9)-P(-q^3)=\frac{8LR}{\pi^2}(\sqrt{l'r'}-\sqrt{lr})^2=0\tag{22}$$ as $lr=l'r'$ and $$3P(-q^{27})-P(-q^9)=\frac{8KL}{\pi^2}(\sqrt{k'l'}-\sqrt{kl})^2\tag{23}$$ Multiplying equation $(23)$ by $9$, $(22)$ by $3$ and then adding these to $(21)$ we get $$27P(-q^{27})-P(-q)=\frac{144KL}{\pi^2}(\sqrt{k'l'}-\sqrt{kl})^2\tag{24}$$ Noting that $k$ is of degree $3$ over $l$ we have using equations $(20b),\dots,(20e)$ $$\sqrt {kl} - \sqrt{k'l'} =\frac{(m-1)(m+3)-(m+1)(3-m)}{4m}=\frac{m^2-3}{2m}$$ and hence $$27P(-q^{27})-P(-q)=\frac{4K^2}{\pi^2}\cdot\frac{9(m^2-3)^2}{m}$$ Next we note that the value $2kk'$ equals $G_{27}^{-12}$ and it is known that $$G_{27}=\frac{\sqrt[12]{2}} {\sqrt[3]{\sqrt[3]{2}-1}} $$ and hence $$2kk'=\frac{(a-1)^4}{2}$$ where $a=2^{1/3}$. We also note that $ll'=1/4$ and use the equation $(20a) $ to get $$m^2=1+2^{7/4}(2kk')^{3/4}=1+2(a-1)^3=1+2(a^3-3a^2+3a-1)=3+6a-6a^2$$ or $$m=\sqrt{3}\sqrt{1+2a-2a^2}$$ and thus $$\frac{9(m^2-3)^2}{m}=\frac{324(a^2+2a-4)}{\sqrt{3}\sqrt{1+2a-2a^2}}$$ and finally we have $$27P(-q^{27})-P(-q)=\frac{4K^2}{\pi^2}\cdot\frac{324(a^2+2a-4)}{\sqrt{3}\sqrt{1+2a-2a^2}}\tag{25}$$ We also have the identity $$27P(-q^{27})+P(-q)=\frac{12\sqrt{n}}{\pi}=\frac{36\sqrt{3}}{\pi}$$ From the above two equations we get $$\sqrt {n} P_n-\frac{6}{\pi}=\sqrt{27}P(-q^{27})-\frac{6}{\pi}=\frac{4K^2}{\pi^2}\cdot\frac{18(a^2+2a-4)}{\sqrt{1+2a-2a^2}}$$ Next we have $$Q_n=\left(\frac{2K}{\pi}\right)^4(1-4G_{27}^{-24})$$ or $$\sqrt{Q_n} =\left(\frac{2K}{\pi}\right)^2\sqrt{1-(a-1)^8}$$ and thus $$\frac{1}{\sqrt{Q_n}}\left(\sqrt {n} P_n-\frac{6}{\pi}\right)=\frac{18(a^2+2a-4)}{\sqrt{(1+2a-2a^2) (1-(a-1)^8) }}\tag{26} $$ The square of right hand side is a rational function of $a$ and one can use the relation $a^3=2$ to reduce numerator and denominator each to quadratic polynomials in $a$. First we deal with square of denominator which equals $$(1+2a-2a^2)(1-(a-1)^8)$$ which is same as $$(1+2a-2a^2)(-a^8 + 8a^7 - 28a^6 + 56a^5 - 70a^4 + 56a^3 - 28a^2 + 8a)$$ Reducing the above using $a^3=2$ we get $$(1+2a-2a^2)(-4a^2+32a-112+112a^2-140a+112-28a^2+8a)$$ or $$20a(1+2a-2a^2)(4a-5)=20a(-8a^3 + 18a^2 - 6a - 5) $$ and the above finally equals $$60a(6a^2-2a-7)=60(6a^3-2a^2-7a)=60(12-7a-2a^2)$$ The square of numerator of right side of $(26)$ is $$324(a^2+2a-4)^2=324(a^4+4a^2+16+4a^3-16a-8a^2)$$ which reduces to $$648(12-7a-2a^2)$$ Thus the square of right hand side of $(26)$ equals $648/60=54/5$ and hence the expression on right hand side of $(26)$ equals $\sqrt{54/5}=3\sqrt{6/5}$.

We postpone the evaluation of $P_n$ for other values of $n$ to next post to keep the length of the current post at a manageable level. In the next post we will also obtain certain nice series for $1/\pi$ based on these ideas.

Print/PDF Version

0 comments :: Ramanujan's take on Chudnovsky series for 1/π(PI): Part 1

Post a Comment