Residue theorem

{{short description|Concept of complex analysis}}

{{Complex analysis sidebar}}

In complex analysis, the residue theorem, sometimes called Cauchy's residue theorem, is a powerful tool to evaluate line integrals of analytic functions over closed curves; it can often be used to compute real integrals and infinite series as well. It generalizes the Cauchy integral theorem and Cauchy's integral formula. The residue theorem should not be confused with special cases of the generalized Stokes' theorem; however, the latter can be used as an ingredient of its proof.

Statement of Cauchy's residue theorem

{{See also|Residue (complex analysis)}}

The statement is as follows:

File:Residue theorem illustration.svg

Residue theorem: Let U be a simply connected open subset of the complex plane containing a finite list of points a_1, \ldots, a_n, U_0 = U \smallsetminus \{a_1, \ldots, a_n\}, and a function f holomorphic on U_0. Letting \gamma be a closed rectifiable curve in U_0, and denoting the residue of f at each point a_k by \operatorname{Res}(f, a_k) and the winding number of \gamma around a_k by \operatorname{I}(\gamma, a_k), the line integral of f around \gamma is equal to 2\pi i times the sum of residues, each counted as many times as \gamma winds around the respective point:

\oint_\gamma f(z)\, dz = 2\pi i \sum_{k=1}^n \operatorname{I}(\gamma, a_k) \operatorname{Res}(f, a_k).

If \gamma is a positively oriented simple closed curve, \operatorname{I}(\gamma, a_k) is 1 if a_k is in the interior of \gamma and 0 if not, therefore

\oint_\gamma f(z)\, dz = 2\pi i \sum \operatorname{Res}(f, a_k)

with the sum over those a_k inside {{nobr|\gamma.{{harvnb|Whittaker|Watson|1920|loc=§6.1|page=112}}.}}

The relationship of the residue theorem to Stokes' theorem is given by the Jordan curve theorem. The general plane curve {{mvar|γ}} must first be reduced to a set of simple closed curves \{\gamma_i\} whose total is equivalent to \gamma for integration purposes; this reduces the problem to finding the integral of f\, dz along a Jordan curve \gamma_i with interior V. The requirement that f be holomorphic on U_0 = U \smallsetminus \{a_k\} is equivalent to the statement that the exterior derivative d(f\, dz) = 0 on U_0. Thus if two planar regions V and W of U enclose the same subset \{a_j\} of \{a_k\}, the regions V \smallsetminus W and W \smallsetminus V lie entirely in U_0, hence

\int_{V \smallsetminus W} d(f \, dz) - \int_{W \smallsetminus V} d(f \, dz)

is well-defined and equal to zero. Consequently, the contour integral of f\, dz along \gamma_j = \partial V is equal to the sum of a set of integrals along paths \gamma_j, each enclosing an arbitrarily small region around a single a_j — the residues of f (up to the conventional factor 2\pi i at \{a_j\}. Summing over \{\gamma_j\}, we recover the final expression of the contour integral in terms of the winding numbers \{\operatorname{I}(\gamma, a_k)\}.

In order to evaluate real integrals, the residue theorem is used in the following manner: the integrand is extended to the complex plane and its residues are computed (which is usually easy), and a part of the real axis is extended to a closed curve by attaching a half-circle in the upper or lower half-plane, forming a semicircle. The integral over this curve can then be computed using the residue theorem. Often, the half-circle part of the integral will tend towards zero as the radius of the half-circle grows, leaving only the real-axis part of the integral, the one we were originally interested in.

Calculation of residues

{{Excerpt|Residue (complex analysis)|Calculation of residues|subsections=yes}}

Examples

=An integral along the real axis=

The integral

\int_{-\infty}^\infty \frac{e^{itx}}{x^2+1}\,dx

Image:Contour example.svg

arises in probability theory when calculating the characteristic function of the Cauchy distribution. It resists the techniques of elementary calculus but can be evaluated by expressing it as a limit of contour integrals.

Suppose {{math|t > 0}} and define the contour {{mvar|C}} that goes along the real line from {{math|−a}} to {{mvar|a}} and then counterclockwise along a semicircle centered at 0 from {{mvar|a}} to {{math|−a}}. Take {{mvar|a}} to be greater than 1, so that the imaginary unit {{mvar|i}} is enclosed within the curve. Now consider the contour integral

\int_C {f(z)}\,dz = \int_C \frac{e^{itz}}{z^2+1}\,dz.

Since {{math|eitz}} is an entire function (having no singularities at any point in the complex plane), this function has singularities only where the denominator {{math|z2 + 1}} is zero. Since {{math|1=z2 + 1 = (z + i)(zi)}}, that happens only where {{math|1=z = i}} or {{math|1=z = −i}}. Only one of those points is in the region bounded by this contour. Because {{math|f(z)}} is

\begin{align}

\frac{e^{itz}}{z^2+1} & =\frac{e^{itz}}{2i}\left(\frac{1}{z-i}-\frac{1}{z+i}\right) \\

& =\frac{e^{itz}}{2i(z-i)} -\frac{e^{itz}}{2i(z+i)} ,

\end{align}

the residue of {{math|f(z)}} at {{math|1=z = i}} is

\operatorname{Res}_{z=i}f(z)=\frac{e^{-t}}{2i}.

According to the residue theorem, then, we have

\int_C f(z)\,dz=2\pi i\cdot\operatorname{Res}\limits_{z=i}f(z)=2\pi i \frac{e^{-t}}{2i} = \pi e^{-t}.

The contour {{mvar|C}} may be split into a straight part and a curved arc, so that

\int_{\mathrm{straight}} f(z)\,dz+\int_{\mathrm{arc}} f(z)\,dz=\pi e^{-t}

and thus

\int_{-a}^a f(z)\,dz =\pi e^{-t}-\int_{\mathrm{arc}} f(z)\,dz.

Using some estimations, we have

\left|\int_{\mathrm{arc}}\frac{e^{itz}}{z^2+1}\,dz\right| \leq \pi a \cdot \sup_{\text{arc}} \left| \frac{e^{itz}}{z^2+1} \right| \leq \pi a \cdot \sup_{\text{arc}} \frac{1}

z^2+1
\leq \frac{\pi a}{a^2 - 1},

and

\lim_{a \to \infty} \frac{\pi a}{a^2-1} = 0.

The estimate on the numerator follows since {{math|t > 0}}, and for complex numbers {{mvar|z}} along the arc (which lies in the upper half-plane), the argument {{mvar|φ}} of {{mvar|z}} lies between 0 and {{pi}}. So,

\left|e^{itz}\right| = \left|e^{it|z|(\cos\varphi + i\sin\varphi)}\right|=\left|e^{-t|z|\sin\varphi + it|z|\cos\varphi}\right|=e^{-t|z| \sin\varphi} \le 1.

Therefore,

\int_{-\infty}^\infty \frac{e^{itz}}{z^2+1}\,dz=\pi e^{-t}.

If {{math|t < 0}} then a similar argument with an arc {{math|{{prime|C}}}} that winds around {{math|−i}} rather than {{math|i}} shows that

Image:Contour example 2.svg

\int_{-\infty}^\infty\frac{e^{itz}}{z^2+1}\,dz=\pi e^t,

and finally we have

\int_{-\infty}^\infty\frac{e^{itz}}{z^2+1}\,dz=\pi e^{-\left|t\right|}.

(If {{math|1=t = 0}} then the integral yields immediately to elementary calculus methods and its value is {{pi}}.)

=Evaluating zeta functions=

The fact that {{math|π cot(πz)}} has simple poles with residue 1 at each integer can be used to compute the sum

\sum_{n=-\infty}^\infty f(n).

Consider, for example, {{math|1=f(z) = z−2}}. Let {{math|ΓN}} be the rectangle that is the boundary of {{math|[−N − {{sfrac|1|2}}, N + {{sfrac|1|2}}]2}} with positive orientation, with an integer {{mvar|N}}. By the residue formula,

\frac{1}{2 \pi i} \int_{\Gamma_N} f(z) \pi \cot(\pi z) \, dz = \operatorname{Res}\limits_{z = 0} + \sum_{n = -N \atop n\ne 0}^N n^{-2}.

The left-hand side goes to zero as {{math|N → ∞}} since |\cot(\pi z)| is uniformly bounded on the contour, thanks to using x = \pm \left(\frac 12 + N\right) on the left and right side of the contour, and so the integrand has order O(N^{-2}) over the entire contour. On the other hand,{{harvnb|Whittaker|Watson|1920|loc=§7.2|page=125}}. Note that the Bernoulli number B_{2n} is denoted by B_{n} in Whittaker & Watson's book.

\frac{z}{2} \cot\left(\frac{z}{2}\right) = 1 - B_2 \frac{z^2}{2!} + \cdots where the Bernoulli number B_2 = \frac{1}{6}.

(In fact, {{math|1={{sfrac|z|2}} cot({{sfrac|z|2}}) = {{sfrac|iz|1 − eiz}} − {{sfrac|iz|2}}}}.) Thus, the residue {{math|Res{{sub|1=z=0}}}} is {{math|−{{sfrac|π2|3}}}}. We conclude:

\sum_{n = 1}^\infty \frac{1}{n^2} = \frac{\pi^2}{6}

which is a proof of the Basel problem.

The same argument works for all f(x) = x^{-2n} where n is a positive integer, giving us \zeta(2n) = \frac{(-1)^{n+1}B_{2n}(2\pi)^{2n}}{2(2n)!}.The trick does not work when f(x) = x^{-2n-1}, since in this case, the residue at zero vanishes, and we obtain the useless identity 0 + \zeta(2n+1) - \zeta(2n+1) = 0.

= Evaluating Eisenstein series =

The same trick can be used to establish the sum of the Eisenstein series:\pi \cot(\pi z) = \lim_{N \to \infty} \sum_{n=-N}^N (z - n)^{-1}.

{{Math proof|title=Proof|proof=

Pick an arbitrary w \in \mathbb C\setminus \Z. As above, define

g(z) := \frac{1}{w-z} \pi \cot(\pi z)

By the Cauchy residue theorem, for all N large enough such that \Gamma_N encircles w,

\frac{1}{2 \pi i} \oint_{\Gamma_N} g(z) dz = -\pi \cot(\pi z) + \sum_{n=-N}^N \frac{1}{z-n}

It remains to prove the integral converges to zero. Since \pi\cot(\pi z) /z is an even function, and \Gamma_N is symmetric about the origin, we have \oint_{\Gamma_N} \pi\cot(\pi z) /z dz = 0, and so

\oint_{\Gamma_N} g(z) dz = \oint_{\Gamma_N} \left(\frac 1z + \frac{1}{w-z}\right) \pi\cot(\pi z)dz = -w \oint_{\Gamma_N} \frac{1}{z(z-w)} \pi\cot(\pi z) dz = O(1/N)

}}

See also

Notes

{{reflist}}

References

  • {{cite book

| last = Ahlfors

| first = Lars

| author-link = Lars Ahlfors

| title = Complex Analysis

| publisher = McGraw Hill

| year = 1979

| isbn = 0-07-085008-9

}}

  • {{cite book

| last = Lindelöf

| first = Ernst L.

| author-link = Ernst Leonard Lindelöf

| title = Le calcul des résidus et ses applications à la théorie des fonctions

| publisher = Editions Jacques Gabay

| year = 1905

| language = fr

| publication-date = 1989

| isbn = 2-87647-060-8

}}

  • {{cite book

| last1 = Mitrinović

| first1 = Dragoslav

| last2 = Kečkić

| first2 = Jovan

| title = The Cauchy method of residues: Theory and applications

| publisher = D. Reidel Publishing Company

| year = 1984

| isbn = 90-277-1623-4

}}

  • {{cite book

| last1 = Whittaker

| first1 = E. T.

| author1-link = E. T. Whittaker

| last2 = Watson

| first2 = G. N.

| author2-link = G. N. Watson

| title = A Course of Modern Analysis

| publisher = Cambridge University Press

| edition = 3rd

| year = 1920

}}