18
$\begingroup$

Hi,

I asked some time ago the following question on math.stackexchange, but I ask it here too since it remains unanswered.

The question concerns a function I encountered during research :

$$f(k):= k K(k) \sinh \left(\frac{\pi}{2} \frac{K(\sqrt{1-k^2})}{K(k)}\right)$$ for $k \in (0,1)$.

Here $K$ is the Complete elliptic integral of the first kind, defined by $$K(k):= \int_{0}^{1} \frac{dt}{\sqrt{1-t^2} \sqrt{1-k^2t^2}}.$$

More specifically, my question is the following :

Is $f$ decreasing on $(0,1)$?

This seems to be true, as the graph below suggests (obtained with Maple) :

graph of $f$

In fact, as remarked by Henry Cohn, much more seems to be true : all the derivatives of $f$ seem to be negative. This can be seen by looking at the Taylor series expansion of $f$ (see the link to math.stackexchange). The Taylor series expansion seems to have all negative coefficients (except the constant term), and the coefficient of $k^{2j}$ seems to be $\pi$ times a rational number with denominator dividing $16^j$...

Any comment or relevant reference is welcome.

Thank you, Malik

EDIT (20-07-2012) It was remarked by J.M. on M.SE that $f$ can be written as $$f(k)=kK(k)\frac{1-q(k)}{2\sqrt{q(k)}},$$ where $q(k)$ is the Elliptic nome. Maybe this is useful...

$\endgroup$
6
  • $\begingroup$ I have a feeling that looking at Laplace transforms and Bernstein's theorem on completely monotonic functions along with some inverse Laplace transforms might help...but translating this feeling into a proof might be non-trivial... $\endgroup$
    – Suvrit
    Jul 19, 2012 at 19:50
  • $\begingroup$ @Suvrit : Good idea, I'll check this out.. Thank you! $\endgroup$ Jul 20, 2012 at 13:44
  • $\begingroup$ Are you sure you typed in the function all right? In Mathematica I get a very different plot. For example, at k=4/5 it tells me K(3/5)=1.94957, K(4/5)=2.25721, so f(4/5)=0.8*2.25721*Sinh[(pi/2)*1.94957/2.25721]=3.27375, while your graph suggests a value below 3. $\endgroup$
    – GH from MO
    Jul 20, 2012 at 15:20
  • 2
    $\begingroup$ I clarified: K(k)=EllipticK[k^2] in Mathematica. $\endgroup$
    – GH from MO
    Jul 20, 2012 at 15:38
  • $\begingroup$ Yes, according to reference.wolfram.com/mathematica/ref/EllipticK.html. Weird... $\endgroup$ Jul 20, 2012 at 15:47

4 Answers 4

16
$\begingroup$

First we put it in the notation of Mathematica $K(k)$ is $K(k^2)$. So our function will be $$f(k)= k K(k^2)\sinh\Bigl(\frac{\pi}{2}\frac{K(1-k^2)}{K(k^2)}\Bigr).$$ Now we change variables (W486, Whittaker, Watson p.~486) . $$k=\frac{\vartheta_2^2(q)}{\vartheta_3^2(q)} \quad (*)$$ Where$\newcommand\Z{\mathbb{Z}}$ $$\vartheta_2(q)=2q^{\frac14}(1+q^2+q^6+\cdots)= \sum_{n\in\Z}q^{(n-\frac12)^2}= 2q^{\frac14}\prod_{n=1}^\infty\{(1-q^{2n})(1+q^{2n})^2\}$$ $$\vartheta_3(q)=1+2q+2q^4+2q^9+\cdots=\sum_{n\in\Z}q^{n^2}= \prod_{n=1}^\infty\{(1-q^{2n})(1+q^{2n-1})^2\}$$ $$\vartheta_4(q)=1-2q+2q^4-2q^9+\cdots=\sum_{n\in\Z}(-1)^nq^{n^2}= \prod_{n=1}^\infty\{(1-q^{2n})(1-q^{2n-1})^2\}$$ The function of $q$ in (*) is differentiable and increasing on $(0,1)$ it is $0$ in $0$ and $1$ in $1$. (we shall write $\vartheta_j$ to denote $\vartheta_j(q)$). Since (W467) $$\vartheta_2^4+\vartheta_4^4=\vartheta_3^4$$ we have $$1-k^2=1-\frac{\vartheta_2^4}{\vartheta_3^4}=\frac{\vartheta_4^4}{\vartheta_3^4}$$ The interesting thing about this change of variables is that $$K(k^2)=K\Bigl(\frac{\vartheta_2^4}{\vartheta_3^4}\Bigr)=\frac{\pi}{2}\vartheta_3^2,\qquad K(1-k^2)=K\Bigl(\frac{\vartheta_4^4}{\vartheta_3^4}\Bigr)=\frac{\log(1/q)}{2}\vartheta_3^2.$$

Now our function is $$k K(k^2)\sinh\Bigl(\frac{\pi}{2}\frac{K(1-k^2)}{K(k^2)}\Bigr) =\frac{\pi}{4}\Bigl(\frac{1}{\sqrt{q}}-\sqrt{q}\Bigr) \vartheta_2^2$$ that must be decreasing in $q$.

We add to the above some comments:

We need to show that $f(q):=(1-q)\vartheta_2^2/4\sqrt{q}$ is decreasing for $0 < q < 1 $. But

$$f(q)=(1-q)\prod_{n=1}^\infty (1-q^{2n})^2(1+q^{2n})^4= (1-q)\prod_{n=1}^\infty (1-q^{4n})^2(1+q^{2n})^2$$ This is the same to prove that the logarithmic derivative is negative $$-\frac{1}{1-q}-\sum_{n=1}^\infty\frac{8nq^{4n-1}}{1-q^{4n}}+\sum_{n=1}^\infty \frac{4nq^{2n-1}}{1+q^{2n}}$$ multiply this by $q>0$ and expand in series $$-\sum_{m=1}^\infty q^m-\sum_{n=1}^\infty \sum_{k=1}^\infty 8nq^{4nk}+ \sum_{n=1}^\infty\sum_{k=1}^\infty (-1)^{k+1}4nq^{2nk}$$

To show that this is negative observe that the only positive terms are those in the third sum with $k=2j+1$ odd, we will pair this term with that in the first sum corresponding to the same $n$ and $k=j$, these two terms are $$-8nq^{4nj}+4nq^{2n(2j+1)}=-4nq^{4nj}(2-q^{2n})<0.$$ This leaves only the terms with $j=0$ without pair. The remaining positive terms adds to $$\sum_{n=1}^\infty 4nq^{2n}=\frac{4q^2}{(1-q^2)^2}.$$ These terms we compensate with the first sum $$-\frac{q}{1-q}+\frac{4q^2}{(1-q^2)^2}=-\frac{q(1+q)(1-q^2)-4q^2}{(1-q^2)^2}$$ This is negative for $ 0 < q < 0.295598 $.

There is an intrinsic difficulty to treat larger values of $q$.

I propose to use the modularity of the theta function:

We have the equality considering $\vartheta_j$ as functions of $q$ $$\vartheta_2(e^{-\frac{\pi}{x}})=\sqrt{x}\vartheta_4(e^{-\pi x})$$ It follows that putting $q=e^{-\frac{\pi}{x}}$ $$\frac{1}{4}\Bigl(\frac{1}{\sqrt{q}}-\sqrt{q}\Bigr)\vartheta_2^2(e^{-\pi/x})= \frac{x}{2}\sinh\frac{\pi}{2x}\, \vartheta_4^2(e^{-\pi x}).$$ We must show this function is decreasing Since $$\vartheta_4(e^{-\pi x}) = \prod_{n=1}^\infty (1-e^{-2\pi n x})(1-e^{-\pi(2n-1)x})^2$$ is almost $1$ for $x$ near infinite, we only have to show that for $x$ large $\frac{x}{2}\sinh\frac{\pi}{2x}$ is decreasing

In this way we are proving that our function decrease when the initial $q$ is near 1. This was the difficult part before.

This strategy must be sufficient.

$\endgroup$
10
  • $\begingroup$ I have a little confusion with increasing decreasing, but this must be the solution $\endgroup$
    – juan
    Jul 22, 2012 at 17:31
  • 1
    $\begingroup$ I have changed the problem in other I think more tractable. But I see that still there is something missing. $\endgroup$
    – juan
    Jul 22, 2012 at 17:47
  • $\begingroup$ Why is the function of $q$ in (*) increasing on $(0,1)$? $\endgroup$
    – GH from MO
    Jul 23, 2012 at 1:07
  • 2
    $\begingroup$ I see: $\vartheta_3(q)/\vartheta_4(q)$ is increasing as it has positive Taylor coefficients, hence $\vartheta_2^2(q)/\vartheta_3^2(q)=(1-\vartheta_4^4(q)/\vartheta_3^4(q))^{1/2}$ is also increasing. $\endgroup$
    – GH from MO
    Jul 23, 2012 at 1:48
  • 1
    $\begingroup$ We have $f(e^{-\pi/x}) = g(x) \vartheta_4^2(e^{-\pi x})$. The factor $\vartheta_4^2(e^{-\pi x})$ is almost constant $=1$. Therefore the character of $f(e^{-\pi/x}) $ will be that of $g(x)$. Of course this need careful bounds. Those of $\vartheta_4^2(e^{-\pi x})$ are easily obtained from the expression given for the $\vartheta$ functions at the start of my answer. Those of $g(x)$ must be elementary. For example the error in $\vartheta_4^2(e^{-\pi x})\approx1$ is of the order of $C e^{-\pi x}$. This will work for $x\gg1$. $\endgroup$
    – juan
    Jul 24, 2012 at 15:48
6
$\begingroup$

The following is a strengthening of juan's conclusion. We need to show that the logarithmic derivative of $$ f(q):=\left(\frac{1}{\sqrt{q}}-\sqrt{q}\right) \vartheta_2^2(q) $$ is negative. Multiplying the logarithmic derivative by $q$, this means that $$ -\frac{q}{1-q}-\sum_{n=1}^\infty\frac{8nq^{4n}}{1-q^{4n}}+\sum_{n=1}^\infty \frac{4nq^{2n}}{1+q^{2n}} < 0. $$ For $0 < q < 1$ we have $$ \sum_{n=1}^\infty\frac{8nq^{4n}}{1-q^{4n}} > \sum_{n=1}^\infty 8nq^{4n} = \frac{8q^4}{(1-q^4)^2} $$ and $$ \sum_{n=1}^\infty\frac{4nq^{2n}}{1+q^{2n}} < \sum_{n=1}^\infty 4nq^{2n} = \frac{4q^2}{(1-q^2)^2}, $$ hence it suffices to show that $$ -\frac{q}{1-q}-\frac{8q^4}{(1-q^4)^2}+\frac{4q^2}{(1-q^2)^2} < 0. $$ This holds for $0 < q < 0.37795$, hence in this range we are done.

Remark: One can generate larger ranges by keeping the first few terms in the sums, and estimating the tail similarly as above. The sums converge uniformly on any interval $[0,1-\epsilon]$, hence with a complementary argument as outlined by juan for $q\in[1-\epsilon,1]$, the above strategy should indeed work. All that is left now is numerical work, namely specifying the $\epsilon>0$ and the number of terms to be kept in the above sums.

$\endgroup$
7
  • 1
    $\begingroup$ For $q$ near $1$ this type of reasoning fails, but then we may use the other form of the function. Adding the two arguments all the result may be proved. $\endgroup$
    – juan
    Jul 23, 2012 at 17:36
  • $\begingroup$ @juan: I agree, I was typing the same as you! See the "Remark" at the end of my post. Of course the credit is yours. $\endgroup$
    – GH from MO
    Jul 23, 2012 at 17:37
  • $\begingroup$ +1, I'm convinced..! For my needs, a proof like that suffices (that is, with some estimates and then a numerical argument as in GH's remark), but I still wonder if there's something more going on. More precisely, the following question remains : Is it true that the coefficients in the Taylor series of f are all negative? (this seems to be true, and in particular would imply that all the derivative of $f$ are negative). Anyway, thank you both for this nice argument! $\endgroup$ Jul 23, 2012 at 20:43
  • $\begingroup$ @Malik: That is an interesting observation, and I have no idea how to explain it. Just to clarify: you talk about your $f(k)$, not juan's $f(q)$ which has positive coefficients as well. $\endgroup$
    – GH from MO
    Jul 23, 2012 at 20:56
  • $\begingroup$ @GH: Yes indeed, I'm talking about the function $f(k)$ in the question. $\endgroup$ Jul 23, 2012 at 21:06
4
$\begingroup$

Without the $\pi$, we have integer coefficients, paired $+$ and $-$ ... $$ \frac{1}{4}\Bigl(\frac{1}{\sqrt{q}}-\sqrt{q}\Bigr) \vartheta_2^2 = 1 - q + 2 q^{2} - 2 q^{3} + q^{4} - q^{5} + 2 q^{6} - 2 q^{7} + 2 q^{8} - 2 q^{9} + 3 q^{12} - 3 q^{13} + 2 q^{14} - 2 q^{15} + 2 q^{18} - 2 q^{19} + 2 q^{20} - 2 q^{21} + 2 q^{22} - 2 q^{23} + q^{24} - q^{25} + 2 q^{26} - 2 q^{27} + 2 q^{30} - 2 q^{31} + 4 q^{32} - 4 q^{33} + 2 q^{36} - 2 q^{37} + \operatorname{O} \bigl(q^{40}\bigr) $$

$\endgroup$
1
  • $\begingroup$ I think a better idea is to take log. But then the coefficients are also positive and negative, and interesting. $\endgroup$
    – juan
    Jul 22, 2012 at 19:55
0
$\begingroup$

This is more of a long comment to the remark in GH's answer. I did some calculations regarding the expression $$-\frac{q}{1-q} - \sum_{n=1}^{\infty} \frac{8nq^{4n}}{1-q^{4n}} + \sum_{n=1}^{\infty}\frac{4nq^{2n}}{1+q^{2n}}.$$

Keeping the first terms in the two series and estimating the remaining tail as in GH's answer, we obtain the following expression (if I didn't make any mistake..) :

$$-\frac{q}{1-q} - \sum_{n=1}^{k-1} \frac{8nq^{4n}}{1-q^{4n}} + \sum_{n=1}^{k-1}\frac{4nq^{2n}}{1+q^{2n}} + \frac{8(k-1)q^{4(k+1)}-8kq^{4k}}{(1-q^4)^2} + \frac{4kq^{2k}-4(k-1)q^{2(k+1)}}{(1-q^2)^2}.$$

Given small $\epsilon$, we want to find $k$ such that the above is negative on $(0,1-\epsilon)$.

With $k=1$, we get GH's result that the function is negative for $q \in (0,0.37795)$. With $k=8$, Maple gives that the above is negative for $q \in (0, 0.78177)$. Large values of $k$ take longer to solve explicitely, but plotting the function with $k=30$ gives a negative function for $q \in (0,\alpha)$ where $\alpha>0.9$.

$\endgroup$
2
  • $\begingroup$ Solving $\exp(-\pi/x)=0.78177$ gives us x=12.7606. The remaining portion is $12.76 <x <+\infty$. In this range the other argument I think will be relatively easy to implement. For these large values of $x$ we have $1-\vartheta_4(e^{-\pi x})=8.48 \cdot 10^{-17}$. $\endgroup$
    – juan
    Jul 24, 2012 at 16:00
  • $\begingroup$ @juan Yes, this works! $\endgroup$ Jul 24, 2012 at 18:15

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.