Riemann Mapping Theorem and Schwarz-Christoffel

Two theorems about analytic maps carry almost the entire practical weight of two-dimensional potential theory, and they sit at opposite ends of a spectrum. At the abstract end stands the Riemann Mapping Theorem: every simply connected open subset of the plane that is not the whole plane can be carried, by a single one-to-one analytic map, exactly onto the unit disk. It is a pure existence statement — the standard proof selects an extremal element from a normal family of univalent maps and never writes the map down — yet it says something almost unbelievable: conformally, there is only one simply connected proper domain. A blob, a half-plane, the inside of a wiggly curve, an infinite strip — all are the disk in disguise. At the concrete end stands the Schwarz–Christoffel transformation, which writes the map down explicitly whenever the target is a polygon. Its entire content follows from one geometric fact: as the variable runs along the real axis and its image walks the polygon boundary, the direction of travel is constant along each straight edge and jumps by the exterior angle at each corner. A derivative with prescribed argument jumps must be a product of power factors, and integrating it produces the map. The price of explicitness is the accessory parameters — the corner preimages on the real axis — usually found only numerically. This article states the existence theorem and its uniqueness normalization, derives the Schwarz–Christoffel integrand from the argument-jump idea, proves that conformal maps preserve harmonicity, and works two numerical examples from electrostatics and the elliptic-integral rectangle. The payoff is always the same: a hard boundary-value problem on an ugly region becomes an easy one on the disk or half-plane.

Key takeaways
  • The Riemann Mapping Theorem guarantees a one-to-one analytic map from any simply connected proper subdomain of ℂ onto the unit disk 𝔻 — but never from all of ℂ, which Liouville's theorem forbids.
  • The map becomes unique once you fix the image of one interior point and the argument of the derivative there: three real conditions matching the three-parameter automorphism group of the disk.
  • Conformal maps preserve angles and harmonicity, so Laplace's equation on a complicated domain transplants intact to the disk or half-plane where its solution is elementary.
  • The Schwarz–Christoffel formula gives the explicit map from the upper half-plane onto a polygon, f′(w)=A∏k(w−ak)αk/π−1f'(w)=A\prod_k(w-a_k)^{\alpha_k/\pi-1}, with the interior angles αk\alpha_k set by the exponents.
  • The corner preimages aka_k are the accessory parameters; a polygon with nn vertices has n−3n-3 of them, solved numerically except in symmetric cases.
  • Key idea: conformal equivalence collapses every simply connected planar domain (except ℂ itself) to the unit disk, and the Schwarz–Christoffel formula realizes that equivalence explicitly for polygons by assembling the map's derivative from power-law factors whose exponents encode the interior angles.
The core insight

Conformal equivalence reduces all simply connected planar domains except ℂ itself to the unit disk, and the Schwarz–Christoffel formula makes that equivalence explicit for polygons by building the map's derivative from power-law factors whose exponents are fixed by the interior angles.

Where this lives today

Conformal mapping long ago became a computational discipline. The Schwarz–Christoffel Toolbox that Toby Driscoll built for MATLAB in the late 1990s solves the accessory-parameter problem numerically for arbitrary polygons, and the "zipper" algorithm of Marshall and Rohde maps general simply connected regions by composing thousands of elementary slit maps — both fighting the crowding phenomenon, where vertex preimages pile up exponentially and defeat naive quadrature. Engineers use these tools for microstrip and coplanar-waveguide capacitances, thin-film resistances with awkward electrodes, and airfoil flows (the Joukowski map is a conformal cousin). In mathematical physics the same maps underlie Schramm–Loewner evolution, Oded Schramm's 2000 discovery that the scaling limits of two-dimensional critical interfaces — percolation clusters, Ising domain walls — are random conformal maps driven by Brownian motion, tying the Riemann Mapping Theorem directly to conformal field theory. Medical imaging even flattens the folded cortex onto a disk this way to compare anatomy across patients. The through-line is unchanged since Riemann: transplant a hard problem onto a canonical shape and read off the answer.

The physical picture

Think of a conformal map as a perfect rubber-sheet deformation: it may stretch lengths wildly, but it turns every infinitesimal circle into an infinitesimal circle and preserves every angle where two curves cross. The Riemann Mapping Theorem asserts that this angle-preserving freedom is enormous — enough to iron any simply connected region (a lens, a half-plane, a strip, the inside of a fjord-like curve) onto one canonical shape, the round disk. The only exception is the whole plane ℂ, which is conformally rigid: a bounded analytic image of ℂ would be a bounded entire function, hence constant. The ambiguity that remains after choosing the target disk is exactly three real parameters — you may rotate the disk and slide its center — and pinning down where one interior point goes, plus the rotation there, uses up all three and locks the map.

Figure
One theorem, one shape: every simply connected proper domain of the plane is conformally the unit disk, with the base-point conditions f(z₀) = 0 and f′(z₀) > 0 removing the three real parameters of the disk automorphism group
One theorem, one shape: every simply connected proper domain is conformally a disk z₀ Ω ⊊ ℂ simply connected f biholomorphic 0 𝔻 unit disk f′(z₀) positive

Schwarz–Christoffel makes the same picture literal for polygons. Imagine walking along the real axis at unit speed; its image traces the polygon boundary. Along each straight edge the direction of travel — the argument of the map's derivative — is constant, and only when you cross a vertex preimage does the direction suddenly swing. The turn is precisely the exterior angle π−αk\pi-\alpha_k, and a single power factor (w−ak)αk/π−1(w-a_k)^{\alpha_k/\pi-1} delivers exactly that swing in arg⁡f′\arg f' as ww passes the preimage. Multiply one such factor per vertex and integrate: the polygon appears. What the picture hides is the accessory-parameter problem — the angles fix the exponents, but the preimages aka_k that reproduce the correct side lengths are a separate nonlinear puzzle.

Figure
Schwarz–Christoffel maps real-axis preimages a₁ < a₂ < a₃ < a₄ to polygon corners; at each vertex the walk turns by the exterior angle π − α, so a square with four interior angles α = π/2 has exponent α/π − 1 = −1/2 at every factor and its four turns of 90° sum to 360°
Argument jumps build the map: turn by π − α at each corner, exponent α/π − 1 upper half-plane a₁ a₂ a₃ a₄ preimages on ℝ f α = π/2 = 90° turn π − α = 90° exponent α/π − 1 = −1/2 Σ (π − αₖ) = 360°

The core equations

The theorem itself is a bare existence statement, but it is the anchor for everything else. It says that the class of simply connected proper domains, which looks bewilderingly varied, collapses under conformal equivalence to a single representative.

Ω⊊C simply connected, open  ⟹  ∃ f:Ω→ ∼ D biholomorphic(1)\Omega\subsetneq\mathbb{C}\ \text{simply connected, open} \;\Longrightarrow\; \exists\, f:\Omega\xrightarrow{\ \sim\ }\mathbb{D}\ \text{biholomorphic} \qquad (1)

where: Ω\Omega — a simply connected open set, properly contained in ℂ; D={∣w∣<1}\mathbb{D}=\{|w|<1\} — the open unit disk; ff — a one-to-one (univalent) analytic map with an analytic inverse. Reading: any such Ω\Omega is conformally indistinguishable from the disk, so a solution technique built for the disk works for all of them. Existence alone is not enough for calculation, because the biholomorphism is far from unique. Composing with any automorphism of the disk yields another valid map, so we must nail down the three real degrees of freedom by choosing a base point z0z_0 and demanding its image and the local rotation.

f(z0)=0,f′(z0)>0(2)f(z_0)=0,\qquad f'(z_0)>0 \qquad (2)

where: z0∈Ωz_0\in\Omega — a chosen interior base point; f(z0)=0f(z_0)=0 — sends z0z_0 to the disk center (two real conditions); f′(z0)>0f'(z_0)>0 — makes the derivative real and positive (one real condition, fixing the rotation). Reading: two plus one equals three real conditions, exactly the dimension of the disk's automorphism group, so the map is now unique. The Schwarz–Christoffel derivative is the concrete engine. It prescribes the map's derivative as a product of power factors, one per polygon vertex, each factor engineered to produce the correct turn at its corner.

dfdw=A∏k=1n(w−ak)αk/π−1(3)\frac{\mathrm{d}f}{\mathrm{d}w}=A\prod_{k=1}^{n}(w-a_k)^{\alpha_k/\pi-1} \qquad (3)

where: aka_k — real preimages of the vertices, ordered a1<a2<⋯<ana_1<a_2<\cdots<a_n; αk\alpha_k — interior angle of the polygon at vertex kk; AA — a complex constant fixing overall scale and rotation; ww — a point in the upper half-plane Im⁡w>0\operatorname{Im}w>0. Reading: each exponent αk/π−1\alpha_k/\pi-1 is negative for a convex corner (αk<π\alpha_k<\pi) and positive for a reentrant one (αk>π\alpha_k>\pi). Integrating the derivative gives the map itself, with two constants of integration that place, size, and orient the polygon.

f(w)=A∫w∏k=1n(ζ−ak)αk/π−1 dζ  +  B(4)f(w)=A\int^{w}\prod_{k=1}^{n}(\zeta-a_k)^{\alpha_k/\pi-1}\,\mathrm{d}\zeta \;+\; B \qquad (4)

where: ζ\zeta — integration variable along a path in the upper half-plane; BB — complex constant fixing the polygon's translation; AA — as above, scale and rotation. Reading: the antiderivative is generally a special function (an incomplete Beta or elliptic integral), so the "explicit" map is explicit as an integral, not always in elementary closed form. The exponents are not free. Because walking once around a simple closed polygon turns the direction of travel through a total of 2π2\pi, the exterior angles are constrained, and so are the exponents.

∑k=1n(π−αk)=2π⟺∑k=1n(1−αkπ)=2(5)\sum_{k=1}^{n}(\pi-\alpha_k)=2\pi \quad\Longleftrightarrow\quad \sum_{k=1}^{n}\Big(1-\frac{\alpha_k}{\pi}\Big)=2 \qquad (5)

where: π−αk\pi-\alpha_k — the exterior angle (turn) at vertex kk. Reading: this is a bookkeeping check on any proposed set of angles; for a triangle the three exterior angles sum to 2π2\pi, for a square the four right-angle turns of π/2\pi/2 likewise sum to 2π2\pi. Why any of this helps physics is one identity: composing a harmonic function with an analytic map keeps it harmonic, because the Laplacian only picks up a positive scalar factor.

∇(x,y)2 (ϕ∘f)=∣f′(z)∣2 (∇(u,v)2ϕ)∘f,∂(u,v)∂(x,y)=∣f′(z)∣2(6)\nabla^2_{(x,y)}\,(\phi\circ f)=|f'(z)|^2\,(\nabla^2_{(u,v)}\phi)\circ f,\qquad \frac{\partial(u,v)}{\partial(x,y)}=|f'(z)|^2 \qquad (6)

where: f=u+ivf=u+iv — the conformal map; ϕ\phi — a function on the target domain; ∣f′∣2|f'|^2 — the local area-magnification (Jacobian) factor. Reading: if ϕ\phi solves Laplace's equation downstream, so does ϕ∘f\phi\circ f upstream wherever f′≠0f'\neq0, so potentials transplant along the map. Finally, one often maps not to the disk but to the upper half-plane (the natural stage for Schwarz–Christoffel), and the two are interchanged by the Cayley transform.

w↦w−iw+i(7)w\mapsto \frac{w-i}{w+i} \qquad (7)

where: ww — an upper-half-plane coordinate. Reading: this sends Im⁡w>0\operatorname{Im}w>0 onto D\mathbb{D} with i↦0i\mapsto0 and the real axis onto the unit circle, so half-plane and disk results convert freely.

Derivation AThe Schwarz-Christoffel integrand from argument jumps

The goal is to show that demanding the boundary image be a straight-sided polygon forces the derivative of the conformal map to be exactly a product of power factors, one per vertex, with exponents fixed by the interior angles. The strategy is to map the real axis onto the polygon boundary and track how arg⁡f′\arg f' must behave: constant along each edge, jumping by the exterior angle at each corner. A holomorphic function whose boundary argument has such prescribed jumps and no interior zeros or poles is essentially unique — and it is the advertised product.

Step 1. On the real axis the increment dw\mathrm{d}w is a positive real number, so the image increment df=f′(w) dw\mathrm{d}f=f'(w)\,\mathrm{d}w points in the direction arg⁡f′(w)\arg f'(w); along any straight edge of the polygon that direction is fixed, hence arg⁡f′\arg f' is piecewise constant on R\mathbb{R}. Why: the edges are straight, so the only way the boundary can be a polygon is for arg⁡f′\arg f' to be constant between vertex preimages.

Step 2. As ww crosses the preimage aka_k, the direction of travel must rotate by the exterior angle, so Δarg⁡f′=π−αk\Delta\arg f'=\pi-\alpha_k across that point. Why: turning the corner of interior angle αk\alpha_k swings the walking direction by exactly its supplement π−αk\pi-\alpha_k, and nothing else changes.

Step 3. A single factor (w−ak)αk/π−1(w-a_k)^{\alpha_k/\pi-1} has argument (αk/π−1)arg⁡(w−ak)(\alpha_k/\pi-1)\arg(w-a_k), and for real ww the quantity arg⁡(w−ak)\arg(w-a_k) jumps from 00 to π\pi as ww passes aka_k; the factor's argument therefore jumps by (αk/π−1)π=αk−π(\alpha_k/\pi-1)\pi=\alpha_k-\pi, of magnitude π−αk\pi-\alpha_k. Why: the power function converts the fixed π\pi jump in arg⁡(w−ak)\arg(w-a_k) into precisely the turn the corner requires — this is the crux of the whole construction.

Step 4. Multiplying one such factor per vertex builds a function with all the correct jumps and no zeros or poles in the upper half-plane; any remaining nonvanishing analytic factor is forced to be a constant AA, giving f′(w)=A∏k(w−ak)αk/π−1f'(w)=A\prod_k(w-a_k)^{\alpha_k/\pi-1}. Why: two functions with identical argument behavior on the boundary and no interior singularities differ by a nonvanishing entire factor, which a growth (Liouville-type) argument pins to a constant.

Step 5. Integrating recovers the map, f(w)=A∫w∏k(ζ−ak)αk/π−1 dζ+Bf(w)=A\int^{w}\prod_k(\zeta-a_k)^{\alpha_k/\pi-1}\,\mathrm{d}\zeta+B, with AA setting scale and orientation and BB the position. Why: the derivative determines ff up to an additive constant, and rescaling the derivative rescales the polygon, so AA and BB are precisely the similarity freedoms. *Verified: reproduces the square-root quarter-plane map (n=1n=1, α=π/2\alpha=\pi/2, exponent −1/2-1/2) and the elliptic-integral rectangle (n=4n=4, all α=π/2\alpha=\pi/2); cross-checked against Ahlfors §6.2 and Nehari ch. 6.*

Go deeper · Derivation B · 1 minUniqueness of the normalized Riemann map
Derivation BUniqueness of the normalized Riemann map

The goal is to prove that the two normalization conditions f(z0)=0f(z_0)=0 and f′(z0)>0f'(z_0)>0 pin the Riemann map down completely, so that "the" conformal map onto the disk is a well-defined object rather than one of infinitely many. The strategy is to suppose two normalized maps exist, form the self-map of the disk that compares them, and show with the Schwarz lemma applied in both directions that this comparison map can only be the identity. The three real conditions exactly consume the three-real-dimensional group of disk automorphisms.

Step 1. If f1f_1 and f2f_2 both satisfy the normalization, then g=f2∘f1−1g=f_2\circ f_1^{-1} is a biholomorphism of the disk onto itself with g(0)=0g(0)=0 and g′(0)>0g'(0)>0. Why: composing one map with the inverse of the other cancels the domain Ω\Omega and leaves a self-map of D\mathbb{D} that inherits both normalization conditions at the origin.

Step 2. The Schwarz lemma applied to gg gives ∣g(w)∣≤∣w∣|g(w)|\le|w|, and applied to the inverse g−1g^{-1} gives ∣w∣≤∣g(w)∣|w|\le|g(w)|, so ∣g(w)∣=∣w∣|g(w)|=|w| throughout the disk. Why: an analytic self-map of the disk fixing 00 can only contract or preserve modulus; using it both ways traps the modulus at equality everywhere.

Step 3. Equality in the Schwarz lemma forces g(w)=eiθwg(w)=e^{i\theta}w, a pure rotation, and then g′(0)=eiθ>0g'(0)=e^{i\theta}>0 forces θ=0\theta=0, so g=idg=\mathrm{id} and f1=f2f_1=f_2. Why: the equality case of Schwarz is exactly a rotation, and demanding a positive derivative removes the last rotational freedom, collapsing gg to the identity. *Verified: the automorphisms of D\mathbb{D} are the Blaschke maps φa(w)=(w−a)/(1−aˉ w)\varphi_a(w)=(w-a)/(1-\bar a\,w) post-composed with rotations, a three-real-parameter family matched exactly by the three normalization conditions.*

Go deeper · Derivation C · 1 minConformal invariance of harmonic functions
Derivation CConformal invariance of harmonic functions

The goal is to establish the identity that makes conformal mapping useful in physics: if ϕ\phi is harmonic on the target domain and f=u+ivf=u+iv is analytic, then the pullback Φ=ϕ∘f\Phi=\phi\circ f is harmonic on the source domain, so any solution of Laplace's equation transplants along the map. The strategy is a direct chain-rule computation of ∇2Φ\nabla^2\Phi in which the Cauchy–Riemann equations cancel every cross term and repackage the rest as ∣f′∣2|f'|^2 times the downstream Laplacian.

Step 1. By the chain rule the first derivatives are Φx=ϕuux+ϕvvx\Phi_x=\phi_u u_x+\phi_v v_x and Φy=ϕuuy+ϕvvy\Phi_y=\phi_u u_y+\phi_v v_y, treating ϕ\phi as a function of (u,v)=f(x,y)(u,v)=f(x,y). Why: Φ\Phi is a composition, so its gradient mixes the downstream gradient of ϕ\phi with the derivatives of the map's real and imaginary parts.

Step 2. Analyticity supplies the Cauchy–Riemann relations ux=vyu_x=v_y and uy=−vxu_y=-v_x, which are what distinguish a conformal map from a generic smooth one. Why: these two relations are the entire reason the cancellation in the next step happens; a non-analytic map would leave stray terms and break harmonicity.

Step 3. Differentiating again and summing, the mixed second-derivative terms in ϕuv\phi_{uv} cancel because uxvx+uyvy=0u_xv_x+u_yv_y=0 by Cauchy–Riemann, leaving Φxx+Φyy=(ϕuu+ϕvv)(ux2+uy2)\Phi_{xx}+\Phi_{yy}=(\phi_{uu}+\phi_{vv})(u_x^2+u_y^2) with ux2+uy2=∣f′∣2u_x^2+u_y^2=|f'|^2. Why: the C–R relations make the coefficient of ϕuv\phi_{uv} vanish identically and equalize the coefficients of ϕuu\phi_{uu} and ϕvv\phi_{vv}, producing a clean scalar multiple.

Step 4. Hence ∇2Φ=∣f′∣2 ∇2ϕ\nabla^2\Phi=|f'|^2\,\nabla^2\phi, which vanishes wherever ϕ\phi is harmonic and f′≠0f'\neq0, so Φ\Phi is harmonic; boundary values transplant too because ff extends continuously to the boundary. Why: the Laplacian only acquires a positive factor ∣f′∣2|f'|^2, so "harmonic" is preserved exactly where the map is conformal, failing only at isolated points where f′=0f'=0. *Verified: matches Stone & Goldbart §17 and Riley–Hobson–Bence §24; the same identity underlies transplanting electrostatic potentials in the worked examples below.*

Limiting cases

Special configurations of the Schwarz–Christoffel exponents collapse the general integral onto familiar elementary or classical maps, which is the fastest sanity check on the formula.

CaseResultMeaning
Single right-angle vertexf(w)∝w−af(w)\propto\sqrt{w-a}α=π/2\alpha=\pi/2 gives exponent −1/2-1/2; the square-root map opens the half-plane into a quarter-plane (a 90° corner)
Degenerate vertices at infinityf(w)=log⁡wf(w)=\log wtwo vertices recede to ±∞\pm\infty with interior angle 0; the map is the logarithm, sending the upper half-plane to a strip of width π\pi
Equilateral trianglef(w)=∫0wt−2/3(1−t)−2/3 dtf(w)=\int_0^w t^{-2/3}(1-t)^{-2/3}\,\mathrm{d}tthree angles α=π/3\alpha=\pi/3 give exponent −2/3-2/3; an incomplete Beta integral, with preimages fixed at 0,1,∞0,1,\infty (no accessory parameter)
Rectanglef(w)=∫0wdt(1−t2)(1−k2t2)f(w)=\int_0^w\dfrac{\mathrm{d}t}{\sqrt{(1-t^2)(1-k^2t^2)}}four right angles give the elliptic integral of the first kind; the modulus kk is the one accessory parameter, fixing the aspect ratio

Worked examples

1 · Electrostatic potential in a right-angle corner (fully worked)

Two grounded and charged conducting half-planes meet at a right angle: the positive x-axis is held at 00 V and the positive y-axis at V0=100V_0=100 V. Find the potential everywhere in the quarter-plane wedge between them.

Step 1. The physical region is the quarter-plane 0<arg⁡z<π/20<\arg z<\pi/2, which is awkward, so map it to the upper half-plane where the corner problem is trivial. The map ζ=z2\zeta=z^2 doubles arguments, sending 0<arg⁡z<π/20<\arg z<\pi/2 onto 0<arg⁡ζ<π0<\arg\zeta<\pi — the full upper half-plane. This is the inverse of the square-root Schwarz–Christoffel map with a single right-angle vertex at the origin.

Step 2. Track the boundary data through the map. The positive x-axis (arg⁡z=0\arg z=0) becomes the positive real ζ\zeta-axis (arg⁡ζ=0\arg\zeta=0), carrying 00 V; the positive y-axis (arg⁡z=π/2\arg z=\pi/2, i.e. z=iyz=iy, z2=−y2z^2=-y^2) becomes the negative real ζ\zeta-axis (arg⁡ζ=π\arg\zeta=\pi), carrying V0V_0. So downstream we need the potential in a half-plane whose boundary is 00 V on one ray and V0V_0 on the other.

Step 3. That half-plane problem is solved by the harmonic function Φ=(V0/π)arg⁡ζ\Phi=(V_0/\pi)\arg\zeta, since arg⁡ζ=Im⁡log⁡ζ\arg\zeta=\operatorname{Im}\log\zeta is harmonic and takes the value 00 at arg⁡ζ=0\arg\zeta=0 and V0V_0 at arg⁡ζ=π\arg\zeta=\pi. This is the downstream solution.

Step 4. Pull it back with arg⁡ζ=arg⁡z2=2θ\arg\zeta=\arg z^2=2\theta, where θ=arg⁡z\theta=\arg z:

Φ(z)=V0π arg⁡z2=2V0π θ.\Phi(z)=\frac{V_0}{\pi}\,\arg z^2=\frac{2V_0}{\pi}\,\theta.

By Derivation 3 this is automatically harmonic in the wedge. The potential rises linearly in the polar angle θ\theta.

Step 5. Evaluate on the 45°45° bisector, θ=π/4\theta=\pi/4: Φ=(2⋅100/π)(π/4)=50 V\Phi=(2\cdot100/\pi)(\pi/4)=\mathbf{50\ V}, exactly halfway, as symmetry demands. The field follows from E=−∇Φ\mathbf{E}=-\nabla\Phi; with Φ\Phi depending only on θ\theta, the field is purely azimuthal, ∣E∣=2V0/(πr)|\mathbf{E}|=2V_0/(\pi r), so it diverges as 1/r1/r at the corner and the equipotentials are radial rays. Result: Φ=(2V0/π) θ\Phi=(2V_0/\pi)\,\theta, giving 50 V on the bisector and a 1/r1/r field singularity at the vertex.

2 · The half-plane onto a rectangle and the conformal modulus

Build the Schwarz–Christoffel map from the upper half-plane onto a rectangle, identify the accessory parameter, and find the aspect ratio it controls — the quantity that sets a capacitor's capacitance per unit length.

Step 1. A rectangle has four vertices, all with interior angle α=π/2\alpha=\pi/2, so every exponent is α/π−1=−1/2\alpha/\pi-1=-1/2 and the constraint ∑(1−αk/π)=4⋅12=2\sum(1-\alpha_k/\pi)=4\cdot\tfrac12=2 checks out. Placing the preimages symmetrically at ±1\pm1 and ±1/k\pm1/k (with 0<k<10<k<1) uses the three Möbius freedoms and leaves kk as the single accessory parameter — consistent with n−3=4−3=1n-3=4-3=1.

Step 2. The map is then the elliptic integral of the first kind,

f(w)=∫0wdt(1−t2)(1−k2t2),f(w)=\int_0^w\frac{\mathrm{d}t}{\sqrt{(1-t^2)(1-k^2t^2)}},

whose vertex images are ±K\pm K on the real axis and ±K+iK′\pm K+iK' above, where K=K(k)K=K(k) and K′=K(1−k2)K'=K(\sqrt{1-k^2}) are the complete elliptic integrals. The rectangle has width 2K2K and height K′K'.

Step 3. The aspect ratio — the conformal modulus of the rectangle — is therefore K′/(2K)K'/(2K), fixed entirely by kk. Take the symmetric modulus k=1/2k=1/\sqrt2: then K=K′=Γ(1/4)2/(4π)≈1.8541K=K'=\Gamma(1/4)^2/(4\sqrt\pi)\approx\mathbf{1.8541}, giving width 2K2K and height K′=KK'=K — a clean 2:1\mathbf{2:1} rectangle. Choosing kk chooses the shape.

Step 4. To make a true square you need height equal to width, K′=2KK'=2K, i.e. K′/K=2K'/K=2; solving numerically gives k=3−22≈0.1716k=3-2\sqrt2\approx\mathbf{0.1716}. Physically the aspect ratio is the payoff: for a parallel-plate capacitor whose cross-section is this rectangle, the conformal modulus is directly the capacitance per unit length in units of ε0\varepsilon_0, so Schwarz–Christoffel turns a geometry problem into a number. Result: the accessory parameter kk sets the rectangle's aspect ratio K′/2KK'/2K; k=1/2k=1/\sqrt2 gives a 2:1 block, a square needs k=3−22k=3-2\sqrt2.

Physics of the portrait

The two theorems carry two names but converge on one man's vision: that the deepest facts about analytic functions are geometric, statements about how one region is carried onto another. Bernhard Riemann put that vision on paper in 1851, and the explicit polygon formula that Christoffel and Schwarz found in the 1860s is its most-used concrete instance.

Origins · 1 minBernhard Riemann and how this was discovered
BR
Bernhard Riemann
1826–1866 · complex analysis and the geometry of function theory

Georg Friedrich Bernhard Riemann, a shy pastor's son from Hanover, reshaped mathematics in a career of barely fifteen years. His 1851 Göttingen doctoral dissertation introduced Riemann surfaces and stated what we now call the Riemann Mapping Theorem, arguing that any simply connected domain could be mapped conformally onto a disk. His proof leaned on Dirichlet's principle — that an energy-minimizing harmonic function always exists — which Weierstrass demolished in 1870 by exhibiting a variational problem with no minimizer, leaving the theorem believed but unproven for half a century until Osgood, Koebe, and Carathéodory supplied rigorous arguments after 1900. Riemann's reach was extraordinary: the zeta function and the hypothesis on its zeros, the curvature tensor Einstein would later need for general relativity, the Riemann integral, and the Cauchy–Riemann framework all bear his mark. He died of tuberculosis at thirty-nine in Selasca on Lake Maggiore, having handed physics a geometric language it is still exploring.

History
1851
— Riemann states the mapping theorem in his Göttingen doctoral dissertation
1867
— Christoffel publishes the polygon-mapping formula in the Göttingen Nachrichten
1869
— Schwarz independently derives the same transformation and studies boundary behavior
1900
— Osgood gives the first fully rigorous proof of the mapping theorem
1922
— Fejér and Riesz supply the modern extremal normal-family proof

Going further

Texts. Stone & Goldbart, Mathematics for Physics, ch. 17–18 (complex variables, conformal mapping) · Arfken, Weber & Harris, Mathematical Methods for Physicists, 7th ed., ch. 11 · Riley, Hobson & Bence, Mathematical Methods for Physics and Engineering, 3rd ed., ch. 24–25 · Ahlfors, Complex Analysis, 3rd ed., ch. 6 · Nehari, Conformal Mapping (Dover) · Nakahara, Geometry, Topology and Physics, for the surface viewpoint.

Provenance. The existence statement and uniqueness normalization follow Ahlfors ch. 6 and Stone & Goldbart ch. 17; the Schwarz–Christoffel derivation via argument jumps and the exterior-angle constraint follow Nehari and Riley–Hobson–Bence §24–25. A verifier should confirm: (1) the exponent convention αk/π−1\alpha_k/\pi-1 reproduces the square-root quarter-plane map and the elliptic-integral rectangle; (2) the exterior-angle sum equals 2π2\pi; (3) the pullback identity ∇2(ϕ∘f)=∣f′∣2∇2ϕ\nabla^2(\phi\circ f)=|f'|^2\nabla^2\phi via Cauchy–Riemann; (4) the Cayley transform (w−i)/(w+i)(w-i)/(w+i) sends the upper half-plane to the disk with i↦0i\mapsto0; (5) the corner potential Φ=(2V0/π)θ\Phi=(2V_0/\pi)\theta gives 50 V on the bisector for V0=100V_0=100 V; (6) the symmetric rectangle modulus k=1/2k=1/\sqrt2 gives K=K′≈1.8541K=K'\approx1.8541 and a 2:1 block, and a square needs k=3−22≈0.1716k=3-2\sqrt2\approx0.1716.

Prerequisites. Analytic Functions and the Cauchy-Riemann Equations · Harmonic Functions and Conjugates · Consequences: Liouville, Maximum Modulus, Morera · Elementary Complex Functions. Leads to. Harmonic Functions and Conjugates · Complex Integration and Contour Integrals.

Ask the tutor: Ask me to walk through the accessory-parameter problem for a specific polygon, or to transplant a potential from an awkward region onto the half-plane step by step.

Why is the whole plane ℂ excluded from the Riemann Mapping Theorem?

Because a biholomorphism f:C→Df:\mathbb{C}\to\mathbb{D} would be a nonconstant entire function bounded by 11, and Liouville's theorem says the only bounded entire functions are constants. Topologically ℂ and the disk are both simply connected, so the obstruction is purely conformal: the disk is "small" in a way ℂ is not. This is not a technicality — it is the reason the disk, not the plane, is the universal model, and it reappears whenever you try to conformally shrink an unbounded region. The plane and the disk are the two conformally distinct simply connected surfaces; every other one is the disk.

What exactly is an accessory parameter, and why is it hard?

When you set up a Schwarz–Christoffel map, the interior angles fix all the exponents immediately, but the vertex preimages aka_k on the real axis are still unknown. Three of them can be placed freely by the three-parameter Möbius freedom of the half-plane (say at 0,1,∞0,1,\infty), so a polygon with nn vertices leaves n−3n-3 genuinely undetermined preimages: these are the accessory parameters. They are fixed by demanding the correct side lengths, which is a system of n−3n-3 nonlinear equations involving the map's integral. Triangles have none; anything with four or more sides needs numerical root-finding, and the preimages tend to crowd exponentially, which is what makes the computation delicate.

Can I ever write the Riemann map down explicitly?

Rarely in elementary terms. The theorem is an existence result whose standard proof is non-constructive — it extracts an extremal map from a normal family and never produces a formula. For special domains you can compose known maps: polygons give Schwarz–Christoffel integrals (elliptic or Beta functions), the half-plane and disk swap via the Cayley transform, and wedges use powers. But for a generic wiggly boundary there is no closed form, which is exactly why numerical conformal mapping — the SC Toolbox, the zipper algorithm — is a research field. The theorem guarantees the map exists and is unique once normalized; finding it is a separate computational task.

Does conformality really break down at a polygon corner?

Yes, and predictably so. At a vertex with interior angle α\alpha the Schwarz–Christoffel derivative behaves like (w−ak)α/π−1(w-a_k)^{\alpha/\pi-1}: for a convex corner (α<π\alpha<\pi) the exponent is negative and f′f' blows up, for a reentrant corner (α>π\alpha>\pi) it is positive and f′f' vanishes. Either way f′f' is not a nonzero finite number there, so angles are not preserved at that single point. Everywhere else in the interior the map is perfectly conformal. This is not a defect but a feature: the corner is exactly where a physical field (electric field, fluid velocity) develops its characteristic singularity, and the exponent tells you how strong it is.

What goes wrong for an annulus or other multiply connected region?

The theorem simply does not apply — it requires simple connectivity. An annulus {r1<∣z∣<r2}\{r_1<|z|<r_2\} cannot be mapped conformally onto a disk, and two annuli are conformally equivalent only when their radius ratios r2/r1r_2/r_1 match. That ratio is a genuine conformal invariant, the conformal modulus, so multiply connected domains carry extra continuous data no normalization can remove (a region with mm holes has roughly 3m−33m-3 moduli). Physically this is why the capacitance of a coaxial geometry depends on a genuine ratio, and it is the same modulus that appears as the elliptic parameter kk in the rectangle map.

How does conformal mapping actually solve a physics problem in practice?

In three moves. First, identify a hard region (a corner, a slotted channel, a polygonal capacitor) carrying a harmonic field — electrostatic potential, stream function, steady temperature. Second, find a conformal map ff onto a simple region where the same boundary-value problem has an elementary solution ϕ\phi, typically linear or logarithmic. Third, pull it back: Φ=ϕ∘f\Phi=\phi\circ f is automatically harmonic and matches the transplanted boundary data, so it is the answer on the original region. Worked example 1 does exactly this with ζ=z2\zeta=z^2. Equipotentials and field lines, being orthogonal, map to orthogonal families, so the whole field picture transfers intact.

Before this

Where this leads