Mathematics Stack Exchange is a question and answer site for people studying math at any level and professionals in related fields. Join them; it only takes a minute:

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

Working with a scalar field in 2 dimensions I've come to the following integral, from which I can extract the proper ultraviolet behavior ($a \ll 1$) of the theory:

$\int_0^\infty e^{-(4+a^2)x}\left[I_0(2x)\right]^2 ds$.

It is obvious to me that this is the Laplace transform of $\left[I_0(2x)\right]^2$ evaluated at $s = (4+a^2)$. From Wikipedia I got the formula

$\int_0^\infty e^{-sx} f(x)g(x) dx = \frac{1}{2\pi i} \lim_{T \to \infty} \int_{c-iT}^{c+it} F(\sigma)G(s-\sigma) d\sigma$,

where $F(\sigma)$ and $G(\sigma)$ are the Laplace transforms of $f(x)$ and $g(x)$, respectively.

I'm encountering some trouble trying to get an analytical result for that integral. What I actually need is its $a \approx 0$ behavior, but a full analytical answer would be great.

Thanks in advance for any help!

share|cite|improve this question
up vote 6 down vote accepted

Given that $I_0(2 x)^2 = 1 + 2 x^2 + \frac{3}{2} x^4 + \frac{5}{9}x^6 + \mathcal{o}(x^6)$ we see that it is a hypergeometric function: $$ I_0(2x)^2 = {}_1F_2\left(\frac{1}{2}; 1,1; 4 x^2\right) = \sum_{n=0}^\infty \frac{\left(\frac{1}{2}\right)_n}{(1)_n (1)_n} \frac{(4 x^2)^n}{n!} = \sum_{n=0}^\infty \left(\frac{x^{n}}{n!} \right)^2 \binom{2n}{n} $$ Now, integrate term-wise: $$ \int_0^\infty \mathrm{e}^{-k x} [ I_0(2x) ]^2 \mathrm{d} x = \sum_{n=0}^\infty \frac{(2n)!}{n!^4} \int_0^\infty x^{2n} \mathrm{e}^{-k x} \mathrm{d} x = \sum_{n=0}^\infty \frac{(2n)!}{n!^4} \frac{(2n)!}{k^{2n+1}} = \frac{1}{k} \sum_{n=0}^\infty \left( \binom{2n}{n} \frac{1}{k^n} \right)^2 $$ The sum is again hypergeometric, since ratio of subsequent summands is $\frac{4 (2n+1)^2}{k^2 (n+1)^2} = \frac{16}{k^2} \frac{\left( n+ \frac{1}{2}\right)^2}{(n+1)^2} $, thus the sum equals: $$ \int_0^\infty \mathrm{e}^{-k x} [ I_0(2x) ]^2 \mathrm{d} x = \frac{1}{k} \cdot {}_2F_1\left( \frac{1}{2}, \frac{1}{2}; 1; \frac{16}{k^2}\right) = \frac{2}{ \pi k} K\left(\frac{16}{k^2}\right) $$ where $K(m)$ is a complete elliptic integral: $$K(m) = \int_0^{\pi/2} \frac{\mathrm{d} \phi}{\sqrt{1-m \cdot \sin^2(\phi)}} $$

share|cite|improve this answer
Whoa! Thanks a lot! – Misora Grilo Jul 27 '12 at 19:47
@FelipeAttanasio: so you should consider accepting this answer. – enzotib Jul 27 '12 at 21:38

I have used Mathematica, I must confess, but this integral is known in a closed form $$ \int_0^\infty dx e^{-kx}[I_0(2x)]^2=\frac{2}{k\pi}K\left(\frac{16}{k^2}\right). $$ The condition $k>4$ must hold. Here $K(m)$ is an elliptic integral given by $$ K(m) = \int_0^{\pi/2} \frac{d\theta}{\sqrt{1-m \sin^2\theta}} = \int_0^1 \frac{dt}{\sqrt{(1-t^2)(1-m t^2)}}. $$

share|cite|improve this answer
That's nice, thanks! But do you have any idea on how to prove it? I don't really need the proof, but it definately wouldn't hurt :D – Misora Grilo Jul 27 '12 at 18:21
I think Sasha's answer fits the bill. – Jon Jul 27 '12 at 19:07
You should edit your definition of the complete elliptic integral. Sasha's answer has the one being used by Mathematica. – J. M. Jul 29 '12 at 5:17
@J.M.: Thanks for your comment. I did it taking into account that I worked this out with Mathematica and held the further definition taken Wikipedia. – Jon Jul 29 '12 at 9:55

Your Answer


By posting your answer, you agree to the privacy policy and terms of service.

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