Notes on Complex Analysis
Export PDFGitHub

§5.4 The Schwarz–Christoffel Transformation

5.4 The Schwarz–Christoffel Transformation

The Riemann Mapping Theorem is elegant in its own simplicity and definitions. However, it is only a theorem that guarantees existence of biholomorphisms. No information whatsoever can be straightforwardly extracted regarding the explicit construction of such biholomorphisms. However, in the explicit case that 𝑈 is the open interior of a polygon, the result is provided by the Schwarz–Christoffel Transformation.

Let 𝑎1<𝑎2<⋯<𝑎𝑛 be 𝑛∈ℕ distinct real numbers and let 𝑎0=−∞ and 𝑎𝑛+1=∞. Suppose 𝛼1,𝛼2,…,𝛼𝑛 are 𝑛 positive real numbers satisfying ∑𝑘=1𝑛𝛼𝑘<𝑛−1. Let

𝛽(𝜁)=(𝜁−𝑎1)𝛼1−1⋯(𝜁−𝑎𝑛)𝛼𝑛−1=∏𝑘=1𝑛(𝜁−𝑎𝑘)𝛼𝑘−1,

where the branch of each factor is selected to be

(𝜁−𝑎𝑘)𝛼𝑘−1=e(𝛼𝑘−1)(log(𝜁−𝑎𝑘)),

where the branch of log is selected such that −π2<ℑ𝔪(log)≤3π2, holomorphic on ℂ∖iℝ≤0 (the lower imaginary axis is known as a branch cut). For 𝜁<𝑎𝑘, the argument of this factor is π(𝛼𝑘−1). For 𝜁<𝑎1,

arg(𝛽(𝜁))=π(−𝑛+∑𝑘=1𝑛𝛼𝑘),

achieved by selecting branches of each factor by the method described earlier.

Let 𝑘 be fixed. If 𝜁∈(𝑎𝑘−1,𝑎𝑘), the branches of all (𝜁−𝑎𝑗)𝛼𝑗−1 where 1≤𝑗≤𝑘−1 have vanishing arguments; hence,

arg(𝛽(𝜁))=π(−𝑛+𝑘−1+∑𝑗=𝑘𝑛𝛼𝑗).

If 𝜁>𝑎𝑛, we have

arg(𝛽(𝜁))=0.

Then, we define 𝑛+2 complex numbers via

𝑤𝑘=∫0𝑎𝑘𝛽(𝜁)d𝜁,(5.4.1)

where we integrate over the real line as an improper Riemann integral.

The absolute integrability of 𝛽(𝜁) along the real axis concerns only the convergence at each singularity 𝜁=𝑎𝑘 and the behavior as 𝜁→±∞. For each fixed 𝑘,

𝛽(𝜁)=ℎ𝑘(𝜁)(𝜁−𝑎𝑘)𝛼𝑘−1forℎ𝑘=∏𝑗=1𝑗≠𝑘𝑛(𝜁−𝑎𝑗)𝛼𝑗−1=𝒪︀(1)

as 𝜁→𝑎𝑘 (where ℎ𝑘 is holomorphic and nonzero in a compact neighborhood of 𝑎𝑘). Since 𝛼𝑘−1>−1,

|𝛽(𝜁)|=𝒪︀(|𝜁−𝑎𝑘|𝛼𝑘−1),(5.4.2)

and therefore 𝜁=𝑎𝑘 is an integrable singularity. Since 𝛽(𝜁)∼𝜁∑𝛼𝑘−𝑛 as 𝜁→±∞ and ∑𝑘=1𝑛𝛼𝑘−𝑛<−1, 𝛽 is integrable on ℝ.

Let

𝑓(𝑧)=∫0𝑧𝛽(𝜁)d𝜁.(5.4.3)

We first aim to show that this integral is path-independent for 𝑧∈ℍ+‾. This would be a direct consequence of Cauchy–Goursat (Theorem 3.1.9Theorem 3.1.9Theorem 3.1.9), if it weren’t for the singularities along the real line. We invoke a method used in Section 4.5Section 4.5Section 4.5, namely the indentation of contours; Let Γ⊂ℍ+‾ be an arbitrary piecewise 𝐶1 Jordan curve. There exists a 𝜀′>0 such that for all 0<𝜀<𝜀′, each 𝐷(𝑎𝑘,𝜀)‾ is disjoint. Then the indented contour

Γ̃=𝜕(intΓ∖⋃𝑘=1𝑛𝐷(𝑎𝑘,𝜀)‾)

is a Jordan curve in ℍ+‾, and can be decomposed as Γ̃=𝛾∪𝐶𝜀 where 𝛾⊆Γ and 𝐶𝜀⊆⋃𝑘=1𝑛𝜕𝐷(𝑎𝑘,𝜀). Then by Cauchy–Goursat (Theorem 3.1.9Theorem 3.1.9Theorem 3.1.9)

∮Γ̃𝛽(𝜁)d𝜁=∫𝛾𝛽(𝜁)d𝜁+∫𝐶𝜀𝛽(𝜁)d𝜁=0.

As 𝜀→0+, we have that

∮Γ𝛽(𝜁)d𝜁=lim𝜀→0+∫𝛾𝛽(𝜁)d𝜁=−lim𝜀→0+∫𝐶𝜀𝛽(𝜁)d𝜁,

where

|∫𝐶𝜀𝛽(𝜁)d𝜁|≤∑𝑘=1𝑛∮𝜕𝐷(𝑎𝑘,𝜀)|𝛽(𝜁)||d𝜁|≤∑𝑘=1𝑛∫02π|𝛽(𝑎𝑘+𝜀ei𝜃)|𝜀d𝜃=𝑛∫02π𝒪︀(𝜀𝛼𝑘)d𝜃=𝒪︀(𝜀𝛼𝑘)→0.(by(5.4.2)(5.4.2)(5.4.2))

Hence, ∮Γ𝛽(𝜁)d𝜁=0 for any piecewise 𝐶1 Jordan curve Γ⊂ℍ+‾, and the integral in (5.4.3)(5.4.3)(5.4.3) is path-independent for 𝑧∈ℍ+‾.

The holomorphy of 𝑓 in ℍ+ is equally evident: choosing the polygonal chain 0→i→𝑧 and applying 𝜕𝜕𝑧‾ gives

𝜕𝑓𝜕𝑧‾=𝜕𝜕𝑧‾∫i𝑧𝛽(𝜁)d𝜁=𝜕𝜕𝑧‾(𝑧−i)∫01𝛽(i+(𝑧−i)𝑡)d𝑡=(𝑧−i)∫01𝜕𝜕𝑧‾𝛽(i+(𝑧−i)𝑡)d𝑡=0,

since 𝛽(𝜁) is bounded on compact subsets of ℍ+ (this is the reason for which we choose 0→i→𝑧, as 𝛽 may be unbounded near 0). Thus 𝑓 is holomorphic on ℍ+; it is also trivial that 𝑓 is continuous on ℍ+‾.

For 𝑥∈ℝ where 𝑥∈(𝑎𝑘−1,𝑎𝑘) (where 1≤𝑘≤𝑛+1) and 𝑎0=−∞ and 𝑎𝑛+1=∞,

𝑓(𝑥)=∫0𝑥𝛽(𝑡)d𝑡=𝑤𝑘−1+∫𝑎𝑘−1𝑥𝛽(𝑡)d𝑡=𝑤𝑘−1+∫𝑎𝑘−1𝑥|𝛽(𝑡)|𝑒iarg(𝛽(𝑡))d𝑡=𝑤𝑘−1+exp[iπ(−𝑛+𝑘−1+∑𝑗=𝑘𝑛𝛼𝑗)]∫𝑎𝑘−1𝑥|𝛽(𝑡)|d𝑡,

where each 𝑤𝑘−1 is defined as in (5.4.1)(5.4.1)(5.4.1) (and observe that 𝑤𝑘−1=𝑓(𝑎𝑘−1)). Therefore,

arg[𝑓(𝑥)−𝑤𝑘−1]=π(−𝑛+𝑘−1+∑𝑗=𝑘𝑛𝛼𝑗)(5.4.4)

and

|𝑓(𝑥)−𝑤𝑘−1|=∫𝑎𝑘−1𝑥|𝛽(𝑡)|d𝑡(5.4.5)

for 𝑥∈(𝑎𝑘−1,𝑎𝑘). Define

𝑙𝑘=∫𝑎𝑘−1𝑎𝑘|𝛽(𝑡)|d𝑡andΔ𝑘=𝑓([𝑎𝑘−1,𝑎𝑘])(5.4.6)

for 1≤𝑘≤𝑛+1. Note that Δ𝑘 is a curve from 𝑤𝑘−1 to 𝑤𝑘. By (5.4.4)(5.4.4)(5.4.4), since arg[𝑓(𝑥)−𝑤𝑘−1] is constant for 𝑎𝑘−1<𝑥<𝑎𝑘, Δ𝑘 is a straight line segment. Furthermore, by (5.4.5)(5.4.5)(5.4.5) with 𝑥=𝑎𝑘, 𝑙𝑘 measures the length of said segment and its slope angle is given by

𝜃𝑘=arg[𝑤𝑘−𝑤𝑘−1]=π(−𝑛+𝑘−1+∑𝑗=𝑘𝑛𝛼𝑗).

Therefore, 𝑓 maps the real axis ℝ to a polygonal chain with vertices at 𝑤0,𝑤1,⋯,𝑤𝑛,𝑤𝑛+1. We now aim to prove that 𝑤0=𝑤𝑛+1, which is equivalent to showing

𝑤𝑛+1−𝑤0=∫−∞∞𝛽(𝑡)d𝑡=0.

Let 𝑅′>max𝑘∈ℕ≤𝑛|𝑎𝑘| be arbitrary; then ∃𝜀′>0 such that

⋃𝑘∈ℕ≤𝑛𝐷(𝑎𝑘,𝜀′)‾⊆𝐷(0,𝑅′)∩ℍ+

and each 𝐷(𝑎𝑘,𝜀′)‾ are disjoint. Moreover, ∀𝑅>𝑅′, ∀𝜀∈(0,𝜀′),

⋃𝑘∈ℕ≤𝑛𝐷(𝑎𝑘,𝜀)‾⊆𝐷(0,𝑅)∩ℍ+.

Therefore, we may define a closed (indented semicircle) contour (oriented counterclockwise)

Γ={𝑅ei𝜃:0≤𝜃≤π}∪([−𝑅,𝑅]∖⋃𝑘∈ℕ≤𝑛(𝑎𝑘−𝜀,𝑎𝑘+𝜀))∪⋃𝑘∈ℕ≤𝑛(𝜕𝐷(𝑎𝑘,𝜀)∩ℍ+‾)=𝐶𝑅∪𝐼𝑅∪𝐶𝜀,

such that 𝛽 is continuous on intΓ‾ and holomorphic on intΓ. By Cauchy–Goursat (Theorem 3.1.9Theorem 3.1.9Theorem 3.1.9),

∳Γ𝛽(𝜁)d𝜁=0,

and moreover, taking 𝑅→∞ and 𝜀→0+ gives

∫−∞∞𝛽(𝑡)d𝑡=−lim𝑅→∞∫𝐶𝑅𝛽(𝜁)d𝜁−lim𝜀→0+∫𝐶𝜀𝛽(𝜁)d𝜁.

By (5.4.2)(5.4.2)(5.4.2), we know that integrals over small indentations vanish; therefore

|𝑤𝑛+1−𝑤0|=lim𝑅→∞|∫𝐶𝑅𝛽(𝜁)d𝜁|≤lim𝑅→∞𝑅∫0π|∏𝑘=1𝑛(𝑅ei𝜃−𝑎𝑘)𝛼𝑘−1|d𝜃=lim𝑅→∞𝑅∫0π∏𝑘=1𝑛|𝑅ei𝜃−𝑎𝑘|𝛼𝑘−1d𝜃≤lim𝑅→∞𝑅∫0π∏𝑘=1𝑛max{(𝑅+|𝑎𝑘|)𝛼𝑘−1,|𝑅−|𝑎𝑘||𝛼𝑘−1}d𝜃≤lim𝑅→∞𝑅π𝒪︀(𝑅)−𝑛+∑𝑗=1𝑛𝛼𝑗=lim𝑅→∞𝒪︀(𝑅1−𝑛+∑𝑗=1𝑛𝛼𝑗)=0,

since ∑𝑘=1𝑛𝛼𝑘<𝑛−1 (by assumption). Here, we developed the upper bound max{(𝑅+|𝑎𝑘|)𝛼𝑘−1,|𝑅−|𝑎𝑘||𝛼𝑘−1}d𝜃, as the first part accounts solely for 𝛼𝑘−1≥0 and the second for 𝛼𝑘−1≤0 (although here it made little difference).

Hence, 𝑤0=𝑤𝑛+1 and 𝑓 maps ℝ̂=ℝ∪{∞} to a closed polygonal chain with 𝑛+1 sides (Δ𝑘). For 𝑘∈ℕ≤𝑛, at the vertex 𝑤𝑘 connecting Δ𝑘 and Δ𝑘+1, the tangent angle changes from 𝜃𝑘 to 𝜃𝑘+1, with an exterior angle of 𝜃𝑘+1−𝜃𝑘=π−π𝛼𝑘. Then at this vertex, the polygon has an interior angle of π𝛼𝑘.

At the vertex 𝑤0=𝑤𝑛+1, the polygon has an interior angle of π(𝑛−1)−∑𝑘=1𝑛π𝛼𝑘=π(𝑛−1−∑𝑘=1𝑛𝛼𝑘)>0, since the interior angles of the polygon add to π(𝑛−1).

If the values {𝛼𝑘} are chosen such that ∑𝑘=1𝑛π𝛼𝑘=𝑛−2, then the interior angle at vertex 𝑤0=𝑤𝑛+1 can be considered to be π; the vertex then is simply a straight line and we have a 𝑛-gon instead.

Now assume that the polygonal chain is simple (therefore, we concern ourselves only when there is a polygonal interior).

It remains to show that 𝑓 maps ℍ+ biholomorphically onto the polygon interior. First, the continuity of 𝑓 at infinity must be established (we have established that the limits exist from only the real line, which is insufficient to prove that a 𝛿-infimum exists in every direction). As |𝑧|→∞ (uniformly in all directions),

𝛽(𝑧)=𝒪︀(|𝑧|∑𝑘=1𝑛𝛼𝑘−𝑛).(5.4.7)

Let 𝑧=𝑟ei𝜃∈ℍ+‾, where 0≤𝜃≤π. Join 𝑧∈ℍ+‾ to 𝑟∈ℝ along the circular arc of radius 𝑟, and then join 𝑟 to infinity along the positive real axis. By (5.4.7)(5.4.7)(5.4.7),

|𝑤𝑛+1−𝑓(𝑧)|≤|𝑓(𝑟)−𝑓(𝑧)|+|𝑤𝑛+1−𝑓(𝑟)|≤∫0𝜃|𝛽(𝑟ei𝑡)|𝑟d𝑡+∫𝑟∞|𝛽(𝑡)|d𝑡≤𝒪︀(𝑟∑𝑘=1𝑛𝛼𝑘−𝑛+1)+∫𝑟∞𝒪︀(𝑡∑𝑘=1𝑛𝛼𝑘−𝑛)d𝑡=𝒪︀(𝑟∑𝑘=1𝑛𝛼𝑘−𝑛+1)→0as𝑟→∞.

Thus, defining 𝑓(∞)=𝑤𝑛+1=𝑤0 gives a continuous extension of 𝑓 to ℍ+‾∪{∞}.

Let 𝜓(𝑧)=𝑧−i𝑧+i be the Cayley transform from ℍ+ onto 𝔻, extended by 𝜓(∞)=1, and define

𝑔=𝑓∘𝜓−1.

The preceding continuity at infinity shows that 𝑔∈𝐶(𝔻‾), while 𝑔 is holomorphic on 𝔻. Furthermore, 𝑔|𝜕𝔻 maps 𝜕𝔻 injectively onto the polygonal boundary ⋃𝑘Δ𝑘. Therefore, by Theorem 5.1.5Theorem 5.1.5Theorem 5.1.5, 𝑔 maps 𝔻 univalently and surjectively onto int⋃𝑘Δ𝑘. Since 𝜓 is a biholomorphism, it follows that 𝑓 maps ℍ+ biholomorphically onto the polygon.

The definition in (5.4.3)(5.4.3)(5.4.3) is known as the Schwarz–Christoffel Transformation, which maps the upper-half plane to a polygon with specified angles. The side lengths can be controlled by numerically computing suitable values for {𝑎𝑘}. Note that for any 𝑐∈ℂ, 𝑐′∈ℂ∗, 𝑧↦𝑐+𝑐′𝑓(𝑧) also suffices as such a transformation.

Example 5.4.1.

Let

𝑄={𝑧:max{|ℜ𝔢𝑧|,|ℑ𝔪𝑧|}<1}
be a square centered at the origin with side length 2. Find a biholomorphism between 𝔻 and 𝑄.
Figure 5.1: The biholomorphism between 𝔻 and 𝑄 via 𝐹.
Solution to Example 5.4.1Example 5.4.1Example 5.4.1.

By the Cayley transform, it suffices to construct a biholomorphism from ℍ+ to 𝑄.

Since the interior angles of a square are π2 and (5.4.3)(5.4.3)(5.4.3) maps to an 𝑛+1-gon, we define the three exponents

𝛼𝑘=12for𝑘=1,2,3.

Then a biholomorphism is given by

𝑓𝑎1,𝑎2,𝑎3(𝑧)=∫0𝑧𝛽(𝜁)d𝜁=∫0𝑧d𝜁𝜁−𝑎1𝜁−𝑎2𝜁−𝑎3,

where 𝑎1<𝑎2<𝑎3,(±∞) are finite prevertices on the real line. Observe that by the restrictions set upon 𝛼𝑘, any such configuration of 𝑎1,𝑎2,𝑎3 will give a rectangle. Therefore, we want to find a configuration of prevertices that ensures equal side lengths. From the rectangularity of the boundary, it is already known that

𝑙1=𝑙3and𝑙2=𝑙4,(5.4.8)

where

𝑙1=∫−∞𝑎1|𝛽(𝑡)|d𝑡,𝑙2=∫𝑎1𝑎2|𝛽(𝑡)|d𝑡,𝑙3=∫𝑎2𝑎3|𝛽(𝑡)|d𝑡,𝑙4=∫𝑎3𝑎4|𝛽(𝑡)|d𝑡

correspond to the side lengths of the rectangle.

Fix 𝑎1,𝑎3. We now aim to find 𝑎2 in terms of 𝑎1 and 𝑎3 such that 𝑙1=𝑙2 (which would then imply 𝑙1=𝑙2=𝑙3=𝑙4 from (5.4.8)(5.4.8)(5.4.8)). Define a linear affine map that sends −1→𝑎1, 1→𝑎3, namely 𝜑(𝑧)=𝑧(𝑎3−𝑎1)2+𝑎3+𝑎12. By a substitution,

𝑙1=∫−∞𝑎1|𝛽(𝑡)|d𝑡=𝑎3−𝑎12∫−∞−1d𝑡|𝜑(𝑡)−𝑎1||𝜑(𝑡)−𝑎2||𝜑(𝑡)−𝑎3|=∫−∞−1d𝑡(𝑡2−1)|𝑡(𝑎3−𝑎1)2+𝑎3+𝑎1−2𝑎22|,=∫1∞d𝑡(𝑡2−1)|𝑡(𝑎1−𝑎3)2+𝑎3+𝑎1−2𝑎22|,(1)𝑙4=∫𝑎3∞|𝛽(𝑡)|d𝑡=𝑎3−𝑎12∫1∞d𝑡|𝜑(𝑡)−𝑎1||𝜑(𝑡)−𝑎2||𝜑(𝑡)−𝑎3|=∫1∞d𝑡(𝑡2−1)|𝑡(𝑎3−𝑎1)2+𝑎3+𝑎1−2𝑎22|.(2)

Now force 𝑙1=𝑙4.

Our aim is to show that 𝑎2 lies at the midpoint between 𝑎1 and 𝑎3. For the sake of contradiction, assume 𝑎3+𝑎1−2𝑎2>0. Then for all 𝑡>1,

𝑡(𝑎1−𝑎3)2+𝑎3+𝑎1−2𝑎22<𝑎1−𝑎2<0,

implying that

|𝑡(𝑎1−𝑎3)2+𝑎3+𝑎1−2𝑎22|=𝑡(𝑎3−𝑎1)2−𝑎3+𝑎1−2𝑎22<𝑡(𝑎3−𝑎1)2,

while

|𝑡(𝑎3−𝑎1)2+𝑎3+𝑎1−2𝑎22|=𝑡(𝑎3−𝑎1)2+𝑎3+𝑎1−2𝑎22>𝑡(𝑎3−𝑎1)2.

Then we obtain that (1) is strictly greater than (2). Similarly, 𝑎3+𝑎1−2𝑎2<0 forces 𝑙1<𝑙4. Therefore, for this equality to hold, we must have 𝑎2=𝑎1+𝑎32.

Then the general biholomorphism between ℍ+ and a square is in the form of

𝑐′𝑓𝑎1,𝑎1+𝑎32,𝑎3(𝑧)+𝑐.(5.4.9)

Without loss of generality, assume 𝑎1=−1,𝑎3=1,𝑎2=0 and write 𝑓(𝑧)=𝑓−1,0,1(𝑧). Let

𝑙=𝑙1=⋯=𝑙4=∫−10|𝛽(𝑡)|d𝑡=∫01d𝑡𝑡(1−𝑡2)=∫01d𝑡2𝑡𝑡(1−𝑡)=12∫01𝑡−34(1−𝑡)−12d𝑡=12Β(14,12)=Γ(14)Γ(12)2Γ(34)=Γ(14)222π

by Definition 12.1.9Definition 12.1.9Definition 12.1.9, Theorem 12.1.5Theorem 12.1.5Theorem 12.1.5 (the Gamma and Beta function identities). Traversing ℝ from −∞ to +∞, by using (5.4.4)(5.4.4)(5.4.4),

∫−∞−1𝛽=i𝑙,∫−10𝛽=−𝑙,∫01𝛽=−i𝑙,∫1∞𝛽=𝑙.

Hence,

𝑓(−1)=𝑙,𝑤(0)=0,𝑤(1)=−i𝑙,𝑓(∞)=𝑙−i𝑙=𝑓(−∞).

The four distinct vertices 0,𝑙,𝑙−i𝑙,−i𝑙 bound

𝑆={𝑥+i𝑦:0<𝑥<𝑙,−𝑙<𝑦<0},

a square of side 𝑙. Thus 𝑓(ℍ+)=𝑆. The center of 𝑆 is 𝑙2−i𝑙2. Thus translate and dilate to 𝑄=(−1,1)×(−1,1) as in (5.4.9)(5.4.9)(5.4.9), to obtain

Φ(𝑧)=2𝑙𝑓(𝑧)−1+i=2𝑙∫0𝑧d𝜁𝜁+1𝜁𝜁−1−1+i.

Then Φ:ℍ+→𝑄 is biholomorphic with an extension to Φ(0)=−1+i, Φ(−1)=1+i, Φ(±∞)=1−i, Φ(1)=−1−i.

Finally 𝜓(𝑧)=𝑧−i𝑧+i maps ℍ+ onto 𝔻 with inverse 𝜓−1(𝑤)=i(1+𝑤)1−𝑤. Hence

𝐹(𝑤)=Φ∘𝜓−1(𝑤)=42𝜋Γ(14)2∫0i(1+𝑤)1−𝑤d𝜁𝜁+1𝜁𝜁−1−1+i(5.4.10)

gives the required biholomorphism, where the square root branches are selected such that their corresponding logarithms have imaginary parts in (−π2,3π2]. A visualization of this conformal map has been provided in Figure 5.1Figure 5.1Figure 5.1. ⁠□