[Math] n expression for the inverse Fourier transform of a log-normal function

Question

i.e. is there a simple solution for the following integral?

$$\int_{-\infty}^{\infty} \exp(-\log^2(|\omega|/\omega_0)) \; \exp(i \omega t) \; d\omega$$

where $\omega_0 > 0$

Failing that, is there a nice series approximation?

Answer 1

Since mention log-normal function, you probably mean to integrate $\omega$ from 0 to $\infty$, instead of using the absolute value.

The answer to your question depends on your acceptable definition of "simple" or, in other words, on the class of functions you would like to represent the integral with. Let's do some rewriting: $$ \mathcal{F}(\omega_0, t) = \omega_0 \mathcal{I}(\omega_0 t) = \omega_0 \int_0^\infty \exp(-\log^2(\omega)) \exp( i \omega \omega_0 t) \mathrm{d} \omega $$ If you insisted on integrating over the real line while using the modulus, notice that $\int_{-\infty}^\infty \exp(-\log^2 \vert \omega/\omega_0 \vert ) \exp(i \omega t) \mathrm{d} \omega = 2 \int_{0}^\infty \exp(-\log^2( \omega/\omega_0 ) ) \cos(\omega t) \mathrm{d} \omega = 2 \omega_0 \mathfrak{Re}\left( \mathcal{I}(\omega_0 t) \right)$.

The integral $\mathcal{I}(t)$ can not be represented in terms of elementary or even hypergeometric functions.

The following formal manipulations give divergent series: $$\begin{eqnarray} \mathcal{I}(t) &=& \int_0^\infty \exp(-\log^2(\omega)) \exp( i \omega t) \mathrm{d} \omega \stackrel{\omega = \exp(u)}{=} \int_{-\infty}^\infty \exp(u-u^2) \exp(i t \mathrm{e}^u ) \mathrm{d} u \\ & \stackrel{\mathrm{formal}}{=}& \sum_{k=0}^\infty \frac{(i t)^k}{k!} \int_{-\infty}^\infty \exp( (k+1) u - u^2) \mathrm{d} u = \sum_{k=0}^\infty \frac{(i t)^k}{k!} \exp\left( \frac{k^2}{4} \right) \end{eqnarray} $$ This formal series, indeed, coincides with the formal Taylor series, since $$ \mathcal{I}^{(k)}(0) = i^k \int_0^\infty \omega^k \exp(-\log^2 \omega) \mathrm{d} \omega \stackrel{\omega = \exp(u)} = i^k \exp\left( \frac{k^2}{4} \right) $$ This means that $\mathcal{I}(t)$ is not analytic at $t=0$.

Gauss-Kronrod quadratures work fine for small $t$, but Levin collocation method might be needed for large $t$. Both methods are implemented in Mathematica: