Inverse Laplace Transform with a coinciding pole and branch point

Recently, the following Laplace transform was asked to be inverted:

$$F(s) = s^{-a-1} e^{-s^a} $$

where $a \in (0,1)$. This is a tough problem for two reasons. One is that there is very little chance of there being an analytical result for arbitrary values of $a$. The other, however, is more subtle: there is a coinciding pole and branch point. These two facts will collide to create a very difficult situation as we figure out just what is an acceptable end result. It turns out, however, that there is such a result: even though it is not in closed form, it at least may be evaluated using some numerical scheme if need be. Further, there is at least one special case that does yield a simple result (although that result is known to Mathematica and the tables).

The goal here is to address the subtlety in formulating a solution to such a problem. We will encounter divergences along the way, and we will have to know that the divergences must cancel if there is to be a solution. Thus, rather than be put off by the divergences, we must let them be a guide to the correct solution.

Because of the branch point, the way to attack a problem like this is via Cauchy’s theorem on a properly distorted Bromwich contour. Here, we want our contour to avoid the branch point at $z=0$. This, we consider

$$\oint_C dz \, z^{-a-1} e^{-z^a} e^{z t} $$

where $a \in (0,1)$ and $C$ is the following contour:


We will define $\text{Arg}{z} \in (-\pi,\pi]$, so the branch is the negative real axis. There are $6$ pieces to this contour, $C_k$, $k \in \{1,2,3,4,5,6\}$, as follows.

$C_1$ is the contour along the line $z \in [c-i R,c+i R]$ for some large value of $R$.

$C_2$ is the contour along a circular arc of radius $R$ from the top of $C_1$ to just above the negative real axis.

$C_3$ is the contour along a line just above the negative real axis between $[-R, -\epsilon]$ for some small $\epsilon$.

$C_4$ is the contour along a circular arc of radius $\epsilon$ about the origin.

$C_5$ is the contour along a line just below the negative real axis between $[-\epsilon,-R]$.

$C_6$ is the contour along the circular arc of radius $R$ from just below the negative real axis to the bottom of $C_1$.

When $t \gt 0$, the integral over the contours $C_2$ and $C_6$ vanish in the limit as $R \to \infty$.

The contour integral is thus equal to, in this limit,

$$\int_{c-i \infty}^{c+i \infty} ds \, s^{-a-1} e^{-s^a} e^{s t} + e^{-i \pi a} \int_{\infty}^{\epsilon} dx \, x^{-a-1} e^{-e^{i \pi a} x^a} e^{-x t} \\ + i \epsilon^{-a} \int_{\pi}^{-\pi} d\phi \, e^{-i a \phi} e^{-\epsilon^a e^{i a \phi}} e^{\epsilon t e^{i \phi}} + e^{i a \pi} \int_{\epsilon}^{\infty} dx \, x^{-a-1} e^{-e^{-i \pi a} x^a} e^{-x t}$$

Note that there is an apparent singularity at $\epsilon = 0$; however, the divergences cancel in the limit as $\epsilon \to 0$.

In this limit, the third integral has the following leading behavior:

$$-i \frac{2}{a} \epsilon^{-a} \sin{\pi a} +i 2 \pi$$

Rescaling and combining the second and fourth integrals, we get for the contour integral:

$$\int_{c-i \infty}^{c+i \infty} ds \, s^{-a-1} e^{-s^a} e^{s t} -i 2 t^a \operatorname{Im}{\left [e^{-i \pi a}\int_{\epsilon t}^{\infty} \frac{du}{u^{1+a}} e^{-u} e^{-e^{i \pi a} t^{-a} u^a} \right]}-i \frac{2}{a} \epsilon^{-a} \sin{\pi a} +i 2 \pi $$

We may Taylor expand the second exponential in the integrand because it is subdominant to the first exponential (at least for the first $n$ terms, where $n$ is that largest integer such that $\lfloor n a \rfloor = 0$). We need only expand to the first two terms to treat the limit as $\epsilon \to 0$. Note that

$$\int_{\epsilon t}^{\infty} \frac{du}{u^{1+a}} e^{-u} = \frac{t^{-a}}{a} \epsilon^{-a} + \Gamma(-a) + O \left ( \epsilon^a \right ) $$

The second term produces

$$-e^{i \pi a} t^{-a}\int_{\epsilon t}^{\infty} \frac{du}{u} e^{-u} $$

Because the exponentials outside the integral cancel, the imaginary part of the term is zero. Thus, we now take the limit as $\epsilon \to 0$. Because the contour integral is zero by Cauchy’s theorem, we get for the ILT,

$$\frac1{i 2 \pi} \int_{c-i \infty}^{c+i \infty} ds \, s^{-a-1} e^{-s^a} e^{s t} =
\frac{t^a}{\pi} \operatorname{Im}{\left [e^{-i \pi a}\int_{0}^{\infty} \frac{du}{u^{1+a}} e^{-u} \left ( e^{-e^{i \pi a} t^{-a} u^a} – 1 + e^{i \pi a} t^{-a} u^a \right ) \right]} + \frac{t^a}{\Gamma(1+a)} – 1$$

Actually, I have been able to go a little further with this, although I am not sure how practical it really will turn out to be. Basically, I expanded out the integrand and took the imaginary part. Then, keeping in mind Parseval’s theorem as I employ below, I split up this integrand and set up Fourier transforms of the individual pieces. Unfortunately, one of those pieces is the one-sided exponential of an arbitrary power, which does not have an analytical transform. The other piece, however, does and, further, has finite support. Thus, we may in fact use such a formulation to produce a simple numerical scheme for evaluating the integral for arbitrary $a$.

Now, even though the above result is suitable for numerical calculation, we can illustrate the above result with an analytical example that also illustrates my above reasoning with regard to Parseval. Consider the case $a=1/2$. Subbing $u=x^2$ and taking the imaginary part of the integral, we end up with

$$\operatorname{Im}{\left [e^{-i \pi a}\int_{0}^{\infty} \frac{du}{u^{1+a}} e^{-u} \left ( e^{-e^{i \pi a} t^{-a} u^a} – 1 + e^{i \pi a} t^{-a} u^a \right ) \right]} = 2 \int_{-\infty}^{\infty} dx \, e^{-x^2} \frac{\sin^2{\beta x}}{x^2}$$

where $\beta = 1/(2 \sqrt{t})$.

The latter integral may be evaluated using Parseval’s theorem, because the individual factors of the integrand are inverse Fourier transforms of simple functions. For example,

$$\int_{-\infty}^{\infty} dx \, e^{-x^2} e^{i k x} = \sqrt{\pi} e^{-k^2/4} $$
$$\int_{-\infty}^{\infty} dx \, \frac{\sin^2{\beta x}}{x^2} e^{i k x} =\begin{cases} \pi \beta \left ( 1-\frac{|k|}{2 \beta} \right ) & |k| \lt 2 \beta \\ 0 & |k| \gt 2 \beta \end{cases}$$

The integral is then equal to

$$2 \frac1{2 \pi} \sqrt{\pi} \pi \beta \int_{-2 \beta}^{2 \beta} dk \, \left ( 1-\frac{|k|}{2 \beta} \right ) e^{-k^2/4} = \sqrt{\pi} \beta \int_0^{2 \beta} dk \, \left ( 1-\frac{k}{2 \beta} \right ) e^{-k^2/4}$$

The evaluation is fairly straightforward using the definition of the error function. The result is, for the integral,

$$2 \pi \beta \operatorname{erf}{\beta} – 2 \sqrt{\pi} \left (1-e^{-\beta^2}\right ) $$

Now plugging this back into the main result above and using $\beta = 1/(2 \sqrt{t})$, we get that

$$\begin{align}\frac1{i 2 \pi} \int_{c-i \infty}^{c+i \infty} ds \, s^{-3/2} e^{-s^{1/2}} e^{s t} &= \operatorname{erf}{\left ( \frac1{2 \sqrt{t}} \right )} + \frac{2}{\sqrt{\pi}} \sqrt{t} e^{-\frac1{4 t}} – 1 \\ &= \frac{2}{\sqrt{\pi}} \sqrt{t} e^{-\frac1{4 t}} – \operatorname{erfc}{\left ( \frac1{2 \sqrt{t}} \right )}\end{align}$$

This is a known result, and the interested reader should play around with the behavior of this ILT near $t=0$. An exercise for the reader is to contemplate higher-order poles for the LT.

The art of using the Residue Theorem in evaluating definite integrals

As many of you know, using the Residue Theorem to evaluate a definite integral involves not only choosing a contour over which to integrate a function, but also choosing a function as the integrand. Many times, this is an easy task when integrating, say, rational functions over the real line. Sometimes there are less trivial […]

A Tale of Two Integrals

Two integrals look almost the same, and even to those fairly well-versed in the art, are the same. But alas, as we shall see. Consider the integral $$I_1 = \int_0^{\pi} dx \frac{x \sin{x}}{1+\cos^2{x}} $$ This may be evaluated by subbing $x \mapsto \pi-x$ as follows: $$\begin{align} I_1 &= \int_0^{\pi} dx \frac{(\pi-x) \sin{x}}{1+\cos^2{x}} \\&= \pi \int_0^{\pi} […]

Weird integral whose simplicity is only apparent in the complex plane

Every so often we come across an integral that seems absolutely impossible on its face, but is easily attacked – in fact, is custom designed – for the residue theorem. I wonder why a first year complex analysis class doesn’t show off this case as Exhibit A in why the residue theorem is so useful. […]

Expansions of $e^x$

A very basic question was asked recently: What is a better approximation to $e^x$, the usual Taylor approximation, or a similar approximation involving $1/e^{-x}$? More precisely, given an integer $m$, which is a better approximation to $e^x$: $$f_1(x) = \sum_{k=0}^m \frac{x^k}{k!} $$ or $$f_2(x) = \frac1{\displaystyle \sum_{k=0}^m \frac{(-1)^k x^k}{k!}} $$ The answer is amazingly simple: […]

Real evaluation of an improper log integral

The problem posed in M.SE concerns real methods of evaluating $$\int_0^{\infty} dx\frac{\log(x)}{\cosh(x) \sec(x)- \tan(x)} $$ The place I started is the nifty result, proven here, that $$\frac{\sin{x}}{\cosh{t} – \cos{x}} = 2 \sum_{k=1}^{\infty} e^{-k t} \sin{k x} $$ Of course, the integral actually looks like $$\int_0^{\infty} dx \frac{\cos{x}}{\cosh{x} – \sin{x}} \log{x} $$ so we need to […]

Generalizing an already tough integral

I did the case $p=1$ here. The generalization to higher $p$ may involve higher-order derivatives as follows: $$\begin{align}K_p &= \int_0^{\pi/2} dx \frac{x^{2 p}}{1+\cos^2{x}} = \frac1{2^{4 p-1}} \int_{-\pi}^{\pi} dy \frac{y^{2 p}}{3+\cos{y}} \end{align}$$ So define, as before, $$J(a) = \int_{-\pi}^{\pi} dy \frac{e^{i a y}}{3+\cos{y}} $$ Then $$K_p = \frac{(-1)^p}{2^{4 p-1}} \left [\frac{d^{2 p}}{da^{2 p}} J(a) \right ]_{a=0}$$ […]

Mathematica v9.0.1 states that this integral does not converge.

The problem is to evaluate the following integral: $$\int_0^{\infty} dx \frac{\log(1+x)}{\left(1+x^2\right)\,\left(1+x^3\right)}$$ This one turned out to be messy but straightforward. What did surprise is the way in which I would need to employ the residue theorem. Clearly, the integral is more amenable to real methods than a contour integration. What happens, though, is that the […]

An unusually alternating sum

Problem: evaluate the following sum… $$\sum_{k=0}^{\infty}\dfrac{(-1)^{\frac{k(k+1)}{2}}}{(2k+1)^2}$$ This is unusual because the $-1$ is raised to the $k (k+1)/2$ power, rather than the usual $k$th power. On the surface, this problem may look hopeless, but really, it is all about determining the pattern of odd and even numbers from the sequence $k (k+1)/2$, which turns out […]

Algebraically difficult integral

Well, some integrals are not all that hard to evaluate in principle. The one I am posting here should be an open and shut application of the residue theorem, using the unit circle as a contour. The form of the integrand, however, should give a little pause. It turns out that actually computing residues on […]

An integral involving a quadratic phase

This one was first posted on the site Integrals and Series and was brought to my attention on M.SE by Cody. My solution involves a contour integration, although the approach is far from trivial. Yet again, the solution boils down to finding a contour and a function to integrate over the contour. The problem involves […]

An integral that illustrates the beauty of contour integration methods

I attack many of the integrals in this blog using contour integration methods, some of which are obvious and some of which take a little more imagination. Here is one that illustrates the beauty an power of such methods as much as any other integral here. The problem is to derive the integral representation $$\sin […]

Funky triple integral

Note: There was a sign error in the original solution, which I have since fixed. Problem: Evaluate $$ \int _0 ^{\infty}\int _0 ^{\infty}\int _0 ^{\infty} \log(x)\log(y)\log(z)\cos(x^2+y^2+z^2)dzdydx$$ Solution: One way to attack this is to exploit the symmetry of the integral. Start by expanding the cosine term into individual pieces, i.e., $$\begin{align}\cos{(x^2+y^2+z^2)} &= \cos{x^2} \cos{(y^2+z^2)} – […]

An odd-looking integral of a square root of trig functions

The integral to evaluate is $$\int_0^{\pi/3} \big((\sqrt{3}\cos x-\sin x)\sin x\big)^{1/2}\cos x \,dx $$ A cursory glance at this specimen leads to exasperation, as the $\pi/3$ in the limit seems arbitrary, and the integrand seems devoid of an antiderivative. It turns out, however, that one way to attack this integral is to transform it into a […]

Integral of function with deceptive triple pole

One integral posted came from Hermite’s integral representation of the Hurwitz zeta function. The integral is not the most difficult that I have ever evaluated, but is interesting from a pedagogical point of view. The problem is to evaluate $$\int_0^{\infty} dx \frac{x}{(e^x-1) (x^2+4 \pi^2)^2}$$ This integral may be done via the residue theorem, by considering […]

Deceptively Tricky Integral II

Evaluate the following integral $$\int_0^{\pi/2} dx \frac{x^2}{1+\cos^2 x}$$ This integral is deviously difficult. It may look like it has the solution I provided for this integral, but, as you will see, there is an additional wrinkle. As in the linked solution, I will express the integral in a form in which I may attack via […]

Improper integral of a high power of log

I feel like I’ve developed something new in evaluating this integral. Previously, when confronted with an integral of log to the n times a function, I considered the integral of log to the n+1 times the function over a keyhole contour. This gave me the original integral, but also all of the other integrals of […]

Two integrals, each easier than it looks

I start with an integral that stymied a bunch of people for the better part of an hour. The following solution had a lot of people slapping their heads. The problem is to evaluate, for any real $\alpha$, the following integral: $$\int_0^{\pi/2} \frac{dx}{1+\tan^{\alpha}{x}}$$ Solution: use the fact that $$\tan{\left (\frac{\pi}{2}-x\right)} = \frac{1}{\tan{x}}$$ i.e., $$\frac1{1+\tan^{\alpha}{x}} = […]

A crazy-ass integral, the evaluation of which got a lot of love at Math.SE

There are a lot of integrals posted at Math.SE. I attempt to evaluate quite a few of them. Many times, I fail spectacularly; you will typically not hear of those because I feel there is nothing to say. Occasionally, I succeed; of course, you will hear about those because I post and then post again […]

Deceptively tricky integral

Someone wanted us to evaluate the following integral: $$\int\limits_0^\pi dx{\frac{{{x^2}}}{{\sqrt 5-2\cos x}}}$$ This looks straightforward, although the factor of $x^2$ does seem to complicate things a little. Actually, a lot. The usual tricks to deal with powers leads us down a road that leads to interesting territory involving dilogarithms. Let’s define $$J(a) = \int_{-\pi}^{\pi} dx […]

Unlikely application of Parseval’s equality

A poster somehow got a numerical value for an integral and wondered if it could be expressed in terms of simple constants. Let $\operatorname{erfi}(x)$ be the imaginary error function $$\operatorname{erfi}(x)=\frac{2}{\sqrt{\pi}}\int_0^xe^{z^2}dz.$$ Consider the integral $$I=\int_0^\infty\frac{\sin(x)\ \operatorname{erfi}\left(\sqrt{x}\right)\ e^{-x\sqrt{2}}}{x}dx.$$ Its numeric value is approximately $0.625773669454426$ The question is: Is it possible to express $I$ in a closed form […]

A nifty integral worked out via contour integration

Problem: evaluate the following integral: $$\int_0^{\infty}dx \frac{e^{-a x^2(x^2-\pi^2)}\cos(2\pi a x^3)}{\cosh x} $$ To evaluate this integral, consider the following integral in the complex plane: $$\oint_C \frac{dz}{\sinh{z}} e^{-a (z^2-\pi^2/4)^2}$$ where $C=C_1+C_2+C_3+C_4+C_5+C_6$ as illustrated below: $$\int_{C_1} \frac{dz}{\sinh{z}} e^{-a (z^2-\pi^2/4)^2} = \int_{i \pi/2}^{R+i \pi/2} \frac{dx}{\sinh{x}} e^{-a (x^2-\pi^2/4)^2}$$ $$\int_{C_2} \frac{dz}{\sinh{z}} e^{-a (z^2-\pi^2/4)^2} = i\int_{\pi/2}^{-\pi/2} \frac{dy}{\sinh{(R+iy)}} e^{-a [(R+i y)^2-\pi^2/4]^2} $$ […]