previous  home  next   PDF

A Note on Evaluating Integrals by Contour Integration: Finding Residues

Michael Fowler

Contours Meet Singularities

Remember that in evaluating an integral of a function along a closed contour in the complex plane, we can always move the contour around, provided it does not encounter a point where the integrand is not analytic.

So, given an integral, usually along the real axis or part of it, we complete the contour in the complex plane, then typically shrink it down to contours around each of the singularities inside.  (For a cut plane, the procedure might be more complicated.)  The small contour around a pole contributes the residue at that pole, we add these together to get the result.

A simple pole gives a residue, an integrand 1/ z 2  doesn’t, as we can see by integrating around the unit circle:

| z |=1 dz z =2πi, | z |=1 dz z 2 =0,,

However, a second order pole will give a residue if there’s more to the function specifically, if the rest of the function has a nonzero derivative at the pole, because in that  case there is effectively a first order pole there, and

| z |=1 f( z )dz z 2 =2πi f ( 0 )

This is because we can write it as

| z |=1 ( f( 0 )+ f ( 0 )z+ )dz z 2 =2πi f ( 0 ).

Simple Examples of Contour Integration

Let’s start with

I= dx x 2 + a 2 .

 Actually this can be done directly, with the substitution x=atanθ , the result is π/a

But it’s pretty simple as a contour integral too: we have to have a closed contour to use the calculus of residues, as it’s called, but we can put a big semicircular contour in the upper half complex plane (or the lower, choose one) it will have length πR , but the integrand is or order 1/ R 2 , so this contribution from the big semicircle goes to zero on taking R

There are two poles at z=±ia , our contour encloses the one in the upper half plane only, and near that pole

1 z 2 + a 2 = 1 ( z+ia )( zia ) 1 2ia( zia )

So the result is

I=2πi( residue )= 2πi 2ia = π a .

Now let’s look at something more difficult:

e imx ( a 2 + x 2 ) 2 dx.

Again, we close the real axis contour with a large semicircle, but now we’re forced to put it in the upper half plane, because the exponential diverges in the lower half plane.

There is one pole, at z=ia .  It’s a double pole, though (I mean it goes as f( z )/ ( zia ) 2  ).  We need to factor it out like this and figure out the first derivative of the rest of the integrand, that is, f( z ) , at that point.

Putting z=ia+ζ , near ζ=0  we have the integrand

e imz ( a 2 + z 2 ) 2 = e im( ia+ζ ) ( ζ( 2ia+ζ ) ) 2 e ma ( 1+imζ ) 4 a 2 ζ 2 ( 1+ ζ 2ia ) 2 e ma 4 a 2 ζ 2 ( 1+imζ ζ ia )

So the residue is the coefficient of 1/ζ  in this expression, which is i e ma ( 1+ma ) 4 a 3 , so

e imx ( a 2 + x 2 ) 2 dx= π e ma ( 1+ma ) 2 a 3 .

Notice that taking the real part we immediately get

0 cosmx ( a 2 + x 2 ) 2 dx= π e ma ( 1+ma ) 4 a 3 .

   We could not have done a direct contour integral on cos, because it blows up in both half planes, cos for large imaginary argument being cosh for large real argument.

It  was safe to ignore the integral around the large semicircle in the above because of the denominator. In fact, there’s a result known as Jordan’s lemma that says if the integrand has the form e iλz f( z )  with λ  real and positive, and f( z )  goes uniformly to zero as z  in the upper half plane, then the large semicircle contribution goes to zero.

Trigonometric Integrals

Trigonometric integrals can often be evaluated by integrating around the unit circle, z=r e iθ , dz=izdθ , cosθ= 1 2 ( z+1/z ),sinθ= 1 2i ( z1/z ) .

For example,

0 2π sin 2 θdθ ( a+bcosθ ) = i 2b | z |=1 ( z 2 1 ) 2 dz z 2 ( z 2 +2az/b+1 ) .

This has a second-order pole at the origin to calculate the residue there we need to find the derivative of the rest of the function at the origin.

It also has two simple poles from the roots of the quadratic in the denominator, but these roots have a product of 1, so exactly one of them is inside the unit circle, we need to find the value of the rest of the function at that pole.  I’m not going to do it here it’s routine.

Many Valued Functions

Suppose we have an integral 0 x a1 Q( x )dx  where a is not an integer.  Now z a1  is a many valued function if we define it as real on the real axis, than on taking it around the origin and back to the starting point on the real axis, it will have picked up a factor e 2πi( a1 ) = e 2πia .

Let’s now take Q( x )  to be a rational function with no poles on the real axis, and such that the integrand goes to zero faster than 1/z  at infinity.

Then we extend the integral along the real axis with a complete circle at infinity, coming back along the real axis, this last contribution being the same as the original integral, but multiplied by e 2πia .

Therefore,

0 x a1 Q( x )dx = 2πi 1 e 2πia ( sum of residues in whole complex plane ).

Poles on the Real Axis

Poles on the real axis require special treatment the contour can’t go through them, that would be meaningless.  These poles come up in scattering theory, and there the physics (such as causality) tells us if the contour is to be shifted slightly above, or slightly below the pole.

The integral in practice is written as a sum of two terms: the principal part, which is defined by stopping the integral just short of the pole then picking it up at the same distance on the other side, in the limit that this distance goes to zero, and a semicircular contour around the pole, which picks up half the residue, with sign determined by which way we take it around the pole.

For example, the integrand e imz /z,m>0  has a simple pole at the origin, and P e imx x dx =πi .

(Putting the little circle around the origin either way gives the same result!  Can you see why?)

Equating imaginary parts, we find 0 sinmx x dx= 1 2 π , since the P no longer matters if the integrand is finite at the origin.

You will see in physics expressions with propagators terms in integrands like 1 x x iε .  The iε  is an infinitesimal term, whose only purpose is to tell you which way the contour goes around the pole.

previous  home  next   PDF