# On different integral representations of $(1+x^2)^{-1/2}$ via Bessel functions

The function from the title of this question has (at least) two different integral representations in terms of Bessel function. First, we have a Fourier expansion$$\frac{1}{\sqrt{1+x^2}}= \frac{2}{\pi} \int_0^\infty K_0(k) \cos(kx) dk,$$with $$K_0$$ being the modified Bessel function of second kind. Second, we have a Laplace expansion$$\frac{1}{\sqrt{1+x^2}}= \int_0^\infty J_0(k) \exp(-k|x|) dk,$$where $$J_0$$ is the Bessel function of the first kind.

I am very comfortable with the first (Fourier) expansion, both sides of the equation look smooth at $$x=0$$ and all is good. On the other hand, the second (Laplace) expansion seems to imply non-differentiability at $$x=0$$, contrary to the LHS of the equation.

My question: Is the second expansion, in fact, valid, and, if so, how to interpret the different behavior at $$x=0$$ on both sides of the equation?

EDIT: It seems that the essence of my question has nothing to do with Bessel functions, rather, it is a question of taking the derivative of a Laplace transform. Is it possible to differentiate under the integral sign, possibly at the cost of introducing additional terms of distributional nature (delta functions) that take care of possible singularities at $$x=0$$?

• I think the usual claim is that the second displayed equation holds for all real $x\ne0$. Unlike $K_0$, the function $J_0$ is not in $L^p(\mathbb R^+)$ for $p=1$ or $p=2$, and when $x=$ the RHS must be understood as a Cauchy Principal Value or something. – kimchi lover Nov 20 '19 at 13:37
• @kimchilover Would it be possible to add some terms to the RHS so that the equation holds for $x=0$ as well? I'm thinking in terms of distributions, delta functions and so on. – Fizikus Nov 22 '19 at 10:10

You've got two questions about the absolute Laplace transform, which is defined as $$\mathcal L_a[f(k)](x)=\int^\infty_{0}f(k)e^{-k|x|}dk$$ for continuous $$f(k) \in o(e^{\delta x}), \forall\delta>0$$.

Question 1: How to interpret the different behaviour on both sides of the equation under differentiation at $$x=0$$?

In $$\mathbb R$$, if there exists a punctured neighbourhood of $$x=0$$ such that $$\mathcal L_a[f](x)$$ and $$\displaystyle{\sum^\infty_{n=0}a_n x^{2n}}$$ coincide , then both expressions have the same behaviour under differentiation at $$x=0$$, in the sense that $$\int^\infty_{0}f(k) \left( \frac{\partial}{\partial x} \right)_{\text{left}} e^{-k|x|}dk\bigg\vert_{x\to0^-} =\int^\infty_{0}f(k) \left( \frac{\partial}{\partial x} \right)_{\text{right}} e^{-k|x|}dk\bigg\vert_{x\to0^+} =\frac{d}{dx}\sum^\infty_{n=0}a_n x^{2n}\bigg\vert_{x=0}=0$$ where the subscripts denote the one-sided derivatives.

Sketch of proof:

The third equality is trivial.

To prove the right-sided derivative is zero, we want to switch the order of integral and differentiation. Here we will utilise the 'extended Leibniz integral rule':

For $$\frac{d}{dx}\int^\infty_c f(x,t)dt=\int^\infty_c \frac{\partial}{\partial x} f(x,t)dt\qquad{x\in(a,b)}$$ to be true, the sufficient conditions are

1. $$f(x,t)$$ and $$\displaystyle{\frac{\partial}{\partial x} f(x,t)}$$ are continuous in the region $$c\le t<\infty$$, $$a\le x\le b$$.

2. $$\displaystyle{\lim_{N\to\infty}\int^N_c \frac{\partial}{\partial x}f(x,t)dt}$$ converges uniformly for $$x\in(a,b)$$.

3. $$\displaystyle{\lim_{N\to\infty}\int^N_c f(x,t)dt}$$ converges for $$x\in(a,b)$$.

It is straight-forward to prove that the three conditions are satisfied for $$0 ($$r$$ is the radius of convergence of the Taylor series). Thus \begin{align} \int^\infty_{0}f(k) \left( \frac{\partial}{\partial x} \right)_{\text{right}} e^{-k|x|}dk\bigg\vert_{x\to0^+} &=\left( \frac{d}{dx} \right)_{\text{right}}\int^\infty_{0}f(k)e^{-kx}dk\bigg\vert_{x\to0^+} \\ &=\left( \frac{d}{dx} \right)_{\text{right}}\sum^\infty_{n=0}a_n x^{2n}\bigg\vert_{x\to0^+} \\ &=0 \end{align}

Similarly, the left-sided derivative is also zero.

Note: It is a little bit more complicated to show condition 2 is satisfied.

We aim to prove that for $$x>0$$, $$\lim_{N\to\infty}\int^N_{0}f(k) \left( \frac{\partial}{\partial x} \right)_{\text{right}} e^{-k|x|}dk=\lim_{N\to \infty}\int^N_{0}-kf(k) e^{-kx}dk\quad\text{converges uniformly.}$$

To this end, we make use of Cauchy criterion:

for sufficiently large $$m>n>N$$, \begin{align} \left|\int^m_{n}-kf(k) e^{-kx}dk\right| &<\int^m_{n}\left|kf(k) e^{-kx}\right|dk \\ &<\int^m_{n}e^{\delta x} e^{-kx}dk \\ &<2\cdot\frac{e^{(\delta-x)n}}{x-\delta} \\ &<2\cdot\frac{e^{(\delta-x)N}}{x-\delta} \\ &<2\cdot\frac{e^{-\Delta N}}{\Delta} \quad \text{for } x>\delta+\Delta, \Delta>0\\ \end{align}

Choosing $$N=\frac1{\Delta}\ln\frac 2{\epsilon\Delta}$$ would show uniform convergence for $$x>\delta+\Delta$$, and hence justifying the exchange of differentiation and integral for $$x>\delta+\Delta$$. By noticing that $$\delta,\Delta$$ can be chosen to be arbitrarily small, we have shown that condition 2 is satisfied for all $$x>0$$.

Question 2: Is it possible to differentiate under the integral sign, possibly at the cost of introducing additional terms of distributional nature?

Yes.

Suppose $$\mathcal L_a[f](x)$$ and $$\displaystyle{\sum^\infty_{n=0}a_n x^{2n}}$$ coincide in a punctured neighbourhood of $$x=0$$.

Then, indeed, in the sense of distribution, $$\int^\infty_0 kf(k)dk=0$$ and thus $$\frac{d}{dx}\mathcal L_a[f(k)](x)\bigg\vert_{x=0}=-\text{sgn}(x)\int^\infty_0 kf(k)e^{-k|x|}dk\bigg\vert_{x=0}=-\text{sgn}(0)\int^\infty_0 kf(k)dk=0$$

Proof:

It is well-known that $$\int^\infty_0 \delta'(s)e^{-sk}ds=k$$

Therefore, \begin{align} \int^\infty_0 kf(k)dk &=\int^\infty_0 \int^\infty_0 \delta'(s)e^{-sk} f(k) \, ds \, dk \\ &=\int^\infty_0 \int^\infty_0 \delta'(s)e^{-sk} f(k) \, dk \, ds \qquad (1)\\ &=\int^\infty_0 \delta'(s)\left(\int^\infty_0 f(k)e^{-sk} dk\right)ds \\ &=\int^\infty_0 \delta'(s)\sum^\infty_{n=0}a_n s^{2n} ds \qquad (2)\\ &=-\int^\infty_0 \delta(s)\left(\sum^\infty_{n=0}a_n s^{2n}\right)' ds \\ &=-\left(\sum^\infty_{n=0}a_n s^{2n}\right)'_{s=0} \\ &=0 \end{align} $$(1)$$: Changing order of integrals is justified by Fubini's theorem.

$$(2)$$: Due to the formula $$\displaystyle{\int^\infty_{-\infty}\delta'(x)\varphi(x)dx=-\int^\infty_{-\infty}\delta(x)\varphi'(x)dx}$$.

• How do you interpret $\mathrm{sign}(0) \times \int_0^\infty k f(k) dk$, when the integral diverges? This is the case with $f(k) = J_0(k)$. – Fizikus Nov 28 '19 at 14:23
• @Fizikus The integral diverges only in the classical sense, and converges in distributional sense, as shown in the proof in Question 2’s answer. – Szeto Nov 28 '19 at 23:24

Let's find the Fourier series for $$e^{-k|x|}$$ on the interval $$[-1,1]$$.

$$e^{-k|x|}$$ is an even function so the Fourier coefficients are $$\int_0^1 e^{-k|x|}\cos(n\pi x)\mathrm{d}x=\Re\left(\frac{e^{-n\pi-k}-1}{in\pi -k}\right)=k\frac{1-e^{-k}(-1)^n}{k^2+n^2\pi^2}$$.

The case $$n=k=0$$ can be handled by noting that $$\lim_{k\rightarrow 0} \frac{1-e^{-k}}{k}=1$$.

On the interval $$[-1,1]$$ write $$\displaystyle \frac{1}{\sqrt{1+x^2}}=\int_0^\infty J_0(k)\left(\sum_{n=0}^\infty k\frac{1-e^{-k}(-1)^n}{k^2+n^2\pi^2} \cos(n\pi x)\right) \ \mathrm{d}k$$.

Finally, let's check whether the RHS is differentiable. The Fourier series of $$e^{-k|x|}$$ can be legitimately differentiated term-by-term because it is continuous and piecewise smooth.

$$\frac{\mathrm{d}}{\mathrm{d}x}\frac{1}{\sqrt{1+x^2}} = \int_0^{\infty} \sum_n -kn\pi J_0(k) \frac{1-e^{-k}(-1)^n}{k^2+n^2\pi^2} \sin(n\pi x) \mathrm{d}k$$.

At $$x=0$$ the RHS is zero.

\begin{align} & f\left( x \right)=\int\limits_{0}^{\infty }{{{e}^{-k\left| x \right|}}{{J}_{0}}\left( k \right)dk}=\frac{1}{\pi }\int\limits_{0}^{\infty }{{{e}^{-k\left| x \right|}}dk}\int\limits_{0}^{\pi }{{{e}^{ik\cos \left( \theta \right)}}d\theta }=\frac{1}{2\pi }\int\limits_{0}^{2\pi }{\frac{d\theta }{\left| x \right|-i\cos \left( \theta \right)}} \\ & =\frac{1}{\pi i}\int\limits_{0}^{2\pi }{\frac{dz}{2\left| x \right|z-i\left( {{z}^{2}}+1 \right)}}=\frac{1}{\pi }\int\limits_{\left| z \right|=1}^{{}}{\frac{dz}{\left( z+i\left| x \right|-i\sqrt{1+{{x}^{2}}} \right)\left( z+i\left| x \right|+i\sqrt{1+{{x}^{2}}} \right)}} \\ & =\frac{1}{\sqrt{1+{{x}^{2}}}} \\ \end{align}

So the residue is essentially removing the absolute value. If you take the two cases i.e. $$x>0$$ or $$x<0$$, and the derivative, using the same procedure (you would have poles of order 2) you would get the same result as simply taking the derivative of $$1/\sqrt{1+{{x}^{2}}}$$ over those two separate domains, and because of that, both expressions appear to behave the same at $$x=0$$ (even though the derivative of the integral involving $$\left| x \right|$$ is not strictly defined there).