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.
A complex number z=x+iy with x,y∈R has modulus∣z∣=x2+y2 and, if z=0, an argument: any real θ with
x=∣z∣cosθ,y=∣z∣sinθ,
so that z=r(cosθ+isinθ) with r=∣z∣. The argument is fixed only up to an integer multiple of 2π; the whole set of its values is written argz, and the one value in a chosen half-open interval of length 2π 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 (−π,π] and writing the principal argument Argz.
In polar form multiplication adds arguments and multiplies moduli: if z1=r1(cosθ1+isinθ1) and z2=r2(cosθ2+isinθ2), then the addition formulas for cos and sin give
z1z2=r1r2(cos(θ1+θ2)+isin(θ1+θ2)),
and for z2=0 the quotient z1/z2 has modulus r1/r2 and argument θ1−θ2. Once the exponential is available this is written z=reiθ, and the two rules become the single identity eiθ1eiθ2=ei(θ1+θ2).
Definition(The exponential)AdaptedAdapted from Ahlfors (1979), Ch. 2, §3.1, Eq. (21), p. 43
For z∈C set
ez=expz=n=0∑∞n!zn.
The series converges absolutely for every z, because ∣z∣n/n!→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. For real z it is the usual real exponential. Ahlfors takes the differential equation f′=f with f(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
Addition theorem.ea+b=eaeb for all a,b∈C.
No zeros.ez=0 for every z, and e−z=1/ez.
Modulus.ez=ezˉ, hence ∣ex+iy∣=ex for real x,y; in particular ∣eiy∣=1.
Periodicity. The periods of ez are exactly the numbers iω with ω a period of y↦eiy, and these ω are exactly the integer multiples of one smallest positive period. That smallest positive period is denoted 2π, so
ez+2πi=ezfor all z,ez=1⟺z∈2πiZ.
Surjectivity onto C∖{0}. For every w=0 the equation ez=w has a solution, and its solutions form a single coset z0+2πiZ.
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 π enters. Nothing yet tells us that ez has any period at all, so we must produce one: truncating the real series for cos gives numerical inequalities sharp enough to show that cos changes sign before 3, and the first zero y0 of cos delivers the period 4y0. Part 5 then splits ez=w into one real equation for ∣w∣ 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):
2. No zeros. Put b=−a in part 1: eae−a=e0=1. A product equal to 1 has no zero factor.
3. Modulus. The coefficients 1/n! are real, so conjugating the series termwise gives ez=ezˉ. Hence ∣eiy∣2=eiyeiy=eiye−iy=1 for real y, and by part 1, ∣ex+iy∣=∣ex∣∣eiy∣=ex, using that ex>0 for real x: for x≥0 every term of the series is non-negative and the first is 1, and for x<0 we have ex=1/e−x>0 by part 2.
4. Periodicity.A period exists. Write C(y)=21(eiy+e−iy) and S(y)=2i1(eiy−e−iy) for real y; from the series these are the real functions cos and sin, with C′=−S, S′=C, C(0)=1, S(0)=0, and C2+S2=1 because C2+S2=eiye−iy=1 by parts 1 and 2. In particular ∣C∣≤1 and ∣S∣≤1. Four successive integrations from 0 now trap C and S between consecutive partial sums of their series. From C≤1 and S(0)=0 we get S(y)≤y for y≥0; integrating C′=−S≥−y gives C(y)≥1−y2/2; integrating S′=C≥1−y2/2 gives S(y)≥y−y3/6; integrating C′=−S≤−y+y3/6 gives
C(y)≤1−2y2+24y4,y≥0.
At y=3 the right-hand side is 1−23+83=−81<0. Since C is continuous with C(0)=1, it has a zero in (0,3); let y0 be the smallest one. From C(y0)=0 and C2+S2=1 we get S(y0)2=1; and S(y0)>0, because S(y)≥y−y3/6>0 for 0<y<6 and y0<3<6. Hence S(y0)=1, so
eiy0=C(y0)+iS(y0)=i,e4iy0=i4=1:
the number 4y0 is a period of y↦eiy.
It is the smallest. For 0<y<y0 we have C(y)>0 (by minimality of y0 and continuity, since C(0)=1), so S′=C>0 there and S increases strictly from S(0)=0. Hence 0<S(y)<S(y0)=1 on (0,y0), and therefore C(y)2=1−S(y)2∈(0,1). So for 0<y<y0 the number eiy=C(y)+iS(y) has both parts nonzero and is not one of 1,−1,i,−i; consequently e4iy=1, because w4=1 forces w∈{1,−1,i,−i}. No positive period is therefore smaller than 4y0.
All periods are multiples. Let ω>0 be a period and ω0=4y0. Choose the integer n with nω0≤ω<(n+1)ω0. If ω=nω0 then ω−nω0 is a positive period smaller than ω0, which we have excluded. Finally, ez+c=ez for all z forces ec=1, so ∣ec∣=1 gives Rec=0 by part 3, and then c=iω with ω a period of y↦eiy. Writing 2π=4y0 gives the displayed statement.
5. Surjectivity. Let w=0 and seek z=x+iy. By part 3, ex+iy=w splits into ex=∣w∣ and eiy=w/∣w∣. The first has the unique real solution x=log∣w∣, since the real exponential is a strictly increasing bijection of R onto (0,∞). For the second, the map y↦eiy sends [0,2π) onto the unit circle: it is continuous; part 4 showed that S>0 on (0,y0), so C′=−S<0 there and C decreases strictly from C(0)=1 to C(y0)=0, which together with C2+S2=1 and S≥0 traces the first quadrant of the circle exactly once, and the four reflections ei(y+kπ/2)=ikeiy carry this arc onto the remaining three quadrants. Two solutions of eiy=w/∣w∣ differ by a period, by part 4, so the solution set of ez=w is a single coset of 2πiZ. □
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
with tanz=sinz/cosz and tanhz=sinhz/coshz where the denominators do not vanish. Substituting the exponential series gives the familiar power series, so on R these are the usual functions; and (cosz)′=−sinz, (sinz)′=cosz. Euler's formula eiz=cosz+isinz and the identity cos2z+sin2z=1 are immediate from the definitions, as are
coshz=cosiz,sinhz=−isiniz,
so the hyperbolic and circular functions are each other's restrictions to perpendicular lines. Separating real and imaginary parts, sin(x+iy)=sinxcoshy+icosxsinhy, whence ∣sin(x+iy)∣2=sin2x+sinh2y: unlike their real restrictions, sin and cos are unbounded on C.
For w=0, a logarithm of w is any solution of ez=w. By part 5 of the statement above there are infinitely many, differing by multiples of 2πi; their imaginary parts are exactly the values of argw, and the whole set is written
Lnw=log∣w∣+iargw=log∣w∣+i(θ+2kπ),k∈Z,
with θ any fixed argument of w and log the real logarithm of the positive number ∣w∣. The number 0 has no logarithm. The principal valuelnw=log∣w∣+iArgw takes θ∈(−π,π].
For a=0 and b∈C the power is defined by
ab=exp(bLna),
again a set of values, which differ by the factors e2πikb. It has a single value exactly when b∈Z, and then agrees with the repeated product; it has q values when b=p/q in lowest terms with q∈N; otherwise it has infinitely many. Its principal value uses lna.
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 and n∈Z,
(cosθ+isinθ)n=cosnθ+isinnθ(de Moivre’s theorem).
For a=0 with a=ρ(cosφ+isinφ) and n∈N, the equation zn=a has exactly n solutions,
zk=ρ1/n(cosnφ+2kπ+isinnφ+2kπ),k=0,1,…,n−1,
where ρ1/n is the positive real nth root. They have equal modulus and equally spaced arguments, so they are the vertices of a regular n-gon centred at 0. Taking a=1 gives the nth roots of unity1,ω,…,ωn−1 with ω=e2πi/n, and every solution of zn=a is ωkz0 for one k.
De Moivre's theorem is the addition theorem in polar clothing: cosθ+isinθ=eiθ by Euler's formula, and (eiθ)n=einθ for n≥0 by part 1 of the statement above and induction, while for n<0 one uses e−iθ=1/eiθ from part 2.
For the roots, write z=reiϑ with r>0; this costs nothing, because zn=a=0 forces z=0, and every nonzero z has such a form by the definition of polar form. Then zn=rneinϑ, and zn=a=ρeiφ holds if and only if the moduli agree and the arguments agree modulo a period:
rn=ρ, which by strict monotonicity of r↦rn on (0,∞) fixes r=ρ1/n uniquely;
einϑ=eiφ, which by part 4 means nϑ−φ∈2πZ, that is ϑ=(φ+2kπ)/n for some k∈Z.
Two integers k and k′ give the same z exactly when (k−k′)/n∈Z, so k=0,1,…,n−1 list the solutions once each: there are exactly n. Their moduli are all ρ1/n and consecutive arguments differ by 2π/n, which is the description of a regular n-gon. Dividing two solutions gives a solution of zn=1, which proves the last sentence. □
RemarkOriginal
Where the hypotheses bite, and which convention is which.
0 is excluded everywhere.ez never vanishes, so 0 has no logarithm and 0b is not defined by the formula above; the nth-root count also needs a=0, since zn=0 has the single root 0 with multiplicity n.
Principal ranges differ between books.Riley, Hobson and Bence (Riley, Hobson & Bence (2006), §3.3, p. 93 and Eq. (24.21), p. 834) use −π<θ≤π, and we follow them throughout. Ahlfors (Ahlfors (1979), Ch. 2, §3.4, p. 46) normalises arg to 0≤y<2π when solving eiy=w/∣w∣. 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 n, and §11.6 says only that the branch n=0 is the one usually called principal — their own Example 11.6.1 (p. 499) runs θ from 0 to 2π and explicitly declines the value −π. All three readings give ln(−1)=iπ, but just below the negative real axis the first two differ by 2πi: at z=−1−iε the first gives an argument near −π and the second one near +π.
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∈Z, which is the form Ahlfors (Ahlfors (1979), Ch. 1, §2.2, Eq. (20), p. 15) gives, because for non-integer n 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=21 and θ=3π, for instance, the left side is (cos3π+isin3π)1/2=(−1)1/2=±i with principal value +i, while the right side is cos23π+isin23π=−i: the two agree only after the other branch is chosen.
The addition law for Ln is a statement about sets.Ln(z1z2)=Lnz1+Lnz2 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πi. For principal values it can fail outright: ln(−1)+ln(−1)=2πi, while ln((−1)(−1))=ln1=0.
Powers are sets too, so the school rules need care: (ab)c=abc 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
A logarithm. For z=−1+i3 we have ∣z∣=2 and Argz=2π/3 (the point lies in the second quadrant), so lnz=log2+32πi≈0.693147+2.094395i, and Lnz=log2+i(2π/3+2kπ).
A power all of whose values are real. Take the base a=21(1+i)=eiπ/4, which has ∣a∣=1, and the purely imaginary exponent b=4i. Then log∣a∣=0, so Lna=i(π/4+2kπ) and
ab=exp(4i⋅i(π/4+2kπ))=e−π−8kπ,k∈Z,
a set of positive real numbers: the principal value is e−π≈0.0432139, and the next value up is e7π≈3.5533×109. The mechanism is visible in general. For a purely imaginary exponent b=it with t real, bLna=itlog∣a∣−t(θ+2kπ), so the values of ab are e−t(θ+2kπ)eitlog∣a∣: all of them are real exactly when tlog∣a∣∈πZ, and all positive exactly when tlog∣a∣∈2πZ. Our base has ∣a∣=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−2i.
A power that is not.(1+i)i has principal value exp(i(21log2+iπ/4))=e−π/4(cos(21log2)+isin(21log2))≈0.428829+0.154872i, since ln(1+i)=21log2+iπ/4≈0.346574+0.785398i.
Cube roots. For a=8i we have ρ=8, φ=π/2, so the roots are 2eiπ/6,2ei5π/6,2ei3π/2, that is 3+i, −3+i and −2i. Check the first: (3+i)3=(3+i)(3+23i−1)=(3+i)(2+23i)=23+6i+2i−23=8i.
A sine off the real axis.sin(1+2i)=sin1cosh2+icos1sinh2≈3.165779+1.959601i, with ∣sin(1+2i)∣2=sin21+sinh22≈13.86219 — far outside the real range [−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−2xcos5x and e−2xsin5x, and then the value of ∫0∞e−2xcos5xdx.
One complex integrand instead of two real ones. Both functions are components of e(−2+5i)x=e−2x(cos5x+isin5x), whose antiderivative is e(−2+5i)x/(−2+5i) because dxdeλx=λeλx for complex λ (differentiate the exponential series term by term). Rationalising with (−2+5i)=−2−5i and ∣−2+5i∣2=29,
∫e(−2+5i)xdx=29(−2−5i)e−2x(cos5x+isin5x)+c.
Separate the parts. The real part of the right-hand side is 291e−2x(−2cos5x+5sin5x) and the imaginary part is 291e−2x(−2sin5x−5cos5x), so
Check by differentiating back. For the first, 291e−2x(−2(5sin5x−2cos5x)+(25cos5x+10sin5x))=291e−2x⋅29cos5x=e−2xcos5x.
A number to check against. The improper integral ∫0∞e−2xcos5xdx converges absolutely, since ∣e(−2+5i)x∣=e−2x by part 3 of the statement above. The antiderivative tends to 0 as x→∞ and equals −2/29 at x=0, so the integral is 2/29≈0.0689655; likewise ∫0∞e−2xsin5xdx=5/29≈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 split into its parts.
ApplicationOriginal
Every linear system driven at a single frequency is handled with these functions. A real signal Acos(ωt+ϕ) is written Re(A^eiωt) with the complex amplitude A^=Aeiϕ, so that amplitude and phase travel together in one number and differentiation becomes multiplication by iω. Impedances in AC circuit theory, the complex refractive index n+iκ of an absorbing medium (where ∣eikz∣=e−κωz/c is the Beer-Lambert decay), and the e−iEt/ℏ phase of a stationary quantum state are all the same device.
The nth roots of unity are the characters of the cyclic group 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=a sit at the vertices of a regular polygon is why normal modes of a ring of n identical masses come in such symmetric families (normal modes). The unboundedness of sin and cos off the real axis, visible in ∣sin(x+iy)∣2=sin2x+sinh2y, is what makes evanescent waves possible: a "wavenumber" forced to be imaginary turns an oscillation into a growing or decaying exponential.
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.
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 f be defined on a neighbourhood of z∈C. Then f is complex differentiable at z if
f′(z)=h→0limhf(z+h)−f(z)
exists, the limit being taken over complex h→0 and required to be the same however h approaches 0. A function that is single-valued and complex differentiable at every point of an open set Ω is analytic (equivalently holomorphic) on Ω; a function analytic on all of C is entire. A point where f is defined on a neighbourhood but not differentiable is a singular point.
Writing f=u+iv with u,v real-valued functions of (x,y), we use subscripts for partial derivatives: ux=∂u/∂x and so on.
Differentiability at a point implies continuity there, since f(z+h)−f(z)=h⋅hf(z+h)−f(z)→0⋅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+iv be defined on an open set Ω⊆C.
Necessity. If f is complex differentiable at z0=x0+iy0, then all four partial derivatives ux,uy,vx,vy exist at z0 and satisfy there the Cauchy-Riemann equations
ux=vy,uy=−vx,
equivalently ∂f/∂x=−i∂f/∂y; and then f′(z0)=ux+ivx.
Sufficiency. If u and v have continuous first partial derivatives on Ω and satisfy the Cauchy-Riemann equations at every point of Ω, then f is analytic on Ω, and f′=ux+ivx 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.)
Harmonic conjugates. If in addition u and v have continuous second partial derivatives, then both are harmonic on Ω:
uxx+uyy=0,vxx+vyy=0.
A pair (u,v) of harmonic functions satisfying the Cauchy-Riemann equations is a pair of conjugate harmonic functions; v is determined by u up to an additive constant on a connected Ω, and u is the conjugate of −v. Connectedness is needed for that uniqueness clause alone: parts 1 and 2, and the harmonicity of u and v, 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 v in the proof. It can be removed: Ahlfors (Ahlfors (1979), Ch. 2, §1.2, p. 25) obtains harmonicity instead from the analyticity of f′, which Cauchy's integral formula supplies for every analytic f 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), 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) is multiplication by the complex number ux+ivx, 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.
Part 1 (necessity). Suppose f′(z0) exists. Restrict h to real values: then y is held at y0 and the difference quotient becomes the difference quotient of x↦f(x+iy0). Its limit exists and equals f′(z0), so ux and vx exist at z0 and
f′(z0)=ux+ivx.
Now restrict h to purely imaginary values h=ik with k real. Since 1/(ik)=−i/k,
so uy and vy exist at z0 as well. The two expressions for f′(z0) are the same complex number; equating real parts gives ux=vy and imaginary parts gives vx=−uy. (The combined form is fx=−ify, which on substituting f=u+iv reproduces the pair.)
Part 2 (sufficiency). Fix z0∈Ω and let h+ik be a small complex increment, h,k real. Because u has continuous first partials on the open set Ω, 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),
with the partials evaluated at z0 and ε1(h,k)/∣h+ik∣→0 as h+ik→0; and likewise
v(x0+h,y0+k)−v(x0,y0)=vxh+vyk+ε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+ivx is continuous. Adding i times the second to the first,
Now apply the hypothesis: by the Cauchy-Riemann equations at z0, uy+ivy=−vx+iux=i(ux+ivx). Hence the two linear terms combine,
f(z0+h+ik)−f(z0)=(ux+ivx)(h+ik)+ε,
and dividing by h+ik=0,
h+ikf(z0+h+ik)−f(z0)=(ux+ivx)+h+ikε.
The last term tends to 0 as h+ik→0, because ∣ε∣/∣h+ik∣≤∣ε1∣/∣h+ik∣+∣ε2∣/∣h+ik∣→0. So the limit exists, independently of the direction of approach, and f′(z0)=ux+ivx. As z0∈Ω was arbitrary, f is analytic on Ω, and f′=ux+ivx is continuous because ux and vx are.
Part 3 (harmonicity). Differentiate ux=vy with respect to x and uy=−vx with respect to y:
uxx=vyx,uyy=−vxy.
The second partial derivatives of v are continuous by hypothesis, so the mixed partials agree, vyx=vxy. Adding the two displayed equations therefore gives uxx+uyy=0. For v, note that −if=v−iu is differentiable wherever f is, with real part v and imaginary part −u; applying the result just proved to −if gives vxx+vyy=0. (Equivalently, differentiate ux=vy with respect to y and uy=−vx with respect to x and subtract.)
Uniqueness of the conjugate. If v and v~ both satisfy the Cauchy-Riemann equations with the same u, then w=v−v~ has wx=wy=0 on Ω, so w is constant on each connected component: constant on Ω if Ω is connected. □
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)=zzˉ2(z=0),f(0)=0.
In polar form f(reiθ)=re−3iθ, so ∣f(z)∣=∣z∣ and f is continuous at 0. On the real axis f(x)=x, giving ux(0)=1, vx(0)=0; on the imaginary axis f(iy)=iy, giving uy(0)=0, vy(0)=1. Both Cauchy-Riemann equations hold at the origin. Yet f(z)/z=e−4iθ, which is 1 along θ=0 and −1 along θ=π/4: no limit, so f is not differentiable at 0. 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θ0 shows the equations fail: moving radially, dz=eiθ0dr and df=e−3iθ0dr, so the quotient is e−4iθ0, whereas moving along the circle ∣z∣=r0, dz=ir0eiθ0dθ and df=−3ir0e−3iθ0dθ, so the quotient is −3e−4iθ0. The two differ, so no derivative exists at z0 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′ 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 Lnz satisfies the Cauchy-Riemann equations on C∖{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 v is the conjugate of u then u is the conjugate of −v, not of v. For (v,u) to be a conjugate pair in that order, the first Cauchy-Riemann equation would have to read vx=uy; what the equations for (u,v) actually give is vx=−uy, so the sign is wrong unless uy=0. Concretely, for z2=(x2−y2)+i(2xy) we have vx=2y and uy=−2y, which is the relation vx=−uy and not vx=uy; and indeed z2 is analytic while 2xy+i(x2−y2)=iz2 is nowhere analytic, its Cauchy-Riemann equations holding at the origin only.
A test for non-analyticity.f(z)=zˉ has u=x, v=−y, so ux=1=−1=vy everywhere: zˉ is continuous on C and differentiable nowhere. The same computation kills ∣z∣2=zzˉ except at the origin, and Rez everywhere.
Example(Constructing a harmonic conjugate)Original
Let u(x,y)=x3−3xy2+2y. First check that a conjugate can exist: uxx=6x and uyy=−6x, so u is harmonic on R2.
Integrate the first equation.vy=ux=3x2−3y2, so integrating in y at fixed x,
v(x,y)=3x2y−y3+φ(x)
with φ an unknown function of x alone.
Impose the second equation.vx=6xy+φ′(x) must equal −uy=6xy−2, so φ′(x)=−2 and φ(x)=−2x+c.
Conclusion.v(x,y)=3x2y−y3−2x+c, and by part 2 of the statement (the partials here are polynomials, hence continuous) f=u+iv is entire.
Identify f. Taking c=0 and comparing with z3=(x3−3xy2)+i(3x2y−y3) and −2iz=2y−2ix, we get f(z)=z3−2iz. Check at z=2+i: the formula gives (2+i)3−2i(2+i)=(2+11i)+(2−4i)=4+7i, while u(2,1)=8−6+2=4 and v(2,1)=12−1−4=7. They agree.
Derivative two ways.f′(z)=3z2−2i gives f′(2+i)=3(3+4i)−2i=9+10i; and ux+ivx at (2,1) is (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, the Cauchy-Riemann equations make the level curves of u and of v meet at right angles, because
∇u⋅∇v=uxvx+uyvy=ux(−uy)+uy(ux)=0.
Take f(z)=z2, so u=x2−y2 and v=2xy. The curves u=c are the rectangular hyperbolas x2−y2=c and the curves v=d are xy=d/2. At the point z0=2+i we have u=3, v=4, and the gradients are ∇u=(2x,−2y)=(4,−2) and ∇v=(2y,2x)=(2,4): their dot product is 8−8=0, and the tangent directions (2,4)/20 and (4,−2)/20 are perpendicular unit vectors. The orthogonality fails exactly where f′(z)=2z vanishes, that is at the origin, where the two hyperbola families degenerate into the pairs of lines y=±x and xy=0 crossing at 45∘ instead of 90∘.
ApplicationOriginal
Part 3 is the bridge between complex analysis and two-dimensional field theory. Any analytic f=u+iv on a plane region supplies, for free, a pair of solutions of Laplace's equation whose level sets are mutually orthogonal, so f can be read as a complex potential: in electrostatics the curves u=const are equipotentials and v=const are field lines; in ideal fluid flow u is the velocity potential and v the stream function, and 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.
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.
Let f be a multivalued function: a rule assigning to each z in a set a finite or infinite set of values, such as Lnz or z1/2 (see the elementary functions).
Fix z0 and one of the values f(z0). Along a path γ starting at z0 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 γ: the continuation of the chosen value along γ. The question is what happens when γ closes up and returns to z0.
A point z∗ is a branch point of f if continuation once around some sufficiently small circle about z∗ returns a value of f 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 f; 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)
Logarithm.Lnz=log∣z∣+i(θ+2kπ), k∈Z, has a single branch point, at z=0, and it is of infinite order: each counterclockwise circuit of the origin adds 2πi. One cut from 0 to ∞, along any simple arc, makes Ln single-valued. With the cut along the negative real axis, the principal branch is
lnz=log∣z∣+iArgz,Argz∈(−π,π),
which is analytic on C∖(−∞,0] with
dzdlnz=z1.
The open interval here is deliberate, and differs by one endpoint from the convention of the elementary functions, where Argz∈(−π,π] so that Arg and ln are defined at every z=0, the negative real axis included. The cut plane C∖(−∞,0] is where ln is analytic, and it omits that axis; on the axis itself ln has the value log∣z∣+iπ 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.
Powers.zp=exp(pLnz) for p∈C is single-valued if p∈Z; it has exactly q branches, and a branch point of order q at z=0, if p=s/q is rational in lowest terms with q∈N; and it has infinitely many branches otherwise. In particular z1/2 has a branch point of order 2 at the origin and two branches, with (z1/2)2=z on each.
No branch survives on an annulus. Let 0≤ρ1<ρ2≤∞. There is no single-valued continuous branch of Ln(z−a) on the annulus ρ1<∣z−a∣<ρ2 — hence none that is analytic, and none on a punctured disk about a either. The reason is that such a branch would deliver a continuous argument along a circle ∣z−a∣=ρ with ρ1<ρ<ρ2, and that circle is exactly the loop on which part 1 shows the argument to increase by 2π.
This is Ahlfors (Ahlfors (1979), Ch. 4, §1.3, p. 108)'s observation. He deduces it from the value 2πi of the integral of 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π)), 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 Lnz, since Ln(1/z)=−Lnz has a branch point at the origin; we follow him in treating ∞ on the same footing when deciding where a cut may end. Note the consequence for part 1: the cut from 0 to ∞ 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], set
fn(reiθ)=logr+i(θ+2nπ),−π<θ<π,r>0,n∈Z.
Each fn is single-valued and analytic on the cut plane, f0=ln is the principal branch, and fn reproduces on the cut plane the values Ln would take for (2n−1)π<θ<(2n+1)π. Gluing the upper edge of the cut in the fn copy to the lower edge in the fn+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 ew is unchanged by w↦w+2πi, inverting it can only ever return an answer modulo 2πi, and which representative continuation hands back depends on how the argument of z has moved. So the bookkeeping reduces to tracking arg along the path: a loop that encircles the origin advances argz by 2π whether we like it or not, and a loop that does not encircle it returns argz to its starting value. For zp the shift 2πi in Lnz becomes the factor e2πip, which is 1 exactly when p is an integer and of finite order q exactly when p has denominator q. Analyticity of the branch, finally, is the Cauchy-Riemann test applied in polar coordinates.
1. The logarithm. Let z0=r0eiθ0 with r0>0, and continue a chosen value logr0+iθ0 of Ln along the circle z(t)=r0ei(θ0+t), 0≤t≤2π. The function t↦logr0+i(θ0+t) is continuous and is at each t a value of Lnz(t), so it is the continuation; at t=2π it equals the starting value plus 2πi. Since 2πi=0 and n⋅2πi=0 for every n∈N, the origin is a branch point of infinite order. For z∗=0, a circle about z∗ of radius less than ∣z∗∣ lies in a half-plane not containing 0, 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 off the cut. On C∖(−∞,0] write u=logr and v=Argz with r=x2+y2. Then
ux=r2x,uy=r2y,vx=−r2y,vy=r2x,
the formulas for v coming from differentiating tanv=y/x (or cotv=x/y near the imaginary axis) and holding on the cut plane because Arg is there a single-valued smooth function. These four functions are continuous on C∖(−∞,0] and satisfy ux=vy and uy=−vx, so by the sufficiency half of the Cauchy-Riemann theoremln is analytic there, with
dzdlnz=ux+ivx=r2x−iy=zzˉzˉ=z1.
2. Powers. By definition the values of zp at a fixed z=0 are exp(pℓ) as ℓ runs over the values of Lnz; two of them differ by a factor e2πikp, k∈Z. One counterclockwise circuit of the origin replaces ℓ by ℓ+2πi, hence multiplies the continued value by e2πip, and n circuits multiply it by e2πinp. This factor is 1 if and only if np∈Z. So:
if p∈Z, already one circuit changes nothing and zp is single-valued;
if p=s/q in lowest terms with q∈N, then np∈Z first at n=q (because gcd(s,q)=1), so the branch point at 0 has order q, and the distinct values at a fixed z are the q numbers e2πiks/qexp(plnz), k=0,…,q−1;
otherwise np∈/Z for every n∈N, the order is infinite, and the values e2πikp are distinct for distinct k.
For p=1/2 this gives two branches, ±exp(21lnz), and squaring either gives exp(lnz)=z.
3. No branch on an annulus. Fix ρ with ρ1<ρ<ρ2, so that the circle ∣z−a∣=ρ lies in the annulus, and suppose a single-valued continuous L existed there with eL(z)=z−a at every point — which is what "a branch of Ln(z−a)" means. Set
φ(t)=ImL(a+ρeit),0≤t≤2π.
φ is continuous, being the imaginary part of a continuous function of t.
Exponentiating the defining property at z=a+ρeit gives eL=ρeit, so eReL=ρ and φ(t) is an argument of eit. Two arguments of the same number differ by an integer multiple of 2π, so φ(t)−t∈2πZ for every t.
The function t↦φ(t)−t is therefore a continuous map of the interval [0,2π] into the set 2πZ. Its image is a connected subset of 2πZ, hence a single point: φ(t)=t+2πk for one fixed k∈Z.
Consequently φ(2π)−φ(0)=2π=0. But a+ρe2πi and a+ρei0 are the same point of the annulus, and L was assumed single-valued, so φ(2π)=φ(0). The contradiction rules out such an L. Only continuity was used, so the analytic case is covered as well, and the argument applies verbatim to a punctured disk (ρ1=0) and to the punctured plane. □
For f(z)=∏j=1m(z−aj)pj, and more generally for any function built from logarithms and powers of polynomials:
List the candidate branch points. These are the points aj where a factor vanishes (and, after the substitution z=1/t, possibly ∞).
Write each factor in polar form about its own branch point:z−aj=rjeiθj, so that
f(z)=(j∏rjpj)exp(ij∑pjθj)
on a chosen branch (for real pj).
Take each closed loop in turn and record how each θj changes: by 2π if the loop encircles aj counterclockwise, by 0 if not. The continued value is multiplied by exp(2πi∑j∈Spj), where S is the set of encircled branch points.
A loop is harmless when that factor is 1, that is when ∑j∈Spj∈Z. Only the harmful loops must be blocked.
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 S of branch points is harmless when encircled together, a single finite cut joining them is enough, and is often the most convenient choice.
Fix the branch by naming the value of f at one point of the cut plane, or equivalently the range of each θ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; the function, the construction in step 4 and all the numbers are ours.
Let f(z)=((z−1)(z−4))1/2, with p1=p2=21 at a1=1, a2=4.
Candidate branch points:z=1 and z=4. At infinity, putting z=1/t gives f=t−1((1−t)(1−4t))1/2, whose second factor is single-valued near t=0; so ∞ is not a branch point.
Loop test. A loop around z=1 alone, or around z=4 alone, multiplies f by e2πi/2=−1: harmful. A loop around both multiplies it by e2πi(1/2+1/2)=1: harmless. A loop around neither does nothing.
Cut. It therefore suffices to block the two single loops while allowing the double loop, and the finite segment [1,4] does exactly that. (Two cuts (−∞,1] and [4,∞) would also work but forbid more than necessary.)
A concrete branch, built so that its analyticity is visible. The naive recipe z−1z−4 with two principal square roots is cut along (−∞,4], which is more than we want, so we build the branch out of a single principal logarithm instead. Let
m(z)=z−4z−1,F(z)=(z−4)exp(21lnm(z))on C∖[1,4],
with ln the principal logarithm, analytic on C∖(−∞,0] by part 1 of the statement. Three observations finish the construction.
m omits the cut of ln.m is a Möbius map of the sphere, so it is injective; it sends 1↦0, 4↦∞, and a real x with 1<x<4 to (x−1)/(x−4)<0. So m carries the segment [1,4] onto (−∞,0]∪{∞} and, being injective, carries nothing else there. For z∈/[1,4] we therefore have m(z)∈C∖(−∞,0], which is where ln lives.
Hence F is analytic on all of C∖[1,4], being a product of the polynomial z−4 with a composition of analytic functions. No question of continuity across anything arises: the formula is analytic where it stands.
F is a branch of f.exp(21lnw)2=exp(lnw)=w, so F(z)2=(z−4)2m(z)=(z−4)(z−1).
And Fis the naive product after all: evaluating both at every point of a grid over −8≤Rez≤12, ∣Imz∣≤6 (about 48000 points, the segment [1,4] excluded) gives F(z)=z−1z−4 to machine precision everywhere, the two sign jumps of the principal roots on (−∞,1) cancelling each other. The point of building F 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.
Numbers.F(0)=(−4)exp(21ln41)=−4⋅21=−2, and indeed F(0)2=4=(0−1)(0−4). F(10)=6exp(21ln23)=63/2=36≈7.348469, consistent with F(z)≈z−25 for large z (at z=10 that estimate gives 7.5). Across the cut, m(2.5+iε) tends to −1 through the lower half-plane and m(2.5−iε) through the upper half-plane, so 21lnm→∓iπ/2 and the two limits of F at z=2.5 are (−1.5)(∓i)=±1.5i: equal and opposite, as a sign-changing pair of branch points requires. Note that this branch is negative for real z<1 and positive for real z>4; a single-valued branch need not have a single sign.
Example(Values and orders)Original
Order three.z1/3 has a branch point of order 3 at the origin: one counterclockwise circuit multiplies a value by e2πi/3, two circuits by e4πi/3, three by 1. Starting from the value 2 at z=8 and going once round, we return with 2e2πi/3=−1+i3, whose cube is 8 as well.
A logarithm with two cuts, or one.g(z)=lnz+1z−1 has branch points at ±1, each of infinite order, since a loop around +1 alone adds 2πi and a loop around −1 alone subtracts it. A loop around both adds 0, so the finite cut [−1,1] suffices; this is the cut taken for Legendre's functions of the second kind Qn, the first of which is Q0(z)=21lnz−1z+1.
Principal values on the cut. With the convention Argz∈(−π,π], the principal ln is defined at z=−4, with ln(−4)=log4+iπ≈1.386294+3.141593i, but it is not continuous there: ln(−4−iε)→log4−iπ. The cut is the half-line where the function must be regarded as having two boundary values, differing by 2πi.
ApplicationOriginal
Branch cuts are where the physics lives, not a technicality. Three examples.
Propagating versus evanescent waves. A wave number k=ω2/c2−κ2, as a function of frequency, is a two-valued function with branch points at ω=±cκ: 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 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πip on part of the plane. The example above, where z−1z−4 turned out to be analytic across (−∞,1) although each factor is not, is the kind of check worth making before trusting such a formula.
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.
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]→C, g=g1+ig2 with g1,g2 real, define
∫abg(t)dt=∫abg1(t)dt+i∫abg2(t)dt.
A path (or contour) γ is a piecewise continuously differentiable map z:[a,b]→C; it is closed if z(a)=z(b), and simple if it does not meet itself otherwise. For f continuous on the image of γ,
∫γf(z)dz=∫abf(z(t))z′(t)dt,
the interval being subdivided at the finitely many points where z′ jumps. Writing f=u+iv and z=x+iy this is the pair of real line integrals
∫γf(z)dz=∫γ(udx−vdy)+i∫γ(udy+vdx),
which is how it is computed in practice. It is also the limit of the Riemann sums
j∑f(ζj)(zj−zj−1)
over refining subdivisions a=t0<⋯<tN=b, with zj=z(tj) and ζj a point of the arc between zj−1 and zj; 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 C1 path that is a small theorem rather than a restatement of the definition, and one line proves it: writing zj−zj−1=∫tj−1tjz′dt turns the difference between the sum and the integral into ∑j∫tj−1tj(f(ζj)−f(z(t)))z′(t)dt, whose modulus is at most L(γ) times the largest oscillation of f∘z over a subinterval, and that tends to 0 with the mesh because f∘z is uniformly continuous on [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 γ are used below. With respect to arc length,
∫γf∣dz∣=∫abf(z(t))∣z′(t)∣dt,soL(γ)=∫γ∣dz∣
is the length of γ. 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 ∮γ 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 f be continuous on the path γ, with z:[a,b]→C.
Linearity and subdivision.∫γ(αf+βg)dz=α∫γfdz+β∫γgdz for α,β∈C, and if γ is split into consecutive subarcs γ1,…,γn then ∫γfdz=∑k∫γkfdz.
Invariance and reversal. If t=φ(τ) is a piecewise C1 increasing bijection of [α,β] onto [a,b], then z∘φ gives the same integral: the contour integral depends on γ as an oriented curve, not on its parametrization. Reversing the orientation changes the sign, ∫−γfdz=−∫γfdz, whereas ∫−γf∣dz∣=∫γf∣dz∣.
The ML estimate. If ∣f(z)∣≤M at every point of γ and L=L(γ), then
∫γf(z)dz≤∫γ∣f∣∣dz∣≤ML.
Fundamental theorem. Let Ω be an open set. The following are equivalent for f continuous on Ω:
(a) ∫γfdz depends only on the endpoints of γ, for all paths γ in Ω;
(b) ∮γfdz=0 for every closed path γ in Ω;
(c) f=F′ for some analytic F on Ω (a primitive, or indefinite integral, of f).
When they hold, ∫γfdz=F(z(b))−F(z(a)). If Ω 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) ⇒ (c) handles this by working one component at a time.
Powers. For every integer n≥0 and every closed path γ, ∮γ(z−a)ndz=0. The same holds for n≤−2 and every closed path avoiding a. For n=−1 it fails: on the positively oriented circle C of radius ρ about a,
∮Cz−adz=2π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∣ 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 "fdz is exact" says that the potential satisfies the Cauchy-Riemann equations, so the potential is analytic with derivative f. Part 5 is then a corollary for n=−1, since (z−a)n has the obvious primitive; and the single exception n=−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 ∫abcgdt=c∫abgdt for a complex constant c=α+iβ: writing g=g1+ig2, both sides equal ∫(αg1−βg2)dt+i∫(αg2+βg1)dt.
2. By the real substitution rule applied to each component,
∫αβf(z(φ(τ)))z′(φ(τ))φ′(τ)dτ=∫abf(z(t))z′(t)dt,
and the chain rule identifies z′(φ(τ))φ′(τ) as the derivative of z∘φ, so the left side is the integral over the reparametrized path. For −γ, parametrized by w(t)=z(−t) on [−b,−a], three factors of −1 appear and their product is what we want:
one sign from the derivative w′(t)=−z′(−t), one from dt=−ds under s=−t, and one from the exchange of the limits b,a into a,b: an odd number, hence ∫−γfdz=−∫γfdz. With ∣dz∣ the first of the three disappears, since ∣w′(t)∣=∣z′(−t)∣ carries no sign, and the remaining two cancel, leaving ∫−γf∣dz∣=∫γf∣dz∣.
3. First, for continuous g on [a,b] with a≤b we show ∫abgdt≤∫ab∣g∣dt. If ∫abgdt=0 there is nothing to prove. Otherwise put θ=arg∫abgdt and c=e−iθ; then by part 1,
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 Rew≤∣w∣ pointwise and ∣c∣=1. Applying this to g(t)=f(z(t))z′(t) gives the first inequality of part 3, since ∣g∣=∣f(z(t))∣∣z′(t)∣ is the integrand of ∫γ∣f∣∣dz∣. For the second, use the hypothesis ∣f∣≤M on γ:
∫γ∣f∣∣dz∣≤M∫γ∣dz∣=ML.
4. (a) ⇔ (b). If the integral depends only on the endpoints and γ is closed, then γ and −γ have the same endpoints, so ∫γfdz=∫−γfdz=−∫γfdz, giving ∮γfdz=0. Conversely, if all closed integrals vanish and γ1,γ2 share their endpoints, then γ1 followed by −γ2 is closed, so by part 1 the two integrals agree.
(c) ⇒ (a). If f=F′ with F analytic, then on each smooth piece t↦F(z(t)) is differentiable with derivative 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 ∫γfdz=F(z(b))−F(z(a)), which depends only on the endpoints.
(a) ⇒ (c). Fix a point z0 in a component of Ω and set F(z)=∫σfdz along any path σ in Ω from z0 to z; by (a) this is well defined. To compute ∂F/∂x at z, take the last piece of σ to be a horizontal segment, which is possible because Ω is open and so contains a disk about z: with x as parameter, F(z)=const+∫xf(s+iy)ds, so ∂F/∂x=f(z) by the fundamental theorem of calculus (using continuity of f). Choosing the last piece vertical instead gives ∂F/∂y=if(z). Hence Fx=−iFy, which is the Cauchy-Riemann condition in combined form, and Fx=f is continuous; so by the sufficiency half of the Cauchy-Riemann theorem, F is analytic with F′=Fx=f.
5. For n≥0, (z−a)n=dzdn+1(z−a)n+1 with the primitive analytic on all of C, so (c) ⇒ (b) of part 4 gives 0. For n≤−2 the same primitive (z−a)n+1/(n+1) is analytic on C∖{a}, which contains γ, so again the integral is 0. For n=−1 the conclusion fails, and the direct computation on z(t)=a+ρeit, 0≤t≤2π, shows it:
∮Cz−adz=∫02πρeitiρeitdt=∫02πidt=2πi,
independent of ρ and of a. Read through part 4 in the direction (b) ⇒ (c) this says something more: 1/(z−a) has no primitive on any open set containing a circle about a, since a primitive would force the integral to vanish. □
RemarkOriginal
Where the hypotheses are used.
Piecewise C1 makes 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 f on γ is what makes the integrand integrable. Nothing in parts 1 to 3 requires f to be analytic, or even to be defined off γ: the ML estimate is a statement about continuous functions.
Analyticity of F, not just of f, is the content of part 4(c). It is the primitive that must be single-valued and analytic on all of Ω. The failure for n=−1 in part 5 is a failure of single-valuedness, not of differentiability: ln(z−a) has derivative 1/(z−a) on any cut plane, but no cut plane contains a circle around a.
Path dependence is the generic case. For f continuous but not analytic the integral usually does depend on the path, as the second example below shows; the discovery that for analytic f 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
Using a primitive. Let γ be the arc of the unit circle from 1 to i traversed counterclockwise. Since z2 has the entire primitive z3/3, part 4 gives
∫γz2dz=3i3−13=3−1−i≈−0.333333−0.333333i,
with no parametrization needed. Direct evaluation confirms it: ∫0π/2e2itieitdt=[e3it/3]0π/2=(e3iπ/2−1)/3=(−i−1)/3.
Path dependence of a non-analytic integrand. Integrate f(z)=zˉ from 0 to 1+i.
Along the straight segment z=t(1+i), 0≤t≤1: zˉ=t(1−i), dz=(1+i)dt, and (1−i)(1+i)=2, so the integral is ∫012tdt=1.
Along the two legs 0→1→1+i: on the first, z=t gives ∫01tdt=21; on the second, z=1+it gives ∫01(1−it)idt=i−i2∫01tdt=i+21. The total is 1+i.
The two values differ, as they must: zˉ satisfies no Cauchy-Riemann equation anywhere. Equivalently, around the closed triangle the integral is (1+i)−1=i=0.
The ML estimate, attained and applied. On the unit circle ∣zˉ∣=1 and L=2π, so ∣∮zˉdz∣≤2π. The exact value is ∮∣z∣=1zˉdz=∫02πe−itieitdt=2πi, of modulus exactly 2π: here the bound is sharp. Used the other way, on the circle ∣z∣=R with R>1 we have ∣z4+1∣≥R4−1 by the triangle inequality, so
∮∣z∣=Rz4+1dz≤R4−12πR,
which is 0.235619 at R=3 and 0.006284 at R=10, and tends to 0 as R→∞. 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) and velocity field v=(u1,u2), the combination w′(z)=u1+iu2 packages the velocity, and for a closed contour γ
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 lnz and whose complex velocity is accordingly a multiple of 1/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 and the flux of 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.
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.
A region here is a non-empty connected open subset of C. Call Ωsimply connected when each closed curve lying in it can be contracted, continuously and without ever leaving Ω, down to a single point of Ω — informally, when Ω 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∣<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 f be analytic on a simply connected region Ω and let γ be any closed path in Ω. Then
∮γf(z)dz=0.
Two special cases carry the proof and are worth stating separately.
Rectangle (Goursat). If R is a closed rectangle with sides parallel to the axes and f is analytic on an open set containing R, then ∮∂Rf(z)dz=0.
Disk. If f is analytic on an open disk Δ, then ∮γf(z)dz=0 for every closed path γ in Δ; equivalently, f has a primitive on Δ.
No hypothesis beyond analyticity of f 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 u and v are continuous (equivalently that 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.
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.
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 u and v 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.
Two paths. If γ1 and γ2 are paths in Ω with the same initial and the same final point, and the region enclosed between them lies in Ω, then ∫γ1fdz=∫γ2fdz. "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 followed by the reverse of γ2 is a simple closed curve. We read the hypothesis in the form the proof actually uses: γ1 followed by the reverse of γ2 can be contracted to a point inside Ω, which is what sweeping across an enclosed region in Ω accomplishes when there is one.
Deformation. If C and σ are positively oriented simple closed contours with σ inside C, and f is analytic on the closed region between them, then
∮Cf(z)dz=∮σf(z)dz.
More generally, the value of ∮γfdz is unchanged by any continuous deformation of γ within the region of analyticity, which precludes moving the contour across a point where f fails to be analytic.
Several holes. If C is a positively oriented simple closed contour and σ1,…,σm are positively oriented simple closed contours inside C, mutually exterior to one another, and f is analytic on the closed region inside C and outside all the σk, then
∮Cf(z)dz=k=1∑m∮σkf(z)dz.
Idea of the proofOriginal
Three ideas, in increasing subtlety.
For the simply connected statement under the extra hypothesis that f′ is continuous, there is nothing to invent: fdz 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∗. Near z∗ the function agrees with its own affine approximation to within ε∣z−z∗∣, for any ε we fix in advance and all z close enough; and the affine part integrates to zero around any closed contour because 1 and z have primitives. The error contributes at most ε⋅(diameter)⋅(perimeter), and under bisection diameter and perimeter each halve, so this bound decays like 4−n — at exactly the rate at which the bisection inflates the integral. The two rates cancel and leave ∣∮∣≤εdL for every ε.
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.
Assume for this proof that u and v have continuous first partial derivatives, and let γ be a positively oriented simple closed contour in Ω bounding the region A. We need A⊆Ω. That holds because Ω 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⊆Ω as a hypothesis: it is satisfied in every application below.
Step 2 (Green's theorem in the plane). Both are real line integrals of the form ∮γ(Pdx+Qdy) with P,Q having continuous first partials on A and its boundary, so Stokes' theorem in the plane applies to each:
∮γ(Pdx+Qdy)=∬A(∂x∂Q−∂y∂P)dxdy.
Step 3 (the integrands vanish). For the first integral take P=u, Q=−v; for the second, P=v, Q=u. The two area integrands are
−vx−uyandux−vy,
and both are identically zero on A by the Cauchy-Riemann equations, which hold at every point of Ω because f is analytic there. Hence ∮γfdz=0.
Steps 1 to 3 prove the theorem for a simple closed contour in Ω, which is the case needed in everything that follows. □
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′ is assumed. For a rectangle S with sides parallel to the axes write η(S)=∮∂Sfdz, the boundary positively oriented.
Step 1 (bisection does not lose everything). Divide S into four congruent rectangles S(1),…,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)).
By the triangle inequality at least one k has ∣η(S(k))∣≥41∣η(S)∣; choose one such (say the first in a fixed ordering) and call it R1.
Step 2 (nested rectangles). Repeating, we get R⊃R1⊃R2⊃⋯ with
∣η(Rn)∣≥4−n∣η(R)∣,
and with diameter dn=2−nd and perimeter Ln=2−nL, where d and L are those of R. The Rn are nested non-empty compact sets with diameters tending to 0, so by compactness they have exactly one common point z∗, and every neighbourhood of z∗ contains all but finitely many Rn.
Step 3 (use differentiability at the one point z∗). Let ε>0. Since f is analytic on an open set containing R, it is differentiable at z∗, so there is δ>0 with
Step 4 (the affine part integrates to zero). The functions 1 and z have the primitives z and z2/2, entire; so by part 5 of the properties of contour integrals, ∮∂Rndz=0 and ∮∂Rnzdz=0. Therefore the constant f(z∗)−z∗f′(z∗) and the linear term zf′(z∗) contribute nothing and
η(Rn)=∮∂Rn[f(z)−f(z∗)−(z−z∗)f′(z∗)]dz.
Step 5 (estimate and compare the two rates). By the ML estimate in its arc-length form, and because ∣z−z∗∣≤dn on ∂Rn (both points lie in Rn),
∣η(Rn)∣≤ε∮∂Rn∣z−z∗∣∣dz∣≤εdnLn=ε4−ndL.
Combining with Step 2, 4−n∣η(R)∣≤ε4−ndL, so ∣η(R)∣≤εdL. As ε>0 was arbitrary and d,L are fixed, η(R)=0. □
Proof(Part 2: from a rectangle to a disk)AdaptedAdapted from Ahlfors (1979), Ch. 4, §1.5, Thm 4, p. 113
Let f be analytic on the open disk Δ with centre z0=x0+iy0, and define, for z=x+iy∈Δ,
F(z)=∫σzf(ζ)dζ,
where σz is the horizontal segment from (x0,y0) to (x,y0) followed by the vertical segment from (x,y0) to (x,y). Both segments lie in Δ, because a disk contains the whole axis-parallel rectangle determined by its centre and any of its points.
Step 1 (∂F/∂y). Only the second segment depends on y, and on it ζ=x+is with dζ=ids, so F(z)=c(x)+i∫y0yf(x+is)ds. Since f is continuous, the fundamental theorem of calculus gives ∂F/∂y=if(z).
Step 2 (∂F/∂x). Let σ~z be the other staircase: vertical from (x0,y0) to (x0,y), then horizontal to (x,y). Together, σz and the reverse of σ~z traverse the boundary of the rectangle with opposite corners z0 and z, which lies in Δ; by part 1 that boundary integral is 0, so ∫σ~zf=∫σzf=F(z). Now only the final horizontal segment of σ~z depends on x, and there ζ=s+iy with dζ=ds, so F(z)=c~(y)+∫x0xf(s+iy)ds and ∂F/∂x=f(z).
Step 3 (F is analytic with F′=f). From Steps 1 and 2, Fx=f=−i(if)=−iFy, which is the Cauchy-Riemann condition in the combined form of the Cauchy-Riemann theorem; and Fx=f is continuous. Hence F is analytic on Δ with F′=f.
Step 4 (conclusion).f has a primitive on Δ, so by the fundamental theorem for contour integrals its integral over every closed path in Δ vanishes. □
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 γ be γ1 followed by the reverse of γ2; since the two paths share both endpoints, γ is a closed path in Ω. By the hypothesis as read in the statement, γ can be contracted to a point within Ω. A closed path contractible in Ω has winding number zero about every point of the complement of Ω, since the winding number is unchanged under a homotopy that never crosses the point and is 0 for a constant path; in Ahlfors (Ahlfors (1979), §4.3, Definition 2, p. 141)'s language γ is then homologous to zero in Ω, 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 ∮γfdz=0, that is ∫γ1fdz−∫γ2fdz=0 by the reversal and subdivision rules for contour integrals.
Note what is and is not assumed here: γ need not be simple and Ω 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 γ happens to be a simple closed curve whose inside lies in Ω, or to lie inside a disk contained in Ω, 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 C and σ along two nearby arcs λ and λ′ joining a point of σ to a point of C. Travelling along C in the positive sense, then inward along λ′, then around σ in the negative sense, then outward along λ, produces a single closed contour Γ which does not enter the hole; the region it encloses is simply connected and contained in the region of analyticity, so ∮Γfdz=0. Now let λ and λ′ approach each other. Their contributions tend to the integrals of f along the same arc in opposite directions, which cancel, because f is continuous across the cut — unlike a branch cut, this barrier is of our own making. What is left is
∮Cfdz−∮σfdz=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 Ω; each such step changes nothing by part 1.
Part 3 (several holes). Insert one cut from C to each σk. The same construction yields one closed contour enclosing a simply connected region of analyticity, and the cuts cancel in pairs, leaving ∮Cfdz−∑k∮σkfdz=0. □
RemarkOriginal
Where the hypotheses are used, and what fails without them.
Simple connectivity is indispensable. On the region Ω={0<∣z∣<2}, which is not simply connected, f(z)=1/z is analytic, yet ∮∣z∣=1dz/z=2πi=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/z is analytic on the disk ∣z−1∣<1, which is simply connected, so every closed path inside that disk integrates to zero. Conversely the circle ∣z∣=1 lies in no simply connected subregion of Ω on which 1/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 although 1/z is perfectly well behaved on the circle itself. The theorem is a statement about what the contour surrounds: if f is analytic on an open set containing a simple closed contour and its inside, the integral does vanish, and the 1/z example escapes only because the enclosed point 0 is excluded.
Continuity of f′ is used only in the Green's theorem proof, at Step 2, where Stokes' theorem in the plane needs continuously differentiable P and Q. 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′ 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 f 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 are traversed clockwise, as they are when one regards Γ as a single boundary, the signs flip and the total is zero.
Example(Deforming a square to a circle)Original
Evaluate I=∮Sdz/z, where S is the boundary of the square with vertices ±2±2i, positively oriented.
Choose the deformation.1/z is analytic on C∖{0}, and the closed region between S and the unit circle ∣z∣=1 contains no singularity, so part 2 applies.
Conclusion.I=2πi≈6.283185i. Direct numerical integration along the four sides of the square gives 6.283185i to six digits, which is the check.
Contrast. For the square S′ with vertices 3±i and 5±i, which does not enclose the origin, 1/z is analytic on a simply connected region containing S′ and its interior, so ∮S′dz/z=0.
Example(A contour with two holes)Original
Let g(z)=z2−11=21(z−11−z+11) and let C be the circle ∣z∣=3, positively oriented.
Isolate the holes.g is analytic except at z=±1. Take σ± to be the positively oriented circles of radius 21 about ±1; they lie inside C and are mutually exterior, and g is analytic on the closed region inside C and outside both.
Apply part 3.∮Cg=∮σ+g+∮σ−g.
Evaluate each. On σ+, the term −21/(z+1) is analytic inside, so by Cauchy-Goursat it integrates to zero, and ∮σ+g=21∮σ+z−1dz=21⋅2πi=πi. Symmetrically ∮σ−g=−21⋅2πi=−πi.
Total.∮Cgdz=πi−πi=0.
Check by the ML estimate. On ∣z∣=R we have ∣z2−1∣≥R2−1, so ∣∮∣z∣=Rg∣≤2πR/(R2−1), which tends to 0 as R→∞; by part 2 the integral is the same for every R>1, so it must be 0 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 κ at the origin is κ/(2πiz). The circulation around any loop enclosing the vortex is κ, 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).
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.
Let C be a positively oriented simple closed contour, let f be analytic on and inside C, and let z0 be a point inside C. Then
f(z0)=2πi1∮Cz−z0f(z)dz,
and for every n∈N the derivative f(n)(z0) exists and
f(n)(z0)=2πin!∮C(z−z0)n+1f(z)dz.
Consequences.
Infinite differentiability. A function analytic on an open set has derivatives of all orders there, and every one of them is again analytic.
Point outside. If instead z0 lies outside C, the whole integrand is analytic on and inside C, so ∮Cf(z)(z−z0)−1dz=0.
Mean value property. If the closed disk ∣z−z0∣≤r lies in the region of analyticity, then
f(z0)=2π1∫02πf(z0+reiθ)dθ.
Morera's theorem (converse of Cauchy's theorem). If f is continuous on a region Ω and ∮γf(z)dz=0 for every closed contour γ in Ω, then f is analytic on Ω.
Idea of the proofOriginal
The formula is Cauchy's theorem applied to the right function. The quotient (f(z)−f(z0))/(z−z0) is analytic except at z0, where it stays bounded, so its integral around C vanishes; splitting the numerator leaves ∮f(z)/(z−z0)dz on one side and f(z0)∮dz/(z−z0)=2πif(z0) on the other. Equivalently, and this is the version easiest to compute with: deform C to a tiny circle about z0, on which f is nearly the constant f(z0), and the integral is then f(z0) times 2π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ζ, with φ merely continuous on the contour, is a differentiable function of z 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 φ. The payoff is disproportionate: since the right-hand side makes sense for all n, so does the left, and analyticity at once implies infinite differentiability — a conclusion with no counterpart for functions of a real variable.
Step 1 (reduce to a small circle). Let σρ be the positively oriented circle ∣z−z0∣=ρ, with ρ small enough that the closed disk it bounds lies inside C. The integrand f(z)/(z−z0) is analytic on the closed region between C and σρ, because f is analytic there and the only zero of the denominator is z0, which is not in that region. By the deformation principle,
∮Cz−z0f(z)dz=∮σρz−z0f(z)dzfor every such ρ.
Step 2 (evaluate on the circle). Parametrize σρ by z=z0+ρeiθ, 0≤θ≤2π, so dz=iρeiθdθ and z−z0=ρeiθ. The factors of ρeiθ cancel:
∮σρz−z0f(z)dz=i∫02πf(z0+ρeiθ)dθ.
Dividing by 2πi gives the mean value of f 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π,
by the modulus inequality for integrals. Since f is continuous at z0 — it is differentiable there — the right-hand side tends to 0 as ρ→0+. But by Steps 1 and 2 the left-hand side is 2πi1∮Cf(z)(z−z0)−1dz−f(z0), a quantity independent of ρ. A constant that is smaller than every positive number is zero, which is the formula. □
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 γ be a path and let φ be continuous on γ. For n∈N and z not on γ put
Fn(z)=∫γ(ζ−z)nφ(ζ)dζ.
Then every Fn is analytic off γ, and Fn′=nFn+1.
Set-up for the proof of the lemma. Fix z1 off γ and pick δ>0 with ∣ζ−z1∣≥δ for all ζ∈γ, which is possible because γ is compact. Throughout, z is restricted to ∣z−z1∣<δ/2, so that ∣ζ−z∣>δ/2 on γ. Write K=∫γ∣φ∣∣dζ∣, a finite constant. All estimates below are the ML estimate of contour integrals.
Step 1 (F1 is continuous at z1). From the identity
and note that ψ is continuous on γ, since ζ=z1 there. The right-hand side is therefore the function "F1 built from ψ", which by Step 1 is continuous at z1; its value at z=z1 is ∫γφ(ζ)(ζ−z1)−2dζ=F2(z1). Hence F1′(z1)=F2(z1).
Step 3 (the induction identity). Let n≥2 and keep ψ(ζ)=φ(ζ)/(ζ−z1). Denote by Gk the integrals built from ψ, that is Gk(z)=∫γψ(ζ)(ζ−z)−kdζ. Dividing (∗) by (ζ−z)n−1 and rearranging gives
(ζ−z)n1=(ζ−z)n−1(ζ−z1)1+(ζ−z)n(ζ−z1)z−z1,
and integrating against φ and subtracting Fn(z1),
Here we used Gn−1(z1)=∫γφ(ζ)(ζ−z1)−ndζ=Fn(z1) and, for the last term, ∫γφ(ζ)(ζ−z)−n(ζ−z1)−1dζ=Gn(z).
Step 4 (induction). Suppose, as the induction hypothesis, that Fn−1 is continuous and differentiable at z1 with Fn−1′=(n−1)Fn, for every continuous density on γ — in particular for ψ, so Gn−1′(z1)=(n−1)Gn(z1). Steps 1 and 2 are this statement for n−1=1.
Continuity of Fn at z1. In (†) the bracket tends to 0 because Gn−1 is continuous at z1, and Gn is bounded on ∣z−z1∣<δ/2, since ∣ψ∣≤∣φ∣/δ on γ gives ∣Gn(z)∣≤(2/δ)nK/δ independently of z; so the second term tends to 0 as well. Since the density in (†) was arbitrary, this proves continuity of Fn at z1 for every continuous density, hence also continuity of Gn at z1.
Differentiability. Divide (†) by z−z1 and let z→z1. The first quotient tends to Gn−1′(z1)=(n−1)Gn(z1) by the induction hypothesis; the second term is Gn(z), which tends to Gn(z1) by the continuity just established. Therefore
where the third equality is the definition of Gn with ψ=φ/(ζ−z1). This closes the induction and proves the lemma.
Application to f. Let Ω be open, f analytic on Ω, and a∈Ω. Choose a positively oriented circle Ca about a whose closed disk lies in Ω. For every z inside Ca, the formula already proved gives
f(z)=2πi1∮Caζ−zf(ζ)dζ=2πi1F1(z),
with the continuous density φ=f on Ca. The lemma then gives, by induction on n,
f(n)(z)=2πin!∮Ca(ζ−z)n+1f(ζ)dζ(n=0,1,2,…),
since the case n=0 is the display above and, assuming the formula for n, the lemma gives f(n+1)(z)=2πin!Fn+1′(z)=2πin!(n+1)Fn+2(z)=2πi(n+1)!Fn+2(z). In particular every f(n) is differentiable, hence analytic, which is consequence 1.
Finally, for z0 inside the original contour C, the integrand f(ζ)(ζ−z0)−n−1 is analytic on the closed region between C and a small circle Cz0 about z0, so the deformation principle replaces Cz0 by C and the formula holds as stated. □
Let f be continuous on the region Ω with all closed contour integrals zero. By the fundamental theorem for contour integrals (contour integrals, part 4), f has a primitive: there is an analytic F on Ω with F′=f. By consequence 1, applied to F, the derivative F′=f is itself analytic on Ω. □
RemarkOriginal
Where the hypotheses are used, and why the order of results matters.
Analyticity inside, not only on, C. The proof uses the deformation principle on the region between C and a small circle, which requires f to be analytic throughout. Dropping it destroys the conclusion: take f(z)=1/z, which is analytic on the unit circle C but not inside it, and z0=21. Partial fractions give z(z−21)1=2(z−211−z1), so 2πi1∮Cz−211/zdz=2(1−1)=0, whereas f(21)=2. Numerical integration confirms the value 0.
z0 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 φ merely continuous, which is what allows the induction to be run on densities such as φ(ζ)/(ζ−z1)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′ is continuous, then consequence 1 — which produces continuity of 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 Ω 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 Ω, so the primitive it produces is single-valued whatever the topology of Ω, 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∣ is differentiable on R with derivative 2∣x∣, which is not differentiable at 0. One complex derivative forces infinitely many.
Example(Evaluating contour integrals with the formulas)Original
Order zero.∮∣z∣=2z−1ezdz=2πie1=2πie≈17.079468i, since ez is entire and z0=1 lies inside the circle. Numerical integration over the parametrized circle reproduces 17.079468i.
Order one.∮∣z∣=2(z−1)2ezdz=1!2πidzdezz=1=2πie as well — the same value here because ez is its own derivative.
Order two.∮∣z∣=1z3coszdz=2!2πidz2d2coszz=0=πi(−cos0)=−πi≈−3.141593i, again confirmed numerically.
Point outside.∮∣z∣=1z−3ezdz=0, by consequence 2.
The mean value property as a check. Take f(z)=1/(2−z), analytic on ∣z∣≤1. Then 2π1∫02π2−eiθdθ=f(0)=21. Numerical quadrature of the left side gives 0.500000, with imaginary part zero to the same accuracy.
Example(A function reconstructed from its boundary values)Original
Let f be analytic on ∣z∣≤1 with boundary values f(eiθ)=e2iθ. The formula fixes f on the whole disk. For ∣z0∣<1,
f(z0)=2πi1∮∣z∣=1z−z0z2dz
because z2 agrees with f on the contour, and the right-hand side equals z02 by the formula applied to the entire function z2. So f(z)=z2: the only analytic function on the disk with those boundary values. Checking at z0=21+21i: f(z0)=(21+21i)2=21i, and the integral 2πi1∫02πeiθ−z0e2iθieiθdθ evaluates numerically to 0.000000+0.500000i.
Note what happens if the boundary data are not the restriction of an analytic function. With f(eiθ)=e−iθ, that is zˉ on the circle, we may write zˉ=1/z there, and partial fractions give
for every z0 with 0<∣z0∣<1 (numerically 0.000000 at z0=21+21i). So if some analytic g on the closed disk had g=zˉ on the circle, the formula would force g≡0, contradicting ∣g∣=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 nth derivative formula turns differentiation into integration. Evaluating f(n)(z0)=2πrnn!∫02πf(z0+reiθ)e−inθdθ by the trapezoidal rule on m equally spaced points converges geometrically in m for analytic f 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) by the resolvent (zI−A)−1 of a matrix or operator A turns the formula into the definition f(A)=2πi1∮Cf(z)(zI−A)−1dz, with C enclosing the spectrum. This is how etA, 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.
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.
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
Cauchy's estimate. Let f be analytic on a region containing the closed disk ∣z−a∣≤r and put M(r)=max∣z−a∣=r∣f(z)∣. Then
f(n)(a)≤rnn!M(r),n=0,1,2,…
Liouville's theorem. A function analytic on the whole plane and bounded there is constant.
Fundamental theorem of algebra. A polynomial P of degree n≥1 with complex coefficients has a root; consequently it factors as P(z)=an∏k=1n(z−λk) and has exactly n roots counted with multiplicity.
Maximum modulus principle. Let f be analytic on a region (a connected open set) Ω.
(a) Global form. If there is z0∈Ω with ∣f(z)∣≤∣f(z0)∣ for every z∈Ω, then f is constant on Ω. Equivalently, if f is non-constant then ∣f∣ has no maximum in Ω.
(b) Local form. If ∣f∣ has a local maximum at z0∈Ω — that is, ∣f(z)∣≤∣f(z0)∣ for all z in some disk about z0 — then f is constant on a disk about z0. In particular ∣f∣ has no strict local maximum anywhere in Ω; 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∣ 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 f has a (non-strict) local maximum of ∣f∣ at every point. Form (b) is the conditional version — a local maximum is allowed, but only if f 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.
Boundary form. If E⊂C is closed and bounded, f is continuous on E and analytic on the interior of E, then maxE∣f∣ is attained at a point of the boundary of E. In particular, if ∣f∣≤M on the boundary then ∣f∣≤M throughout E.
Idea of the proofOriginal
Everything here comes from one inequality applied to one formula. Cauchy's integral formula writes f(n)(a) as an average of f over a circle against the kernel (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 r do not cancel — they leave r−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=1 with r→∞: if M(r) never grows, 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 z0 contains a disk around f(z0), which contains points of larger modulus. The analytic one, which is self-contained here: the integral formula on a circle says f(z0) is the average of f 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)∣ were maximal, ∣f∣ would be constant on every small circle about z0, hence near z0, and then f itself would be constant.
On C the integrand has modulus ∣f(z)∣/rn+1≤M(r)/rn+1, since ∣z−a∣=r exactly. The ML estimate (contour integrals) gives
f(n)(a)≤2πn!⋅rn+1M(r)⋅2πr=rnn!M(r).
Part 2. Let f be entire with ∣f∣≤M everywhere, and let a∈C. Then M(r)≤M for every r>0, so part 1 with n=1 gives
∣f′(a)∣≤rM(r)≤rMfor all r>0.
Here the hypothesis that f is entire is used twice: to know that everyr is allowed, and that the bound M holds on every circle. Letting r→∞ gives f′(a)=0. As a was arbitrary, f′≡0 on C, so ux=vx=0 and, by the Cauchy-Riemann equations, uy=vy=0 as well; a function with vanishing partial derivatives on the connected set C is constant.
Part 3. Let P(z)=anzn+⋯+a0 with n≥1 and an=0, and suppose P had no root. Then g=1/P is entire. For ∣z∣=R,
by the triangle inequality, so g(z)→0 as ∣z∣→∞. Choose R0 with ∣g∣≤1 for ∣z∣≥R0; on the compact disk ∣z∣≤R0 the continuous function ∣g∣ is bounded by compactness. So g is a bounded entire function, hence constant by part 2 — and the constant is 0 because g→0, which is impossible since g=1/P never vanishes. Therefore P has a root λ1.
Dividing, P(z)=(z−λ1)Q(z) with degQ=n−1: indeed P(z)−P(λ1)=P(z) is a linear combination of the terms zk−λ1k, each divisible by z−λ1. Repeating n times reduces the degree to 0 and gives the stated factorization, with the leading coefficient an; a product of linear factors vanishes exactly at the λk, so P has no further roots. □
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 f be analytic on the region Ω and suppose first only that ∣f∣ has a local maximum at z0∈Ω: there is r0>0 with ∣f(z)∣≤∣f(z0)∣ whenever ∣z−z0∣<r0, and we may shrink r0 so that this disk lies in Ω.
Step 1 (the mean value property). Choose 0<r<r0, so that the closed disk ∣z−z0∣≤r lies in Ω and in the disk where the maximality holds. By the mean value consequence of Cauchy's formula,
f(z0)=2π1∫02πf(z0+ρeiθ)dθfor every 0<ρ≤r,
and therefore, taking moduli,
∣f(z0)∣≤2π1∫02πf(z0+ρeiθ)dθ.(∗)
Step 2 (∣f∣ is constant on small circles). Put h(θ)=∣f(z0)∣−∣f(z0+ρeiθ)∣ for 0<ρ≤r. By the assumed local maximality, h≥0; and h is continuous, since f is. Rearranging (∗) gives 2π1∫02πh(θ)dθ≤0. A continuous non-negative function with non-positive integral vanishes identically, so ∣f∣=∣f(z0)∣ on every circle ∣z−z0∣=ρ with ρ≤r, hence on the whole disk D={∣z−z0∣≤r}.
Step 3 (constant modulus forces constant f; this proves form (b)). On D we have u2+v2=c with c=∣f(z0)∣2. If c=0 then f≡0 on D. Otherwise differentiate: uux+vvx=0 and uuy+vvy=0. Using the Cauchy-Riemann equations uy=−vx, vy=ux, the second becomes −uvx+vux=0. The two linear equations
uux+vvxvux−uvx=0,=0
in the unknowns (ux,vx) have determinant −(u2+v2)=−c=0, so ux=vx=0 on D, and then uy=vy=0 too. Hence f is constant on D.
Step 4 (form (a): spread the constancy over Ω). Assume now the global hypothesis: ∣f(z)∣≤∣f(z0)∣ for everyz∈Ω. Put c0=f(z0) and
A={z∈Ω:f(z)=c0}.
A=∅, since z0∈A.
A is closed in Ω: it is f−1({c0}) and f is continuous.
A is open. Let z∈A. Then ∣f(z)∣=∣c0∣=∣f(z0)∣, and by the global hypothesis ∣f∣≤∣f(z)∣ on all of Ω, so ∣f∣ certainly has a local maximum at z. By Steps 1 to 3 applied at z in place of z0, f is constant on a disk about z, and that constant is f(z)=c0; so this disk lies in A.
A non-empty subset of a connected set that is both open and closed in it is the whole set (connectedness), so A=Ω and f≡c0 on Ω.
Part 5.E is closed and bounded and ∣f∣ is continuous on it, so by compactness the maximum is attained, say at z0∈E. If z0∈∂E, we are done. Otherwise z0 lies in the interior of E; let Ω0 be the connected component of intE containing z0, a region. Since ∣f∣≤∣f(z0)∣ on all of E, in particular near z0, part 4 applies on Ω0 and f is constant there, of modulus ∣f(z0)∣.
Now ∂Ω0 is non-empty, because Ω0 is a non-empty bounded open set, and ∂Ω0⊆∂E: if some p∈∂Ω0 lay in intE, take a disk N⊆intE about p; as p is a limit of points of Ω0, the connected set N meets Ω0, so N∪Ω0 is a connected subset of intE and hence lies in the component Ω0; then p∈Ω0, which is open, contradicting p∈∂Ω0. And p∈/E is impossible because ∂Ω0⊆Ω0⊆E, E being closed. So p∈∂E.
By continuity of f on E, ∣f∣=∣f(z0)∣ on Ω0 and in particular at any point of ∂Ω0⊆∂E: the maximum is attained on the boundary of E after all. The final sentence of part 5 is this statement rearranged. □
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) is analytic and bounded on the disk ∣z∣<21 but not constant: the proof needs arbitrarily large circles, which only an entire function provides. Likewise eiz is entire and bounded on the closed upper half-plane, and is not constant.
Boundedness cannot be weakened to "bounded on the real axis".sinz is entire with ∣sinx∣≤1 for real x, yet ∣sin(10i)∣=sinh10≈11013.23. Part 2 demands a bound on the whole plane.
Polynomial growth gives polynomials. If f is entire with ∣f(z)∣≤A+B∣z∣p for some p≥0, then part 1 with n>p gives ∣f(n)(0)∣≤n!(A+Brp)/rn→0, so f(n)≡0 for every n>p (the same argument at every centre) and f is a polynomial of degree at most ⌊p⌋ — the floor, because n runs over integers and the smallest integer exceeding p is ⌊p⌋+1. Liouville is the case p=0.
Cauchy's estimate is rarely sharp, and the choice of r matters. The advice to choose r so as to minimise M(r)r−n, and the example of a function whose maximum modulus blows up like 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) on ∣z∣<2 we have M(r)=1/(2−r) and f(n)(0)=n!/2n+1 exactly. The bound is n!/(rn(2−r)), minimised at r=2n/(n+1): for n=1 this gives 1 against the true value 0.25, and for n=3 it gives 3.556 against 0.375. The estimate is an order-of-magnitude tool.
The maximum principle says nothing about the minimum unless f has no zero: ∣z∣ on the closed unit disk has its minimum at the interior point 0. Applying part 4 to 1/f gives a minimum principle wherever f=0.
"Region" means connected. On a disconnected open set f can be one constant on one piece and a larger constant on another, and then ∣f∣ 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 ef, which is analytic with ∣ef∣=eRef: since t↦et is strictly increasing, a maximum of Ref is a maximum of ∣ef∣, so there is none unless f is constant. For the imaginary part the right auxiliary function is e−if, which has ∣e−if∣=eImf; ef says nothing about Imf. Either way this is the maximum principle for harmonic functions (harmonic functions), and minima are covered by e−f and eif.
Example(Cauchy's estimate, Liouville's theorem and a growth bound)Original
A numerical estimate. Let f(z)=1/(2−z) and a=0. On ∣z∣=1, ∣2−z∣≥1, so M(1)=1 (attained at z=1) and part 1 gives ∣f′(0)∣≤1. The true value is f′(z)=(2−z)−2, so f′(0)=0.25: the bound holds with a factor of 4 to spare, as it must since ∣f∣ is far below M(1) on most of the circle.
A growth bound applied. Suppose f is entire with ∣f(z)∣≤5+2∣z∣3/2 for all z. Take n=2 in part 1 at any centre a: ∣f′′(a)∣≤2(5+2(∣a∣+r)3/2)/r2, and letting r→∞ with a fixed gives f′′(a)=0. So f′′≡0 and f(z)=αz+β. The conclusion "degree at most 1" cannot be improved to "constant", because one non-constant affine function does satisfy the hypothesis: f(z)=z has ∣z∣≤5+2∣z∣3/2 for every z — for ∣z∣≤5 because the left side is then at most 5, and for ∣z∣>5 because 2∣z∣3/2>25∣z∣>∣z∣. (Most affine functions do not satisfy it: f≡100 already fails at z=0, where the bound is 5. The hypothesis has fixed constants, so it is a genuine restriction on α and β; all that sharpness needs is one non-constant example.) Note what the growth exponent 23 buys: not a 23-power function, which is not entire, but a polynomial of degree ⌊23⌋=1.
An entire function cannot omit a disk. Suppose f is entire and its values avoid some disk: ∣f(z)−w0∣≥δ>0 for all z, with δ>0 fixed. Then g=1/(f−w0) is entire, since the denominator never vanishes, and ∣g∣≤1/δ everywhere. By part 2, g is constant, hence so is f. So the image of a non-constant entire function is dense in C. The hypothesis really is about omitting a disk and not a point: ez is entire and non-constant and its image is exactly C∖{0}, which omits one point but no disk. Concretely, ez=10−6 is solved by z=−6log10≈−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)
An exact maximum. Let f(z)=z2−2z on the closed unit disk E={∣z∣≤1}. By part 5 the maximum of ∣f∣ is on ∣z∣=1, where ∣f(z)∣=∣z∣∣z−2∣=∣z−2∣. For z=eiθ, ∣z−2∣2=5−4cosθ, maximal at θ=π. So maxE∣f∣=3, attained at z=−1 and nowhere else; a scan of 1000 points of the circle confirms the maximum 3.000000 at θ=π. At the interior point 0, ∣f(0)∣=0.
A bound obtained without computing anything. Let f(z)=ez and E=[0,1]×[0,1], the closed unit square. Then ∣ez∣=eRez≤e on E, with equality exactly on the boundary edge Rez=1: the maximum e≈2.718282 is on ∂E, as part 5 requires, and no interior point comes close to it.
A Schwarz-type estimate. Suppose g is analytic on an open set containing the closed disk ∣z∣≤2, with ∣g∣≤10 on the circle ∣z∣=2 and g(0)=0. The conclusion is that
∣g(z)∣≤5∣z∣ for ∣z∣≤2,∣g′(0)∣≤5,
which is Ahlfors (Ahlfors (1979), Thm 13 and the paragraph after it, p. 135)'s Schwarz lemma with his bound M set to 10 and his radius R to 2, so that M/R=5. His proof puts h(z)=g(z)/z with 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 0 never has to be discussed.
The derivative bound first. Part 1 with n=1, a=0, r=2 and M(2)≤10 gives ∣g′(0)∣≤1!⋅10/2=5 outright.
Boundary values of h=g/z on an annulus. Fix 0<ε<2 and let Eε={ε≤∣z∣≤2}. On Eε the function h(z)=g(z)/z is continuous and, on the interior, analytic — the denominator does not vanish anywhere on Eε, so no extension is needed. The boundary of Eε is the two circles ∣z∣=2, where ∣h∣=∣g∣/2≤5, and ∣z∣=ε, where ∣h∣≤Mg(ε)/ε with Mg(ε)=max∣z∣=ε∣g∣.
The inner circle is harmless as ε→0. Since g(0)=0, the difference quotient g(z)/z tends to g′(0) as z→0, so Mg(ε)/ε=max∣z∣=ε∣g(z)/z∣→∣g′(0)∣≤5. Hence for any η>0 there is εη with Mg(ε)/ε≤5+η for all ε<εη.
Conclude. By part 5 applied to h on Eε, ∣h∣≤max{5,Mg(ε)/ε}≤5+η on all of Eε, for every ε<εη. Any z with 0<∣z∣≤2 lies in Eε for all small enough ε, so ∣h(z)∣≤5+η; letting η→0 gives ∣h(z)∣≤5, that is ∣g(z)∣≤5∣z∣, for every such z, and at z=0 the inequality reads 0≤0.
Both bounds are attained: g(z)=5z has ∣g∣=10 on ∣z∣=2, g(0)=0, ∣g(z)∣=5∣z∣ everywhere and g′(0)=5.
ApplicationOriginal
Earnshaw's theorem. An electrostatic potential Φ 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Φ 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×n matrix has n 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 etA (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.