Mathematical Methods Handbook

Analytic functions and Cauchy's theorem

Differentiability in the complex plane and contour integration.

In this chapter 7

Tick a unit to add it to your PDF basket.

  1. Complex functions
  2. Contour integration

Complex functions

3

Definition

Complex numbers and elementary functions

Adapted

Polar form and de Moivre's theorem, the exponential as a power series with its addition theorem and period 2 pi i, the trigonometric and hyperbolic functions built from it, and the logarithm, argument and complex powers as its multivalued inverses.

Prerequisites: Power series, Absolute and conditional convergence

Definition(Polar form)AdaptedAdapted from Riley, Hobson & Bence (2006), §3.3, pp. 93–94; Ahlfors (1979), Ch. 1, §2.1, pp. 13–14

A complex number z=x+iyz = x + iy with x,y∈Rx, y \in \mathbb{R} has modulus ∣z∣=x2+y2\lvert z\rvert = \sqrt{x^2 + y^2} and, if z≠0z \ne 0, an argument: any real θ\theta with

x=∣z∣cos⁡θ,y=∣z∣sin⁡θ,x = \lvert z\rvert \cos\theta, \qquad y = \lvert z\rvert \sin\theta ,

so that z=r(cos⁡θ+isin⁡θ)z = r(\cos\theta + i\sin\theta) with r=∣z∣r = \lvert z\rvert. The argument is fixed only up to an integer multiple of 2π2\pi; the whole set of its values is written arg⁡z\arg z, and the one value in a chosen half-open interval of length 2π2\pi is a principal value. We follow Riley, Hobson and Bence (Riley, Hobson & Bence (2006), Eq. (24.21), p. 834) in taking that interval to be (−π,π](-\pi, \pi] and writing the principal argument Arg⁡z\operatorname{Arg} z.

In polar form multiplication adds arguments and multiplies moduli: if z1=r1(cos⁡θ1+isin⁡θ1)z_1 = r_1(\cos\theta_1 + i \sin\theta_1) and z2=r2(cos⁡θ2+isin⁡θ2)z_2 = r_2(\cos\theta_2 + i\sin\theta_2), then the addition formulas for cos⁡\cos and sin⁡\sin give

z1z2=r1r2(cos⁡(θ1+θ2)+isin⁡(θ1+θ2)),z_1 z_2 = r_1 r_2 \bigl(\cos(\theta_1 + \theta_2) + i\sin(\theta_1 + \theta_2)\bigr),

and for z2≠0z_2 \ne 0 the quotient z1/z2z_1/z_2 has modulus r1/r2r_1/r_2 and argument θ1−θ2\theta_1 - \theta_2. Once the exponential is available this is written z=reiθz = re^{i\theta}, and the two rules become the single identity eiθ1eiθ2=ei(θ1+θ2)e^{i\theta_1}e^{i\theta_2} = e^{i(\theta_1+\theta_2)}.

Definition(The exponential)AdaptedAdapted from Ahlfors (1979), Ch. 2, §3.1, Eq. (21), p. 43

For z∈Cz \in \mathbb{C} set

ez=exp⁡z=∑n=0∞znn!.e^z = \exp z = \sum_{n=0}^{\infty} \frac{z^n}{n!} .

The series converges absolutely for every zz, because ∣z∣n/n!→0\lvert z\rvert^n/n! \to 0 faster than any geometric sequence; by the theory of power series its radius of convergence is infinite, the sum is differentiable term by term, and (ez)′=ez(e^z)' = e^z. For real zz it is the usual real exponential. Ahlfors takes the differential equation f′=ff' = f with f(0)=1f(0) = 1 as the starting point and derives the series from it; we take the series as the definition and read the differential equation off it.

Definition(Properties of the exponential)AdaptedAdapted from Ahlfors (1979), Ch. 2, §3.1 and §3.3, pp. 43, 45
  1. Addition theorem. ea+b=eaebe^{a+b} = e^a e^b for all a,b∈Ca, b \in \mathbb{C}.

  2. No zeros. ez≠0e^z \ne 0 for every zz, and e−z=1/eze^{-z} = 1/e^z.

  3. Modulus. ez‾=ezˉ\overline{e^z} = e^{\bar z}, hence ∣ex+iy∣=ex\lvert e^{x+iy}\rvert = e^x for real x,yx, y; in particular ∣eiy∣=1\lvert e^{iy}\rvert = 1.

  4. Periodicity. The periods of eze^z are exactly the numbers iωi\omega with ω\omega a period of y↦eiyy \mapsto e^{iy}, and these ω\omega are exactly the integer multiples of one smallest positive period. That smallest positive period is denoted 2π2\pi, so

    ez+2πi=ezfor all z,ez=1  ⟺  z∈2πi Z.e^{z + 2\pi i} = e^z \quad\text{for all } z, \qquad e^{z} = 1 \iff z \in 2\pi i\,\mathbb{Z} .
  5. Surjectivity onto C∖{0}\mathbb{C}\setminus\{0\}. For every w≠0w \ne 0 the equation ez=we^z = w has a solution, and its solutions form a single coset z0+2πi Zz_0 + 2\pi i \,\mathbb{Z}.

Idea of the proofOriginal

Parts 1 to 3 are algebra on the series: the Cauchy product of two absolutely convergent exponential series is, by the binomial theorem, the exponential series of the sum. Part 4 is the only part with real content, and it is where the number π\pi enters. Nothing yet tells us that eze^z has any period at all, so we must produce one: truncating the real series for cos⁡\cos gives numerical inequalities sharp enough to show that cos⁡\cos changes sign before 3\sqrt3, and the first zero y0y_0 of cos⁡\cos delivers the period 4y04y_0. Part 5 then splits ez=we^z = w into one real equation for ∣w∣\lvert w\rvert and one equation on the unit circle.

ProofAdaptedAdapted from Ahlfors (1979), Ch. 2, §3.1, p. 43 and §3.3, p. 45

1. Addition theorem. Both series converge absolutely, so their product may be summed by diagonals (absolute convergence):

eaeb=∑N=0∞∑k=0Nakk!bN−k(N−k)!=∑N=0∞1N!∑k=0N(Nk)akbN−k=∑N=0∞(a+b)NN!=ea+b,e^a e^b = \sum_{N=0}^\infty \sum_{k=0}^{N} \frac{a^k}{k!}\frac{b^{N-k}}{(N-k)!} = \sum_{N=0}^\infty \frac{1}{N!}\sum_{k=0}^N \binom{N}{k} a^k b^{N-k} = \sum_{N=0}^\infty \frac{(a+b)^N}{N!} = e^{a+b},

the middle step being the binomial theorem.

2. No zeros. Put b=−ab = -a in part 1: eae−a=e0=1e^a e^{-a} = e^0 = 1. A product equal to 11 has no zero factor.

3. Modulus. The coefficients 1/n!1/n! are real, so conjugating the series termwise gives ez‾=ezˉ\overline{e^z} = e^{\bar z}. Hence ∣eiy∣2=eiyeiy‾=eiye−iy=1\lvert e^{iy}\rvert^2 = e^{iy}\overline{e^{iy}} = e^{iy}e^{-iy} = 1 for real yy, and by part 1, ∣ex+iy∣=∣ex∣ ∣eiy∣=ex\lvert e^{x+iy}\rvert = \lvert e^x\rvert\,\lvert e^{iy}\rvert = e^x, using that ex>0e^x > 0 for real xx: for x≥0x \ge 0 every term of the series is non-negative and the first is 11, and for x<0x < 0 we have ex=1/e−x>0e^x = 1/e^{-x} > 0 by part 2.

4. Periodicity. A period exists. Write C(y)=12(eiy+e−iy)C(y) = \tfrac12(e^{iy} + e^{-iy}) and S(y)=12i(eiy−e−iy)S(y) = \tfrac{1}{2i}(e^{iy} - e^{-iy}) for real yy; from the series these are the real functions cos⁡\cos and sin⁡\sin, with C′=−SC' = -S, S′=CS' = C, C(0)=1C(0)=1, S(0)=0S(0)=0, and C2+S2=1C^2 + S^2 = 1 because C2+S2=eiye−iy=1C^2 + S^2 = e^{iy}e^{-iy} = 1 by parts 1 and 2. In particular ∣C∣≤1\lvert C\rvert \le 1 and ∣S∣≤1\lvert S \rvert \le 1. Four successive integrations from 00 now trap CC and SS between consecutive partial sums of their series. From C≤1C \le 1 and S(0)=0S(0)=0 we get S(y)≤yS(y) \le y for y≥0y \ge 0; integrating C′=−S≥−yC' = -S \ge -y gives C(y)≥1−y2/2C(y) \ge 1 - y^2/2; integrating S′=C≥1−y2/2S' = C \ge 1 - y^2/2 gives S(y)≥y−y3/6S(y) \ge y - y^3/6; integrating C′=−S≤−y+y3/6C' = -S \le -y + y^3/6 gives

C(y)≤1−y22+y424,y≥0.C(y) \le 1 - \frac{y^2}{2} + \frac{y^4}{24}, \qquad y \ge 0 .

At y=3y = \sqrt3 the right-hand side is 1−32+38=−18<01 - \tfrac32 + \tfrac38 = -\tfrac18 < 0. Since CC is continuous with C(0)=1C(0) = 1, it has a zero in (0,3)(0, \sqrt 3); let y0y_0 be the smallest one. From C(y0)=0C(y_0) = 0 and C2+S2=1C^2 + S^2 = 1 we get S(y0)2=1S(y_0)^2 = 1; and S(y0)>0S(y_0) > 0, because S(y)≥y−y3/6>0S(y) \ge y - y^3/6 > 0 for 0<y<60 < y < \sqrt 6 and y0<3<6y_0 < \sqrt 3 < \sqrt 6. Hence S(y0)=1S(y_0) = 1, so

eiy0=C(y0)+iS(y0)=i,e4iy0=i4=1:e^{iy_0} = C(y_0) + iS(y_0) = i, \qquad e^{4iy_0} = i^4 = 1 :

the number 4y04y_0 is a period of y↦eiyy \mapsto e^{iy}.

It is the smallest. For 0<y<y00 < y < y_0 we have C(y)>0C(y) > 0 (by minimality of y0y_0 and continuity, since C(0)=1C(0)=1), so S′=C>0S' = C > 0 there and SS increases strictly from S(0)=0S(0) = 0. Hence 0<S(y)<S(y0)=10 < S(y) < S(y_0) = 1 on (0,y0)(0,y_0), and therefore C(y)2=1−S(y)2∈(0,1)C(y)^2 = 1 - S(y)^2 \in (0,1). So for 0<y<y00 < y < y_0 the number eiy=C(y)+iS(y)e^{iy} = C(y) + iS(y) has both parts nonzero and is not one of 1,−1,i,−i1, -1, i, -i; consequently e4iy≠1e^{4iy} \ne 1, because w4=1w^4 = 1 forces w∈{1,−1,i,−i}w \in \{1,-1,i,-i\}. No positive period is therefore smaller than 4y04y_0.

All periods are multiples. Let ω>0\omega > 0 be a period and ω0=4y0\omega_0 = 4y_0. Choose the integer nn with nω0≤ω<(n+1)ω0n\omega_0 \le \omega < (n+1)\omega_0. If ω≠nω0\omega \ne n\omega_0 then ω−nω0\omega - n\omega_0 is a positive period smaller than ω0\omega_0, which we have excluded. Finally, ez+c=eze^{z+c} = e^z for all zz forces ec=1e^c = 1, so ∣ec∣=1\lvert e^c \rvert = 1 gives Re⁡c=0\operatorname{Re} c = 0 by part 3, and then c=iωc = i\omega with ω\omega a period of y↦eiyy \mapsto e^{iy}. Writing 2π=4y02\pi = 4y_0 gives the displayed statement.

5. Surjectivity. Let w≠0w \ne 0 and seek z=x+iyz = x+iy. By part 3, ex+iy=we^{x+iy} = w splits into ex=∣w∣e^x = \lvert w\rvert and eiy=w/∣w∣e^{iy} = w/\lvert w\rvert. The first has the unique real solution x=log⁡∣w∣x = \log\lvert w\rvert, since the real exponential is a strictly increasing bijection of R\mathbb{R} onto (0,∞)(0,\infty). For the second, the map y↦eiyy \mapsto e^{iy} sends [0,2π)[0, 2\pi) onto the unit circle: it is continuous; part 4 showed that S>0S > 0 on (0,y0)(0,y_0), so C′=−S<0C' = -S < 0 there and CC decreases strictly from C(0)=1C(0) = 1 to C(y0)=0C(y_0) = 0, which together with C2+S2=1C^2 + S^2 = 1 and S≥0S \ge 0 traces the first quadrant of the circle exactly once, and the four reflections ei(y+kπ/2)=ikeiye^{i(y + k\pi/2)} = i^k e^{iy} carry this arc onto the remaining three quadrants. Two solutions of eiy=w/∣w∣e^{iy} = w/\lvert w\rvert differ by a period, by part 4, so the solution set of ez=we^z = w is a single coset of 2πi Z2\pi i\,\mathbb{Z}. □\square

Definition(Trigonometric and hyperbolic functions)AdaptedAdapted from Ahlfors (1979), Ch. 2, §3.2, Eq. (23), p. 43; Riley, Hobson & Bence (2006), §3.7.1–3.7.2, pp. 102–103

For z∈Cz \in \mathbb{C},

cos⁡z=eiz+e−iz2,sin⁡z=eiz−e−iz2i,cosh⁡z=ez+e−z2,sinh⁡z=ez−e−z2,\cos z = \frac{e^{iz} + e^{-iz}}{2}, \quad \sin z = \frac{e^{iz} - e^{-iz}}{2i}, \quad \cosh z = \frac{e^{z} + e^{-z}}{2}, \quad \sinh z = \frac{e^{z} - e^{-z}}{2} ,

with tan⁡z=sin⁡z/cos⁡z\tan z = \sin z/\cos z and tanh⁡z=sinh⁡z/cosh⁡z\tanh z = \sinh z/\cosh z where the denominators do not vanish. Substituting the exponential series gives the familiar power series, so on R\mathbb{R} these are the usual functions; and (cos⁡z)′=−sin⁡z(\cos z)' = -\sin z, (sin⁡z)′=cos⁡z(\sin z)' = \cos z. Euler's formula eiz=cos⁡z+isin⁡ze^{iz} = \cos z + i \sin z and the identity cos⁡2z+sin⁡2z=1\cos^2 z + \sin^2 z = 1 are immediate from the definitions, as are

cosh⁡z=cos⁡iz,sinh⁡z=−isin⁡iz,\cosh z = \cos iz, \qquad \sinh z = -i \sin iz ,

so the hyperbolic and circular functions are each other's restrictions to perpendicular lines. Separating real and imaginary parts, sin⁡(x+iy)=sin⁡xcosh⁡y+icos⁡xsinh⁡y\sin(x+iy) = \sin x\cosh y + i\cos x \sinh y, whence ∣sin⁡(x+iy)∣2=sin⁡2x+sinh⁡2y\lvert \sin(x+iy)\rvert^2 = \sin^2 x + \sinh^2 y: unlike their real restrictions, sin⁡\sin and cos⁡\cos are unbounded on C\mathbb{C}.

Definition(Logarithm, argument and powers)AdaptedAdapted from Ahlfors (1979), Ch. 2, §3.4, p. 46; Riley, Hobson & Bence (2006), §24.4, Eqs. (24.20)–(24.22), p. 834

For w≠0w \ne 0, a logarithm of ww is any solution of ez=we^z = w. By part 5 of the statement above there are infinitely many, differing by multiples of 2πi2\pi i; their imaginary parts are exactly the values of arg⁡w\arg w, and the whole set is written

Ln⁡w=log⁡∣w∣+iarg⁡w=log⁡∣w∣+i(θ+2kπ),k∈Z,\operatorname{Ln} w = \log\lvert w\rvert + i\arg w = \log \lvert w \rvert + i(\theta + 2k\pi), \quad k \in \mathbb{Z},

with θ\theta any fixed argument of ww and log⁡\log the real logarithm of the positive number ∣w∣\lvert w\rvert. The number 00 has no logarithm. The principal value ln⁡w=log⁡∣w∣+iArg⁡w\ln w = \log\lvert w\rvert + i\operatorname{Arg} w takes θ∈(−π,π]\theta \in (-\pi, \pi].

For a≠0a \ne 0 and b∈Cb \in \mathbb{C} the power is defined by

ab=exp⁡(bLn⁡a),a^{b} = \exp\bigl(b \operatorname{Ln} a\bigr),

again a set of values, which differ by the factors e2πikbe^{2\pi i k b}. It has a single value exactly when b∈Zb \in \mathbb{Z}, and then agrees with the repeated product; it has qq values when b=p/qb = p/q in lowest terms with q∈Nq \in \mathbb{N}; otherwise it has infinitely many. Its principal value uses ln⁡a\ln a.

Definition(de Moivre and the nth roots)AdaptedAdapted from Riley, Hobson & Bence (2006), §3.4, Eq. (3.27), p. 95 and §3.4.2, p. 97; Ahlfors (1979), Ch. 1, §2.2, Eqs. (19)–(20), pp. 15–16

For θ∈R\theta \in \mathbb{R} and n∈Zn \in \mathbb{Z},

(cos⁡θ+isin⁡θ)n=cos⁡nθ+isin⁡nθ(de Moivre’s theorem).(\cos\theta + i \sin\theta)^n = \cos n\theta + i \sin n\theta \qquad\text{(de Moivre's theorem).}

For a≠0a \ne 0 with a=ρ(cos⁡φ+isin⁡φ)a = \rho(\cos\varphi + i\sin\varphi) and n∈Nn \in \mathbb{N}, the equation zn=az^n = a has exactly nn solutions,

zk=ρ1/n(cos⁡φ+2kπn+isin⁡φ+2kπn),k=0,1,…,n−1,z_k = \rho^{1/n}\Bigl(\cos\frac{\varphi + 2k\pi}{n} + i \sin\frac{\varphi + 2k\pi}{n}\Bigr), \qquad k = 0, 1, \dots, n-1 ,

where ρ1/n\rho^{1/n} is the positive real nnth root. They have equal modulus and equally spaced arguments, so they are the vertices of a regular nn-gon centred at 00. Taking a=1a = 1 gives the nnth roots of unity 1,ω,…,ωn−11, \omega, \dots, \omega^{n-1} with ω=e2πi/n\omega = e^{2\pi i/n}, and every solution of zn=az^n = a is ωkz0\omega^k z_0 for one kk.

ProofAdaptedAdapted from Ahlfors (1979), Ch. 1, §2.2, pp. 15–16; Riley, Hobson & Bence (2006), §3.4.2, p. 97

De Moivre's theorem is the addition theorem in polar clothing: cos⁡θ+isin⁡θ=eiθ\cos\theta + i\sin\theta = e^{i\theta} by Euler's formula, and (eiθ)n=einθ(e^{i\theta})^n = e^{in\theta} for n≥0n \ge 0 by part 1 of the statement above and induction, while for n<0n < 0 one uses e−iθ=1/eiθe^{-i\theta} = 1/e^{i\theta} from part 2.

For the roots, write z=reiϑz = r e^{i\vartheta} with r>0r > 0; this costs nothing, because zn=a≠0z^n = a \ne 0 forces z≠0z \ne 0, and every nonzero zz has such a form by the definition of polar form. Then zn=rneinϑz^n = r^n e^{in\vartheta}, and zn=a=ρeiφz^n = a = \rho e^{i\varphi} holds if and only if the moduli agree and the arguments agree modulo a period:

  • rn=ρr^n = \rho, which by strict monotonicity of r↦rnr \mapsto r^n on (0,∞)(0,\infty) fixes r=ρ1/nr = \rho^{1/n} uniquely;
  • einϑ=eiφe^{in\vartheta} = e^{i\varphi}, which by part 4 means nϑ−φ∈2πZn\vartheta - \varphi \in 2\pi\mathbb{Z}, that is ϑ=(φ+2kπ)/n\vartheta = (\varphi + 2k\pi)/n for some k∈Zk \in \mathbb{Z}.

Two integers kk and k′k' give the same zz exactly when (k−k′)/n∈Z(k - k')/n \in \mathbb{Z}, so k=0,1,…,n−1k = 0, 1, \dots, n-1 list the solutions once each: there are exactly nn. Their moduli are all ρ1/n\rho^{1/n} and consecutive arguments differ by 2π/n2\pi/n, which is the description of a regular nn-gon. Dividing two solutions gives a solution of zn=1z^n = 1, which proves the last sentence. □\square

RemarkOriginal

Where the hypotheses bite, and which convention is which.

  • 00 is excluded everywhere. eze^z never vanishes, so 00 has no logarithm and 0b0^b is not defined by the formula above; the nnth-root count also needs a≠0a \ne 0, since zn=0z^n = 0 has the single root 00 with multiplicity nn.
  • Principal ranges differ between books. Riley, Hobson and Bence (Riley, Hobson & Bence (2006), §3.3, p. 93 and Eq. (24.21), p. 834) use −π<θ≤π-\pi < \theta \le \pi, and we follow them throughout. Ahlfors (Ahlfors (1979), Ch. 2, §3.4, p. 46) normalises arg⁡\arg to 0≤y<2π0 \le y < 2\pi when solving eiy=w/∣w∣e^{iy} = w/\lvert w\rvert. Arfken, Weber and Harris (Arfken, Weber & Harris (2013), §1.8, Eq. (1.138), p. 60 and §11.6, Eq. (11.55), p. 501) fix no range at all: Eq. (1.138) carries the free integer nn, and §11.6 says only that the branch n=0n = 0 is the one usually called principal — their own Example 11.6.1 (p. 499) runs θ\theta from 00 to 2π2\pi and explicitly declines the value −π-\pi. All three readings give ln⁡(−1)=iπ\ln(-1) = i\pi, but just below the negative real axis the first two differ by 2πi2\pi i: at z=−1−iεz = -1 - i\varepsilon the first gives an argument near −π-\pi and the second one near +π+\pi.
  • De Moivre's theorem, for which exponents. Riley, Hobson and Bence (Riley, Hobson & Bence (2006), §3.4, Eq. (3.27), p. 95) assert the identity for every exponent, real or complex alike. Our statement restricts it to n∈Zn \in \mathbb{Z}, which is the form Ahlfors (Ahlfors (1979), Ch. 1, §2.2, Eq. (20), p. 15) gives, because for non-integer nn the left-hand side is a set of values and the right-hand side a single one, so the equation has to be read as "some value of the left side equals the right side". At n=12n = \tfrac12 and θ=3π\theta = 3\pi, for instance, the left side is (cos⁡3π+isin⁡3π)1/2=(−1)1/2=±i(\cos 3\pi + i \sin 3\pi)^{1/2} = (-1)^{1/2} = \pm i with principal value +i+i, while the right side is cos⁡3π2+isin⁡3π2=−i\cos\tfrac{3\pi}{2} + i\sin\tfrac{3\pi}{2} = -i: the two agree only after the other branch is chosen.
  • The addition law for Ln⁡\operatorname{Ln} is a statement about sets. Ln⁡(z1z2)=Ln⁡z1+Ln⁡z2\operatorname{Ln}(z_1z_2) = \operatorname{Ln} z_1 + \operatorname{Ln} z_2 holds as an equality of sets of complex numbers, but a chosen value on one side need only agree with one on the other up to 2πi2\pi i. For principal values it can fail outright: ln⁡(−1)+ln⁡(−1)=2πi\ln(-1) + \ln(-1) = 2\pi i, while ln⁡((−1)(−1))=ln⁡1=0\ln\bigl((-1)(-1)\bigr) = \ln 1 = 0.
  • Powers are sets too, so the school rules need care: (ab)c=abc(a^b)^c = a^{bc} is false in general. The multivaluedness is not a defect to be suppressed but the beginning of branch points and branch cuts.
Example(Values, roots and a real-valued power)Original
  1. A logarithm. For z=−1+i3z = -1 + i\sqrt 3 we have ∣z∣=2\lvert z\rvert = 2 and Arg⁡z=2π/3\operatorname{Arg} z = 2\pi/3 (the point lies in the second quadrant), so ln⁡z=log⁡2+2πi3≈0.693147+2.094395 i\ln z = \log 2 + \tfrac{2\pi i}{3} \approx 0.693147 + 2.094395\,i, and Ln⁡z=log⁡2+i(2π/3+2kπ)\operatorname{Ln} z = \log 2 + i(2\pi/3 + 2k\pi).

  2. A power all of whose values are real. Take the base a=12(1+i)=eiπ/4a = \tfrac{1}{\sqrt2}(1+i) = e^{i\pi/4}, which has ∣a∣=1\lvert a\rvert = 1, and the purely imaginary exponent b=4ib = 4i. Then log⁡∣a∣=0\log\lvert a\rvert = 0, so Ln⁡a=i(π/4+2kπ)\operatorname{Ln} a = i(\pi/4 + 2k\pi) and

    ab=exp⁡(4i⋅i(π/4+2kπ))=e−π−8kπ,k∈Z,a^{b} = \exp\bigl(4i\cdot i(\pi/4 + 2k\pi)\bigr) = e^{-\pi - 8k\pi}, \qquad k \in \mathbb{Z},

    a set of positive real numbers: the principal value is e−π≈0.0432139e^{-\pi} \approx 0.0432139, and the next value up is e7π≈3.5533×109e^{7\pi}\approx 3.5533 \times 10^{9}. The mechanism is visible in general. For a purely imaginary exponent b=itb = it with tt real, bLn⁡a=itlog⁡∣a∣−t(θ+2kπ)b\operatorname{Ln}a = it\log\lvert a\rvert - t(\theta + 2k\pi), so the values of aba^{b} are e−t(θ+2kπ) eitlog⁡∣a∣e^{-t(\theta + 2k\pi)}\,e^{it\log\lvert a\rvert}: all of them are real exactly when tlog⁡∣a∣∈πZt\log\lvert a\rvert \in \pi\mathbb{Z}, and all positive exactly when tlog⁡∣a∣∈2πZt\log\lvert a\rvert \in 2\pi\mathbb{Z}. Our base has ∣a∣=1\lvert a \rvert = 1, which is the easiest way to arrange it. Riley, Hobson and Bence (Riley, Hobson & Bence (2006), §3.5, p. 100) make the same observation with i−2ii^{-2i}.

  3. A power that is not. (1+i)i(1+i)^{i} has principal value exp⁡(i(12log⁡2+iπ/4))=e−π/4(cos⁡(12log⁡2)+isin⁡(12log⁡2))≈0.428829+0.154872 i\exp\bigl(i(\tfrac12\log 2 + i\pi/4)\bigr) = e^{-\pi/4}\bigl(\cos(\tfrac12 \log 2) + i \sin(\tfrac12\log 2)\bigr) \approx 0.428829 + 0.154872\,i, since ln⁡(1+i)=12log⁡2+iπ/4≈0.346574+0.785398 i\ln(1+i) = \tfrac12\log 2 + i\pi/4 \approx 0.346574 + 0.785398\,i.

  4. Cube roots. For a=8ia = 8i we have ρ=8\rho = 8, φ=π/2\varphi = \pi/2, so the roots are 2eiπ/6,2ei5π/6,2ei3π/22e^{i\pi/6}, 2e^{i5\pi/6}, 2e^{i3\pi/2}, that is 3+i\sqrt3 + i, −3+i-\sqrt3 + i and −2i-2i. Check the first: (3+i)3=(3+i)(3+23i−1)=(3+i)(2+23i)=23+6i+2i−23=8i(\sqrt3+i)^3 = (\sqrt3+i)(3 + 2\sqrt3 i - 1) = (\sqrt3+i)(2 + 2\sqrt3 i) = 2\sqrt3 + 6i + 2i - 2\sqrt3 = 8i.

  5. A sine off the real axis. sin⁡(1+2i)=sin⁡1cosh⁡2+icos⁡1sinh⁡2≈3.165779+1.959601 i\sin(1 + 2i) = \sin 1\cosh 2 + i \cos 1 \sinh 2 \approx 3.165779 + 1.959601\,i, with ∣sin⁡(1+2i)∣2=sin⁡21+sinh⁡22≈13.86219\lvert \sin(1+2i)\rvert^2 = \sin^2 1 + \sinh^2 2 \approx 13.86219 — far outside the real range [−1,1][-1,1].

Example(A real integral through the complex exponential)AdaptedAdapted from Riley, Hobson & Bence (2006), §3.6, Eq. (3.37), p. 101

The route is Riley, Hobson and Bence (Riley, Hobson & Bence (2006), §3.6, Eq. (3.37), p. 101)'s; the function and the numbers are ours. We want antiderivatives of e−2xcos⁡5xe^{-2x}\cos 5x and e−2xsin⁡5xe^{-2x}\sin 5x, and then the value of ∫0∞e−2xcos⁡5x dx\int_0^\infty e^{-2x}\cos 5x\,dx.

  1. One complex integrand instead of two real ones. Both functions are components of e(−2+5i)x=e−2x(cos⁡5x+isin⁡5x)e^{(-2+5i)x} = e^{-2x}(\cos 5x + i \sin 5x), whose antiderivative is e(−2+5i)x/(−2+5i)e^{(-2+5i)x}/(-2+5i) because ddxeλx=λeλx\frac{d}{dx}e^{\lambda x} = \lambda e^{\lambda x} for complex λ\lambda (differentiate the exponential series term by term). Rationalising with (−2+5i)‾=−2−5i\overline{(-2+5i)} = -2-5i and ∣−2+5i∣2=29\lvert -2+5i\rvert^2 = 29,

    ∫e(−2+5i)x dx=(−2−5i)29 e−2x(cos⁡5x+isin⁡5x)+c.\int e^{(-2+5i)x}\,dx = \frac{(-2-5i)}{29}\,e^{-2x}\bigl(\cos 5x + i \sin 5x\bigr) + c .
  2. Separate the parts. The real part of the right-hand side is 129e−2x(−2cos⁡5x+5sin⁡5x)\tfrac{1}{29}e^{-2x}(-2\cos 5x + 5 \sin 5x) and the imaginary part is 129e−2x(−2sin⁡5x−5cos⁡5x)\tfrac{1}{29}e^{-2x}(-2\sin 5x - 5\cos 5x), so

    ∫e−2xcos⁡5x dx=e−2x29(5sin⁡5x−2cos⁡5x)+c1,∫e−2xsin⁡5x dx=−e−2x29(5cos⁡5x+2sin⁡5x)+c2.\int e^{-2x}\cos 5x\,dx = \frac{e^{-2x}}{29}\bigl(5\sin 5x - 2\cos 5x\bigr) + c_1, \qquad \int e^{-2x}\sin 5x\,dx = -\frac{e^{-2x}}{29}\bigl(5\cos 5x + 2\sin 5x\bigr) + c_2 .

    Two antiderivatives for the price of one.

  3. Check by differentiating back. For the first, 129e−2x(−2(5sin⁡5x−2cos⁡5x)+(25cos⁡5x+10sin⁡5x))=129e−2x⋅29cos⁡5x=e−2xcos⁡5x\tfrac{1}{29}e^{-2x}\bigl(-2(5\sin 5x - 2\cos 5x) + (25\cos 5x + 10 \sin 5x)\bigr) = \tfrac{1}{29}e^{-2x}\cdot 29\cos 5x = e^{-2x}\cos 5x.

  4. A number to check against. The improper integral ∫0∞e−2xcos⁡5x dx\int_0^\infty e^{-2x}\cos 5x\,dx converges absolutely, since ∣e(−2+5i)x∣=e−2x\lvert e^{(-2+5i)x}\rvert = e^{-2x} by part 3 of the statement above. The antiderivative tends to 00 as x→∞x \to \infty and equals −2/29-2/29 at x=0x = 0, so the integral is 2/29≈0.06896552/29 \approx 0.0689655; likewise ∫0∞e−2xsin⁡5x dx=5/29≈0.1724138\int_0^\infty e^{-2x}\sin 5x\,dx = 5/29 \approx 0.1724138. Both agree with numerical quadrature to the digits shown. The pair is just ∫0∞e(−2+5i)xdx=1/(2−5i)=(2+5i)/29\int_0^\infty e^{(-2+5i)x}dx = 1/(2-5i) = (2+5i)/29 split into its parts.

ApplicationOriginal

Every linear system driven at a single frequency is handled with these functions. A real signal Acos⁡(ωt+ϕ)A\cos(\omega t + \phi) is written Re⁡(A^eiωt)\operatorname{Re}\bigl(\hat A e^{i\omega t}\bigr) with the complex amplitude A^=Aeiϕ\hat A = Ae^{i\phi}, so that amplitude and phase travel together in one number and differentiation becomes multiplication by iωi\omega. Impedances in AC circuit theory, the complex refractive index n+iκn + i\kappa of an absorbing medium (where ∣eikz∣=e−κωz/c\lvert e^{ikz}\rvert = e^{-\kappa \omega z/c} is the Beer-Lambert decay), and the e−iEt/ℏe^{-iEt/\hbar} phase of a stationary quantum state are all the same device.

The nnth roots of unity are the characters of the cyclic group Z/n\mathbb{Z}/n and so the entries of the discrete Fourier matrix (the DFT and FFT); in crystallography they are the allowed phase factors of a rotation axis, and the fact that the roots of zn=az^n = a sit at the vertices of a regular polygon is why normal modes of a ring of nn identical masses come in such symmetric families (normal modes). The unboundedness of sin⁡\sin and cos⁡\cos off the real axis, visible in ∣sin⁡(x+iy)∣2=sin⁡2x+sinh⁡2y\lvert \sin(x+iy)\rvert^2 = \sin^2 x + \sinh^2 y, is what makes evanescent waves possible: a "wavenumber" forced to be imaginary turns an oscillation into a growing or decaying exponential.

Theorem

Analyticity and the Cauchy-Riemann equations

Adapted

A complex function is differentiable at a point only if its real and imaginary parts satisfy the Cauchy-Riemann equations there, and it is analytic on an open set if those equations hold with continuous first partial derivatives; the two parts are then harmonic conjugates.

Prerequisites: Complex numbers and elementary functions, Partial derivatives and the chain rule

Definition(Complex differentiability and analyticity)AdaptedAdapted from Ahlfors (1979), Ch. 2, §1.1–1.2, pp. 23–24; Arfken, Weber & Harris (2013), §11.2, Eq. (11.3) and the subsection on analytic functions, pp. 471–472; Riley, Hobson & Bence (2006), §24.1, Eq. (24.1), p. 825, and p. 826

Let ff be defined on a neighbourhood of z∈Cz \in \mathbb{C}. Then ff is complex differentiable at zz if

f′(z)=lim⁡h→0f(z+h)−f(z)hf'(z) = \lim_{h \to 0} \frac{f(z+h) - f(z)}{h}

exists, the limit being taken over complex h→0h \to 0 and required to be the same however hh approaches 00. A function that is single-valued and complex differentiable at every point of an open set Ω\Omega is analytic (equivalently holomorphic) on Ω\Omega; a function analytic on all of C\mathbb{C} is entire. A point where ff is defined on a neighbourhood but not differentiable is a singular point.

Writing f=u+ivf = u + iv with u,vu, v real-valued functions of (x,y)(x,y), we use subscripts for partial derivatives: ux=∂u/∂xu_x = \partial u/\partial x and so on.

Differentiability at a point implies continuity there, since f(z+h)−f(z)=h⋅f(z+h)−f(z)h→0⋅f′(z)=0f(z+h) - f(z) = h\cdot\frac{f(z+h)-f(z)}{h} \to 0 \cdot f'(z) = 0. The sum, product and (where the denominator does not vanish) quotient of differentiable functions are differentiable, with the usual rules; the proofs are the real ones verbatim, because they use only field algebra and the limit laws.

Theorem(Cauchy-Riemann equations)AdaptedAdapted from Ahlfors (1979), Ch. 2, §1.2, Eqs. (5)–(6) and the displayed theorem, pp. 25–26; Arfken, Weber & Harris (2013), §11.2, Eqs. (11.9)–(11.12), pp. 471–472

Let f=u+ivf = u + iv be defined on an open set Ω⊆C\Omega \subseteq \mathbb{C}.

  1. Necessity. If ff is complex differentiable at z0=x0+iy0z_0 = x_0 + iy_0, then all four partial derivatives ux,uy,vx,vyu_x, u_y, v_x, v_y exist at z0z_0 and satisfy there the Cauchy-Riemann equations

    ux=vy,uy=−vx,u_x = v_y, \qquad u_y = -v_x ,

    equivalently ∂f/∂x=−i ∂f/∂y\partial f/\partial x = -i\,\partial f/\partial y; and then f′(z0)=ux+ivxf'(z_0) = u_x + i v_x.

  2. Sufficiency. If uu and vv have continuous first partial derivatives on Ω\Omega and satisfy the Cauchy-Riemann equations at every point of Ω\Omega, then ff is analytic on Ω\Omega, and f′=ux+ivxf'= u_x + iv_x is continuous. (Riley, Hobson and Bence (Riley, Hobson & Bence (2006), p. 828) record the same three conditions — existence, continuity, and the relations — as sufficient, and note that it is continuity that separates them from the necessary conditions.)

  3. Harmonic conjugates. If in addition uu and vv have continuous second partial derivatives, then both are harmonic on Ω\Omega:

    uxx+uyy=0,vxx+vyy=0.u_{xx} + u_{yy} = 0, \qquad v_{xx} + v_{yy} = 0 .

    A pair (u,v)(u,v) of harmonic functions satisfying the Cauchy-Riemann equations is a pair of conjugate harmonic functions; vv is determined by uu up to an additive constant on a connected Ω\Omega, and uu is the conjugate of −v-v. Connectedness is needed for that uniqueness clause alone: parts 1 and 2, and the harmonicity of uu and vv, hold on any open set.

    The hypothesis of continuous second partials is ours, not the source's, and it is there only to let us interchange the mixed partials of vv in the proof. It can be removed: Ahlfors (Ahlfors (1979), Ch. 2, §1.2, p. 25) obtains harmonicity instead from the analyticity of f′f', which Cauchy's integral formula supplies for every analytic ff with no extra assumption, so an analytic function automatically has partial derivatives of all orders. We assume the hypothesis here because that route is not available until later in the chapter.

Idea of the proofOriginal

Necessity costs nothing: the limit must not depend on the direction of approach, so computing it along the two coordinate directions gives two expressions for the same number f′(z0)f'(z_0), and equating real and imaginary parts gives the two equations. Sufficiency is where work is needed, and the work is all in one step: the Cauchy-Riemann equations are algebra, but to turn the existence of partial derivatives into a statement about every direction of approach we need the first-order Taylor expansion of a function of two variables, which is exactly what continuity of the partials buys. Once that expansion is available, the Cauchy-Riemann equations say precisely that the real-linear map (h,k)↦(uxh+uyk,vxh+vyk)(h,k) \mapsto (u_xh + u_yk, v_xh + v_yk) is multiplication by the complex number ux+ivxu_x + iv_x, and a difference quotient whose numerator is multiplication by a fixed complex number has a limit. Harmonicity is then one more differentiation of each equation.

ProofAdaptedAdapted from Ahlfors (1979), Ch. 2, §1.2, pp. 24–26; Arfken, Weber & Harris (2013), §11.2, Eqs. (11.7)–(11.13), pp. 471–473

Part 1 (necessity). Suppose f′(z0)f'(z_0) exists. Restrict hh to real values: then yy is held at y0y_0 and the difference quotient becomes the difference quotient of x↦f(x+iy0)x \mapsto f(x + iy_0). Its limit exists and equals f′(z0)f'(z_0), so uxu_x and vxv_x exist at z0z_0 and

f′(z0)=ux+ivx.f'(z_0) = u_x + i v_x .

Now restrict hh to purely imaginary values h=ikh = ik with kk real. Since 1/(ik)=−i/k1/(ik) = -i/k,

f′(z0)=lim⁡k→0f(z0+ik)−f(z0)ik=−ilim⁡k→0f(z0+ik)−f(z0)k=−i(uy+ivy)=vy−iuy,f'(z_0) = \lim_{k\to 0} \frac{f(z_0 + ik) - f(z_0)}{ik} = -i \lim_{k \to 0}\frac{f(z_0 + ik) - f(z_0)}{k} = -i(u_y + iv_y) = v_y - iu_y ,

so uyu_y and vyv_y exist at z0z_0 as well. The two expressions for f′(z0)f'(z_0) are the same complex number; equating real parts gives ux=vyu_x = v_y and imaginary parts gives vx=−uyv_x = -u_y. (The combined form is fx=−ifyf_x = -if_y, which on substituting f=u+ivf = u+iv reproduces the pair.)

Part 2 (sufficiency). Fix z0∈Ωz_0 \in \Omega and let h+ikh + ik be a small complex increment, h,kh, k real. Because uu has continuous first partials on the open set Ω\Omega, it is differentiable there in the sense of the calculus of several variables:

u(x0+h,y0+k)−u(x0,y0)=uxh+uyk+ε1(h,k),u(x_0 + h, y_0 + k) - u(x_0,y_0) = u_x h + u_y k + \varepsilon_1(h,k),

with the partials evaluated at z0z_0 and ε1(h,k)/∣h+ik∣→0\varepsilon_1(h,k)/\lvert h + ik\rvert \to 0 as h+ik→0h + ik \to 0; and likewise

v(x0+h,y0+k)−v(x0,y0)=vxh+vyk+ε2(h,k).v(x_0 + h, y_0 + k) - v(x_0,y_0) = v_x h + v_y k + \varepsilon_2(h,k) .

This is the one place where continuity of the partials does essential work; it is used once more, trivially, at the end of the part, to see that f′=ux+ivxf' = u_x + iv_x is continuous. Adding ii times the second to the first,

f(z0+h+ik)−f(z0)=(ux+ivx)h+(uy+ivy)k+ε,ε=ε1+iε2.f(z_0 + h + ik) - f(z_0) = (u_x + iv_x)h + (u_y + iv_y)k + \varepsilon, \qquad \varepsilon = \varepsilon_1 + i\varepsilon_2 .

Now apply the hypothesis: by the Cauchy-Riemann equations at z0z_0, uy+ivy=−vx+iux=i(ux+ivx)u_y + iv_y = -v_x + iu_x = i(u_x + i v_x). Hence the two linear terms combine,

f(z0+h+ik)−f(z0)=(ux+ivx)(h+ik)+ε,f(z_0 + h + ik) - f(z_0) = (u_x + iv_x)(h + ik) + \varepsilon ,

and dividing by h+ik≠0h + ik \ne 0,

f(z0+h+ik)−f(z0)h+ik=(ux+ivx)+εh+ik.\frac{f(z_0 + h + ik) - f(z_0)}{h+ik} = (u_x + iv_x) + \frac{\varepsilon}{h+ik} .

The last term tends to 00 as h+ik→0h + ik \to 0, because ∣ε∣/∣h+ik∣≤∣ε1∣/∣h+ik∣+∣ε2∣/∣h+ik∣→0\lvert \varepsilon\rvert/\lvert h+ik\rvert \le \lvert\varepsilon_1\rvert/\lvert h+ik\rvert + \lvert \varepsilon_2\rvert/\lvert h+ik\rvert \to 0. So the limit exists, independently of the direction of approach, and f′(z0)=ux+ivxf'(z_0) = u_x + iv_x. As z0∈Ωz_0 \in \Omega was arbitrary, ff is analytic on Ω\Omega, and f′=ux+ivxf' = u_x + iv_x is continuous because uxu_x and vxv_x are.

Part 3 (harmonicity). Differentiate ux=vyu_x = v_y with respect to xx and uy=−vxu_y = -v_x with respect to yy:

uxx=vyx,uyy=−vxy.u_{xx} = v_{yx}, \qquad u_{yy} = -v_{xy} .

The second partial derivatives of vv are continuous by hypothesis, so the mixed partials agree, vyx=vxyv_{yx} = v_{xy}. Adding the two displayed equations therefore gives uxx+uyy=0u_{xx} + u_{yy} = 0. For vv, note that −if=v−iu-if = v - iu is differentiable wherever ff is, with real part vv and imaginary part −u-u; applying the result just proved to −if-if gives vxx+vyy=0v_{xx} + v_{yy} = 0. (Equivalently, differentiate ux=vyu_x = v_y with respect to yy and uy=−vxu_y = -v_x with respect to xx and subtract.)

Uniqueness of the conjugate. If vv and v~\tilde v both satisfy the Cauchy-Riemann equations with the same uu, then w=v−v~w = v - \tilde v has wx=wy=0w_x = w_y = 0 on Ω\Omega, so ww is constant on each connected component: constant on Ω\Omega if Ω\Omega is connected. □\square

RemarkOriginal

Where the hypotheses are used, and what fails without them.

  • The Cauchy-Riemann equations at a single point are not enough. Let

    f(z)=zˉ 2z (z≠0),f(0)=0.f(z) = \frac{\bar z^{\,2}}{z}\ (z \ne 0), \qquad f(0) = 0 .

    In polar form f(reiθ)=re−3iθf(re^{i\theta}) = r e^{-3i\theta}, so ∣f(z)∣=∣z∣\lvert f(z)\rvert = \lvert z\rvert and ff is continuous at 00. On the real axis f(x)=xf(x) = x, giving ux(0)=1u_x(0) = 1, vx(0)=0v_x(0) = 0; on the imaginary axis f(iy)=iyf(iy) = iy, giving uy(0)=0u_y(0) = 0, vy(0)=1v_y(0) = 1. Both Cauchy-Riemann equations hold at the origin. Yet f(z)/z=e−4iθf(z)/z = e^{-4i\theta}, which is 11 along θ=0\theta = 0 and −1-1 along θ=π/4\theta = \pi/4: no limit, so ff is not differentiable at 00. The hypothesis in part 2 is that the equations hold throughout an open set, and here they hold at the one point only. Off the origin, comparing the radial and tangential difference quotients at z0=r0eiθ0z_0 = r_0e^{i\theta_0} shows the equations fail: moving radially, dz=eiθ0drdz = e^{i\theta_0}dr and df=e−3iθ0drdf = e^{-3i\theta_0}dr, so the quotient is e−4iθ0e^{-4i\theta_0}, whereas moving along the circle ∣z∣=r0\lvert z\rvert = r_0, dz=ir0eiθ0dθdz = ir_0e^{i\theta_0}d\theta and df=−3ir0e−3iθ0dθdf = -3ir_0e^{-3i\theta_0}d\theta, so the quotient is −3e−4iθ0-3e^{-4i\theta_0}. The two differ, so no derivative exists at z0z_0 either.

  • Continuity of the partials does its real work once, in getting the first-order expansions in part 2 (it is used again only to read off that f′f' is continuous, which is immediate). It cannot simply be deleted from that argument: mere existence of the four partials at a point does not give the expansions. It is, however, not the weakest possible hypothesis; it is the one Ahlfors (Ahlfors (1979), Ch. 2, §1.2, p. 26) adopts explicitly, and the one we prove, and it is satisfied by every function met in practice.

  • Single-valuedness is part of the definition of analytic, as Arfken, Weber and Harris (Arfken, Weber & Harris (2013), §11.2, p. 472) stress. The multivalued Ln⁡z\operatorname{Ln} z satisfies the Cauchy-Riemann equations on C∖{0}\mathbb{C}\setminus\{0\} branch by branch but is analytic only after a cut is made (branch points and branch cuts).

  • Conjugacy is a relation, not a symmetry. If vv is the conjugate of uu then uu is the conjugate of −v-v, not of vv. For (v,u)(v,u) to be a conjugate pair in that order, the first Cauchy-Riemann equation would have to read vx=uyv_x = u_y; what the equations for (u,v)(u,v) actually give is vx=−uyv_x = -u_y, so the sign is wrong unless uy=0u_y = 0. Concretely, for z2=(x2−y2)+i(2xy)z^2 = (x^2-y^2) + i(2xy) we have vx=2yv_x = 2y and uy=−2yu_y = -2y, which is the relation vx=−uyv_x = -u_y and not vx=uyv_x = u_y; and indeed z2z^2 is analytic while 2xy+i(x2−y2)=iz2‾2xy + i(x^2-y^2) = i\overline{z^2} is nowhere analytic, its Cauchy-Riemann equations holding at the origin only.

  • A test for non-analyticity. f(z)=zˉf(z) = \bar z has u=xu = x, v=−yv = -y, so ux=1≠−1=vyu_x = 1 \ne -1 = v_y everywhere: zˉ\bar z is continuous on C\mathbb{C} and differentiable nowhere. The same computation kills ∣z∣2=zzˉ\lvert z\rvert^2 = z\bar z except at the origin, and Re⁡z\operatorname{Re} z everywhere.

Example(Constructing a harmonic conjugate)Original

Let u(x,y)=x3−3xy2+2yu(x,y) = x^3 - 3xy^2 + 2y. First check that a conjugate can exist: uxx=6xu_{xx} = 6x and uyy=−6xu_{yy} = -6x, so uu is harmonic on R2\mathbb{R}^2.

  1. Integrate the first equation. vy=ux=3x2−3y2v_y = u_x = 3x^2 - 3y^2, so integrating in yy at fixed xx,

    v(x,y)=3x2y−y3+φ(x)v(x,y) = 3x^2 y - y^3 + \varphi(x)

    with φ\varphi an unknown function of xx alone.

  2. Impose the second equation. vx=6xy+φ′(x)v_x = 6xy + \varphi'(x) must equal −uy=6xy−2-u_y = 6xy - 2, so φ′(x)=−2\varphi'(x) = -2 and φ(x)=−2x+c\varphi(x) = -2x + c.

  3. Conclusion. v(x,y)=3x2y−y3−2x+cv(x,y) = 3x^2y - y^3 - 2x + c, and by part 2 of the statement (the partials here are polynomials, hence continuous) f=u+ivf = u + iv is entire.

  4. Identify ff. Taking c=0c = 0 and comparing with z3=(x3−3xy2)+i(3x2y−y3)z^3 = (x^3 - 3xy^2) + i(3x^2y - y^3) and −2iz=2y−2ix-2iz = 2y - 2ix, we get f(z)=z3−2izf(z) = z^3 - 2iz. Check at z=2+iz = 2+i: the formula gives (2+i)3−2i(2+i)=(2+11i)+(2−4i)=4+7i(2+i)^3 - 2i(2+i) = (2 + 11i) + (2 - 4i) = 4 + 7i, while u(2,1)=8−6+2=4u(2,1) = 8 - 6 + 2 = 4 and v(2,1)=12−1−4=7v(2,1) = 12 - 1 - 4 = 7. They agree.

  5. Derivative two ways. f′(z)=3z2−2if'(z) = 3z^2 - 2i gives f′(2+i)=3(3+4i)−2i=9+10if'(2+i) = 3(3+4i) - 2i = 9 + 10i; and ux+ivxu_x + iv_x at (2,1)(2,1) is (12−3)+i(12−2)=9+10i(12 - 3) + i(12 - 2) = 9 + 10i.

Example(Orthogonal level curves)AdaptedAdapted from Arfken, Weber & Harris (2013), §11.2, Eqs. (11.14)–(11.15), pp. 473–474

Where ∇u≠0\nabla u \ne 0, the Cauchy-Riemann equations make the level curves of uu and of vv meet at right angles, because

∇u⋅∇v=uxvx+uyvy=ux(−uy)+uy(ux)=0.\nabla u \cdot \nabla v = u_x v_x + u_y v_y = u_x(-u_y) + u_y(u_x) = 0 .

Take f(z)=z2f(z) = z^2, so u=x2−y2u = x^2 - y^2 and v=2xyv = 2xy. The curves u=cu = c are the rectangular hyperbolas x2−y2=cx^2 - y^2 = c and the curves v=dv = d are xy=d/2xy = d/2. At the point z0=2+iz_0 = 2+i we have u=3u = 3, v=4v = 4, and the gradients are ∇u=(2x,−2y)=(4,−2)\nabla u = (2x, -2y) = (4,-2) and ∇v=(2y,2x)=(2,4)\nabla v = (2y, 2x) = (2,4): their dot product is 8−8=08 - 8 = 0, and the tangent directions (2,4)/20(2,4)/\sqrt{20} and (4,−2)/20(4,-2)/\sqrt{20} are perpendicular unit vectors. The orthogonality fails exactly where f′(z)=2zf'(z) = 2z vanishes, that is at the origin, where the two hyperbola families degenerate into the pairs of lines y=±xy = \pm x and xy=0xy = 0 crossing at 45∘45^{\circ} instead of 90∘90^{\circ}.

ApplicationOriginal

Part 3 is the bridge between complex analysis and two-dimensional field theory. Any analytic f=u+ivf = u + iv on a plane region supplies, for free, a pair of solutions of Laplace's equation whose level sets are mutually orthogonal, so ff can be read as a complex potential: in electrostatics the curves u=constu = \text{const} are equipotentials and v=constv = \text{const} are field lines; in ideal fluid flow uu is the velocity potential and vv the stream function, and f′(z)‾\overline{f'(z)} is the velocity field. Problems that are awkward in real variables — the field near an edge, flow past a cylinder, the capacitance of an odd-shaped pair of conductors — become a search for the right analytic function, which is the subject of conformal mapping and connects to harmonic functions.

Read the other way, the equations are a tight constraint and so a strong tool. In dispersion theory a response function is analytic in a half-plane because of causality; the Cauchy-Riemann equations then tie its real part (dispersion) to its imaginary part (absorption) point by point, which is the local form of the Kramers-Kronig relations (dispersion relations). In elasticity, the Airy stress function and its conjugate play the same role, and in aerodynamics the Joukowski aerofoil is literally the image of a circle under an analytic map.

Definition

Branch points and branch cuts

Adapted

How a multivalued function such as a logarithm or a fractional power is made single-valued and analytic: locate the branch points, where continuation around a loop changes the value, and cut the plane so that no such loop is left.

Prerequisites: Complex numbers and elementary functions, Analyticity and the Cauchy-Riemann equations

Definition(Continuation along a path, branch point, cut, branch)AdaptedAdapted from Arfken, Weber & Harris (2013), §11.6, pp. 499–501; Riley, Hobson & Bence (2006), §24.5, pp. 835–836

Let ff be a multivalued function: a rule assigning to each zz in a set a finite or infinite set of values, such as Ln⁡z\operatorname{Ln} z or z1/2z^{1/2} (see the elementary functions).

Fix z0z_0 and one of the values f(z0)f(z_0). Along a path γ\gamma starting at z0z_0 we may try to choose, at each point, the value that varies continuously with position. On a path that never revisits a point this prescription determines one value at each point of γ\gamma: the continuation of the chosen value along γ\gamma. The question is what happens when γ\gamma closes up and returns to z0z_0.

  • A point z∗z_* is a branch point of ff if continuation once around some sufficiently small circle about z∗z_* returns a value of ff at the starting point different from the one it started with. The order of a branch point is the number of circuits needed before the original value returns; if no number of circuits restores it, the order is infinite.
  • A branch cut (or branch line) is a curve in the plane, declared off limits to continuation paths, which starts at a branch point and ends either at infinity or at another branch point. It is chosen so that no admissible closed path encircles a single branch point around which the value would change.
  • With the cuts in place, a choice of value at one point determines, by continuation, one value at every point of the cut plane. The resulting single-valued function is a branch of ff; a conventionally preferred one is the principal branch and its values are principal values.

A branch is discontinuous across its cut: a cut is a genuine boundary of the region of analyticity, not a bookkeeping device, so integrals along the two sides of a cut in opposite directions do not cancel. This is unlike the auxiliary barriers used to extend Cauchy's theorem to multiply connected regions, across which the integrand is continuous.

Definition(The branch structure of the logarithm and of powers)AdaptedAdapted from Arfken, Weber & Harris (2013), §11.6, Example 11.6.3 and Eqs. (11.55)–(11.56), p. 501; Riley, Hobson & Bence (2006), §24.4, Eqs. (24.20)–(24.22), p. 834; Ahlfors (1979), Ch. 2, §3.4, p. 46 and Ch. 4, §1.3, p. 108 (part 3)
  1. Logarithm. Ln⁡z=log⁡∣z∣+i(θ+2kπ)\operatorname{Ln} z = \log\lvert z\rvert + i(\theta + 2k\pi), k∈Zk \in \mathbb{Z}, has a single branch point, at z=0z = 0, and it is of infinite order: each counterclockwise circuit of the origin adds 2πi2\pi i. One cut from 00 to ∞\infty, along any simple arc, makes Ln⁡\operatorname{Ln} single-valued. With the cut along the negative real axis, the principal branch is

    ln⁡z=log⁡∣z∣+iArg⁡z,Arg⁡z∈(−π,π),\ln z = \log\lvert z\rvert + i\operatorname{Arg} z, \qquad \operatorname{Arg} z \in (-\pi,\pi),

    which is analytic on C∖(−∞,0]\mathbb{C}\setminus(-\infty,0] with

    ddzln⁡z=1z.\frac{d}{dz}\ln z = \frac{1}{z} .

    The open interval here is deliberate, and differs by one endpoint from the convention of the elementary functions, where Arg⁡z∈(−π,π]\operatorname{Arg} z \in (-\pi,\pi] so that Arg⁡\operatorname{Arg} and ln⁡\ln are defined at every z≠0z \ne 0, the negative real axis included. The cut plane C∖(−∞,0]\mathbb{C}\setminus(-\infty,0] is where ln⁡\ln is analytic, and it omits that axis; on the axis itself ln⁡\ln has the value log⁡∣z∣+iπ\log\lvert z\rvert + i\pi but is not continuous, as item 3 of the examples below shows. Both ranges are in use in this unit, for those two different purposes.

  2. Powers. zp=exp⁡(pLn⁡z)z^{p} = \exp(p \operatorname{Ln} z) for p∈Cp \in \mathbb{C} is single-valued if p∈Zp \in \mathbb{Z}; it has exactly qq branches, and a branch point of order qq at z=0z = 0, if p=s/qp = s/q is rational in lowest terms with q∈Nq \in \mathbb{N}; and it has infinitely many branches otherwise. In particular z1/2z^{1/2} has a branch point of order 22 at the origin and two branches, with (z1/2)2=z(z^{1/2})^2 = z on each.

  3. No branch survives on an annulus. Let 0≤ρ1<ρ2≤∞0 \le \rho_1 < \rho_2 \le \infty. There is no single-valued continuous branch of Ln⁡(z−a)\operatorname{Ln}(z-a) on the annulus ρ1<∣z−a∣<ρ2\rho_1 < \lvert z-a\rvert < \rho_2 — hence none that is analytic, and none on a punctured disk about aa either. The reason is that such a branch would deliver a continuous argument along a circle ∣z−a∣=ρ\lvert z - a\rvert = \rho with ρ1<ρ<ρ2\rho_1 < \rho < \rho_2, and that circle is exactly the loop on which part 1 shows the argument to increase by 2π2\pi.

    This is Ahlfors (Ahlfors (1979), Ch. 4, §1.3, p. 108)'s observation. He deduces it from the value 2πi2\pi i of the integral of 1/(z−a)1/(z-a) around such a circle, a route that needs contour integrals and so is not open to us here; the proof below argues from the argument instead, and in exchange needs no analyticity.

Note(Two conventions, and the sheets)AdaptedAdapted from Hassani (2013), §12.2, Definition 12.2.1, p. 365 and §12.2.1, p. 366

Hassani (Hassani (2013), Definition 12.2.1, p. 365) defines a branch point directly by the loop condition f(z0+reiθ)≠f(z0+rei(θ+2π))f(z_0 + re^{i\theta}) \ne f(z_0 + re^{i(\theta + 2\pi)}), which is the one-circuit case of the definition above. On that reading he also counts the point at infinity as a branch point of Ln⁡z\operatorname{Ln} z, since Ln⁡(1/z)=−Ln⁡z\operatorname{Ln}(1/z) = -\operatorname{Ln} z has a branch point at the origin; we follow him in treating ∞\infty on the same footing when deciding where a cut may end. Note the consequence for part 1: the cut from 00 to ∞\infty joins two branch points, exactly as the general rule asks.

He also records the bookkeeping that gives the branches names. With the cut along (−∞,0](-\infty,0], set

fn(reiθ)=log⁡r+i(θ+2nπ),−π<θ<π, r>0, n∈Z.f_n(re^{i\theta}) = \log r + i(\theta + 2n\pi), \qquad -\pi < \theta < \pi,\ r > 0,\ n \in \mathbb{Z} .

Each fnf_n is single-valued and analytic on the cut plane, f0=ln⁡f_0 = \ln is the principal branch, and fnf_n reproduces on the cut plane the values Ln⁡\operatorname{Ln} would take for (2n−1)π<θ<(2n+1)π(2n-1)\pi < \theta < (2n+1)\pi. Gluing the upper edge of the cut in the fnf_n copy to the lower edge in the fn+1f_{n+1} copy builds the Riemann surface of the logarithm, an infinite spiral of sheets on which the logarithm is genuinely single-valued; see Riemann surfaces.

Idea of the proofOriginal

Everything follows from one fact about the exponential: its period. Because ewe^{w} is unchanged by w↦w+2πiw \mapsto w + 2\pi i, inverting it can only ever return an answer modulo 2πi2\pi i, and which representative continuation hands back depends on how the argument of zz has moved. So the bookkeeping reduces to tracking arg⁡\arg along the path: a loop that encircles the origin advances arg⁡z\arg z by 2π2\pi whether we like it or not, and a loop that does not encircle it returns arg⁡z\arg z to its starting value. For zpz^p the shift 2πi2\pi i in Ln⁡z\operatorname{Ln} z becomes the factor e2πipe^{2\pi i p}, which is 11 exactly when pp is an integer and of finite order qq exactly when pp has denominator qq. Analyticity of the branch, finally, is the Cauchy-Riemann test applied in polar coordinates.

ProofAdaptedAdapted from Arfken, Weber & Harris (2013), §11.6, Examples 11.6.1 and 11.6.3, pp. 499–501; Riley, Hobson & Bence (2006), §24.5, pp. 835–836; Arfken, Weber & Harris (2013), §11.2, Example 11.2.3, p. 475; Ahlfors (1979), Ch. 4, §1.3, p. 108 (part 3, statement only: the argument below is our own)

1. The logarithm. Let z0=r0eiθ0z_0 = r_0 e^{i\theta_0} with r0>0r_0 > 0, and continue a chosen value log⁡r0+iθ0\log r_0 + i\theta_0 of Ln⁡\operatorname{Ln} along the circle z(t)=r0ei(θ0+t)z(t) = r_0e^{i(\theta_0 + t)}, 0≤t≤2π0 \le t \le 2\pi. The function t↦log⁡r0+i(θ0+t)t \mapsto \log r_0 + i(\theta_0 + t) is continuous and is at each tt a value of Ln⁡z(t)\operatorname{Ln} z(t), so it is the continuation; at t=2πt = 2\pi it equals the starting value plus 2πi2\pi i. Since 2πi≠02\pi i \ne 0 and n⋅2πi≠0n \cdot 2\pi i \ne 0 for every n∈Nn \in \mathbb{N}, the origin is a branch point of infinite order. For z∗≠0z_* \ne 0, a circle about z∗z_* of radius less than ∣z∗∣\lvert z_*\rvert lies in a half-plane not containing 00, on which a continuous argument exists; continuation around it returns the starting value, so no other point is a branch point.

Next, analyticity of ln⁡\ln off the cut. On C∖(−∞,0]\mathbb{C}\setminus(-\infty,0] write u=log⁡ru = \log r and v=Arg⁡zv = \operatorname{Arg} z with r=x2+y2r = \sqrt{x^2+y^2}. Then

ux=xr2,uy=yr2,vx=−yr2,vy=xr2,u_x = \frac{x}{r^2}, \quad u_y = \frac{y}{r^2}, \qquad v_x = -\frac{y}{r^2}, \quad v_y = \frac{x}{r^2},

the formulas for vv coming from differentiating tan⁡v=y/x\tan v = y/x (or cot⁡v=x/y\cot v = x/y near the imaginary axis) and holding on the cut plane because Arg⁡\operatorname{Arg} is there a single-valued smooth function. These four functions are continuous on C∖(−∞,0]\mathbb{C}\setminus(-\infty,0] and satisfy ux=vyu_x = v_y and uy=−vxu_y = -v_x, so by the sufficiency half of the Cauchy-Riemann theorem ln⁡\ln is analytic there, with

ddzln⁡z=ux+ivx=x−iyr2=zˉzzˉ=1z.\frac{d}{dz}\ln z = u_x + i v_x = \frac{x - iy}{r^2} = \frac{\bar z}{z \bar z} = \frac1z .

2. Powers. By definition the values of zpz^p at a fixed z≠0z \ne 0 are exp⁡(p ℓ)\exp(p\,\ell) as ℓ\ell runs over the values of Ln⁡z\operatorname{Ln} z; two of them differ by a factor e2πikpe^{2\pi i k p}, k∈Zk \in \mathbb{Z}. One counterclockwise circuit of the origin replaces ℓ\ell by ℓ+2πi\ell + 2\pi i, hence multiplies the continued value by e2πipe^{2\pi i p}, and nn circuits multiply it by e2πinpe^{2\pi i n p}. This factor is 11 if and only if np∈Znp \in \mathbb{Z}. So:

  • if p∈Zp \in \mathbb{Z}, already one circuit changes nothing and zpz^p is single-valued;
  • if p=s/qp = s/q in lowest terms with q∈Nq \in \mathbb{N}, then np∈Znp \in \mathbb{Z} first at n=qn = q (because gcd⁡(s,q)=1\gcd(s,q) = 1), so the branch point at 00 has order qq, and the distinct values at a fixed zz are the qq numbers e2πiks/qexp⁡(pln⁡z)e^{2\pi i k s/q}\exp(p\ln z), k=0,…,q−1k = 0,\dots,q-1;
  • otherwise np∉Znp \notin \mathbb{Z} for every n∈Nn \in \mathbb{N}, the order is infinite, and the values e2πikpe^{2\pi i kp} are distinct for distinct kk.

For p=1/2p = 1/2 this gives two branches, ±exp⁡(12ln⁡z)\pm\exp(\tfrac12 \ln z), and squaring either gives exp⁡(ln⁡z)=z\exp(\ln z) = z.

3. No branch on an annulus. Fix ρ\rho with ρ1<ρ<ρ2\rho_1 < \rho < \rho_2, so that the circle ∣z−a∣=ρ\lvert z - a\rvert = \rho lies in the annulus, and suppose a single-valued continuous LL existed there with eL(z)=z−ae^{L(z)} = z - a at every point — which is what "a branch of Ln⁡(z−a)\operatorname{Ln}(z-a)" means. Set

φ(t)=Im⁡L(a+ρeit),0≤t≤2π.\varphi(t) = \operatorname{Im} L\bigl(a + \rho e^{it}\bigr), \qquad 0 \le t \le 2\pi .
  • φ\varphi is continuous, being the imaginary part of a continuous function of tt.
  • Exponentiating the defining property at z=a+ρeitz = a + \rho e^{it} gives eL=ρeite^{L} = \rho e^{it}, so eRe⁡L=ρe^{\operatorname{Re}L} = \rho and φ(t)\varphi(t) is an argument of eite^{it}. Two arguments of the same number differ by an integer multiple of 2π2\pi, so φ(t)−t∈2πZ\varphi(t) - t \in 2\pi\mathbb{Z} for every tt.
  • The function t↦φ(t)−tt \mapsto \varphi(t) - t is therefore a continuous map of the interval [0,2π][0,2\pi] into the set 2πZ2\pi\mathbb{Z}. Its image is a connected subset of 2πZ2\pi\mathbb{Z}, hence a single point: φ(t)=t+2πk\varphi(t) = t + 2\pi k for one fixed k∈Zk \in \mathbb{Z}.

Consequently φ(2π)−φ(0)=2π≠0\varphi(2\pi) - \varphi(0) = 2\pi \ne 0. But a+ρe2πia + \rho e^{2\pi i} and a+ρei0a + \rho e^{i0} are the same point of the annulus, and LL was assumed single-valued, so φ(2π)=φ(0)\varphi(2\pi) = \varphi(0). The contradiction rules out such an LL. Only continuity was used, so the analytic case is covered as well, and the argument applies verbatim to a punctured disk (ρ1=0\rho_1 = 0) and to the punctured plane. □\square

Method(Choosing cuts for a product of powers)AdaptedAdapted from Riley, Hobson & Bence (2006), §24.5, pp. 836–837; Arfken, Weber & Harris (2013), §11.6, Example 11.6.4, pp. 502–503

For f(z)=∏j=1m(z−aj)pjf(z) = \prod_{j=1}^{m}(z - a_j)^{p_j}, and more generally for any function built from logarithms and powers of polynomials:

  1. List the candidate branch points. These are the points aja_j where a factor vanishes (and, after the substitution z=1/tz = 1/t, possibly ∞\infty).

  2. Write each factor in polar form about its own branch point: z−aj=rjeiθjz - a_j = r_j e^{i\theta_j}, so that

    f(z)=(∏jrj pj)exp⁡(i∑jpjθj)f(z) = \Bigl(\prod_j r_j^{\,p_j}\Bigr)\exp\Bigl(i\sum_j p_j\theta_j\Bigr)

    on a chosen branch (for real pjp_j).

  3. Take each closed loop in turn and record how each θj\theta_j changes: by 2π2\pi if the loop encircles aja_j counterclockwise, by 00 if not. The continued value is multiplied by exp⁡(2πi∑j∈Spj)\exp(2\pi i \sum_{j \in S} p_j), where SS is the set of encircled branch points.

  4. A loop is harmless when that factor is 11, that is when ∑j∈Spj∈Z\sum_{j\in S}p_j \in \mathbb{Z}. Only the harmful loops must be blocked.

  5. Draw cuts that block exactly the harmful loops. Each cut runs from a branch point to another branch point or to infinity; its shape is free, only its endpoints matter. If a set SS of branch points is harmless when encircled together, a single finite cut joining them is enough, and is often the most convenient choice.

  6. Fix the branch by naming the value of ff at one point of the cut plane, or equivalently the range of each θj\theta_j.

Example(A square root with two branch points and a finite cut)AdaptedAdapted from Riley, Hobson & Bence (2006), §24.5, pp. 836–837

The loop analysis in steps 1 to 3 follows Riley, Hobson and Bence (Riley, Hobson & Bence (2006), §24.5, pp. 836–837)'s route for (z2+1)1/2(z^2+1)^{1/2}; the function, the construction in step 4 and all the numbers are ours.

Let f(z)=((z−1)(z−4))1/2f(z) = \bigl((z-1)(z-4)\bigr)^{1/2}, with p1=p2=12p_1 = p_2 = \tfrac12 at a1=1a_1 = 1, a2=4a_2 = 4.

  1. Candidate branch points: z=1z = 1 and z=4z = 4. At infinity, putting z=1/tz = 1/t gives f=t−1((1−t)(1−4t))1/2f = t^{-1}\bigl((1-t)(1-4t)\bigr)^{1/2}, whose second factor is single-valued near t=0t = 0; so ∞\infty is not a branch point.

  2. Loop test. A loop around z=1z=1 alone, or around z=4z=4 alone, multiplies ff by e2πi/2=−1e^{2\pi i/2} = -1: harmful. A loop around both multiplies it by e2πi(1/2+1/2)=1e^{2\pi i(1/2 + 1/2)} = 1: harmless. A loop around neither does nothing.

  3. Cut. It therefore suffices to block the two single loops while allowing the double loop, and the finite segment [1,4][1,4] does exactly that. (Two cuts (−∞,1](-\infty,1] and [4,∞)[4,\infty) would also work but forbid more than necessary.)

  4. A concrete branch, built so that its analyticity is visible. The naive recipe z−1 z−4\sqrt{z-1}\,\sqrt{z-4} with two principal square roots is cut along (−∞,4](-\infty,4], which is more than we want, so we build the branch out of a single principal logarithm instead. Let

    m(z)=z−1z−4,F(z)=(z−4)exp⁡(12ln⁡m(z))on C∖[1,4],m(z) = \frac{z-1}{z-4}, \qquad F(z) = (z-4)\exp\Bigl(\tfrac12 \ln m(z)\Bigr) \quad\text{on } \mathbb{C}\setminus[1,4] ,

    with ln⁡\ln the principal logarithm, analytic on C∖(−∞,0]\mathbb{C}\setminus(-\infty,0] by part 1 of the statement. Three observations finish the construction.

    • mm omits the cut of ln⁡\ln. mm is a Möbius map of the sphere, so it is injective; it sends 1↦01 \mapsto 0, 4↦∞4\mapsto\infty, and a real xx with 1<x<41 < x < 4 to (x−1)/(x−4)<0(x-1)/(x-4) < 0. So mm carries the segment [1,4][1,4] onto (−∞,0]∪{∞}(-\infty,0]\cup\{\infty\} and, being injective, carries nothing else there. For z∉[1,4]z \notin [1,4] we therefore have m(z)∈C∖(−∞,0]m(z) \in \mathbb{C}\setminus(-\infty,0], which is where ln⁡\ln lives.
    • Hence FF is analytic on all of C∖[1,4]\mathbb{C}\setminus[1,4], being a product of the polynomial z−4z-4 with a composition of analytic functions. No question of continuity across anything arises: the formula is analytic where it stands.
    • FF is a branch of ff. exp⁡(12ln⁡w)2=exp⁡(ln⁡w)=w\exp(\tfrac12\ln w)^2 = \exp(\ln w) = w, so F(z)2=(z−4)2m(z)=(z−4)(z−1)F(z)^2 = (z-4)^2 m(z) = (z-4)(z-1).

    And FF is the naive product after all: evaluating both at every point of a grid over −8≤Re⁡z≤12-8 \le \operatorname{Re}z \le 12, ∣Im⁡z∣≤6\lvert\operatorname{Im}z\rvert \le 6 (about 48 00048\,000 points, the segment [1,4][1,4] excluded) gives F(z)=z−1z−4F(z) = \sqrt{z-1}\sqrt{z-4} to machine precision everywhere, the two sign jumps of the principal roots on (−∞,1)(-\infty,1) cancelling each other. The point of building FF the other way is that the construction proves the analyticity that the cancellation only suggests: continuity across a line, together with analyticity off it, is not by itself a reason for analyticity across the line.

  5. Numbers. F(0)=(−4)exp⁡(12ln⁡14)=−4⋅12=−2F(0) = (-4)\exp\bigl(\tfrac12\ln\tfrac14\bigr) = -4\cdot\tfrac12 = -2, and indeed F(0)2=4=(0−1)(0−4)F(0)^2 = 4 = (0-1)(0-4). F(10)=6exp⁡(12ln⁡32)=63/2=36≈7.348469F(10) = 6\exp\bigl(\tfrac12\ln\tfrac32\bigr) = 6\sqrt{3/2} = 3\sqrt 6 \approx 7.348469, consistent with F(z)≈z−52F(z) \approx z - \tfrac52 for large zz (at z=10z=10 that estimate gives 7.57.5). Across the cut, m(2.5+iε)m(2.5 + i\varepsilon) tends to −1-1 through the lower half-plane and m(2.5−iε)m(2.5 - i\varepsilon) through the upper half-plane, so 12ln⁡m→∓iπ/2\tfrac12\ln m \to \mp i\pi/2 and the two limits of FF at z=2.5z = 2.5 are (−1.5)(∓i)=±1.5 i(-1.5)(\mp i) = \pm 1.5\,i: equal and opposite, as a sign-changing pair of branch points requires. Note that this branch is negative for real z<1z<1 and positive for real z>4z>4; a single-valued branch need not have a single sign.

Example(Values and orders)Original
  1. Order three. z1/3z^{1/3} has a branch point of order 33 at the origin: one counterclockwise circuit multiplies a value by e2πi/3e^{2\pi i/3}, two circuits by e4πi/3e^{4\pi i/3}, three by 11. Starting from the value 22 at z=8z = 8 and going once round, we return with 2e2πi/3=−1+i32e^{2\pi i/3} = -1 + i\sqrt3, whose cube is 88 as well.
  2. A logarithm with two cuts, or one. g(z)=ln⁡z−1z+1g(z) = \ln\frac{z-1}{z+1} has branch points at ±1\pm 1, each of infinite order, since a loop around +1+1 alone adds 2πi2\pi i and a loop around −1-1 alone subtracts it. A loop around both adds 00, so the finite cut [−1,1][-1,1] suffices; this is the cut taken for Legendre's functions of the second kind QnQ_n, the first of which is Q0(z)=12ln⁡z+1z−1Q_0(z) = \tfrac12\ln\frac{z+1}{z-1}.
  3. Principal values on the cut. With the convention Arg⁡z∈(−π,π]\operatorname{Arg} z \in (-\pi,\pi], the principal ln⁡\ln is defined at z=−4z = -4, with ln⁡(−4)=log⁡4+iπ≈1.386294+3.141593 i\ln(-4) = \log 4 + i\pi \approx 1.386294 + 3.141593\,i, but it is not continuous there: ln⁡(−4−iε)→log⁡4−iπ\ln(-4 - i\varepsilon) \to \log 4 - i\pi. The cut is the half-line where the function must be regarded as having two boundary values, differing by 2πi2\pi i.
ApplicationOriginal

Branch cuts are where the physics lives, not a technicality. Three examples.

  • Propagating versus evanescent waves. A wave number k=ω2/c2−κ2k = \sqrt{\omega^2/c^2 - \kappa^2}, as a function of frequency, is a two-valued function with branch points at ω=±cκ\omega = \pm c\kappa: the cut separates the propagating band from the evanescent one, and which branch of the square root is taken is fixed by the physical requirement of decay at infinity. The same choice appears in waveguides, in the Sommerfeld half-space problem and in the index of refraction of a plasma.
  • Thresholds in scattering. A scattering amplitude, continued in the energy, has a branch point at each threshold where a new channel opens; the discontinuity across the cut is the production cross-section, by the optical theorem, and the sheet reached by crossing the cut carries the resonance poles. This is the setting of dispersion relations, and the multi-sheeted picture is made precise by Riemann surfaces.
  • Integrals around a cut. Because a branch is discontinuous across its cut, the two sides contribute with opposite signs rather than cancelling, and the resulting contour integral evaluates real integrals that elementary methods miss — for instance ∫0∞x−s/(1+x) dx\int_0^\infty x^{-s}/(1+x)\,dx via a keyhole contour (integrals with branch cuts).

In numerical work the practical point is sharper still: library functions implement principal branches, so a formula assembled from them can be wrong by a sign or a factor e2πipe^{2\pi i p} on part of the plane. The example above, where z−1z−4\sqrt{z-1}\sqrt{z-4} turned out to be analytic across (−∞,1)(-\infty,1) although each factor is not, is the kind of check worth making before trusting such a formula.

Contour integration

4

Definition

Contour integrals

Adapted

The integral of a continuous complex function along a piecewise smooth path, its independence of the parametrization, the ML bound on its size, and the fact that it depends only on the endpoints exactly when the integrand is the derivative of an analytic function.

Prerequisites: Complex numbers and elementary functions, Analyticity and the Cauchy-Riemann equations, Line integrals and conservative fields

Definition(Paths and contour integrals)AdaptedAdapted from Ahlfors (1979), Ch. 4, §1.1, Eqs. (1) and (4), pp. 101–102; Arfken, Weber & Harris (2013), §11.3, Eqs. (11.16)–(11.17), pp. 477–478

For a continuous g:[a,b]→Cg:[a,b] \to \mathbb{C}, g=g1+ig2g = g_1 + ig_2 with g1,g2g_1, g_2 real, define

∫abg(t) dt=∫abg1(t) dt+i∫abg2(t) dt.\int_a^b g(t)\,dt = \int_a^b g_1(t)\,dt + i \int_a^b g_2(t)\,dt .

A path (or contour) γ\gamma is a piecewise continuously differentiable map z:[a,b]→Cz:[a,b] \to \mathbb{C}; it is closed if z(a)=z(b)z(a) = z(b), and simple if it does not meet itself otherwise. For ff continuous on the image of γ\gamma,

∫γf(z) dz=∫abf(z(t)) z′(t) dt,\int_\gamma f(z)\,dz = \int_a^b f\bigl(z(t)\bigr)\,z'(t)\,dt ,

the interval being subdivided at the finitely many points where z′z' jumps. Writing f=u+ivf = u+iv and z=x+iyz = x+iy this is the pair of real line integrals

∫γf(z) dz=∫γ(u dx−v dy)+i∫γ(u dy+v dx),\int_\gamma f(z)\,dz = \int_\gamma (u\,dx - v\,dy) + i \int_\gamma (u\,dy + v\,dx) ,

which is how it is computed in practice. It is also the limit of the Riemann sums

∑jf(ζj)(zj−zj−1)\sum_j f(\zeta_j)(z_j - z_{j-1})

over refining subdivisions a=t0<⋯<tN=ba = t_0 < \dots < t_N = b, with zj=z(tj)z_j = z(t_j) and ζj\zeta_j a point of the arc between zj−1z_{j-1} and zjz_j; this is the form in which Arfken, Weber and Harris (Arfken, Weber & Harris (2013), §11.3, Eq. (11.16), p. 477) introduce the contour integral. For a piecewise C1C^1 path that is a small theorem rather than a restatement of the definition, and one line proves it: writing zj−zj−1=∫tj−1tjz′ dtz_j - z_{j-1} = \int_{t_{j-1}}^{t_j}z'\,dt turns the difference between the sum and the integral into ∑j∫tj−1tj(f(ζj)−f(z(t)))z′(t) dt\sum_j \int_{t_{j-1}}^{t_j}\bigl(f(\zeta_j) - f(z(t))\bigr)z'(t)\,dt, whose modulus is at most L(γ)L(\gamma) times the largest oscillation of f∘zf \circ z over a subinterval, and that tends to 00 with the mesh because f∘zf\circ z is uniformly continuous on [a,b][a,b]. Nothing below uses the Riemann-sum form; it is recorded because it is how the integral is often defined.

Two further integrals along γ\gamma are used below. With respect to arc length,

∫γf ∣dz∣=∫abf(z(t)) ∣z′(t)∣ dt,soL(γ)=∫γ∣dz∣\int_\gamma f\,\lvert dz\rvert = \int_a^b f\bigl(z(t)\bigr)\,\lvert z'(t)\rvert\,dt, \qquad\text{so}\qquad L(\gamma) = \int_\gamma \lvert dz\rvert

is the length of γ\gamma. A closed contour is traversed in the positive sense when the region it encloses lies to its left, which for a circle about an interior point means counterclockwise, the direction of increasing polar angle. We write ∮γ\oint_\gamma for an integral over a closed contour.

Definition(Properties of the contour integral)AdaptedAdapted from Ahlfors (1979), Ch. 4, §1.1, Eqs. (2), (3), (5), (6), (9), pp. 102–104, and §1.3, Thm 1 and the statement following it, pp. 106–107; Riley, Hobson & Bence (2006), §24.9, Eq. (24.39), p. 849

Let ff be continuous on the path γ\gamma, with z:[a,b]→Cz:[a,b]\to\mathbb{C}.

  1. Linearity and subdivision. ∫γ(αf+βg) dz=α∫γf dz+β∫γg dz\int_\gamma (\alpha f + \beta g)\,dz = \alpha\int_\gamma f\,dz + \beta \int_\gamma g\,dz for α,β∈C\alpha,\beta \in \mathbb{C}, and if γ\gamma is split into consecutive subarcs γ1,…,γn\gamma_1,\dots,\gamma_n then ∫γf dz=∑k∫γkf dz\int_\gamma f\,dz = \sum_k \int_{\gamma_k} f \, dz.

  2. Invariance and reversal. If t=φ(τ)t = \varphi(\tau) is a piecewise C1C^1 increasing bijection of [α,β][\alpha,\beta] onto [a,b][a,b], then z∘φz\circ\varphi gives the same integral: the contour integral depends on γ\gamma as an oriented curve, not on its parametrization. Reversing the orientation changes the sign, ∫−γf dz=−∫γf dz\int_{-\gamma} f\,dz = -\int_\gamma f\,dz, whereas ∫−γf ∣dz∣=∫γf ∣dz∣\int_{-\gamma} f \,\lvert dz\rvert = \int_{\gamma} f\,\lvert dz\rvert.

  3. The ML estimate. If ∣f(z)∣≤M\lvert f(z)\rvert \le M at every point of γ\gamma and L=L(γ)L = L(\gamma), then

    ∣∫γf(z) dz∣  ≤  ∫γ∣f∣ ∣dz∣  ≤  ML.\Bigl\lvert \int_\gamma f(z)\,dz \Bigr\rvert \;\le\; \int_\gamma \lvert f\rvert\,\lvert dz\rvert \;\le\; ML .
  4. Fundamental theorem. Let Ω\Omega be an open set. The following are equivalent for ff continuous on Ω\Omega:

    (a) ∫γf dz\int_\gamma f\,dz depends only on the endpoints of γ\gamma, for all paths γ\gamma in Ω\Omega;

    (b) ∮γf dz=0\oint_\gamma f \,dz = 0 for every closed path γ\gamma in Ω\Omega;

    (c) f=F′f = F' for some analytic FF on Ω\Omega (a primitive, or indefinite integral, of ff).

    When they hold, ∫γf dz=F(z(b))−F(z(a))\int_\gamma f\,dz = F(z(b)) - F(z(a)). If Ω\Omega is disconnected, read (a) with care: no path joins two different components, so (a) imposes a condition within each component and nothing between them, and the primitive of (c) is fixed only up to a separate constant on each component. The proof of (a) ⇒\Rightarrow (c) handles this by working one component at a time.

  5. Powers. For every integer n≥0n \ge 0 and every closed path γ\gamma, ∮γ(z−a)ndz=0\oint_\gamma (z-a)^n dz = 0. The same holds for n≤−2n \le -2 and every closed path avoiding aa. For n=−1n = -1 it fails: on the positively oriented circle CC of radius ρ\rho about aa,

    ∮Cdzz−a=2πi.\oint_{C} \frac{dz}{z-a} = 2\pi i .
Idea of the proofOriginal

Parts 1 and 2 are the chain rule and the substitution rule for real integrals, imported one real and one imaginary component at a time. Part 3 has one non-obvious ingredient: for a complex-valued integral over a real interval, ∣∫g∣≤∫∣g∣\lvert \int g\rvert \le \int \lvert g\rvert is not the triangle inequality for real integrals, because the modulus of a complex integral is not the integral of anything; the trick is to rotate the integral onto the positive real axis by a constant unimodular factor, after which taking real parts does the work. Part 4 is the two-variable statement that a line integral is path independent exactly when the integrand is an exact differential; the complex content is that "f dzf\,dz is exact" says that the potential satisfies the Cauchy-Riemann equations, so the potential is analytic with derivative ff. Part 5 is then a corollary for n≠−1n \ne -1, since (z−a)n(z-a)^n has the obvious primitive; and the single exception n=−1n = -1 is a one-line computation that will turn out to carry almost all the later theory.

ProofAdaptedAdapted from Ahlfors (1979), Ch. 4, §1.1, pp. 102–104 and §1.3, Thm 1 and pp. 107–108

1. Both claims are the corresponding statements for real integrals applied to real and imaginary parts, together with the fact that ∫abc g dt=c∫abg dt\int_a^b c\,g\,dt = c\int_a^b g\,dt for a complex constant c=α+iβc = \alpha + i\beta: writing g=g1+ig2g = g_1 + ig_2, both sides equal ∫(αg1−βg2)dt+i∫(αg2+βg1)dt\int(\alpha g_1 - \beta g_2)dt + i\int(\alpha g_2 + \beta g_1)dt.

2. By the real substitution rule applied to each component,

∫αβf(z(φ(τ)))z′(φ(τ))φ′(τ) dτ=∫abf(z(t))z′(t) dt,\int_\alpha^\beta f\bigl(z(\varphi(\tau))\bigr)z'(\varphi(\tau))\varphi'(\tau)\,d\tau = \int_a^b f(z(t))z'(t)\,dt ,

and the chain rule identifies z′(φ(τ))φ′(τ)z'(\varphi(\tau))\varphi'(\tau) as the derivative of z∘φz\circ\varphi, so the left side is the integral over the reparametrized path. For −γ-\gamma, parametrized by w(t)=z(−t)w(t) = z(-t) on [−b,−a][-b,-a], three factors of −1-1 appear and their product is what we want:

∫−γf dz=∫−b−af(z(−t))(−z′(−t)) dt=∫baf(z(s))(−z′(s))(−ds)=−∫abf(z(s))z′(s) ds,\int_{-\gamma}f\,dz = \int_{-b}^{-a}f\bigl(z(-t)\bigr)\bigl(-z'(-t)\bigr)\,dt = \int_{b}^{a} f\bigl(z(s)\bigr)\bigl(-z'(s)\bigr)(-ds) = -\int_a^b f\bigl(z(s)\bigr)z'(s)\,ds ,

one sign from the derivative w′(t)=−z′(−t)w'(t) = -z'(-t), one from dt=−dsdt = -ds under s=−ts = -t, and one from the exchange of the limits b,ab, a into a,ba, b: an odd number, hence ∫−γf dz=−∫γf dz\int_{-\gamma}f\,dz = -\int_\gamma f\,dz. With ∣dz∣\lvert dz\rvert the first of the three disappears, since ∣w′(t)∣=∣z′(−t)∣\lvert w'(t)\rvert = \lvert z'(-t)\rvert carries no sign, and the remaining two cancel, leaving ∫−γf ∣dz∣=∫γf ∣dz∣\int_{-\gamma}f\,\lvert dz\rvert = \int_\gamma f\,\lvert dz\rvert.

3. First, for continuous gg on [a,b][a,b] with a≤ba \le b we show ∣∫abg dt∣≤∫ab∣g∣ dt\bigl\lvert\int_a^b g\,dt\bigr\rvert \le \int_a^b \lvert g\rvert\,dt. If ∫abg dt=0\int_a^b g\,dt = 0 there is nothing to prove. Otherwise put θ=arg⁡∫abg dt\theta = \arg\int_a^b g\,dt and c=e−iθc = e^{-i\theta}; then by part 1,

∣∫abg dt∣=c∫abg dt=∫abcg dt=Re⁡∫abcg dt=∫abRe⁡(cg) dt≤∫ab∣cg∣ dt=∫ab∣g∣ dt,\Bigl\lvert \int_a^b g\,dt\Bigr\rvert = c\int_a^b g\,dt = \int_a^b cg\,dt = \operatorname{Re}\int_a^b cg\,dt = \int_a^b \operatorname{Re}(cg)\,dt \le \int_a^b \lvert cg\rvert \,dt = \int_a^b \lvert g\rvert\,dt,

where the third equality holds because the left-hand side is real, the fourth because taking real parts commutes with a real integral, and the inequality because Re⁡w≤∣w∣\operatorname{Re} w \le \lvert w\rvert pointwise and ∣c∣=1\lvert c\rvert = 1. Applying this to g(t)=f(z(t))z′(t)g(t) = f(z(t))z'(t) gives the first inequality of part 3, since ∣g∣=∣f(z(t))∣∣z′(t)∣\lvert g \rvert = \lvert f(z(t))\rvert\lvert z'(t)\rvert is the integrand of ∫γ∣f∣∣dz∣\int_\gamma \lvert f\rvert \lvert dz\rvert. For the second, use the hypothesis ∣f∣≤M\lvert f\rvert \le M on γ\gamma:

∫γ∣f∣ ∣dz∣≤M∫γ∣dz∣=ML.\int_\gamma \lvert f\rvert\,\lvert dz\rvert \le M\int_\gamma \lvert dz\rvert = ML .

4. (a) ⇔\Leftrightarrow (b). If the integral depends only on the endpoints and γ\gamma is closed, then γ\gamma and −γ-\gamma have the same endpoints, so ∫γf dz=∫−γf dz=−∫γf dz\int_\gamma f\,dz = \int_{-\gamma}f\,dz = -\int_\gamma f\,dz, giving ∮γf dz=0\oint_\gamma f\,dz=0. Conversely, if all closed integrals vanish and γ1,γ2\gamma_1,\gamma_2 share their endpoints, then γ1\gamma_1 followed by −γ2-\gamma_2 is closed, so by part 1 the two integrals agree.

(c) ⇒\Rightarrow (a). If f=F′f = F' with FF analytic, then on each smooth piece t↦F(z(t))t \mapsto F(z(t)) is differentiable with derivative F′(z(t))z′(t)=f(z(t))z′(t)F'(z(t))z'(t) = f(z(t))z'(t) by the chain rule; summing over the pieces and applying, to the real and to the imaginary part separately, the fundamental theorem of one-variable calculus gives ∫γf dz=F(z(b))−F(z(a))\int_\gamma f\,dz = F(z(b)) - F(z(a)), which depends only on the endpoints.

(a) ⇒\Rightarrow (c). Fix a point z0z_0 in a component of Ω\Omega and set F(z)=∫σf dzF(z) = \int_{\sigma} f\,dz along any path σ\sigma in Ω\Omega from z0z_0 to zz; by (a) this is well defined. To compute ∂F/∂x\partial F/\partial x at zz, take the last piece of σ\sigma to be a horizontal segment, which is possible because Ω\Omega is open and so contains a disk about zz: with xx as parameter, F(z)=const+∫xf(s+iy) dsF(z) = \text{const} + \int^x f(s + iy)\,ds, so ∂F/∂x=f(z)\partial F/\partial x = f(z) by the fundamental theorem of calculus (using continuity of ff). Choosing the last piece vertical instead gives ∂F/∂y=if(z)\partial F/\partial y = i f(z). Hence Fx=−iFyF_x = -iF_y, which is the Cauchy-Riemann condition in combined form, and Fx=fF_x = f is continuous; so by the sufficiency half of the Cauchy-Riemann theorem, FF is analytic with F′=Fx=fF' = F_x = f.

5. For n≥0n \ge 0, (z−a)n=ddz(z−a)n+1n+1(z-a)^n = \frac{d}{dz}\frac{(z-a)^{n+1}}{n+1} with the primitive analytic on all of C\mathbb{C}, so (c) ⇒\Rightarrow (b) of part 4 gives 00. For n≤−2n \le -2 the same primitive (z−a)n+1/(n+1)(z-a)^{n+1}/(n+1) is analytic on C∖{a}\mathbb{C}\setminus\{a\}, which contains γ\gamma, so again the integral is 00. For n=−1n = -1 the conclusion fails, and the direct computation on z(t)=a+ρeitz(t) = a + \rho e^{it}, 0≤t≤2π0 \le t \le 2\pi, shows it:

∮Cdzz−a=∫02πiρeitρeit dt=∫02πi dt=2πi,\oint_C \frac{dz}{z-a} = \int_0^{2\pi}\frac{i\rho e^{it}}{\rho e^{it}}\,dt = \int_0^{2\pi} i\,dt = 2\pi i,

independent of ρ\rho and of aa. Read through part 4 in the direction (b) ⇒\Rightarrow (c) this says something more: 1/(z−a)1/(z-a) has no primitive on any open set containing a circle about aa, since a primitive would force the integral to vanish. □\square

RemarkOriginal

Where the hypotheses are used.

  • Piecewise C1C^1 makes z′z' integrable, so the defining formula means something, and gives the length as an integral. Continuous paths of infinite length (nowhere differentiable curves) have no contour integral in this sense; the natural weakening is to rectifiable paths, which Ahlfors (Ahlfors (1979), Ch. 4, §1.2, pp. 104–105) treats and which the elementary theory does not need.
  • Continuity of ff on γ\gamma is what makes the integrand integrable. Nothing in parts 1 to 3 requires ff to be analytic, or even to be defined off γ\gamma: the ML estimate is a statement about continuous functions.
  • Analyticity of FF, not just of ff, is the content of part 4(c). It is the primitive that must be single-valued and analytic on all of Ω\Omega. The failure for n=−1n=-1 in part 5 is a failure of single-valuedness, not of differentiability: ln⁡(z−a)\ln(z-a) has derivative 1/(z−a)1/(z-a) on any cut plane, but no cut plane contains a circle around aa.
  • Path dependence is the generic case. For ff continuous but not analytic the integral usually does depend on the path, as the second example below shows; the discovery that for analytic ff it does not is Cauchy's theorem.
  • The ML bound can be attained, so it is not merely a crude estimate; see part 3 of the examples.
Example(A primitive, and path dependence)Original
  1. Using a primitive. Let γ\gamma be the arc of the unit circle from 11 to ii traversed counterclockwise. Since z2z^2 has the entire primitive z3/3z^3/3, part 4 gives

    ∫γz2 dz=i3−133=−1−i3≈−0.333333−0.333333 i,\int_\gamma z^2\,dz = \frac{i^3 - 1^3}{3} = \frac{-1 - i}{3} \approx -0.333333 - 0.333333\,i ,

    with no parametrization needed. Direct evaluation confirms it: ∫0π/2e2it ieit dt=[e3it/3]0π/2=(e3iπ/2−1)/3=(−i−1)/3\int_0^{\pi/2} e^{2it}\,ie^{it}\,dt = \bigl[e^{3it}/3\bigr]_0^{\pi/2} = (e^{3i\pi/2} - 1)/3 = (-i-1)/3.

  2. Path dependence of a non-analytic integrand. Integrate f(z)=zˉf(z) = \bar z from 00 to 1+i1+i.

    • Along the straight segment z=t(1+i)z = t(1+i), 0≤t≤10\le t\le 1: zˉ=t(1−i)\bar z = t(1-i), dz=(1+i) dtdz = (1+i)\,dt, and (1−i)(1+i)=2(1-i)(1+i) = 2, so the integral is ∫012t dt=1\int_0^1 2t\,dt = 1.
    • Along the two legs 0→1→1+i0 \to 1 \to 1+i: on the first, z=tz = t gives ∫01t dt=12\int_0^1 t\,dt = \tfrac12; on the second, z=1+itz = 1+it gives ∫01(1−it) i dt=i−i2 ⁣∫01t dt=i+12\int_0^1 (1-it)\,i\,dt = i - i^2\!\int_0^1 t\,dt = i + \tfrac12. The total is 1+i1 + i.

    The two values differ, as they must: zˉ\bar z satisfies no Cauchy-Riemann equation anywhere. Equivalently, around the closed triangle the integral is (1+i)−1=i≠0(1+i) - 1 = i \ne 0.

  3. The ML estimate, attained and applied. On the unit circle ∣zˉ∣=1\lvert \bar z\rvert = 1 and L=2πL = 2\pi, so ∣∮zˉ dz∣≤2π\lvert\oint \bar z\,dz\rvert \le 2\pi. The exact value is ∮∣z∣=1zˉ dz=∫02πe−itieit dt=2πi\oint_{\lvert z\rvert = 1}\bar z\,dz = \int_0^{2\pi} e^{-it}ie^{it}\,dt = 2\pi i, of modulus exactly 2π2\pi: here the bound is sharp. Used the other way, on the circle ∣z∣=R\lvert z \rvert = R with R>1R>1 we have ∣z4+1∣≥R4−1\lvert z^4+1\rvert \ge R^4 - 1 by the triangle inequality, so

    ∣∮∣z∣=Rdzz4+1∣≤2πRR4−1,\Bigl\lvert \oint_{\lvert z\rvert = R}\frac{dz}{z^4+1}\Bigr\rvert \le \frac{2\pi R}{R^4-1},

    which is 0.2356190.235619 at R=3R = 3 and 0.0062840.006284 at R=10R = 10, and tends to 00 as R→∞R \to \infty. Estimates of exactly this shape are what let a large circular arc be discarded when a real integral is computed by contour methods (infinite integrals).

ApplicationOriginal

In two-dimensional ideal flow with complex potential w(z)w(z) and velocity field v=(u1,u2)\mathbf{v} = (u_1,u_2), the combination w′(z)‾=u1+iu2\overline{w'(z)} = u_1 + iu_2 packages the velocity, and for a closed contour γ\gamma

∮γw′(z) dz=∮γ(u1 dx+u2 dy)⏟circulation+  i∮γ(u1 dy−u2 dx)⏟flux,\oint_\gamma w'(z)\,dz = \underbrace{\oint_\gamma (u_1\,dx + u_2\,dy)}_{\text{circulation}} + \; i \underbrace{\oint_\gamma (u_1\,dy - u_2\,dx)}_{\text{flux}} ,

so one complex integral delivers both the circulation around a wing section and the net flux through the curve. The circulation is what the Kutta-Joukowski theorem converts into lift. For a point vortex at the origin, whose complex potential is a multiple of ln⁡z\ln z and whose complex velocity is accordingly a multiple of 1/z1/z, part 5 evaluates the circulation around a circle centred on the vortex and shows it to be nonzero and independent of the radius. That the same value is obtained around every contour enclosing the vortex is a further step, needing either Cauchy's theorem in its deformation form or the winding number; part 5 on its own does not give it.

The same splitting appears in electromagnetism, where ∮E⋅dl\oint \mathbf{E}\cdot d\mathbf{l} and the flux of E\mathbf{E} across a curve are the real and imaginary parts of one contour integral of the complex field, and in the calculation of forces on conductors. Beyond that, the contour integral is the workhorse of the whole chapter: Cauchy's integral formula represents a function by one, the residue theorem evaluates them, and the ML estimate is the tool that makes the limiting arguments work.

Theorem

Cauchy's theorem

Adapted

The integral of an analytic function around a closed contour in a simply connected region is zero; consequently a contour may be deformed freely within a region of analyticity, and in a multiply connected region the integral reduces to one loop around each hole.

Prerequisites: Contour integrals, Analyticity and the Cauchy-Riemann equations, Connectedness and simple connectivity, Stokes' theorem and Green's theorem in the plane, Compactness

Definition(Simply connected regions and orientation)AdaptedAdapted from Arfken, Weber & Harris (2013), §11.3, p. 479; Riley, Hobson & Bence (2006), §24.9, p. 850

A region here is a non-empty connected open subset of C\mathbb{C}. Call Ω\Omega simply connected when each closed curve lying in it can be contracted, continuously and without ever leaving Ω\Omega, down to a single point of Ω\Omega — informally, when Ω\Omega has no holes. (For the precise notion of such a shrinking, see connectedness and homotopy.) A region that is not simply connected is multiply connected; the punctured disk 0<∣z∣<10 < \lvert z\rvert < 1 and any annulus are the standard examples.

A simple closed contour is positively oriented when the bounded region it encloses lies to its left. For a contour encircling an interior point this is the counterclockwise sense, that of increasing polar angle.

Theorem(Cauchy-Goursat theorem)AdaptedAdapted from Arfken, Weber & Harris (2013), §11.3, Eq. (11.18), p. 478; Ahlfors (1979), Ch. 4, §1.4, Thm 2, p. 109 and §1.5, Thm 4, p. 113

Let ff be analytic on a simply connected region Ω\Omega and let γ\gamma be any closed path in Ω\Omega. Then

∮γf(z) dz=0.\oint_\gamma f(z)\,dz = 0 .

Two special cases carry the proof and are worth stating separately.

  1. Rectangle (Goursat). If RR is a closed rectangle with sides parallel to the axes and ff is analytic on an open set containing RR, then ∮∂Rf(z) dz=0\oint_{\partial R} f(z)\,dz = 0.
  2. Disk. If ff is analytic on an open disk Δ\Delta, then ∮γf(z) dz=0\oint_\gamma f(z)\,dz = 0 for every closed path γ\gamma in Δ\Delta; equivalently, ff has a primitive on Δ\Delta.

No hypothesis beyond analyticity of ff is needed. Arfken, Weber and Harris (Arfken, Weber & Harris (2013), §11.3, pp. 481–482) and Riley, Hobson and Bence (Riley, Hobson & Bence (2006), §24.9, p. 849) both prove the theorem under the extra assumption that the first partial derivatives of uu and vv are continuous (equivalently that f′f' is continuous), and both record that the extra assumption is not necessary; part 1 is where it is removed.

Note(What is proved here, and what is quoted)AdaptedAdapted from Ahlfors (1979), Ch. 4, §4.2, Definition 1 and Thm 14, p. 139, and §4.4, Thm 15 with Corollary 1, p. 141; Arfken, Weber & Harris (2013), §11.3, p. 479

Two points of honesty about the statement above.

  1. Which definition of simple connectivity. We use the contractibility definition of Arfken, Weber and Harris (Arfken, Weber & Harris (2013), §11.3, p. 479). Ahlfors (Ahlfors (1979), Definition 1, p. 139) instead calls a region simply connected when its complement in the extended plane is connected, and then proves (his Theorem 14, p. 139) that this is equivalent to every closed curve in the region having winding number zero about every point outside it. For plane regions the three descriptions single out the same regions; we do not prove that equivalence here, and the winding number is where the third one is developed.

  2. Which cases we prove. The proofs below establish the rectangle and disk cases in full generality (no hypothesis beyond analyticity), and the simply connected case for a simple closed contour under the additional hypothesis that the first partials of uu and vv are continuous. The general statement — every closed path, every simply connected region, analyticity alone — is Ahlfors (Ahlfors (1979), Thm 15 and Corollary 1, p. 141)'s homology form of Cauchy's theorem, proved there by reducing an arbitrary cycle homologous to zero to a combination of rectangle boundaries. That reduction is a topological argument beyond the scope of this unit, and we quote the result rather than prove it. Nothing later in this chapter needs more than the disk case plus the deformation principle.

Theorem(Deformation of contours and multiply connected regions)AdaptedAdapted from Arfken, Weber & Harris (2013), §11.3, Eqs. (11.26)–(11.28), pp. 483–484; Riley, Hobson & Bence (2006), §24.9, Eqs. (24.43)–(24.45), pp. 850–851

Let ff be analytic on a region Ω\Omega.

  1. Two paths. If γ1\gamma_1 and γ2\gamma_2 are paths in Ω\Omega with the same initial and the same final point, and the region enclosed between them lies in Ω\Omega, then ∫γ1f dz=∫γ2f dz\int_{\gamma_1}f\,dz = \int_{\gamma_2}f\,dz. "The region enclosed between them" is the wording of Riley, Hobson and Bence (Riley, Hobson & Bence (2006), §24.9, p. 850), and it names nothing definite unless γ1\gamma_1 followed by the reverse of γ2\gamma_2 is a simple closed curve. We read the hypothesis in the form the proof actually uses: γ1\gamma_1 followed by the reverse of γ2\gamma_2 can be contracted to a point inside Ω\Omega, which is what sweeping across an enclosed region in Ω\Omega accomplishes when there is one.

  2. Deformation. If CC and σ\sigma are positively oriented simple closed contours with σ\sigma inside CC, and ff is analytic on the closed region between them, then

    ∮Cf(z) dz=∮σf(z) dz.\oint_C f(z)\,dz = \oint_\sigma f(z)\,dz .

    More generally, the value of ∮γf dz\oint_\gamma f\,dz is unchanged by any continuous deformation of γ\gamma within the region of analyticity, which precludes moving the contour across a point where ff fails to be analytic.

  3. Several holes. If CC is a positively oriented simple closed contour and σ1,…,σm\sigma_1,\dots,\sigma_m are positively oriented simple closed contours inside CC, mutually exterior to one another, and ff is analytic on the closed region inside CC and outside all the σk\sigma_k, then

    ∮Cf(z) dz=∑k=1m∮σkf(z) dz.\oint_C f(z)\,dz = \sum_{k=1}^{m}\oint_{\sigma_k} f(z)\,dz .
Idea of the proofOriginal

Three ideas, in increasing subtlety.

For the simply connected statement under the extra hypothesis that f′f' is continuous, there is nothing to invent: f dzf\,dz is a pair of real one-forms, Green's theorem in the plane converts the loop integral into an area integral, and the Cauchy-Riemann equations say the two integrands vanish identically.

Goursat's proof of part 1 dispenses with that hypothesis by using differentiability where it is strongest, namely at a single point. Bisect the rectangle repeatedly, always keeping a quarter whose integral is at least a quarter of the current one; the rectangles shrink to a point z∗z_*. Near z∗z_* the function agrees with its own affine approximation to within ε∣z−z∗∣\varepsilon\lvert z - z_*\rvert, for any ε\varepsilon we fix in advance and all zz close enough; and the affine part integrates to zero around any closed contour because 11 and zz have primitives. The error contributes at most ε⋅(diameter)⋅(perimeter)\varepsilon \cdot (\text{diameter})\cdot(\text{perimeter}), and under bisection diameter and perimeter each halve, so this bound decays like 4−n4^{-n} — at exactly the rate at which the bisection inflates the integral. The two rates cancel and leave ∣∮∣≤ε dL\lvert\oint\rvert \le \varepsilon\,dL for every ε\varepsilon.

Part 2 then upgrades rectangles to all closed paths in a disk by building a primitive: integrate along a horizontal-then-vertical staircase from the centre. Part 1 says the vertical-then-horizontal staircase gives the same value, and comparing the two ways of reaching a nearby point yields the two partial derivatives of the primitive; the fundamental theorem for primitives finishes it. The deformation results are a bookkeeping trick: cut the region between the two contours with a thin bridge to make it simply connected, apply the theorem, and let the two edges of the bridge cancel.

Proof(Simply connected regions, with f′f' continuous)AdaptedAdapted from Arfken, Weber & Harris (2013), §11.3, Eqs. (11.21)–(11.25), pp. 481–482; Riley, Hobson & Bence (2006), §24.9, Eqs. (24.41)–(24.42), p. 849

Assume for this proof that uu and vv have continuous first partial derivatives, and let γ\gamma be a positively oriented simple closed contour in Ω\Omega bounding the region AA. We need A⊆ΩA \subseteq \Omega. That holds because Ω\Omega is simply connected, by a fact of plane topology — a simply connected plane region contains the inside of every simple closed curve lying in it, equivalently a region with a hole has closed curves it cannot fill in — which we use here without proof; it belongs to the circle of ideas of connectedness and simple connectivity and of the equivalences recorded in the note above. Readers who prefer no appeal to topology may take A⊆ΩA \subseteq \Omega as a hypothesis: it is satisfied in every application below.

Step 1 (split into real line integrals). With f=u+ivf = u+iv and dz=dx+i dydz = dx + i\,dy, as in the definition of a contour integral,

∮γf dz=∮γ(u dx−v dy)+i∮γ(v dx+u dy).\oint_\gamma f\,dz = \oint_\gamma (u\,dx - v\,dy) + i\oint_\gamma (v\,dx + u\,dy) .

Step 2 (Green's theorem in the plane). Both are real line integrals of the form ∮γ(P dx+Q dy)\oint_\gamma (P\,dx + Q\,dy) with P,QP,Q having continuous first partials on AA and its boundary, so Stokes' theorem in the plane applies to each:

∮γ(P dx+Q dy)=∬A(∂Q∂x−∂P∂y) dx dy.\oint_\gamma (P\,dx + Q\,dy) = \iint_A \Bigl(\frac{\partial Q}{\partial x} - \frac{\partial P}{\partial y}\Bigr)\,dx\,dy .

Step 3 (the integrands vanish). For the first integral take P=uP = u, Q=−vQ = -v; for the second, P=vP = v, Q=uQ = u. The two area integrands are

−vx−uyandux−vy,-v_x - u_y \qquad\text{and}\qquad u_x - v_y ,

and both are identically zero on AA by the Cauchy-Riemann equations, which hold at every point of Ω\Omega because ff is analytic there. Hence ∮γf dz=0\oint_\gamma f\,dz = 0.

Steps 1 to 3 prove the theorem for a simple closed contour in Ω\Omega, which is the case needed in everything that follows. □\square

Proof(Part 1: Goursat's theorem for a rectangle)AdaptedAdapted from Ahlfors (1979), Ch. 4, §1.4, Thm 2, pp. 109–111

No continuity of f′f' is assumed. For a rectangle SS with sides parallel to the axes write η(S)=∮∂Sf dz\eta(S) = \oint_{\partial S}f\,dz, the boundary positively oriented.

Step 1 (bisection does not lose everything). Divide SS into four congruent rectangles S(1),…,S(4)S^{(1)},\dots,S^{(4)} by halving each side. Each interior edge of the subdivision occurs in exactly two of the four boundaries, with opposite orientations, so those contributions cancel and

η(S)=η(S(1))+η(S(2))+η(S(3))+η(S(4)).\eta(S) = \eta(S^{(1)}) + \eta(S^{(2)}) + \eta(S^{(3)}) + \eta(S^{(4)}) .

By the triangle inequality at least one kk has ∣η(S(k))∣≥14∣η(S)∣\lvert \eta(S^{(k)})\rvert \ge \tfrac14\lvert\eta(S)\rvert; choose one such (say the first in a fixed ordering) and call it R1R_1.

Step 2 (nested rectangles). Repeating, we get R⊃R1⊃R2⊃⋯R \supset R_1 \supset R_2 \supset \cdots with

∣η(Rn)∣≥4−n∣η(R)∣,\lvert \eta(R_n)\rvert \ge 4^{-n}\lvert\eta(R)\rvert ,

and with diameter dn=2−ndd_n = 2^{-n}d and perimeter Ln=2−nLL_n = 2^{-n}L, where dd and LL are those of RR. The RnR_n are nested non-empty compact sets with diameters tending to 00, so by compactness they have exactly one common point z∗z_*, and every neighbourhood of z∗z_* contains all but finitely many RnR_n.

Step 3 (use differentiability at the one point z∗z_*). Let ε>0\varepsilon > 0. Since ff is analytic on an open set containing RR, it is differentiable at z∗z_*, so there is δ>0\delta > 0 with

∣f(z)−f(z∗)−(z−z∗)f′(z∗)∣≤ε ∣z−z∗∣whenever ∣z−z∗∣<δ.\bigl\lvert f(z) - f(z_*) - (z-z_*)f'(z_*)\bigr\rvert \le \varepsilon\,\lvert z - z_*\rvert \qquad\text{whenever } \lvert z - z_*\rvert < \delta .

Fix nn with Rn⊂{∣z−z∗∣<δ}R_n \subset \{\lvert z - z_*\rvert < \delta\}.

Step 4 (the affine part integrates to zero). The functions 11 and zz have the primitives zz and z2/2z^2/2, entire; so by part 5 of the properties of contour integrals, ∮∂Rndz=0\oint_{\partial R_n} dz = 0 and ∮∂Rnz dz=0\oint_{\partial R_n} z\,dz = 0. Therefore the constant f(z∗)−z∗f′(z∗)f(z_*) - z_*f'(z_*) and the linear term zf′(z∗)zf'(z_*) contribute nothing and

η(Rn)=∮∂Rn[f(z)−f(z∗)−(z−z∗)f′(z∗)]dz.\eta(R_n) = \oint_{\partial R_n}\bigl[f(z) - f(z_*) - (z - z_*)f'(z_*)\bigr]dz .

Step 5 (estimate and compare the two rates). By the ML estimate in its arc-length form, and because ∣z−z∗∣≤dn\lvert z - z_*\rvert \le d_n on ∂Rn\partial R_n (both points lie in RnR_n),

∣η(Rn)∣≤ε∮∂Rn∣z−z∗∣ ∣dz∣≤ε dnLn=ε 4−ndL.\lvert\eta(R_n)\rvert \le \varepsilon \oint_{\partial R_n} \lvert z - z_*\rvert\,\lvert dz\rvert \le \varepsilon\, d_n L_n = \varepsilon\,4^{-n} d L .

Combining with Step 2, 4−n∣η(R)∣≤ε 4−ndL4^{-n}\lvert\eta(R)\rvert \le \varepsilon\,4^{-n}dL, so ∣η(R)∣≤ε dL\lvert \eta(R)\rvert \le \varepsilon\,dL. As ε>0\varepsilon>0 was arbitrary and d,Ld, L are fixed, η(R)=0\eta(R) = 0. □\square

Proof(Part 2: from a rectangle to a disk)AdaptedAdapted from Ahlfors (1979), Ch. 4, §1.5, Thm 4, p. 113

Let ff be analytic on the open disk Δ\Delta with centre z0=x0+iy0z_0 = x_0 + iy_0, and define, for z=x+iy∈Δz = x+iy \in \Delta,

F(z)=∫σzf(ζ) dζ,F(z) = \int_{\sigma_z} f(\zeta)\,d\zeta ,

where σz\sigma_z is the horizontal segment from (x0,y0)(x_0,y_0) to (x,y0)(x,y_0) followed by the vertical segment from (x,y0)(x,y_0) to (x,y)(x,y). Both segments lie in Δ\Delta, because a disk contains the whole axis-parallel rectangle determined by its centre and any of its points.

Step 1 (∂F/∂y\partial F/\partial y). Only the second segment depends on yy, and on it ζ=x+is\zeta = x + is with dζ=i dsd\zeta = i\,ds, so F(z)=c(x)+i∫y0yf(x+is) dsF(z) = c(x) + i\int_{y_0}^{y}f(x+is)\,ds. Since ff is continuous, the fundamental theorem of calculus gives ∂F/∂y=if(z)\partial F/\partial y = if(z).

Step 2 (∂F/∂x\partial F/\partial x). Let σ~z\tilde\sigma_z be the other staircase: vertical from (x0,y0)(x_0,y_0) to (x0,y)(x_0,y), then horizontal to (x,y)(x,y). Together, σz\sigma_z and the reverse of σ~z\tilde\sigma_z traverse the boundary of the rectangle with opposite corners z0z_0 and zz, which lies in Δ\Delta; by part 1 that boundary integral is 00, so ∫σ~zf=∫σzf=F(z)\int_{\tilde\sigma_z}f = \int_{\sigma_z}f = F(z). Now only the final horizontal segment of σ~z\tilde\sigma_z depends on xx, and there ζ=s+iy\zeta = s + iy with dζ=dsd\zeta = ds, so F(z)=c~(y)+∫x0xf(s+iy) dsF(z) = \tilde c(y) + \int_{x_0}^{x}f(s+iy)\,ds and ∂F/∂x=f(z)\partial F/\partial x = f(z).

Step 3 (FF is analytic with F′=fF' = f). From Steps 1 and 2, Fx=f=−i(if)=−iFyF_x = f = -i(if) = -iF_y, which is the Cauchy-Riemann condition in the combined form of the Cauchy-Riemann theorem; and Fx=fF_x = f is continuous. Hence FF is analytic on Δ\Delta with F′=fF' = f.

Step 4 (conclusion). ff has a primitive on Δ\Delta, so by the fundamental theorem for contour integrals its integral over every closed path in Δ\Delta vanishes. □\square

Proof(Deformation and multiply connected regions)AdaptedAdapted from Arfken, Weber & Harris (2013), §11.3, Eqs. (11.26)–(11.28), pp. 483–484; Riley, Hobson & Bence (2006), §24.9, Eqs. (24.43)–(24.45), pp. 850–851; Ahlfors (1979), Ch. 4, §4.4, Thm 15 and Corollary 1, p. 141 (quoted in part 1)

Part 1 (two paths). Let γ\gamma be γ1\gamma_1 followed by the reverse of γ2\gamma_2; since the two paths share both endpoints, γ\gamma is a closed path in Ω\Omega. By the hypothesis as read in the statement, γ\gamma can be contracted to a point within Ω\Omega. A closed path contractible in Ω\Omega has winding number zero about every point of the complement of Ω\Omega, since the winding number is unchanged under a homotopy that never crosses the point and is 00 for a constant path; in Ahlfors (Ahlfors (1979), §4.3, Definition 2, p. 141)'s language γ\gamma is then homologous to zero in Ω\Omega, and his Theorem 15 — the homology form of Cauchy's theorem, which the note above records as quoted and not proved here — gives the conclusion. Both the winding-number fact and the theorem are taken on quotation; they are developed in the winding number. Hence ∮γf dz=0\oint_\gamma f\,dz = 0, that is ∫γ1f dz−∫γ2f dz=0\int_{\gamma_1}f\,dz - \int_{\gamma_2}f\,dz = 0 by the reversal and subdivision rules for contour integrals.

Note what is and is not assumed here: γ\gamma need not be simple and Ω\Omega need not be simply connected, but the price is that this part leans on the homology form of the theorem, which the note above says we quote rather than prove. When γ\gamma happens to be a simple closed curve whose inside lies in Ω\Omega, or to lie inside a disk contained in Ω\Omega, the conclusion follows instead from the Green's theorem proof or from part 2 of the Cauchy-Goursat statement, with nothing quoted.

Part 2 (deformation). Cut the annular region between CC and σ\sigma along two nearby arcs λ\lambda and λ′\lambda' joining a point of σ\sigma to a point of CC. Travelling along CC in the positive sense, then inward along λ′\lambda', then around σ\sigma in the negative sense, then outward along λ\lambda, produces a single closed contour Γ\Gamma which does not enter the hole; the region it encloses is simply connected and contained in the region of analyticity, so ∮Γf dz=0\oint_\Gamma f\,dz = 0. Now let λ\lambda and λ′\lambda' approach each other. Their contributions tend to the integrals of ff along the same arc in opposite directions, which cancel, because ff is continuous across the cut — unlike a branch cut, this barrier is of our own making. What is left is

∮Cf dz−∮σf dz=0.\oint_C f\,dz - \oint_\sigma f\,dz = 0 .

For the general deformation statement, a continuous deformation that stays in the region of analyticity can be performed in finitely many steps each of which replaces an arc by another arc with the same endpoints, the region between them lying in Ω\Omega; each such step changes nothing by part 1.

Part 3 (several holes). Insert one cut from CC to each σk\sigma_k. The same construction yields one closed contour enclosing a simply connected region of analyticity, and the cuts cancel in pairs, leaving ∮Cf dz−∑k∮σkf dz=0\oint_C f\,dz - \sum_k \oint_{\sigma_k} f\,dz = 0. □\square

RemarkOriginal

Where the hypotheses are used, and what fails without them.

  • Simple connectivity is indispensable. On the region Ω={0<∣z∣<2}\Omega = \{0 < \lvert z\rvert < 2\}, which is not simply connected, f(z)=1/zf(z) = 1/z is analytic, yet ∮∣z∣=1dz/z=2πi≠0\oint_{\lvert z\rvert = 1} dz/z = 2\pi i \ne 0. Every later appearance of a nonzero contour integral traces back to this example.
  • What the hypothesis asks for is a simply connected region of analyticity containing the contour — and it is enough that some such region exist, even a small one. The same 1/z1/z is analytic on the disk ∣z−1∣<1\lvert z - 1\rvert < 1, which is simply connected, so every closed path inside that disk integrates to zero. Conversely the circle ∣z∣=1\lvert z\rvert = 1 lies in no simply connected subregion of Ω\Omega on which 1/z1/z is analytic, which is why the theorem has nothing to say about it. In practice one therefore looks for the largest convenient simply connected subregion, rather than asking whether the region one started with is simply connected.
  • Analyticity on the enclosed region, not just on the contour. ∮∣z∣=1dz/z≠0\oint_{\lvert z\rvert=1}dz/z \ne 0 although 1/z1/z is perfectly well behaved on the circle itself. The theorem is a statement about what the contour surrounds: if ff is analytic on an open set containing a simple closed contour and its inside, the integral does vanish, and the 1/z1/z example escapes only because the enclosed point 00 is excluded.
  • Continuity of f′f' is used only in the Green's theorem proof, at Step 2, where Stokes' theorem in the plane needs continuously differentiable PP and QQ. Goursat's proof shows that the conclusion survives without it; this matters because the whole point of Cauchy's integral formula is to deduce smoothness of f′f' rather than assume it, and assuming it here would make the later argument circular. Both Arfken, Weber and Harris (Arfken, Weber & Harris (2013), §11.4, footnote 3, p. 488) and Riley, Hobson and Bence (Riley, Hobson & Bence (2006), §24.9, p. 850) flag the point and omit the stronger proof; we give it.
  • The cancelling bridges in the deformation proof need ff continuous across the bridge. That is why the same picture must not be used across a branch cut, where the two edges carry different values of the function and their contributions add instead of cancelling.
  • Orientation matters in part 3. With all contours positively oriented, the inner integrals are added; if the σk\sigma_k are traversed clockwise, as they are when one regards Γ\Gamma as a single boundary, the signs flip and the total is zero.
Example(Deforming a square to a circle)Original

Evaluate I=∮Sdz/zI = \oint_{S} dz/z, where SS is the boundary of the square with vertices ±2±2i\pm 2 \pm 2i, positively oriented.

  1. Choose the deformation. 1/z1/z is analytic on C∖{0}\mathbb{C}\setminus\{0\}, and the closed region between SS and the unit circle ∣z∣=1\lvert z\rvert = 1 contains no singularity, so part 2 applies.
  2. Evaluate on the easy contour. ∮∣z∣=1dz/z=2πi\oint_{\lvert z\rvert = 1}dz/z = 2\pi i by part 5 of the properties of contour integrals.
  3. Conclusion. I=2πi≈6.283185 iI = 2\pi i \approx 6.283185\,i. Direct numerical integration along the four sides of the square gives 6.283185 i6.283185\,i to six digits, which is the check.
  4. Contrast. For the square S′S' with vertices 3±i3\pm i and 5±i5 \pm i, which does not enclose the origin, 1/z1/z is analytic on a simply connected region containing S′S' and its interior, so ∮S′dz/z=0\oint_{S'}dz/z = 0.
Example(A contour with two holes)Original

Let g(z)=1z2−1=12(1z−1−1z+1)g(z) = \dfrac{1}{z^2 - 1} = \dfrac{1}{2}\Bigl(\dfrac{1}{z-1} - \dfrac{1}{z+1}\Bigr) and let CC be the circle ∣z∣=3\lvert z \rvert = 3, positively oriented.

  1. Isolate the holes. gg is analytic except at z=±1z = \pm1. Take σ±\sigma_{\pm} to be the positively oriented circles of radius 12\tfrac12 about ±1\pm 1; they lie inside CC and are mutually exterior, and gg is analytic on the closed region inside CC and outside both.
  2. Apply part 3. ∮Cg=∮σ+g+∮σ−g\oint_C g = \oint_{\sigma_+} g + \oint_{\sigma_-} g.
  3. Evaluate each. On σ+\sigma_+, the term −12/(z+1)-\tfrac12/(z+1) is analytic inside, so by Cauchy-Goursat it integrates to zero, and ∮σ+g=12∮σ+dzz−1=12⋅2πi=πi\oint_{\sigma_+}g = \tfrac12\oint_{\sigma_+}\frac{dz}{z-1} = \tfrac12 \cdot 2\pi i = \pi i. Symmetrically ∮σ−g=−12⋅2πi=−πi\oint_{\sigma_-}g = -\tfrac12\cdot 2\pi i = -\pi i.
  4. Total. ∮Cg dz=πi−πi=0\oint_C g\,dz = \pi i - \pi i = 0.
  5. Check by the ML estimate. On ∣z∣=R\lvert z\rvert = R we have ∣z2−1∣≥R2−1\lvert z^2-1\rvert \ge R^2 - 1, so ∣∮∣z∣=Rg∣≤2πR/(R2−1)\lvert\oint_{\lvert z\rvert = R} g\rvert \le 2\pi R/(R^2-1), which tends to 00 as R→∞R \to \infty; by part 2 the integral is the same for every R>1R>1, so it must be 00 for all of them — consistent with the computation, and an instance of the standard way to show that a sum of residues vanishes.
ApplicationOriginal

Cauchy's theorem is the reason so many physical quantities are path independent, and the deformation principle is the reason the ones that are not depend only on what is enclosed.

  • Vortices and circulation. In two-dimensional ideal flow the complex velocity of a point vortex of strength κ\kappa at the origin is κ/(2πiz)\kappa/(2\pi i z). The circulation around any loop enclosing the vortex is κ\kappa, by deformation to a small circle, and zero around any loop that does not: the "hole" in the region of analyticity is the vortex core. The same statement in electromagnetism is Ampère's law for a thin wire, and in superfluids and type-II superconductors it is the quantization of circulation and of magnetic flux.
  • Rotating the integration path. Because the integral over a closed contour vanishes, a real integral along the axis may be replaced by one along any deformed path with the same endpoints, provided no singularity is crossed. This is how steepest-descent and saddle-point evaluations are justified (the method of steepest descent and stationary phase), how a Fourier integral is turned into a convergent one by moving into the upper half-plane, and how the Matsubara and Wick rotations of quantum field theory are carried out.
  • Causality and dispersion. A response function analytic in a half-plane has, by Cauchy's theorem, a vanishing integral around a contour closing in that half-plane; moving the contour onto the real axis converts this into a relation between the real and imaginary parts of the response, which is the Kramers-Kronig relation (dispersion relations).

Theorem

Cauchy's integral formula

Adapted

An analytic function and all of its derivatives at an interior point are given by integrals of the function over the surrounding contour; consequently an analytic function is infinitely differentiable, and conversely a continuous function all of whose closed contour integrals vanish is analytic.

Prerequisites: Cauchy's theorem

Theorem(Cauchy's integral formula)AdaptedAdapted from Riley, Hobson & Bence (2006), §24.10, Eqs. (24.46)–(24.48), pp. 851–853; Ahlfors (1979), Ch. 4, §2.2, Eq. (22), p. 119 and §2.3, Eq. (24), p. 120; Arfken, Weber & Harris (2013), §11.4, Eqs. (11.30)–(11.33), pp. 486–488

Let CC be a positively oriented simple closed contour, let ff be analytic on and inside CC, and let z0z_0 be a point inside CC. Then

f(z0)=12πi∮Cf(z)z−z0 dz,f(z_0) = \frac{1}{2\pi i}\oint_C \frac{f(z)}{z - z_0}\,dz ,

and for every n∈Nn \in \mathbb{N} the derivative f(n)(z0)f^{(n)}(z_0) exists and

f(n)(z0)=n!2πi∮Cf(z)(z−z0)n+1 dz.f^{(n)}(z_0) = \frac{n!}{2\pi i}\oint_C \frac{f(z)}{(z-z_0)^{n+1}}\,dz .

Consequences.

  1. Infinite differentiability. A function analytic on an open set has derivatives of all orders there, and every one of them is again analytic.

  2. Point outside. If instead z0z_0 lies outside CC, the whole integrand is analytic on and inside CC, so ∮Cf(z)(z−z0)−1dz=0\oint_C f(z)(z-z_0)^{-1}dz = 0.

  3. Mean value property. If the closed disk ∣z−z0∣≤r\lvert z - z_0\rvert \le r lies in the region of analyticity, then

    f(z0)=12π∫02πf(z0+reiθ) dθ.f(z_0) = \frac{1}{2\pi}\int_0^{2\pi} f\bigl(z_0 + re^{i\theta}\bigr)\,d\theta .
  4. Morera's theorem (converse of Cauchy's theorem). If ff is continuous on a region Ω\Omega and ∮γf(z) dz=0\oint_\gamma f(z)\,dz = 0 for every closed contour γ\gamma in Ω\Omega, then ff is analytic on Ω\Omega.

Idea of the proofOriginal

The formula is Cauchy's theorem applied to the right function. The quotient (f(z)−f(z0))/(z−z0)\bigl(f(z) - f(z_0)\bigr)/(z-z_0) is analytic except at z0z_0, where it stays bounded, so its integral around CC vanishes; splitting the numerator leaves ∮f(z)/(z−z0) dz\oint f(z)/(z-z_0)\,dz on one side and f(z0)∮dz/(z−z0)=2πif(z0)f(z_0)\oint dz/(z-z_0) = 2\pi i f(z_0) on the other. Equivalently, and this is the version easiest to compute with: deform CC to a tiny circle about z0z_0, on which ff is nearly the constant f(z0)f(z_0), and the integral is then f(z0)f(z_0) times 2πi2\pi i.

For the derivatives the temptation is to differentiate under the integral sign, which gives the right answer but needs justification. The honest route is to prove once and for all that an integral of the form ∮φ(ζ)(ζ−z)−ndζ\oint \varphi(\zeta)(\zeta - z)^{-n}d\zeta, with φ\varphi merely continuous on the contour, is a differentiable function of zz off the contour, with the derivative one gets by differentiating the kernel. Each step of the induction only needs the previous step applied to a slightly different φ\varphi. The payoff is disproportionate: since the right-hand side makes sense for all nn, so does the left, and analyticity at once implies infinite differentiability — a conclusion with no counterpart for functions of a real variable.

Proof(The formula)AdaptedAdapted from Ahlfors (1979), Ch. 4, §2.2, pp. 118–119; Riley, Hobson & Bence (2006), §24.10, p. 852

Step 1 (reduce to a small circle). Let σρ\sigma_\rho be the positively oriented circle ∣z−z0∣=ρ\lvert z - z_0\rvert = \rho, with ρ\rho small enough that the closed disk it bounds lies inside CC. The integrand f(z)/(z−z0)f(z)/(z-z_0) is analytic on the closed region between CC and σρ\sigma_\rho, because ff is analytic there and the only zero of the denominator is z0z_0, which is not in that region. By the deformation principle,

∮Cf(z)z−z0 dz=∮σρf(z)z−z0 dzfor every such ρ.\oint_C \frac{f(z)}{z-z_0}\,dz = \oint_{\sigma_\rho}\frac{f(z)}{z - z_0}\,dz \qquad\text{for every such } \rho .

Step 2 (evaluate on the circle). Parametrize σρ\sigma_\rho by z=z0+ρeiθz = z_0 + \rho e^{i\theta}, 0≤θ≤2π0 \le \theta \le 2\pi, so dz=iρeiθdθdz = i\rho e^{i\theta}d\theta and z−z0=ρeiθz - z_0 = \rho e^{i\theta}. The factors of ρeiθ\rho e^{i\theta} cancel:

∮σρf(z)z−z0 dz=i∫02πf(z0+ρeiθ) dθ.\oint_{\sigma_\rho}\frac{f(z)}{z-z_0}\,dz = i\int_0^{2\pi} f\bigl(z_0 + \rho e^{i\theta}\bigr)\,d\theta .

Dividing by 2πi2\pi i gives the mean value of ff over the circle, which is consequence 3 as soon as the formula itself is proved.

Step 3 (let the circle shrink). Subtract the constant: using ∫02πdθ=2π\int_0^{2\pi}d\theta = 2\pi,

∣12π∫02πf(z0+ρeiθ) dθ−f(z0)∣=∣12π∫02π[f(z0+ρeiθ)−f(z0)]dθ∣≤max⁡∣z−z0∣=ρ∣f(z)−f(z0)∣,\Bigl\lvert \frac{1}{2\pi}\int_0^{2\pi}f(z_0+\rho e^{i\theta})\,d\theta - f(z_0)\Bigr\rvert = \Bigl\lvert \frac{1}{2\pi}\int_0^{2\pi}\bigl[f(z_0 + \rho e^{i\theta}) - f(z_0)\bigr]d\theta\Bigr\rvert \le \max_{\lvert z - z_0\rvert = \rho}\lvert f(z) - f(z_0)\rvert ,

by the modulus inequality for integrals. Since ff is continuous at z0z_0 — it is differentiable there — the right-hand side tends to 00 as ρ→0+\rho \to 0^+. But by Steps 1 and 2 the left-hand side is ∣12πi∮Cf(z)(z−z0)−1dz−f(z0)∣\bigl\lvert \frac{1}{2\pi i}\oint_C f(z)(z-z_0)^{-1}dz - f(z_0)\bigr\rvert, a quantity independent of ρ\rho. A constant that is smaller than every positive number is zero, which is the formula. □\square

Proof(The derivative formulas, and infinite differentiability)AdaptedAdapted from Ahlfors (1979), Ch. 4, §2.3, Lemma 3 and its induction, pp. 120–122; Riley, Hobson & Bence (2006), §24.10, Eq. (24.47), p. 852

Lemma. Let γ\gamma be a path and let φ\varphi be continuous on γ\gamma. For n∈Nn \in \mathbb{N} and zz not on γ\gamma put

Fn(z)=∫γφ(ζ)(ζ−z)n dζ.F_n(z) = \int_\gamma \frac{\varphi(\zeta)}{(\zeta - z)^n}\,d\zeta .

Then every FnF_n is analytic off γ\gamma, and Fn′=n Fn+1F_n' = n\,F_{n+1}.

Set-up for the proof of the lemma. Fix z1z_1 off γ\gamma and pick δ>0\delta > 0 with ∣ζ−z1∣≥δ\lvert \zeta - z_1\rvert \ge \delta for all ζ∈γ\zeta \in \gamma, which is possible because γ\gamma is compact. Throughout, zz is restricted to ∣z−z1∣<δ/2\lvert z - z_1\rvert < \delta/2, so that ∣ζ−z∣>δ/2\lvert \zeta - z\rvert > \delta/2 on γ\gamma. Write K=∫γ∣φ∣ ∣dζ∣K = \int_\gamma \lvert\varphi\rvert\,\lvert d\zeta\rvert, a finite constant. All estimates below are the ML estimate of contour integrals.

Step 1 (F1F_1 is continuous at z1z_1). From the identity

1ζ−z−1ζ−z1=z−z1(ζ−z)(ζ−z1)(∗)\frac{1}{\zeta - z} - \frac{1}{\zeta - z_1} = \frac{z - z_1}{(\zeta - z)(\zeta - z_1)} \tag{$\ast$}

we get

F1(z)−F1(z1)=(z−z1)∫γφ(ζ) dζ(ζ−z)(ζ−z1),∣F1(z)−F1(z1)∣≤∣z−z1∣ 2δ⋅1δ⋅K=2Kδ2∣z−z1∣,F_1(z) - F_1(z_1) = (z-z_1)\int_\gamma \frac{\varphi(\zeta)\,d\zeta}{(\zeta - z)(\zeta - z_1)}, \qquad \lvert F_1(z) - F_1(z_1)\rvert \le \lvert z - z_1\rvert\,\frac{2}{\delta}\cdot\frac{1}{\delta}\cdot K = \frac{2K}{\delta^2}\lvert z - z_1\rvert ,

which tends to 00 as z→z1z \to z_1.

Step 2 (F1′=F2F_1' = F_2). Divide the first display of Step 1 by z−z1z - z_1:

F1(z)−F1(z1)z−z1=∫γψ(ζ)ζ−z dζ,ψ(ζ)=φ(ζ)ζ−z1,\frac{F_1(z) - F_1(z_1)}{z-z_1} = \int_\gamma \frac{\psi(\zeta)}{\zeta - z}\,d\zeta, \qquad \psi(\zeta) = \frac{\varphi(\zeta)}{\zeta - z_1},

and note that ψ\psi is continuous on γ\gamma, since ζ≠z1\zeta \ne z_1 there. The right-hand side is therefore the function "F1F_1 built from ψ\psi", which by Step 1 is continuous at z1z_1; its value at z=z1z = z_1 is ∫γφ(ζ)(ζ−z1)−2dζ=F2(z1)\int_\gamma \varphi(\zeta)(\zeta-z_1)^{-2}d\zeta = F_2(z_1). Hence F1′(z1)=F2(z1)F_1'(z_1) = F_2(z_1).

Step 3 (the induction identity). Let n≥2n \ge 2 and keep ψ(ζ)=φ(ζ)/(ζ−z1)\psi(\zeta) = \varphi(\zeta)/(\zeta - z_1). Denote by GkG_k the integrals built from ψ\psi, that is Gk(z)=∫γψ(ζ)(ζ−z)−kdζG_k(z) = \int_\gamma \psi(\zeta)(\zeta - z)^{-k}d\zeta. Dividing (∗)(\ast) by (ζ−z)n−1(\zeta - z)^{n-1} and rearranging gives

1(ζ−z)n=1(ζ−z)n−1(ζ−z1)+z−z1(ζ−z)n(ζ−z1),\frac{1}{(\zeta-z)^n} = \frac{1}{(\zeta - z)^{n-1}(\zeta - z_1)} + \frac{z - z_1}{(\zeta-z)^n(\zeta - z_1)} ,

and integrating against φ\varphi and subtracting Fn(z1)F_n(z_1),

Fn(z)−Fn(z1)=[Gn−1(z)−Gn−1(z1)]+(z−z1) Gn(z).(†)F_n(z) - F_n(z_1) = \bigl[G_{n-1}(z) - G_{n-1}(z_1)\bigr] + (z - z_1)\,G_n(z) . \tag{$\dagger$}

Here we used Gn−1(z1)=∫γφ(ζ)(ζ−z1)−ndζ=Fn(z1)G_{n-1}(z_1) = \int_\gamma \varphi(\zeta)(\zeta - z_1)^{-n}d\zeta = F_n(z_1) and, for the last term, ∫γφ(ζ)(ζ−z)−n(ζ−z1)−1dζ=Gn(z)\int_\gamma \varphi(\zeta)(\zeta-z)^{-n}(\zeta-z_1)^{-1}d\zeta = G_n(z).

Step 4 (induction). Suppose, as the induction hypothesis, that Fn−1F_{n-1} is continuous and differentiable at z1z_1 with Fn−1′=(n−1)FnF_{n-1}' = (n-1)F_n, for every continuous density on γ\gamma — in particular for ψ\psi, so Gn−1′(z1)=(n−1)Gn(z1)G_{n-1}'(z_1) = (n-1)G_n(z_1). Steps 1 and 2 are this statement for n−1=1n-1=1.

  • Continuity of FnF_n at z1z_1. In (†)(\dagger) the bracket tends to 00 because Gn−1G_{n-1} is continuous at z1z_1, and GnG_n is bounded on ∣z−z1∣<δ/2\lvert z - z_1\rvert < \delta/2, since ∣ψ∣≤∣φ∣/δ\lvert \psi\rvert \le \lvert\varphi\rvert/\delta on γ\gamma gives ∣Gn(z)∣≤(2/δ)nK/δ\lvert G_n(z)\rvert \le (2/\delta)^n K/\delta independently of zz; so the second term tends to 00 as well. Since the density in (†)(\dagger) was arbitrary, this proves continuity of FnF_n at z1z_1 for every continuous density, hence also continuity of GnG_n at z1z_1.

  • Differentiability. Divide (†)(\dagger) by z−z1z - z_1 and let z→z1z \to z_1. The first quotient tends to Gn−1′(z1)=(n−1)Gn(z1)G_{n-1}'(z_1) = (n-1)G_n(z_1) by the induction hypothesis; the second term is Gn(z)G_n(z), which tends to Gn(z1)G_n(z_1) by the continuity just established. Therefore

    Fn′(z1)=(n−1)Gn(z1)+Gn(z1)=n Gn(z1)=n∫γφ(ζ) dζ(ζ−z1)n+1=n Fn+1(z1),F_n'(z_1) = (n-1)G_n(z_1) + G_n(z_1) = n\,G_n(z_1) = n\int_\gamma \frac{\varphi(\zeta)\,d\zeta}{(\zeta - z_1)^{n+1}} = n\,F_{n+1}(z_1) ,

    where the third equality is the definition of GnG_n with ψ=φ/(ζ−z1)\psi = \varphi/(\zeta - z_1). This closes the induction and proves the lemma.

Application to ff. Let Ω\Omega be open, ff analytic on Ω\Omega, and a∈Ωa \in \Omega. Choose a positively oriented circle CaC_a about aa whose closed disk lies in Ω\Omega. For every zz inside CaC_a, the formula already proved gives

f(z)=12πi∮Caf(ζ)ζ−z dζ=12πiF1(z),f(z) = \frac{1}{2\pi i}\oint_{C_a}\frac{f(\zeta)}{\zeta - z}\,d\zeta = \frac{1}{2\pi i}F_1(z),

with the continuous density φ=f\varphi = f on CaC_a. The lemma then gives, by induction on nn,

f(n)(z)=n!2πi∮Caf(ζ)(ζ−z)n+1 dζ(n=0,1,2,… ),f^{(n)}(z) = \frac{n!}{2\pi i}\oint_{C_a}\frac{f(\zeta)}{(\zeta - z)^{n+1}}\,d\zeta \qquad (n = 0, 1, 2, \dots),

since the case n=0n=0 is the display above and, assuming the formula for nn, the lemma gives f(n+1)(z)=n!2πiFn+1′(z)=n!(n+1)2πiFn+2(z)=(n+1)!2πiFn+2(z)f^{(n+1)}(z) = \frac{n!}{2\pi i}F_{n+1}'(z) = \frac{n!(n+1)}{2\pi i}F_{n+2}(z) = \frac{(n+1)!}{2\pi i}F_{n+2}(z). In particular every f(n)f^{(n)} is differentiable, hence analytic, which is consequence 1.

Finally, for z0z_0 inside the original contour CC, the integrand f(ζ)(ζ−z0)−n−1f(\zeta)(\zeta - z_0)^{-n-1} is analytic on the closed region between CC and a small circle Cz0C_{z_0} about z0z_0, so the deformation principle replaces Cz0C_{z_0} by CC and the formula holds as stated. □\square

Proof(Morera's theorem)AdaptedAdapted from Ahlfors (1979), Ch. 4, §2.3, p. 122; Arfken, Weber & Harris (2013), §11.4, Eqs. (11.34)–(11.36), pp. 489–490

Let ff be continuous on the region Ω\Omega with all closed contour integrals zero. By the fundamental theorem for contour integrals (contour integrals, part 4), ff has a primitive: there is an analytic FF on Ω\Omega with F′=fF' = f. By consequence 1, applied to FF, the derivative F′=fF' = f is itself analytic on Ω\Omega. □\square

RemarkOriginal

Where the hypotheses are used, and why the order of results matters.

  • Analyticity inside, not only on, CC. The proof uses the deformation principle on the region between CC and a small circle, which requires ff to be analytic throughout. Dropping it destroys the conclusion: take f(z)=1/zf(z) = 1/z, which is analytic on the unit circle CC but not inside it, and z0=12z_0 = \tfrac12. Partial fractions give 1z(z−12)=2(1z−12−1z)\frac{1}{z(z-\frac12)} = 2\bigl(\frac{1}{z - \frac12} - \frac1z\bigr), so 12πi∮C1/zz−12 dz=2(1−1)=0\frac{1}{2\pi i}\oint_C\frac{1/z}{z-\frac12}\,dz = 2(1 - 1) = 0, whereas f(12)=2f(\tfrac12) = 2. Numerical integration confirms the value 00.
  • z0z_0 strictly inside. On the contour the integral is not even absolutely convergent as written; the boundary case is the subject of principal value integrals.
  • The lemma is proved with φ\varphi merely continuous, which is what allows the induction to be run on densities such as φ(ζ)/(ζ−z1)k\varphi(\zeta)/(\zeta - z_1)^k that are not of the original form. Differentiating under the integral sign, as Riley, Hobson and Bence (Riley, Hobson & Bence (2006), §24.10, p. 853) and Hassani (Hassani (2013), Prop. 10.5.1, p. 315) do, reaches the same formula and is the quickest derivation; it needs a theorem on differentiation under the integral sign that the lemma above supplies in exactly the case at hand.
  • Goursat's form of Cauchy's theorem is needed, not Cauchy's. Arfken, Weber and Harris (Arfken, Weber & Harris (2013), §11.4, footnote 3, p. 488) are explicit about this. If Cauchy's theorem had been proved only under the assumption that f′f' is continuous, then consequence 1 — which produces continuity of f′f' — would be circular. This is why Cauchy's theorem is proved by bisection there.
  • Morera needs all closed contours. It is not enough that the integral vanish over the contours of one family. Note that our sources state the theorem with different hypotheses on the region: Ahlfors (Ahlfors (1979), Ch. 4, §2.3, p. 122) asks only that Ω\Omega be a region, while Arfken, Weber and Harris (Arfken, Weber & Harris (2013), §11.4, p. 489) and Hassani (Hassani (2013), Thm 10.5.8, p. 319) require it to be simply connected. We follow Ahlfors, whose version is the stronger one and is what the proof above gives: the hypothesis is about every closed contour in Ω\Omega, so the primitive it produces is single-valued whatever the topology of Ω\Omega, and simple connectivity is never used. On a simply connected region the three statements coincide.
  • Nothing like consequence 1 is true for real functions. x↦x∣x∣x\mapsto x\lvert x\rvert is differentiable on R\mathbb{R} with derivative 2∣x∣2\lvert x\rvert, which is not differentiable at 00. One complex derivative forces infinitely many.
Example(Evaluating contour integrals with the formulas)Original
  1. Order zero. ∮∣z∣=2ezz−1 dz=2πi e1=2πie≈17.079468 i\displaystyle\oint_{\lvert z\rvert = 2}\frac{e^z}{z - 1}\,dz = 2\pi i\,e^1 = 2\pi i e \approx 17.079468\,i, since eze^z is entire and z0=1z_0 = 1 lies inside the circle. Numerical integration over the parametrized circle reproduces 17.079468 i17.079468\,i.
  2. Order one. ∮∣z∣=2ez(z−1)2 dz=2πi1!ddzez∣z=1=2πie\displaystyle\oint_{\lvert z\rvert = 2}\frac{e^z}{(z-1)^2}\,dz = \frac{2\pi i}{1!}\frac{d}{dz}e^z\Big|_{z=1} = 2\pi i e as well — the same value here because eze^z is its own derivative.
  3. Order two. ∮∣z∣=1cos⁡zz3 dz=2πi2!d2dz2cos⁡z∣z=0=πi (−cos⁡0)=−πi≈−3.141593 i\displaystyle\oint_{\lvert z\rvert = 1}\frac{\cos z}{z^3}\,dz = \frac{2\pi i}{2!}\frac{d^2}{dz^2}\cos z\Big|_{z=0} = \pi i\,(-\cos 0) = -\pi i \approx -3.141593\,i, again confirmed numerically.
  4. Point outside. ∮∣z∣=1ezz−3 dz=0\displaystyle\oint_{\lvert z\rvert = 1}\frac{e^z}{z-3}\,dz = 0, by consequence 2.
  5. The mean value property as a check. Take f(z)=1/(2−z)f(z) = 1/(2-z), analytic on ∣z∣≤1\lvert z\rvert \le 1. Then 12π∫02πdθ2−eiθ=f(0)=12\frac{1}{2\pi}\int_0^{2\pi}\frac{d\theta}{2 - e^{i\theta}} = f(0) = \tfrac12. Numerical quadrature of the left side gives 0.5000000.500000, with imaginary part zero to the same accuracy.
Example(A function reconstructed from its boundary values)Original

Let ff be analytic on ∣z∣≤1\lvert z\rvert \le 1 with boundary values f(eiθ)=e2iθf(e^{i\theta}) = e^{2i\theta}. The formula fixes ff on the whole disk. For ∣z0∣<1\lvert z_0\rvert < 1,

f(z0)=12πi∮∣z∣=1z2z−z0 dzf(z_0) = \frac{1}{2\pi i}\oint_{\lvert z\rvert = 1}\frac{z^2}{z - z_0}\,dz

because z2z^2 agrees with ff on the contour, and the right-hand side equals z02z_0^2 by the formula applied to the entire function z2z^2. So f(z)=z2f(z) = z^2: the only analytic function on the disk with those boundary values. Checking at z0=12+12iz_0 = \tfrac12 + \tfrac12 i: f(z0)=(12+12i)2=12if(z_0) = (\tfrac12+\tfrac12 i)^2 = \tfrac12 i, and the integral 12πi∫02πe2iθeiθ−z0ieiθ dθ\frac{1}{2\pi i}\int_0^{2\pi}\frac{e^{2i\theta}}{e^{i\theta} - z_0}ie^{i\theta}\,d\theta evaluates numerically to 0.000000+0.500000 i0.000000 + 0.500000\,i.

Note what happens if the boundary data are not the restriction of an analytic function. With f(eiθ)=e−iθf(e^{i\theta}) = e^{-i\theta}, that is zˉ\bar z on the circle, we may write zˉ=1/z\bar z = 1/z there, and partial fractions give

12πi∮∣z∣=11/zz−z0 dz=1z0⋅12πi∮(1z−z0−1z)dz=1−1z0=0\frac{1}{2\pi i}\oint_{\lvert z\rvert = 1}\frac{1/z}{z - z_0}\,dz = \frac{1}{z_0}\cdot\frac{1}{2\pi i}\oint\Bigl(\frac{1}{z - z_0} - \frac1z\Bigr)dz = \frac{1 - 1}{z_0} = 0

for every z0z_0 with 0<∣z0∣<10 < \lvert z_0\rvert < 1 (numerically 0.0000000.000000 at z0=12+12iz_0 = \tfrac12 + \tfrac12 i). So if some analytic gg on the closed disk had g=zˉg = \bar z on the circle, the formula would force g≡0g \equiv 0, contradicting ∣g∣=1\lvert g\rvert = 1 on the boundary. No analytic function has those boundary values — which is why the Cauchy integral, applied to arbitrary boundary data, produces an analytic function in the disk whose boundary values are generally not the data one started from.

ApplicationOriginal
  • Numerical differentiation without subtractive cancellation. The nnth derivative formula turns differentiation into integration. Evaluating f(n)(z0)=n!2πrn∫02πf(z0+reiθ)e−inθ dθf^{(n)}(z_0) = \frac{n!}{2\pi r^n}\int_0^{2\pi}f(z_0+re^{i\theta})e^{-in\theta}\,d\theta by the trapezoidal rule on mm equally spaced points converges geometrically in mm for analytic ff and avoids the catastrophic cancellation of finite differences; this is the basis of the Lyness-Moler method and of the "complex-step derivative" used in sensitivity analysis and automatic differentiation.
  • Functions of operators. Replacing 1/(z−z0)1/(z - z_0) by the resolvent (zI−A)−1(zI - A)^{-1} of a matrix or operator AA turns the formula into the definition f(A)=12πi∮Cf(z)(zI−A)−1dzf(A) = \frac{1}{2\pi i}\oint_C f(z)(zI - A)^{-1}dz, with CC enclosing the spectrum. This is how etAe^{tA}, A\sqrt{A} and spectral projections are constructed and computed (the matrix exponential, bounded operators), and it is the contour integral that underlies the spectral theory of unbounded operators in quantum mechanics.
  • Boundary data determine the interior. The second example is a template: a response function analytic in a half-plane, or a potential analytic in a region, is fixed everywhere by its values on the boundary. In electrostatics this is the uniqueness theorem for the Dirichlet problem (the Poisson integral formula), and in scattering theory it is why measuring an amplitude on the physical region constrains it off the real axis, which is what makes analytic continuation a physical tool and not just a formal one.

Theorem

Liouville's theorem and the maximum modulus principle

Adapted

Cauchy's estimates bound the derivatives of an analytic function by the maximum of its modulus on a circle; a bounded entire function is therefore constant, every non-constant polynomial has a root, and the modulus of a non-constant analytic function attains its maximum only on the boundary.

Prerequisites: Cauchy's integral formula, Compactness, Connectedness and simple connectivity

Theorem(Cauchy's estimates, Liouville, and the maximum modulus principle)AdaptedAdapted from Ahlfors (1979), Ch. 4, §2.3, Eq. (25) and p. 122, and §3.4, Thms 12 and 12', p. 134; Riley, Hobson & Bence (2006), §24.10, Eq. (24.49), p. 853; Arfken, Weber & Harris (2013), §11.4, Eq. (11.37), p. 490; Hassani (2013), §10.5, Props. 10.5.4–10.5.5 and Thm 10.5.6, pp. 317–318
  1. Cauchy's estimate. Let ff be analytic on a region containing the closed disk ∣z−a∣≤r\lvert z - a\rvert \le r and put M(r)=max⁡∣z−a∣=r∣f(z)∣M(r) = \max_{\lvert z - a\rvert = r}\lvert f(z)\rvert. Then

    ∣f(n)(a)∣≤n! M(r)rn,n=0,1,2,…\bigl\lvert f^{(n)}(a)\bigr\rvert \le \frac{n!\,M(r)}{r^{n}}, \qquad n = 0, 1, 2, \dots
  2. Liouville's theorem. A function analytic on the whole plane and bounded there is constant.

  3. Fundamental theorem of algebra. A polynomial PP of degree n≥1n \ge 1 with complex coefficients has a root; consequently it factors as P(z)=an∏k=1n(z−λk)P(z) = a_n\prod_{k=1}^{n}(z - \lambda_k) and has exactly nn roots counted with multiplicity.

  4. Maximum modulus principle. Let ff be analytic on a region (a connected open set) Ω\Omega.

    (a) Global form. If there is z0∈Ωz_0 \in \Omega with ∣f(z)∣≤∣f(z0)∣\lvert f(z)\rvert \le \lvert f(z_0)\rvert for every z∈Ωz \in \Omega, then ff is constant on Ω\Omega. Equivalently, if ff is non-constant then ∣f∣\lvert f\rvert has no maximum in Ω\Omega.

    (b) Local form. If ∣f∣\lvert f\rvert has a local maximum at z0∈Ωz_0 \in \Omega — that is, ∣f(z)∣≤∣f(z0)∣\lvert f(z)\rvert \le \lvert f(z_0)\rvert for all zz in some disk about z0z_0 — then ff is constant on a disk about z0z_0. In particular ∣f∣\lvert f\rvert has no strict local maximum anywhere in Ω\Omega; the word "strict" is ours, and the next paragraph says why it is needed.

    Form (a) is how Ahlfors (Ahlfors (1979), Thm 12, p. 134) states the principle, with "non-constant" among the hypotheses. Hassani (Hassani (2013), Prop. 10.5.4, p. 317) states instead that ∣f∣\lvert f\rvert cannot have a local maximum in the region of analyticity, with no exception made for constants. Read literally that is not the same as form (b) and it is not available to us: a constant ff has a (non-strict) local maximum of ∣f∣\lvert f\rvert at every point. Form (b) is the conditional version — a local maximum is allowed, but only if ff is constant nearby — and it is what we prove and use. The proof below gives both forms, and (b) yields (a) once one knows that a function constant on a disk is constant on the whole region, which is the identity theorem (Taylor series); we do not use that here and prove (a) directly.

  5. Boundary form. If E⊂CE \subset \mathbb{C} is closed and bounded, ff is continuous on EE and analytic on the interior of EE, then max⁡E∣f∣\max_{E}\lvert f\rvert is attained at a point of the boundary of EE. In particular, if ∣f∣≤M\lvert f\rvert \le M on the boundary then ∣f∣≤M\lvert f\rvert \le M throughout EE.

Idea of the proofOriginal

Everything here comes from one inequality applied to one formula. Cauchy's integral formula writes f(n)(a)f^{(n)}(a) as an average of ff over a circle against the kernel (z−a)−n−1(z-a)^{-n-1}; the ML estimate turns that average into the product of the largest value of the integrand and the length of the circle, and the powers of rr do not cancel — they leave r−nr^{-n}. So the derivatives of an analytic function cannot be large unless the function is large nearby, and the larger the circle one can use, the smaller the bound. Liouville is this read at n=1n=1 with r→∞r \to \infty: if M(r)M(r) never grows, f′f' must vanish identically. The fundamental theorem of algebra follows because a polynomial without roots would have a bounded reciprocal.

The maximum principle is of a different flavour and has two proofs. The topological one: a non-constant analytic function maps open sets to open sets, so the image of a neighbourhood of z0z_0 contains a disk around f(z0)f(z_0), which contains points of larger modulus. The analytic one, which is self-contained here: the integral formula on a circle says f(z0)f(z_0) is the average of ff on the circle, and an average cannot have modulus larger than the largest value averaged unless all the values line up — so if ∣f(z0)∣\lvert f(z_0)\rvert were maximal, ∣f∣\lvert f\rvert would be constant on every small circle about z0z_0, hence near z0z_0, and then ff itself would be constant.

Proof(Parts 1 to 3)AdaptedAdapted from Ahlfors (1979), Ch. 4, §2.3, p. 122; Riley, Hobson & Bence (2006), §24.10, p. 853; Arfken, Weber & Harris (2013), §11.4, p. 490

Part 1. Let CC be the positively oriented circle ∣z−a∣=r\lvert z - a\rvert = r, of length 2πr2\pi r. By Cauchy's integral formula for the derivatives,

f(n)(a)=n!2πi∮Cf(z)(z−a)n+1 dz.f^{(n)}(a) = \frac{n!}{2\pi i}\oint_C \frac{f(z)}{(z-a)^{n+1}}\,dz .

On CC the integrand has modulus ∣f(z)∣/rn+1≤M(r)/rn+1\lvert f(z)\rvert / r^{n+1} \le M(r)/r^{n+1}, since ∣z−a∣=r\lvert z - a\rvert = r exactly. The ML estimate (contour integrals) gives

∣f(n)(a)∣≤n!2π⋅M(r)rn+1⋅2πr=n! M(r)rn.\bigl\lvert f^{(n)}(a)\bigr\rvert \le \frac{n!}{2\pi}\cdot\frac{M(r)}{r^{n+1}}\cdot 2\pi r = \frac{n!\,M(r)}{r^{n}} .

Part 2. Let ff be entire with ∣f∣≤M\lvert f\rvert \le M everywhere, and let a∈Ca \in \mathbb{C}. Then M(r)≤MM(r) \le M for every r>0r > 0, so part 1 with n=1n = 1 gives

∣f′(a)∣≤M(r)r≤Mrfor all r>0.\lvert f'(a)\rvert \le \frac{M(r)}{r} \le \frac{M}{r} \qquad\text{for all } r > 0 .

Here the hypothesis that ff is entire is used twice: to know that every rr is allowed, and that the bound MM holds on every circle. Letting r→∞r \to \infty gives f′(a)=0f'(a) = 0. As aa was arbitrary, f′≡0f' \equiv 0 on C\mathbb{C}, so ux=vx=0u_x = v_x = 0 and, by the Cauchy-Riemann equations, uy=vy=0u_y = v_y = 0 as well; a function with vanishing partial derivatives on the connected set C\mathbb{C} is constant.

Part 3. Let P(z)=anzn+⋯+a0P(z) = a_nz^n + \dots + a_0 with n≥1n \ge 1 and an≠0a_n \ne 0, and suppose PP had no root. Then g=1/Pg = 1/P is entire. For ∣z∣=R\lvert z\rvert = R,

∣P(z)∣≥∣an∣Rn−∑k=0n−1∣ak∣Rk≥Rn(∣an∣−∑k=0n−1∣ak∣Rn−k)→R→∞∞,\lvert P(z)\rvert \ge \lvert a_n\rvert R^n - \sum_{k=0}^{n-1}\lvert a_k\rvert R^k \ge R^n\Bigl(\lvert a_n\rvert - \sum_{k=0}^{n-1}\frac{\lvert a_k\rvert}{R^{n-k}}\Bigr) \xrightarrow[R\to\infty]{} \infty ,

by the triangle inequality, so g(z)→0g(z) \to 0 as ∣z∣→∞\lvert z\rvert \to \infty. Choose R0R_0 with ∣g∣≤1\lvert g\rvert \le 1 for ∣z∣≥R0\lvert z\rvert \ge R_0; on the compact disk ∣z∣≤R0\lvert z\rvert \le R_0 the continuous function ∣g∣\lvert g\rvert is bounded by compactness. So gg is a bounded entire function, hence constant by part 2 — and the constant is 00 because g→0g \to 0, which is impossible since g=1/Pg = 1/P never vanishes. Therefore PP has a root λ1\lambda_1.

Dividing, P(z)=(z−λ1)Q(z)P(z) = (z - \lambda_1)Q(z) with deg⁡Q=n−1\deg Q = n - 1: indeed P(z)−P(λ1)=P(z)P(z) - P(\lambda_1) = P(z) is a linear combination of the terms zk−λ1kz^k - \lambda_1^k, each divisible by z−λ1z - \lambda_1. Repeating nn times reduces the degree to 00 and gives the stated factorization, with the leading coefficient ana_n; a product of linear factors vanishes exactly at the λk\lambda_k, so PP has no further roots. □\square

Proof(Parts 4 and 5)AdaptedAdapted from Ahlfors (1979), Ch. 4, §3.4, Thms 12 and 12' and Eqs. (34)–(35), pp. 133–135

Part 4. Steps 1 to 3 prove form (b); Step 4 then proves form (a). Let ff be analytic on the region Ω\Omega and suppose first only that ∣f∣\lvert f\rvert has a local maximum at z0∈Ωz_0 \in \Omega: there is r0>0r_0 > 0 with ∣f(z)∣≤∣f(z0)∣\lvert f(z)\rvert \le \lvert f(z_0)\rvert whenever ∣z−z0∣<r0\lvert z - z_0\rvert < r_0, and we may shrink r0r_0 so that this disk lies in Ω\Omega.

Step 1 (the mean value property). Choose 0<r<r00 < r < r_0, so that the closed disk ∣z−z0∣≤r\lvert z - z_0\rvert \le r lies in Ω\Omega and in the disk where the maximality holds. By the mean value consequence of Cauchy's formula,

f(z0)=12π∫02πf(z0+ρeiθ) dθfor every 0<ρ≤r,f(z_0) = \frac{1}{2\pi}\int_0^{2\pi}f\bigl(z_0 + \rho e^{i\theta}\bigr)\,d\theta \qquad\text{for every } 0 < \rho \le r ,

and therefore, taking moduli,

∣f(z0)∣≤12π∫02π∣f(z0+ρeiθ)∣ dθ.(∗)\lvert f(z_0)\rvert \le \frac{1}{2\pi}\int_0^{2\pi}\bigl\lvert f(z_0 + \rho e^{i\theta})\bigr\rvert\,d\theta . \tag{$\ast$}

Step 2 (∣f∣\lvert f\rvert is constant on small circles). Put h(θ)=∣f(z0)∣−∣f(z0+ρeiθ)∣h(\theta) = \lvert f(z_0)\rvert - \lvert f(z_0 + \rho e^{i\theta})\rvert for 0<ρ≤r0 < \rho \le r. By the assumed local maximality, h≥0h \ge 0; and hh is continuous, since ff is. Rearranging (∗)(\ast) gives 12π∫02πh(θ) dθ≤0\frac{1}{2\pi}\int_0^{2\pi}h(\theta)\,d\theta \le 0. A continuous non-negative function with non-positive integral vanishes identically, so ∣f∣=∣f(z0)∣\lvert f\rvert = \lvert f(z_0)\rvert on every circle ∣z−z0∣=ρ\lvert z - z_0\rvert = \rho with ρ≤r\rho \le r, hence on the whole disk D={∣z−z0∣≤r}D = \{\lvert z - z_0\rvert \le r\}.

Step 3 (constant modulus forces constant ff; this proves form (b)). On DD we have u2+v2=cu^2 + v^2 = c with c=∣f(z0)∣2c = \lvert f(z_0)\rvert^2. If c=0c = 0 then f≡0f \equiv 0 on DD. Otherwise differentiate: uux+vvx=0uu_x + vv_x = 0 and uuy+vvy=0uu_y + vv_y = 0. Using the Cauchy-Riemann equations uy=−vxu_y = -v_x, vy=uxv_y = u_x, the second becomes −uvx+vux=0-uv_x + vu_x = 0. The two linear equations

u ux+v vx=0,v ux−u vx=0\begin{aligned} u\,u_x + v\,v_x &= 0, \\ v\,u_x - u\,v_x &= 0 \end{aligned}

in the unknowns (ux,vx)(u_x, v_x) have determinant −(u2+v2)=−c≠0-(u^2+v^2) = -c \ne 0, so ux=vx=0u_x = v_x = 0 on DD, and then uy=vy=0u_y = v_y = 0 too. Hence ff is constant on DD.

Step 4 (form (a): spread the constancy over Ω\Omega). Assume now the global hypothesis: ∣f(z)∣≤∣f(z0)∣\lvert f(z)\rvert \le \lvert f(z_0)\rvert for every z∈Ωz \in \Omega. Put c0=f(z0)c_0 = f(z_0) and

A={z∈Ω:f(z)=c0}.A = \{z \in \Omega : f(z) = c_0\} .
  • A≠∅A \ne \emptyset, since z0∈Az_0 \in A.
  • AA is closed in Ω\Omega: it is f−1({c0})f^{-1}(\{c_0\}) and ff is continuous.
  • AA is open. Let z∈Az \in A. Then ∣f(z)∣=∣c0∣=∣f(z0)∣\lvert f(z)\rvert = \lvert c_0\rvert = \lvert f(z_0)\rvert, and by the global hypothesis ∣f∣≤∣f(z)∣\lvert f\rvert \le \lvert f(z)\rvert on all of Ω\Omega, so ∣f∣\lvert f\rvert certainly has a local maximum at zz. By Steps 1 to 3 applied at zz in place of z0z_0, ff is constant on a disk about zz, and that constant is f(z)=c0f(z) = c_0; so this disk lies in AA.

A non-empty subset of a connected set that is both open and closed in it is the whole set (connectedness), so A=ΩA = \Omega and f≡c0f \equiv c_0 on Ω\Omega.

Part 5. EE is closed and bounded and ∣f∣\lvert f\rvert is continuous on it, so by compactness the maximum is attained, say at z0∈Ez_0 \in E. If z0∈∂Ez_0 \in \partial E, we are done. Otherwise z0z_0 lies in the interior of EE; let Ω0\Omega_0 be the connected component of int⁡E\operatorname{int}E containing z0z_0, a region. Since ∣f∣≤∣f(z0)∣\lvert f\rvert \le \lvert f(z_0)\rvert on all of EE, in particular near z0z_0, part 4 applies on Ω0\Omega_0 and ff is constant there, of modulus ∣f(z0)∣\lvert f(z_0)\rvert.

Now ∂Ω0\partial \Omega_0 is non-empty, because Ω0\Omega_0 is a non-empty bounded open set, and ∂Ω0⊆∂E\partial\Omega_0 \subseteq \partial E: if some p∈∂Ω0p \in \partial\Omega_0 lay in int⁡E\operatorname{int}E, take a disk N⊆int⁡EN \subseteq \operatorname{int}E about pp; as pp is a limit of points of Ω0\Omega_0, the connected set NN meets Ω0\Omega_0, so N∪Ω0N \cup \Omega_0 is a connected subset of int⁡E\operatorname{int}E and hence lies in the component Ω0\Omega_0; then p∈Ω0p \in \Omega_0, which is open, contradicting p∈∂Ω0p \in \partial\Omega_0. And p∉Ep \notin E is impossible because ∂Ω0⊆Ω0‾⊆E\partial\Omega_0 \subseteq \overline{\Omega_0} \subseteq E, EE being closed. So p∈∂Ep \in \partial E.

By continuity of ff on EE, ∣f∣=∣f(z0)∣\lvert f\rvert = \lvert f(z_0)\rvert on Ω0‾\overline{\Omega_0} and in particular at any point of ∂Ω0⊆∂E\partial\Omega_0 \subseteq \partial E: the maximum is attained on the boundary of EE after all. The final sentence of part 5 is this statement rearranged. □\square

RemarkOriginal

Where the hypotheses are used, and how sharp the results are.

  • "Entire" is essential in Liouville, not merely "analytic and bounded". The function f(z)=1/(1−z)f(z) = 1/(1-z) is analytic and bounded on the disk ∣z∣<12\lvert z\rvert < \tfrac12 but not constant: the proof needs arbitrarily large circles, which only an entire function provides. Likewise eize^{iz} is entire and bounded on the closed upper half-plane, and is not constant.
  • Boundedness cannot be weakened to "bounded on the real axis". sin⁡z\sin z is entire with ∣sin⁡x∣≤1\lvert \sin x\rvert \le 1 for real xx, yet ∣sin⁡(10i)∣=sinh⁡10≈11013.23\lvert \sin(10i)\rvert = \sinh 10 \approx 11013.23. Part 2 demands a bound on the whole plane.
  • Polynomial growth gives polynomials. If ff is entire with ∣f(z)∣≤A+B∣z∣p\lvert f(z)\rvert \le A + B\lvert z\rvert^{p} for some p≥0p \ge 0, then part 1 with n>pn > p gives ∣f(n)(0)∣≤n!(A+Brp)/rn→0\lvert f^{(n)}(0)\rvert \le n!(A + Br^{p})/r^{n} \to 0, so f(n)≡0f^{(n)} \equiv 0 for every n>pn > p (the same argument at every centre) and ff is a polynomial of degree at most ⌊p⌋\lfloor p\rfloor — the floor, because nn runs over integers and the smallest integer exceeding pp is ⌊p⌋+1\lfloor p \rfloor + 1. Liouville is the case p=0p = 0.
  • Cauchy's estimate is rarely sharp, and the choice of rr matters. The advice to choose rr so as to minimise M(r)r−nM(r)r^{-n}, and the example of a function whose maximum modulus blows up like 1/(R−r)1/(R-r) at the edge of its disk of analyticity, are Ahlfors (Ahlfors (1979), Ch. 4, §2.3, p. 123, with Exercise 4)'s; the following numbers are ours. For f(z)=1/(2−z)f(z) = 1/(2-z) on ∣z∣<2\lvert z\rvert < 2 we have M(r)=1/(2−r)M(r) = 1/(2-r) and f(n)(0)=n!/2n+1f^{(n)}(0) = n!/2^{n+1} exactly. The bound is n!/(rn(2−r))n!/\bigl(r^{n}(2-r)\bigr), minimised at r=2n/(n+1)r = 2n/(n+1): for n=1n = 1 this gives 11 against the true value 0.250.25, and for n=3n=3 it gives 3.5563.556 against 0.3750.375. The estimate is an order-of-magnitude tool.
  • The maximum principle says nothing about the minimum unless ff has no zero: ∣z∣\lvert z\rvert on the closed unit disk has its minimum at the interior point 00. Applying part 4 to 1/f1/f gives a minimum principle wherever f≠0f \ne 0.
  • "Region" means connected. On a disconnected open set ff can be one constant on one piece and a larger constant on another, and then ∣f∣\lvert f\rvert does attain a maximum. Step 4 of the proof is exactly where connectedness enters.
  • Real and imaginary parts obey the same principle. For the real part, apply part 4 to efe^{f}, which is analytic with ∣ef∣=eRe⁡f\lvert e^f\rvert = e^{\operatorname{Re}f}: since t↦ett \mapsto e^t is strictly increasing, a maximum of Re⁡f\operatorname{Re}f is a maximum of ∣ef∣\lvert e^f\rvert, so there is none unless ff is constant. For the imaginary part the right auxiliary function is e−ife^{-if}, which has ∣e−if∣=eIm⁡f\lvert e^{-if}\rvert = e^{\operatorname{Im}f}; efe^{f} says nothing about Im⁡f\operatorname{Im}f. Either way this is the maximum principle for harmonic functions (harmonic functions), and minima are covered by e−fe^{-f} and eife^{if}.
Example(Cauchy's estimate, Liouville's theorem and a growth bound)Original
  1. A numerical estimate. Let f(z)=1/(2−z)f(z) = 1/(2-z) and a=0a = 0. On ∣z∣=1\lvert z\rvert = 1, ∣2−z∣≥1\lvert 2 - z\rvert \ge 1, so M(1)=1M(1) = 1 (attained at z=1z = 1) and part 1 gives ∣f′(0)∣≤1\lvert f'(0)\rvert \le 1. The true value is f′(z)=(2−z)−2f'(z) = (2-z)^{-2}, so f′(0)=0.25f'(0) = 0.25: the bound holds with a factor of 44 to spare, as it must since ∣f∣\lvert f\rvert is far below M(1)M(1) on most of the circle.
  2. A growth bound applied. Suppose ff is entire with ∣f(z)∣≤5+2∣z∣3/2\lvert f(z)\rvert \le 5 + 2\lvert z\rvert^{3/2} for all zz. Take n=2n = 2 in part 1 at any centre aa: ∣f′′(a)∣≤2(5+2(∣a∣+r)3/2)/r2\lvert f''(a)\rvert \le 2(5 + 2(\lvert a\rvert + r)^{3/2})/r^{2}, and letting r→∞r \to \infty with aa fixed gives f′′(a)=0f''(a) = 0. So f′′≡0f'' \equiv 0 and f(z)=αz+βf(z) = \alpha z + \beta. The conclusion "degree at most 11" cannot be improved to "constant", because one non-constant affine function does satisfy the hypothesis: f(z)=zf(z) = z has ∣z∣≤5+2∣z∣3/2\lvert z\rvert \le 5 + 2\lvert z\rvert^{3/2} for every zz — for ∣z∣≤5\lvert z\rvert \le 5 because the left side is then at most 55, and for ∣z∣>5\lvert z\rvert > 5 because 2∣z∣3/2>25 ∣z∣>∣z∣2\lvert z\rvert^{3/2} > 2\sqrt5\,\lvert z\rvert > \lvert z\rvert. (Most affine functions do not satisfy it: f≡100f \equiv 100 already fails at z=0z = 0, where the bound is 55. The hypothesis has fixed constants, so it is a genuine restriction on α\alpha and β\beta; all that sharpness needs is one non-constant example.) Note what the growth exponent 32\tfrac32 buys: not a 32\tfrac32-power function, which is not entire, but a polynomial of degree ⌊32⌋=1\lfloor \tfrac32\rfloor = 1.
  3. An entire function cannot omit a disk. Suppose ff is entire and its values avoid some disk: ∣f(z)−w0∣≥δ>0\lvert f(z) - w_0\rvert \ge \delta > 0 for all zz, with δ>0\delta>0 fixed. Then g=1/(f−w0)g = 1/(f - w_0) is entire, since the denominator never vanishes, and ∣g∣≤1/δ\lvert g\rvert \le 1/\delta everywhere. By part 2, gg is constant, hence so is ff. So the image of a non-constant entire function is dense in C\mathbb{C}. The hypothesis really is about omitting a disk and not a point: eze^{z} is entire and non-constant and its image is exactly C∖{0}\mathbb{C}\setminus\{0\}, which omits one point but no disk. Concretely, ez=10−6e^{z} = 10^{-6} is solved by z=−6log⁡10≈−13.815511z = -6\log 10 \approx -13.815511, so values arbitrarily close to the omitted point do occur.
Example(Maximum modulus on a square and a disk)AdaptedAdapted from Ahlfors (1979), Ch. 4, §3.4, Thm 13 and the rescaling that follows it, p. 135 (item 3 only)
  1. An exact maximum. Let f(z)=z2−2zf(z) = z^2 - 2z on the closed unit disk E={∣z∣≤1}E = \{\lvert z\rvert \le 1\}. By part 5 the maximum of ∣f∣\lvert f\rvert is on ∣z∣=1\lvert z\rvert = 1, where ∣f(z)∣=∣z∣ ∣z−2∣=∣z−2∣\lvert f(z)\rvert = \lvert z\rvert\,\lvert z - 2\rvert = \lvert z-2\rvert. For z=eiθz = e^{i\theta}, ∣z−2∣2=5−4cos⁡θ\lvert z - 2\rvert^2 = 5 - 4\cos\theta, maximal at θ=π\theta = \pi. So max⁡E∣f∣=3\max_E \lvert f\rvert = 3, attained at z=−1z = -1 and nowhere else; a scan of 10001000 points of the circle confirms the maximum 3.0000003.000000 at θ=π\theta = \pi. At the interior point 00, ∣f(0)∣=0\lvert f(0)\rvert = 0.

  2. A bound obtained without computing anything. Let f(z)=ezf(z) = e^{z} and E=[0,1]×[0,1]E = [0,1]\times[0,1], the closed unit square. Then ∣ez∣=eRe⁡z≤e\lvert e^z\rvert = e^{\operatorname{Re}z} \le e on EE, with equality exactly on the boundary edge Re⁡z=1\operatorname{Re}z = 1: the maximum e≈2.718282e \approx 2.718282 is on ∂E\partial E, as part 5 requires, and no interior point comes close to it.

  3. A Schwarz-type estimate. Suppose gg is analytic on an open set containing the closed disk ∣z∣≤2\lvert z\rvert \le 2, with ∣g∣≤10\lvert g\rvert \le 10 on the circle ∣z∣=2\lvert z\rvert = 2 and g(0)=0g(0) = 0. The conclusion is that

    ∣g(z)∣≤5∣z∣  for ∣z∣≤2,∣g′(0)∣≤5,\lvert g(z)\rvert \le 5\lvert z\rvert \ \text{ for } \lvert z\rvert \le 2, \qquad \lvert g'(0)\rvert \le 5 ,

    which is Ahlfors (Ahlfors (1979), Thm 13 and the paragraph after it, p. 135)'s Schwarz lemma with his bound MM set to 1010 and his radius RR to 22, so that M/R=5M/R = 5. His proof puts h(z)=g(z)/zh(z) = g(z)/z with h(0)=g′(0)h(0) = g'(0) and calls that extension analytic; justifying the extension needs the removable-singularity theorem (isolated singularities), which this unit does not have. The following route uses only parts 1 and 5 above, and works on annuli so that the point 00 never has to be discussed.

    • The derivative bound first. Part 1 with n=1n = 1, a=0a = 0, r=2r = 2 and M(2)≤10M(2) \le 10 gives ∣g′(0)∣≤1!⋅10/2=5\lvert g'(0)\rvert \le 1!\cdot 10/2 = 5 outright.
    • Boundary values of h=g/zh = g/z on an annulus. Fix 0<ε<20 < \varepsilon < 2 and let Eε={ε≤∣z∣≤2}E_\varepsilon = \{\varepsilon \le \lvert z\rvert \le 2\}. On EεE_\varepsilon the function h(z)=g(z)/zh(z) = g(z)/z is continuous and, on the interior, analytic — the denominator does not vanish anywhere on EεE_\varepsilon, so no extension is needed. The boundary of EεE_\varepsilon is the two circles ∣z∣=2\lvert z\rvert = 2, where ∣h∣=∣g∣/2≤5\lvert h\rvert = \lvert g\rvert/2 \le 5, and ∣z∣=ε\lvert z\rvert = \varepsilon, where ∣h∣≤Mg(ε)/ε\lvert h\rvert \le M_g(\varepsilon)/\varepsilon with Mg(ε)=max⁡∣z∣=ε∣g∣M_g(\varepsilon) = \max_{\lvert z\rvert = \varepsilon}\lvert g\rvert.
    • The inner circle is harmless as ε→0\varepsilon \to 0. Since g(0)=0g(0) = 0, the difference quotient g(z)/zg(z)/z tends to g′(0)g'(0) as z→0z \to 0, so Mg(ε)/ε=max⁡∣z∣=ε∣g(z)/z∣→∣g′(0)∣≤5M_g(\varepsilon)/\varepsilon = \max_{\lvert z\rvert = \varepsilon}\lvert g(z)/z\rvert \to \lvert g'(0)\rvert \le 5. Hence for any η>0\eta > 0 there is εη\varepsilon_\eta with Mg(ε)/ε≤5+ηM_g(\varepsilon)/\varepsilon \le 5 + \eta for all ε<εη\varepsilon < \varepsilon_\eta.
    • Conclude. By part 5 applied to hh on EεE_\varepsilon, ∣h∣≤max⁡{5, Mg(ε)/ε}≤5+η\lvert h\rvert \le \max\{5,\ M_g(\varepsilon)/\varepsilon\} \le 5 + \eta on all of EεE_\varepsilon, for every ε<εη\varepsilon < \varepsilon_\eta. Any zz with 0<∣z∣≤20 < \lvert z\rvert \le 2 lies in EεE_\varepsilon for all small enough ε\varepsilon, so ∣h(z)∣≤5+η\lvert h(z)\rvert \le 5 + \eta; letting η→0\eta \to 0 gives ∣h(z)∣≤5\lvert h(z)\rvert \le 5, that is ∣g(z)∣≤5∣z∣\lvert g(z)\rvert \le 5\lvert z\rvert, for every such zz, and at z=0z = 0 the inequality reads 0≤00 \le 0.

    Both bounds are attained: g(z)=5zg(z) = 5z has ∣g∣=10\lvert g\rvert = 10 on ∣z∣=2\lvert z\rvert = 2, g(0)=0g(0) = 0, ∣g(z)∣=5∣z∣\lvert g(z)\rvert = 5\lvert z\rvert everywhere and g′(0)=5g'(0) = 5.

ApplicationOriginal
  • Earnshaw's theorem. An electrostatic potential Φ\Phi in a charge-free region is harmonic, and by the last bullet of the remark it has no interior maximum or minimum. A stable equilibrium of a point charge would require a strict minimum of its potential energy qΦq\Phi at an interior point, which is therefore impossible: no arrangement of fixed charges can trap a charged particle in stable static equilibrium. The same argument forbids static magnetic trapping of a paramagnetic particle and is why ion traps use oscillating fields (Paul trap) and why magnetic levitation needs diamagnets, superconductors or feedback.
  • No growth without singularities. Liouville's theorem says a non-constant analytic function must either be unbounded or have a singularity somewhere, including possibly at infinity. In dispersion theory this is used in reverse: a response function that is analytic in a half-plane and tends to a constant at infinity is pinned down by its singularities, which are the resonances of the system (dispersion relations). The statement that the only bounded entire functions are constants is also why a nontrivial scattering amplitude must have cuts or poles.
  • Stability and spectra. The fundamental theorem of algebra guarantees that an n×nn\times n matrix has nn eigenvalues counted with multiplicity, which is the starting point of eigenvalues and eigenvectors and hence of normal-mode analysis and linear stability theory. Cauchy's estimate, applied to the resolvent, bounds how fast a matrix function can vary and underlies error bounds for contour-integral algorithms for etAe^{tA} (the matrix exponential).
  • A maximum principle for diffusion. The argument of part 4 is the complex-analytic ancestor of the maximum principles for elliptic and parabolic equations: the temperature inside a body with no internal sources never exceeds the maximum on its surface and at the initial time (the maximum principle for the heat equation), which is what makes numerical schemes for diffusion stable when they respect it.