View source for Beltrami equation

{{Short description|Partial differential equation}}
{{Distinguish|Laplace–Beltrami equation}}
In [[mathematics]], the '''Beltrami equation''', named after [[Eugenio Beltrami]], is the [[partial differential equation]]

:<math> {\partial w\over \partial \overline{z}} = \mu {\partial w\over \partial z}.</math>

for ''w'' a complex distribution of the [[complex variable]] ''z'' in some open set ''U'', with derivatives that are locally [[L2 function|''L''<sup>2</sup>]], and where ''μ'' is a given complex function in [[L∞ function|''L''<sup>∞</sup>]](''U'') of norm less than 1, called the '''Beltrami coefficient''', and where <math>\partial / \partial z</math> and <math>\partial / \partial \overline{z}</math> are [[Wirtinger derivatives]]. Classically this differential equation was used by [[Gauss]] to prove the existence locally of [[isothermal coordinates]] on a surface with analytic [[Riemannian metric]]. Various techniques have been developed for solving the equation. The most powerful, developed in the 1950s, provides global solutions of the equation on '''C''' and relies on the [[Lp space|L<sup>''p''</sup>]] theory of the [[Beurling transform]], a [[singular integral operator]] defined on L<sup>''p''</sup>('''C''') for all 1 < ''p'' < ∞. The same method applies equally well on the [[unit disk]] and [[upper half plane]] and plays a fundamental role in [[Teichmüller theory]] and the theory of [[quasiconformal mapping]]s. Various [[uniformization theorem]]s can be proved using the equation, including the [[measurable Riemann mapping theorem]] and the [[simultaneous uniformization theorem]]. The existence of [[conformal welding]]s can also be derived using the Beltrami equation. One of the simplest applications is to the [[Riemann mapping theorem]] for simply connected bounded open domains  in the complex plane. When the domain has smooth boundary, [[elliptic regularity]] for the equation can be used to show that the uniformizing map from the unit disk to the domain extends to a C<sup>∞</sup> function from the closed disk to the closure of the domain.

==Metrics on planar domains==

Consider a 2-dimensional [[Riemannian manifold]], say with an (''x'', ''y'') coordinate system on it.  The curves of constant ''x'' on that surface typically don't intersect the curves of constant ''y'' orthogonally.  A new coordinate system (''u'', ''v'') is called [[Isothermal coordinates|isothermal]] when the curves of constant ''u'' do intersect the curves of constant ''v'' orthogonally and, in addition, the parameter spacing is the same — that is, for small enough ''h'', the little region with <math>a\le u\le a+h</math> and <math>b\le v\le b+h</math> is nearly square, not just nearly rectangular.  The Beltrami equation is the equation that has to be solved in order to construct isothermal coordinate systems.

To see how this works, let ''S'' be an open set in '''C''' and let

:<math> \displaystyle ds^2=E\,dx^2 + 2F\, dx\,dy + G\, dy^2 </math>

be a smooth metric ''g'' on ''S''.  The [[first fundamental form]] of ''g''

:<math>\displaystyle g(x,y)=\begin{pmatrix} E & F \\ F & G\end{pmatrix} </math>

is a positive real matrix (''E'' > 0, ''G'' > 0, ''EG'' − ''F''<sup>2</sup> > 0) that varies smoothly with ''x'' and ''y''.

The '''Beltrami coefficient''' of the metric ''g'' is defined to be

:<math>\displaystyle \mu(x,y)={E-G +2iF\over E+G +2\sqrt{EG-F^2}} </math>

This coefficient has modulus strictly less than one since the identity

:<math>\displaystyle (E-G+2iF)(E-G-2iF)=(E+G+2\sqrt{EG-F^2})(E+G-2\sqrt{EG-F^2})</math>

implies that

:<math>\displaystyle |\mu|^2={E+G - 2\sqrt{EG -F^2}\over E+G + 2\sqrt{EG-F^2}}<1. </math>

Let ''f''(''x'',''y'') =(''u''(''x'',''y''),''v''(''x'',''y'')) be a smooth diffeomorphism of ''S'' onto another open set ''T'' in '''C'''.  The map ''f'' preserves orientation just when its [[Jacobian matrix and determinant|Jacobian]] is positive:

:<math>\displaystyle u_x v_y-v_x u_y>0.</math>

And using ''f'' to pull back to ''S'' the standard Euclidean metric ''ds''<sup>2</sup> = ''du''<sup>2</sup> + ''dv''<sup>2</sup> on ''T'' induces a metric on ''S'' given by

:<math>\displaystyle ds^2=du^2 + dv^2=(u_x^2+v_x^2)\, dx^2 + 2 (u_xu_y + v_xv_y)\, dx \, dy + (u_y^2 + v_y^2)\, dy^2, </math>

a metric whose first fundamental form is

:<math>\displaystyle \begin{pmatrix} u_x^2+v_x^2 & u_x u_y+v_x v_y \\ u_x u_y+v_x v_y & u_y^2+v_y^2\end{pmatrix}.</math>

When ''f'' both preserves orientation and induces a metric that differs from the original metric ''g'' only by a positive, smoothly varying scale factor ''r''(''x'', ''y''), the new coordinates ''u'' and ''v'' defined on ''S'' by ''f'' are called '''[[isothermal coordinates]]'''.

To determine when this happens, we reinterpret ''f'' as a complex-valued function of a complex variable ''f''(''x''+i''y'') = ''u''(''x''+i''y'') + i''v''(''x''+i''y'') so that we can apply the [[Wirtinger derivatives]]:

:<math>\displaystyle \partial_z= {1\over 2} (\partial_x-i\partial_y),\,\, \partial_{\overline{z}}={1\over 2} (\partial_x +i\partial_y);\,\,\,\,dz=dx+i\,dy,\,\,d\overline{z}=dx-i\,dy. </math>

Since

:<math>\displaystyle f_z=((u_x+v_y)+i(v_x-u_y))/2</math>
:<math>\displaystyle f_{\overline{z}}=((u_x-v_y)+i(v_x+u_y))/2,</math>

the metric induced by ''f'' is given by

:<math>\displaystyle ds^2=|f_z \, dz + f_{\overline{z}}d\overline{z}|^2= |f_z|^2\, \left|dz + {f_{\overline{z}}\over f_z} \, d\overline{z}\right|^2. </math>

The '''Beltrami quotient''' of this induced metric is defined to be <math>f_{\overline{z}}/f_z</math>.

The Beltrami quotient <math>f_{\overline{z}}/f_z</math> of <math>f</math> equals the Beltrami coefficient <math>\mu(z)</math> of the original metric ''g'' just when

:<math>\displaystyle ((u_x+v_y)+i(v_x-u_y))\bigl(E-G+2iF\bigr)</math>
:::<math>=((u_x-v_y)+i(v_x+u_y))\bigl(E+G+2\sqrt{EG-F^2}\bigr).</math>

The real and imaginary parts of this identity linearly relate <math>u_x,</math> <math>u_y,</math> <math>v_x,</math> and <math>v_y,</math> and solving for <math>u_y</math> and <math>v_y</math> gives

:<math>\displaystyle u_y=\frac{F u_x-\sqrt{EG-F^2}\,v_x}{E}\quad\text{and}
\quad v_y=\frac{\sqrt{EG-F^2}\,u_x+F v_x}{E}.</math>

It follows that the metric induced by ''f'' is then ''r''(''x'', ''y'') ''g''(''x'',''y''), where <math>r=(u_x^2+v_x^2)/E,</math> which is positive, while the Jacobian of ''f'' is then <math>r\sqrt{EG-F^2},</math> which is also positive.  So, when <math>f_{\overline{z}}/f_z=\mu(z),</math> the new coordinate system given by ''f'' is isothermal.

Conversely, consider a diffeomorphism ''f'' that does give us isothermal coordinates.  We then have

:<math>\displaystyle \mu(z)=
\frac{(u_x^2+v_x^2)-(u_y+v_y)^2+2i(u_xu_y+v_xv_y)}
{(u_x^2+v_x^2)+(u_y^2+v_y)^2+2\sqrt{(u_x^2+v_x^2)(u_y^2+v_y^2)-(u_xu_y+v_xv_y)^2}},</math>

where the scale factor ''r''(''x'', ''y'') has dropped out and the expression inside the square root is the perfect square <math>u_x^2v_y^2-2u_xv_xu_yv_y+v_x^2u_y^2.</math>  Since ''f'' must preserve orientation to give isothermal coordinates, the Jacobian <math>u_xv_y-v_xu_y</math> is the positive square root; so we have

:<math>\displaystyle \mu(z)=
\frac{\bigl((u_x+iu_y)+i(v_x+iv_y)\bigr)\bigl((u_x+iu_y)-i(v_x+iv_y)\bigr)}
{\bigl((u_x+v_y)+i(v_x-u_y)\bigr)\bigl((u_x+v_y)-i(v_x-u_y)\bigr)}.</math>

The right-hand factors in the numerator and denominator are equal and, since the Jacobian is positive, their common value can't be zero; so <math>\mu(z)=f_{\overline{z}}/f_z.</math>

Thus, the local coordinate system given by a diffeomorphism ''f'' is isothermal just when ''f'' solves the Beltrami equation for <math>\mu(z).</math>

== Isothermal coordinates for analytic metrics ==

[[Gauss]] proved the existence of isothermal coordinates locally in the analytic case by reducing the Beltrami to an ordinary differential equation in the complex domain.<ref>{{harvnb|Spivak|1999|pp=314–317, which is pp. 455–460 in the first or second edition; but note that there is a typo in equation (**) on page 315 or 457.  The right-hand side, given as −β/α, should be −α/β.}}</ref>  Here is a cookbook presentation of Gauss's technique.

An isothermal coordinate system, say in a neighborhood of the origin (''x'', ''y'') = (0, 0), is given by the real and imaginary parts of a complex-valued function ''f''(''x'', ''y'') that satisfies

:<math>\displaystyle {f_{\overline{z}}\over f_z}=\mu(z)=\frac{E-G+2iF}{E+G+2\sqrt{EG-F^2}}. </math>

Let <math>f</math> be such a function, and let <math>\psi</math> be a complex-valued function of a complex variable that is [[holomorphic function|holomorphic]] and whose derivative is nowhere zero.  Since any holomorphic function <math>\psi</math> has <math>\psi_{\overline{z}}</math> identically zero, we have

:<math>
\begin{align}
\frac{(\psi\circ f)_{\overline{z}}}{(\psi\circ f)_z} & = \frac{(\psi_z\circ f)f_{\overline{z}}+(\psi_{\overline{z}}\circ f)\overline{f_{\overline{z}}}}{(\psi_z\circ f)f_z+(\psi_{\overline{z}}\circ f)\overline{f_z}}
=\frac{f_{\overline{z}}}{f_z}=\mu(z).
\end{align}
</math>

Thus, the coordinate system given by the real and imaginary parts of <math>\psi\circ f</math> is also isothermal.  Indeed, if we fix <math>f</math> to give one isothermal coordinate system, then all of the possible isothermal coordinate systems are given by <math>\psi\circ f</math> for the various holomorphic <math>\psi</math> with nonzero derivative.

When ''E'', ''F'', and ''G'' are real analytic, Gauss constructed a particular isothermal coordinate system <math>h(x,y)=u(x,y)+iv(x,y),</math> the one that he chose being the one with  <math>h(x,0)=x</math> for all ''x''.  So the ''u'' axis of his isothermal coordinate system coincides with the ''x'' axis of the original coordinates and is parameterized in the same way.  All other isothermal coordinate systems are then of the form <math>\psi\circ h</math> for a holomorphic <math>\psi</math> with nonzero derivative.

Gauss lets ''q''(''t'') be some complex-valued function of a real variable ''t'' that satisfies the following ordinary differential equation:

:<math>\displaystyle q'(t)=\frac{-F-i\sqrt{EG-F^2}}{E}(q(t),t),</math>

where ''E'', ''F'', and ''G'' are here evaluated at ''y'' = ''t'' and ''x'' = ''q''(''t'').  If we specify the value of ''q''(''s'') for some start value ''s'', this differential equation determines the values of ''q''(''t'') for ''t'' either less than or greater than ''s''.  Gauss then defines his isothermal coordinate system ''h'' by setting ''h''(''x'', ''y'') to be <math>q(0)</math> along the solution path of that differential equation that passes through the point (''x'', ''y''), and thus has ''q''(''y'') = ''x''.

This rule sets ''h''(''x'', 0) to be <math>x</math>, since the starting condition is then ''q''(0)=''x''.  More generally, suppose that we move by an infinitesimal vector (''dx'', ''dy'') away from some point (''x'', ''y''), where ''dx'' and ''dy'' satisfy

:<math>\displaystyle E\,dx+(F+i\sqrt{EG-F^2}\,)\,dy=0.</math>

Since <math>q'(t)=dx/dy</math>, the vector (''dx'', ''dy'') is then tangent to the solution curve of the differential equation that passes through the point (''x'', ''y'').  Because we are assuming the metric to be analytic, it follows that

:<math>\displaystyle dh =c(x,y)\bigl(E\,dx+(F+i\sqrt{EG-F^2}\,)\,dy\bigr)</math>

for some smooth, complex-valued function <math>c(x,y)=a(x,y)+ib(x,y).</math>  We thus have

:<math>\displaystyle h_z=((aE+bF+a\sqrt{EG-F^2}\,)+i(bE-aF+b\sqrt{EG-F^2}\,))/2</math>
:<math>\displaystyle h_{\overline{z}}=((aE-bF-a\sqrt{EG-F^2}\,)+i(bE+aF-b\sqrt{EG-F^2}\,))/2.</math>

We form the quotient <math>h_{\overline{z}}/h_z</math> and then multiply numerator and denominator by <math>\overline{h_z}</math>, which is the complex conjugate of the denominator.  Simplifying the result, we find that

:<math>\displaystyle
{h_{\overline{z}}\over h_z} = \frac{(a^2+b^2)E\;(E-G+2iF)}
{(a^2+b^2)E\;(E+G+2\sqrt{EG-F^2})}=\frac{E-G+2iF}{E+G+2\sqrt{EG-F^2}} =\mu(z).
</math>

Gauss's function ''h'' thus gives the desired isothermal coordinates.

==Solution in ''L''<sup>2</sup> for smooth Beltrami coefficients==
In the simplest cases the Beltrami equation can be solved using only Hilbert space techniques and the Fourier transform. The method of proof is the prototype for the general solution using L<sup>''p''</sup> spaces, although [[Adrien Douady]] has indicated a method for handling the general case using only Hilbert spaces: the method relies on the classical theory of [[quasiconformal mapping]]s to establish Hölder estimates that are automatic in the L<sup>''p''</sup> theory for ''p'' > 2.<ref>See:
*{{harvnb|Astala|Iwaniec|Martin|2009}}
*{{harvnb|Bers|John|Schechter|1979}}
*{{harvnb|Ahlfors|1966}}
*{{harvnb|Glutsyuk|2008}}
*{{harvnb|Douady|Buff|2000}}</ref>
Let ''T'' be the [[Beurling transform]] on L<sup>2</sup>('''C''') defined on the Fourier transform of an L<sup>2</sup> function ''f'' as a multiplication operator:

:<math>\displaystyle \widehat{Tf}(z)= {\overline{z}\over z} \widehat{f}(z). </math>

It is a unitary operator and if ''f'' is a tempered distribution on '''C''' with partial derivatives in
L<sup>2</sup> then

:<math>\displaystyle T(f_{\overline{z}})=f_z, </math>

where the subscripts denote complex partial derivatives.

The [[fundamental solution]] of the operator

:<math>D=\partial_{\overline{z}}</math>

is given by the distribution

:<math>\displaystyle{E(z) ={1\over \pi z},}</math>

a locally integrable function on '''C'''. Thus on [[Schwartz function]]s ''f''

:<math>\displaystyle{\partial_{\overline{z}}(E\star f) = f.}</math>

The same holds for distributions of compact support on '''C'''. In particular if ''f'' is an L<sup>2</sup> function with compact support, then its '''Cauchy transform''', defined as

:<math>\displaystyle{Cf=E\star f,}</math>

is locally square integrable. The above equation can be written

:<math>\displaystyle{(Cf)_{\overline{z}}=f.}</math>

Moreover, still regarding ''f'' and ''Cf'' as distributions,

:<math>\displaystyle (Cf)_z= Tf. </math>

Indeed, the operator ''D'' is given on Fourier transforms as multiplication by ''iz''/2 and ''C'' as multiplication by its inverse.

Now in the Beltrami equation

:<math>\displaystyle f_{\overline{z}} = \mu f_z, </math>

with ''μ'' a smooth function of compact support, set

:<math>\displaystyle g(z) = f(z) -z </math>

and assume that the first derivatives of ''g'' are L<sup>2</sup>. Let ''h'' = ''g''<sub>''z''</sub> = ''f''<sub>''z''</sub> – 1. Then

:<math>\displaystyle h=g_z=T(g_{\overline{z}})= T(f_{\overline{z}})=T(\mu f_z)=T(\mu h) + T\mu </math>

If ''A'' and ''B'' are the operators defined by

:<math>\displaystyle{AF=T\mu F,\,\,\,\, BF =\mu TF}</math>

then their operator norms are strictly less than 1 and

:<math>\displaystyle{(I-A)h = T\mu.}</math>

Hence

:<math>\displaystyle{h=(I-A)^{-1}T\mu,\,\,\, T^*h=(I-B)^{-1}\mu}</math>

where the right hand sides can be expanded as [[Neumann series]]. It follows that

:<math>\displaystyle{g_{\overline{z}}=T^*h,}</math>

has the same support as ''μ'' and ''g''. Hence ''f'' is given by

:<math>\displaystyle{f(z) = CT^*h(z) + z.}</math>

[[Elliptic regularity]] can now be used to deduce that ''f'' is smooth.

In fact, off the support of ''μ'',

:<math>\displaystyle \partial_{\overline{z}} f =0, </math>

so by [[Weyl's lemma (Laplace equation)|Weyl's lemma]] ''f'' is even holomorphic for |''z''| > ''R''. Since ''f'' = ''CT*h'' + ''z'', it follows that
''f'' tends to 0 uniformly as |''z''| tends to ∞.

The elliptic regularity argument to prove smoothness, however, is the same everywhere and uses the theory of L<sup>2</sup> Sobolev spaces on the torus.<ref>See:
*{{harvnb|Bers|John|Schechter|1979}}
*{{harvnb|Glutsyuk|2008}}</ref> Let ψ be a smooth function of compact support on '''C''', identically equal to 1 on a neighbourhood of the support of ''μ'' and set ''F'' = ''ψ'' ''f''. The support of ''F'' lies in a large square |''x''|, |''y''| ≤ ''R'', so, identifying opposite sides of the square, ''F'' and ''μ'' can be regarded as a distribution and smooth function on a torus '''T'''<sup>2</sup>. By construction ''F'' is in ''L''<sup>2</sup>('''T'''<sup>2</sup>). As a distribution on '''T'''<sup>2</sup> it satisfies

:<math>\displaystyle F_{\overline{z}}= \mu F_{z} + GF, </math>

where ''G'' is smooth. On the canonical basis ''e''<sub>''m''</sub> of L<sup>2</sup>('''T'''<sup>2</sup>) with ''m'' in '''Z''' + ''i'' '''Z''', define

:<math>\displaystyle{Ue_m={\overline{m}\over m}e_m \,\, (m\ne 0), Ue_0=e_0.}</math>

Thus ''U'' is a unitary and on trigonometric polynomials or smooth functions ''P''

:<math>\displaystyle{U(P_{\overline{z}})=P_z.}</math>

Similarly it extends to a unitary on each [[Sobolev space]] H<sub>''k''</sub>('''T'''<sup>2</sup>) with the same property. It is the counterpart on the torus of the Beurling transform. The standard theory of [[Fredholm operator]]s shows that the operators corresponding to ''I'' – ''μ'' ''U'' and ''I'' – ''U'' ''μ'' are invertible on each Sobolev space. On the other hand,

:<math>\displaystyle{\partial_{z}(I-U\mu)F= UG.}</math>

Since ''UG'' is smooth, so too is (''I'' – ''μU'')''F'' and hence also ''F''.

Thus the original function ''f'' is smooth. Regarded as a map of '''C''' = '''R'''<sup>2</sup> into itself, the Jacobian is given by

:<math>\displaystyle{J(f)=|f_z|^2-|f_{\overline{z}}|^2=|f_z|^2(1-|\mu|^2).}</math>

This Jacobian is nowhere vanishing by a classical argument of {{harvtxt|Ahlfors|1966}}. In fact formally writing
''f''<sub>''z''</sub> = ''e''<sup>''k''</sup>, it follows that

:<math>\displaystyle{k_{\overline{z}} = \mu k_{z} + \mu_{z}.}</math>

This equation for ''k'' can be solved by the same methods as above giving a solution tending to 0 at ∞.
By uniqueness ''h'' + 1 = ''e''<sup>''k''</sup> so that

:<math>\displaystyle{J(f)=(1-|\mu|^2)|e^{2 k}|}</math>

is nowhere vanishing. Since ''f'' induces a smooth map of the Riemann sphere '''C''' ∪ ∞ into itself which is locally a diffeomorphism, ''f'' must be a diffeomorphism. In fact ''f'' must be onto by connectedness of the sphere, since its image is an open and closed subset; but then, as a [[covering map]], ''f'' must cover each point of the sphere the same number of times. Since only ∞ is sent to ∞, it follows that ''f'' is one-to-one.

The solution ''f'' is a quasiconformal conformal diffeomorphism. These form a group and their Beltrami coefficients can be computed according to the following rule:<ref>See:
*{{harvnb|Ahlfors|1966|p=9}}
*{{harvnb|Imayoshi|Taniguchi|1992|p=88}}</ref>

:<math>\displaystyle{\mu_{g\circ f^{-1}}\circ f={f_z\over \overline{f_z}} {\mu_g-\mu_f\over 1 -\overline{\mu_f}\mu_g}.}</math>

Moreover, if ''f''(0) = 0 and

:<math>\displaystyle{g(z)=f(z^{-1})^{-1},}</math>

then<ref>{{harvnb|Ahlfors|1966|p=98}}</ref>

:<math>\displaystyle{\mu_g(z)={z^2\over \overline{z}^2} \mu_f(z^{-1}).}</math>

This formula reflects the fact that on a [[Riemann surface]], a Beltrami coefficient is not a function.
Under a holomorphic change of coordinate ''w'' = ''w''(''z''), the coefficient is transformed to

:<math>\displaystyle{\widetilde{\mu}(w)= (w^\prime/ \overline{w^\prime})\cdot\mu(z).}</math>

Defining a smooth Beltrami coefficient on the sphere in this way, if ''μ'' is such a coefficient then, taking a smooth [[bump function]] ψ equal to 0 near 0, equal 1 for |''z''| > 1 and satisfying 0 ≤ ''ψ'' ≤ 1, ''μ'' can be written as a sum of two Beltrami coefficients:

:<math>\displaystyle{\mu=\psi\mu + (1-\psi) \mu=\mu_0 + \mu_\infty.}</math>

Let ''g'' be the quasiconformal diffeomorphism of the sphere fixing 0 and ∞ with coefficient 
''μ''<sub>∞</sub>. Let λ be the Beltrami coefficient of compact support on '''C''' defined by

:<math>\displaystyle \lambda(z)= \left[{\mu_0 \over 1 -\mu\overline{\mu_\infty}}     {g_z\over\overline{g_z}}\right] \circ g^{-1}(z). </math>

If ''f'' is the quasiconformal diffeomorphism of the sphere fixing 0 and ∞ with coefficient λ, then
the transformation formulas above show that ''f'' ∘ ''g''<sup>−1</sup> is a quasiconformal diffeomorphism of the sphere fixing 0 and ∞ with coefficient ''μ''.

The solutions of Beltrami's equation restrict to diffeomorphisms of the upper halfplane or unit disk if the coefficient ''μ'' has extra symmetry properties;<ref>See
*{{harvnb|Ahlfors|1966|p=99}}
*{{harvnb|Bers|John|Schechter|1979|p=277}}
</ref> since the two regions are related by a Möbius transformation (the Cayley transform), the two cases are essentially the same.

For the upper halfplane Im ''z'' > 0, if ''μ'' satisfies

:<math> \displaystyle \mu(z)=\overline{\mu(\overline{z})}, </math>

then by uniqueness the solution ''f'' of the Beltrami equation satisfies

:<math>\displaystyle f(z)=\overline{f(\overline{z})}, </math>

so leaves the real axis and hence the upper halfplane invariant.

Similarly for the unit disc |''z''| < 1, if ''μ'' satisfies

:<math>\displaystyle \mu(z)={z^2\over \overline{z}^2}\overline{\mu(\overline{z}^{-1})}, </math>

then by uniqueness the solution ''f'' of the Beltrami equation satisfies

:<math>\displaystyle{f(z)=\overline{f(\overline{z}^{-1})}^{-1},}</math>

so leaves the unit circle and hence the unit disk invariant.

Conversely Beltrami coefficients defined on the closures of the upper halfplane or unit disk which satisfy these conditions on the boundary can be "reflected" using the formulas above. If the extended functions are smooth the preceding theory can be applied. Otherwise the extensions will be continuous but with a jump in the derivatives at the boundary. In that case the more general theory for measurable coefficients ''μ'' is required and is most directly handled within the L<sup>''p''</sup> theory.

==Smooth Riemann mapping theorem==
{{See also|Riemann mapping theorem}}
Let ''U'' be an open simply connected domain in the complex plane with smooth boundary containing 0 in its interior and let ''F'' be a diffeomorphism of the unit disk ''D'' onto ''U'' extending smoothly to the boundary and the identity on a neighbourhood of 0. Suppose that in addition the induced metric on the closure of the unit disk can be reflected in the unit circle to define a smooth metric on '''C'''. The corresponding Beltrami coefficient is then a smooth function on '''C''' vanishing near 0 and ∞ and satisfying

:<math>\displaystyle \mu(z)={z^2\over \overline{z}^2}\overline{\mu(\overline{z}^{-1})}^{-1}. </math>

The quasiconformal diffeomorphism ''h'' of '''C''' satisfying

:<math>\displaystyle{h_{\overline{z}}=\mu(z) h_z}</math>

preserves the unit circle together with its interior and exterior. From the composition formulas for Beltrami coefficients

:<math>\displaystyle{\mu_{F\circ h^{-1}}=0,}</math>

so that ''f'' = ''F''∘ ''h''<sup>−1</sup> is a smooth diffeomorphism between the closures of ''D'' and ''U'' which is holomorphic on the interior. Thus, if a suitable diffeomorphism ''F'' can be constructed, the mapping ''f'' proves the smooth [[Riemann mapping theorem]] for the domain ''U''.

To produce a diffeomorphism ''F'' with the properties above, it can be assumed after an affine transformation that the boundary of ''U'' has length 2π and that 0 lies in ''U''. The smooth version of the [[Schoenflies theorem]] produces a smooth diffeomorphism ''G'' from the closure of ''D'' onto the closure of ''u'' equal to the identity on a neighbourhood of 0 and with an explicit form on a tubular neighbourhood of the unit circle. In fact taking polar coordinates (''r'',''θ'') in ''R''<sup>2</sup> and letting (''x''(''θ''),''y''(''θ'')) (''θ'' in [0,2{{pi}}]) be a parametrization of ∂''U'' by arclength, ''G'' has the form

:<math>\displaystyle G(r,\theta)=(x(\theta)-(1-r)y^\prime(\theta),y(\theta) +(1-r)x^\prime(\theta)). </math>

Taking ''t'' = 1 − ''r'' as parameter, the induced metric near the unit circle is given by

:<math> \displaystyle ds^2=(1+t\kappa(\theta))^2 d\theta^2 + dt^2, </math>

where

:<math>\displaystyle \kappa(\theta)=y^{\prime\prime}(\theta)x^\prime(\theta) - x^{\prime\prime}(\theta) y^\prime(\theta) </math>

is the [[radius of curvature (mathematics)|curvature]] of the [[plane curve]] (''x''(''θ''),''y''(''θ'')).

Let

: <math>\displaystyle \alpha={1\over 2\pi}\int_0^{2\pi} \kappa(\theta)\, d\theta,\,\,\, h(\theta)=\kappa(\theta) -\alpha. </math>

After a change of variable in the ''t'' coordinate and a conformal change in the metric, the metric takes the form

:<math>\displaystyle ds^2 =(1 + t\psi(t)h(\theta))^2 d\theta^2 + dt^2, </math>

where ψ is an analytic real-valued function of ''t'':

:<math>\displaystyle \psi(t)=1 +a_1t +a_2 t^2 + \cdots </math>

A formal diffeomorphism sending (''θ'',''t'') to (''f''(''θ'',''t''),''t'') can be defined as a formal [[power series]] in ''t'':

:<math>\displaystyle f(\theta,t)=\theta + f_1(\theta)t + f_2(\theta)t^2 + \cdots=\theta + g(\theta,t) </math>

where the coefficients ''f''<sub>''n''</sub> are smooth functions on the circle. These coefficients can be defined by recurrence so that the transformed metric only has even powers of ''t'' in the coefficients. This condition is imposed by demanding that no odd powers of ''t'' appear in the formal power series expansion:

:<math> \left[1+t\psi(t)\left(\sum_{n\ge 0} h^{(n)} g^n/n!\right)\right] \cdot \left[(1+g^\prime) \, d\theta + \left(\sum_{n\ge 1} n f_n t^{n-1}\right) \, dt\right]. </math>

By [[Borel's lemma]], there is a diffeomorphism defined in a neighbourhood of the unit circle, ''t'' = 0, for which the formal expression ''f''(''θ'',''t'') is the [[Taylor series]] expansion in the ''t'' variable. It follows that, after composing with this diffeomorphism, the extension of the metric obtained by reflecting in the line ''t'' = 0 is smooth.

==Hölder continuity of solutions==
Douady and others have indicated ways to extend the ''L''<sup>2</sup> theory to prove the existence and uniqueness of solutions when the Beltrami coefficient ''μ'' is bounded and measurable with ''L''<sup>∞</sup> norm ''k'' strictly less than one. Their approach involved the theory of quasiconformal mappings to establish directly the solutions of Beltrami's equation when ''μ'' is smooth with fixed compact support are uniformly [[Hölder continuous]].<ref>See:
*{{harvnb|Douady|Buff|2000}}
*{{harvnb|Glutsyuk|2008}}
*{{harvnb|Ahlfors|1966}}</ref> In the L<sup>''p''</sup> approach Hölder continuity follows automatically from operator theory.

The ''L''<sup>''p''</sup> theory when ''μ'' is smooth of compact support proceeds as in the L<sup>2</sup> case. By the [[Singular integral operators of convolution type|Calderón–Zygmund theory]] the Beurling transform and its inverse are known to be continuous for the L<sup>''p''</sup> norm. The [[Riesz–Thorin theorem|Riesz–Thorin convexity theorem]] implies that the norms ''C''<sup>''p''</sup> are continuous functions of ''p''. In particular ''C''<sub>''p''</sub> tends to 1 when ''p'' tends to 2.

In the Beltrami equation

:<math>\displaystyle{f_{\overline{z}} = \mu f_z,}</math>

with ''μ'' a smooth function of compact support, set

:<math>\displaystyle g(z) = f(z) -z </math>

and assume that the first derivatives of ''g'' are L<sup>''p''</sup>. Let ''h'' = ''g''<sub>''z''</sub> = ''f''<sub>''z''</sub> – 1. Then

:<math>\displaystyle{h=g_z=T(g_{\overline{z}})= T(f_{\overline{z}})=T(\mu f_z)=T(\mu h) + T\mu.}</math>

If ''A'' and ''B'' are the operators defined by ''AF'' = ''TμF'' and ''BF'' = ''μTF'', then their operator norms are strictly less than 1 and  (''I'' − ''A'')''h'' = ''T''μ. Hence

:<math>\displaystyle h=(I-A)^{-1}T\mu,\,\,\, T^{-1}=(I-B)^{-1}\mu </math>

where the right hand sides can be expanded as [[Neumann series]]. It follows that

:<math>\displaystyle g_{\overline{z}}=T^{-1}h, </math>

has the same support as ''μ'' and ''g''. Hence, up to the addition of a constant, ''f'' is given by

:<math>\displaystyle f(z) = CT^{-1}h(z) + z. </math>

Convergence of functions with fixed compact support in the L<sup>''p''</sup> norm for ''p'' > 2 implies convergence in
L<sup>2</sup>, so these formulas are compatible with the L<sup>2</sup> theory if ''p'' > 2.

The Cauchy transform ''C'' is not continuous on L<sup>2</sup> except as a map into functions of [[vanishing mean oscillation]].
<ref>{{harvnb|Astala|Iwaniec|Martin|2009}}</ref> On L<sup>''p''</sup> its image is contained in Hölder continuous functions with Hölder exponent 1 − 2''p''<sup>−1</sup> once a suitable constant is added. In fact for a function ''f'' of compact support define

:<math>
\begin{align}
Pf(w) & =Cf(w) +\pi ^{-1} \iint_{\mathbf{C}} {f(z)\over z} \, dx \, dy \\[4pt]
& = -{1\over \pi} \iint_{\mathbf{C}} f(z)\left({1\over z-w} -{1\over z}\right)\, dx \, dy= -{1\over \pi}\iint_{\mathbf{C}} {f(z)w\over z(z-w)}\, dx \, dy.
\end{align}
</math>

Note that the constant is added so that ''Pf''(0)&nbsp;=&nbsp;0. Since ''Pf'' only differs from ''Cf'' by a constant, it follows exactly as in the ''L''<sup>2</sup> theory that

:<math>\displaystyle (Pf)_{\overline{z}} = f,\,\,\, (Pf)_z = Tf.</math>

Moreover, ''P'' can be used instead of ''C'' to produce a solution:

:<math>\displaystyle f(z)=PT^{-1}h(z) +z.</math>

On the other hand, the integrand defining ''Pf'' is in L<sup>''q''</sup> if ''q''<sup>−1</sup> = 1 − ''p''<sup>−1</sup>. The [[Hölder inequality]] implies that ''Pf'' is [[Hölder continuous]] with an explicit estimate:

:<math>\displaystyle |Pf(w_1) - Pf(w_2)| \le K_p \|f\|_p|w_1-w_2|^{1-2/p},</math>

where

:<math>\displaystyle K_p= \|z^{-1}(z-1)^{-1}\|_q/\pi.</math>

For any ''p'' > 2 sufficiently close to 2, ''C''<sub>''p''</sub>''k'' <1. Hence the Neumann series for (''I'' − ''A'')<sup>−1</sup> and (''I'' − ''B'')<sup>−1</sup> converge. The Hölder estimates for ''P'' yield the following uniform estimates for the normalized solution of the Beltrami equation:

:<math>\displaystyle |f(w_1)-f(w_2)|\le {K_p\over 1 -\|\mu\|_\infty C_p} \|\mu\|_p|w_1-w_2|^{1-2/p} + |w_1 -w_2|. </math>

If ''μ'' is supported in |''z''| ≤ ''R'', then

:<math>\displaystyle \|\mu\|_p\le (\pi R^2)^{1/p} \|\mu\|_\infty. </math>

Setting ''w''<sub>1</sub> = ''z'' and ''w''<sub>2</sub> = 0, it follows that for |''z''|  ≤ ''R''

:<math>\displaystyle |f(z)| \le   \left[ 1+ {K_p\pi^{1/p} \|\mu\|_\infty\over 1 -\|\mu\|_\infty C_p}\right]\cdot R = CR, </math>

where the constant ''C'' > 0 depends only on the L<sup>∞</sup> norm of ''μ''. So the Beltrami coefficient of ''f''<sup>−1</sup> is smooth and supported in
''z''| ≤ ''CR''. It has the same L<sup>∞</sup> norm as that of ''f''. So the inverse diffeomorphisms also satisfy uniform Hölder estimates.

==Solution for measurable Beltrami coefficients==

===Existence===
The theory of the Beltrami equation can be extended to measurable Beltrami coefficients ''μ''. For simplicity only a special class of ''μ'' will be considered—adequate for most applications—namely those functions which are smooth an open set Ω (the regular set) with complement Λ a closed set of measure zero (the singular set). Thus Λ is a closed set that is contained in open sets of arbitrarily small area. For measurable Beltrami coefficients ''μ'' with compact support in |''z''| < ''R'', the solution of the Beltrami equation can be obtained as a limit of solutions for smooth Beltrami coefficients.<ref>See:
*{{harvnb|Ahlfors|1966}}
*{{harvnb|Douady|Buff|2000}}</ref>

In fact in this case the singular set Λ is compact. Take smooth functions φ<sub>''n''</sub> of compact support with 0 ≤ φ<sub>''n''</sub> ≤ 1, equal to 1 on a neighborhood of Λ and 0 off a slightly larger neighbourhood, shrinking to Λ as ''n'' increases. Set

:<math>\displaystyle{\mu_n=(1-\varphi_n)\cdot\mu.}</math>

The ''μ''<sub>''n''</sub> are smooth with compact support in |''z''| < ''R'' and

:<math>\displaystyle{\|\mu_n\|_\infty \le \|\mu\|_\infty.}</math>

The ''μ''<sub>''n''</sub> tend to ''μ'' in any ''L''<sup>''p''</sup> norm with ''p'' < ∞.

The corresponding normalised solutions ''f''<sub>''n''</sub> of the Beltrami equations and their inverses ''g''<sub>''n''</sub> satisfy uniform Hölder estimates. They are therefore [[equicontinuous]] on any compact subset of '''C'''; they are even holomorphic for |''z''| > ''R''. So by the [[Arzelà–Ascoli theorem]], passing to a subsequence if necessary, it can be assumed that both ''f''<sub>''n''</sub> and ''g''<sub>''n''</sub> converge uniformly on compacta to ''f'' and ''g''. The limits will satisfy the same Hölder estimates and be holomorphic for |''z''| > ''R''. The relations ''f''<sub>''n''</sub>∘''g''<sub>''n''</sub> = id = ''g''<sub>''n''</sub>∘''f''<sub>''n''</sub> imply that in the limit ''f''∘''g'' = id = ''g''∘''f'', so that ''f'' and ''g'' are homeomorphisms.

*The limits ''f'' and ''g'' are weakly differentiable.<ref>{{harvnb|Douady|Buff|2000|pp=319–320}}</ref> In fact let

::<math>\displaystyle{u_n=\partial_{\overline{z}} f_n,\,\, v_n =\partial_z f_n-1.}</math>

:These lie in L<sup>p</sup> and are uniformly bounded:

::<math>\displaystyle{\|u_n\|_p,\,\, \|v_n\|_p \le {\|\mu_n\|_p\over 1 - \|\mu_n\|_\infty}.}</math>

:Passing to a subsequence if necessary, it can be assumed that the sequences have weak limits ''u'' and ''v'' in L<sup>p</sup>. These are the distributional derivatives of ''f''(''z'') – ''z'', since if ψ is smooth of compact support

::<math>\displaystyle{\iint u \cdot \psi =\lim\iint u_n\cdot \psi =-\lim \iint f_n \cdot \partial_{\overline {z}} \psi = -\iint f\cdot   \partial_{\overline {z}} \psi,}</math>

:and similarly for ''v''. A similar argument applies for the ''g'' using the fact that Beltrami coefficients of the ''g''<sub>''n''</sub> are supported in a fixed closed disk.

*''f'' satisfies the Beltrami equation with Beltrami coefficient ''μ''.<ref>{{harvnb|Douady|Buff|2000|pp=319–320}}</ref> In fact the relation ''u'' = ''μ'' ⋅ ''v'' + ''μ'' follows by continuity from the relation ''u''<sub>''n''</sub> = ''μ''<sub>''n''</sub> ⋅ ''v''<sub>''n''</sub> + ''μ''<sub>''n''</sub>. It suffices to show that ''μ''<sub>''n''</sub> ⋅ ''v''<sub>''n''</sub> tends weakly to ''μ'' ⋅ ''v''. The difference can be written

::<math>\displaystyle{\mu v -\mu_n v_n= \mu(v-v_n) +(\mu-\mu_n)v_n.}</math>

:The first term tends weakly to 0, while the second term equals ''μ φ''<sub>''n''</sub> ''v''<sub>''n''</sub>. The terms are uniformly bounded in ''L''<sup>p</sup>, so to check weak convergence to 0 it enough to check inner products with a dense subset of ''L''<sup>2</sup>. The inner products with functions of  compact support in Ω are zero for ''n'' sufficiently large.

*''f'' carries closed sets of measure zero onto closed sets of measure zero.<ref>{{harvnb|Ahlfors|1966|pp=97–98}}</ref> It suffices to check this for a compact set ''K'' of measure zero. If ''U'' is a bounded open set containing ''K'' and ''J'' denotes the Jacobian of a function, then

:: <math>
\begin{align}
A(f_n(U)) & = \iint_U J(f_n)\, dx\, dy =\iint_U |\partial_z f_n|^2 -|\partial_{\overline{z}}f_n|^2\, dx\,dy \\[4pt]
& \le \iint_U  |\partial_z f_n|^2\, dx\, dy\le \|\partial_z f_n|_U\|_p^2 \,A(U)^{1-2/p}.
\end{align}
</math>

:Thus if ''A''(''U'') is small, so is ''A''(''f''<sub>''n''</sub>(''U'')). On the other hand  ''f''<sub>''n''</sub>(''U'') eventually contains ''f''(''K''), for applying the inverse ''g''<sub>''n''</sub>, ''U'' eventually contains ''g''<sub>''n''</sub> ∘''f'' (''K'') since ''g''<sub>''n''</sub> ∘''f'' tends uniformly to the identity on compacta. Hence ''f''(''K'') has measure zero.
 
* ''f'' is smooth on the regular set Ω of ''μ''. This follows from the elliptic regularity results in the ''L''<sup>2</sup> theory.
* ''f'' has non-vanishing Jacobian there. In particular ''f''<sub>''z''</sub> ≠ 0 on Ω.<ref>{{harvnb|Douady|Buff|2000|p=321}}</ref> In fact for ''z''<sub>0</sub> in Ω, if ''n'' is large enough

::<math>\displaystyle{\partial_{\overline{z}}  (f\circ g_n) = 0}</math>

:near ''z''<sub>1</sub> = ''f''<sub>''n''</sub>(''z''<sub>0</sub>). So ''h'' = ''f'' ∘ ''g''<sub>''n''</sub> is holomorphic near ''z''<sub>1</sub>. Since it is locally a homeomorphism, ''h'' ' (''z''<sub>1</sub>) ≠ 0. Since ''f'' =''h'' ∘ ''f''<sub>''n''</sub>. it follows that the Jacobian of ''f'' is non-zero at ''z''<sub>0</sub>. On the other hand ''J''(''f'') = |''f''<sub>''z''</sub>|<sup>2</sup> (1 − |μ|<sup>2</sup>), so ''f''<sub>''z''</sub> ≠ 0 at ''z''<sub>0</sub>.

* ''g'' satisfies the Beltrami equation with Beltrami coefficient

::<math>\displaystyle{\mu^\prime(f(z))=-{f_z\over \overline{f_z}}\cdot \mu(z)}</math>

:or equivalently

::<math>\displaystyle{\mu^\prime(w) =-{\overline{g_w}\over g_w} \cdot \mu(g(w))}</math>

:on the regular set Ω ' = ''f''(Ω), with corresponding singular set Λ ' = ''f''(Λ).

*''g'' satisfies the Beltrami equation for ''μ''&prime;. In fact ''g'' has weak distributional derivatives in 1 + L<sup>''p''</sup> and L<sup>''p''</sup>. Pairing with smooth functions of compact support in Ω, these derivatives coincide with the actual derivatives at points of Ω. Since Λ has measure zero, the distributional derivatives equal the actual derivatives in ''L''<sup>''p''</sup>. Thus ''g'' satisfies Beltrami's equation since the actual derivatives do.
*If ''f''* and ''f'' are solutions constructed as above for ''μ''* and ''μ'' then ''f''* ∘ ''f''<sup>−1</sup> satisfies the Beltrami equation for

::<math>\displaystyle \nu(f(z))={f_z\over \overline{f_z}} \, {\mu^* - \mu\over 1-\overline{\mu} \mu^*}, </math>

:defined on Ω ∩ Ω*. The weak derivatives of ''f''* ∘ ''f''<sup>−1</sup> are given by the actual derivatives on Ω ∩ Ω*. In fact this follows by approximating ''f''* and ''g'' = ''f''<sup>−1</sup> by ''f''*<sub>''n''</sub> and ''g''<sub>''n''</sub>. The derivatives are uniformly bounded in 1 + L<sup>''p''</sup> and L<sup>''p''</sup>, so as before weak limits give the distributional derivatives of ''f''* ∘ ''f''<sup>−1</sup>. Pairing with smooth functions of compact support in  Ω ∩ Ω*, these agree with the usual derivatives. So the distributional derivatives are given by the usual derivatives off Λ ∪ Λ*, a set of measure zero.

This establishes the '''existence''' of homeomorphic solutions of Beltrami's equation in the case of Beltrami coefficients of compact support. It also shows that the inverse homeomorphisms and composed homeomorphisms satisfy Beltrami equations and that all computations can be performed by restricting to regular sets.

If the support is not compact the same trick used in the smooth case can be used to construct a solution in terms of two homeomorphisms associated to compactly supported Beltrami coefficients. Note that, because of the assumptions on the Beltrami coefficient, a Möbius transformation of the extended complex plane can be applied to make the singular set of the Beltrami coefficient compact. In that case one of the homeomorphisms can be chosen to be a diffeomorphism.

===Uniqueness===
There are several proofs of the uniqueness of solutions of the Beltrami equation with a given Beltrami coefficient.<ref>See:
*{{harvnb|Ahlfors|1966}}
*{{harvnb|Imayoshi|Taniguchi|1992}}
*{{harvnb|Lehto|1987}}
*{{harvnb|Lehto|Virtanen|1973}}
*{{harvnb|Douady|Buff|2000|pp=321–322}}</ref> Since applying a Möbius transformation of the complex plane to any solution gives another solution, solutions can be normalised so that they fix 0, 1 and ∞. The method of solution of the Beltrami equation using the Beurling transform also provides a proof of uniqueness for coefficients of compact support ''μ'' and for which the distributional derivatives are in 1 + L<sup>''p''</sup> and L<sup>''p''</sup>. The relations

:<math>\displaystyle (P\psi)_{\overline{z}} = \psi, \,\, \, (P\psi)_{z} = T\psi </math>

for smooth functions ψ of compact support are also valid in the distributional sense for L<sup>''p''</sup> functions ''h'' since they can be written as L<sup>''p''</sup> of  ψ<sub>''n''</sub>'s. If ''f'' is a solution of the Beltrami equation with ''f''(0) = 0 and ''f''<sub>''z''</sub> - 1 in L<sup>''p''</sup> then

:<math>\displaystyle{F=f-P(f_{\overline{z}})}</math>

satisfies

:<math>\displaystyle{\partial_{\overline{z}}F =0.}</math>

So ''F'' is weakly holomorphic. Applying Weyl's lemma <ref>*{{harvnb|Astala|Iwaniec|Martin|2009}}</ref> it is possible to conclude that there exists a holomorphic function ''G'' that is equal to ''F'' almost everywhere. Abusing notation redefine ''F:=G''. The conditions ''F'' '(z) − 1 lies in L<sup>''p''</sup> and ''F''(0) = 0 force ''F''(''z'') = ''z''. Hence

:<math>\displaystyle{P(f_{\overline{z}}) = f(z) + z}</math>

and so differentiating

:<math>\displaystyle{f_z=T(\mu f_z) + 1.}</math>

If ''g'' is another solution then

:<math>\displaystyle{f_z-g_z=T(\mu(f_z-g_z)).}</math>

Since ''T''μ has operator norm on L<sup>''p''</sup> less than 1, this forces

:<math>\displaystyle{f_z=g_z.}</math>

But then from the Beltrami equation

:<math>\displaystyle{f_{\overline{z}}=g_{\overline{z}}.}</math>

Hence ''f'' − ''g'' is both holomorphic and antiholomorphic, so a constant. Since ''f''(0) = 0 = ''g''(0), it follows that ''f'' = ''g''. Note that since ''f'' is holomorphic off the support of ''μ'' and ''f''(∞) = ∞, the conditions that the derivatives are locally in L<sup>''p''</sup> force

:<math>\displaystyle{f_z -1,\,\,f_{\overline{z}}\in L^p(\mathbf{C}).}</math>

For a general ''f'' satisfying Beltrami's equation and with distributional derivatives locally in L<sup>''p''</sup>, it can be assumed after applying a Möbius transformation that 0 is not in the singular set of the Beltrami coefficient ''μ''. If ''g'' is a smooth diffeomorphism ''g'' with Beltrami coefficient λ supported near 0, the Beltrami coefficient ''ν'' for ''f'' ∘ ''g''<sup>−1</sup> can be calculated directly using the change of variables formula for distributional derivatives:

:<math>\displaystyle \nu(g(z))= {g_z\over \overline{g_z}} \,{\mu -\lambda\over 1-\overline{\lambda}\mu}. </math>

''λ'' can be chosen so that ν vanishes near zero. Applying the map ''z''<sup>−1</sup> results in a solution of Beltrami's equation with a Beltrami coefficient of compact support. The directional derivatives are still locally in L<sup>''p''</sup>. The coefficient ν depends only on ''μ'', ''λ'' and ''g'', so any two solutions of the original equation will produce solutions near 0 with distributional derivatives locally in ''L''<sup>''p''</sup> and the same Beltrami coefficient. They are therefore equal. Hence the solutions of the original equation are equal.

==Uniformization of multiply connected planar domains==

The method used to prove the smooth Riemann mapping theorem can be generalized to multiply connected planar regions with smooth boundary. The Beltrami coefficient in these cases is smooth on an open set, the complement of which has measure zero. The theory of the Beltrami equation with measurable coefficients is therefore required.<ref name="Bers 1961">{{harvnb|Bers|1961}}</ref><ref name="Sibner 1965">{{harvnb|Sibner|1965}}</ref>

'''Doubly connected domains.''' If Ω is a doubly connected planar region, then there is a diffeomorphism ''F'' of an annulus ''r'' ≤ |z| ≤ 1 onto the closure of Ω, such that after a conformal change the induced metric on the annulus can be continued smoothly by reflection in both boundaries. The annulus is a fundamental domain for the group generated by the two reflections, which reverse orientation. The images of the fundamental domain under the group fill out '''C''' with 0 removed and the Beltrami coefficient is smooth there. The canonical solution ''h'' of the Beltrami equation on '''C''', by the L<sup>''p''</sup> theory is a homeomorphism. It is smooth on away from 0 by elliptic regularity. By uniqueness it preserves the unit circle, together with its interior and exterior. Uniqueness of the solution also implies that reflection there is a conjugate  Möbius transformation ''g'' such that ''h'' ∘ ''R'' = ''g'' ∘ ''h'' where ''R'' denotes reflection in |''z''| = ''r''.  Composing with a Möbius transformation that fixes the unit circle it can be assumed that ''g'' is a reflection in a circle |''z''| = ''s'' with ''s'' < 1.  It follows that ''F'' ∘ ''h''<sub>−1</sub> is a smooth diffeomorphism of the annulus ''s'' ≤ |''z''| ≤ 1 onto the closure of Ω, holomorphic in the interior.<ref name="Sibner 1965"/>

'''Multiply connected domains.''' For regions with a higher degree of connectivity ''k'' + 1, the result is essentially [[Lipman Bers|Bers']] generalization of the '''retrosection theorem'''.<ref name="Bers 1961"/><ref name="Sibner 1965"/>  There is a smooth diffeomorphism ''F'' of the  region  Ω<sub>1</sub>, given by the unit disk with ''k'' open disks removed, onto the closure of Ω.  It can be assumed that 0 lies in the interior of the domain. Again after a modification of the diffeomorphism and conformal change near the boundary, the metric can be assumed to be compatible with reflection. Let ''G'' be the group generated by reflections in the boundary circles of Ω<sub>1</sub>. The interior of Ω<sub>1</sub> iz a fundamental domain for ''G''. Moreover, the index two normal subgroup ''G''<sub>0</sub> consisting of orientation-preserving mappings is a classical [[Schottky group]]. Its fundamental domain consists of the original fundamental domain with its reflection in the unit circle added. If the reflection is ''R''<sub>0</sub>, it is a [[free group]] with generators ''R''<sub>''i''</sub>∘''R''<sub>0</sub> where ''R''<sub>''i''</sub> are the reflections in the interior circles in the original domain. The images of the original domain by the ''G'', or equivalently the reflected domain by the Schottky group, fill out the regular set for the Schottky group. It acts properly discontinuously there. The complement is the [[limit set]] of ''G''<sub>0</sub>. It has measure zero. The induced metric on Ω<sub>1</sub> extends by reflection to the regular set. The corresponding Beltrami coefficient is invariant for the reflection group generated by the reflections ''R''<sub>''i''</sub> for ''i'' ≥ 0. Since the limit set has measure zero, the Beltrami coefficient extends uniquely to a bounded measurable function on '''C'''. smooth  on the regular set. The normalised solution of the Beltrami equation ''h'' is a smooth diffeomorphism of the closure of  Ω<sub>1</sub> onto itself preserving the unit circle, its exterior and interior. Necessarily ''h'' ∘ ''R''<sub>''i''</sub> = ''S''<sub>''i''</sub> ∘ ''h''. where ''S''<sub>''i''</sub> is the reflection in another circle in the unit disk. Looking at fixed points, the circles arising this way for different ''i'' must be disjoint. It follows that ''F'' ∘ ''h''<sup>−1</sup> defines a smooth diffeomorphism of the unit disc with the interior of these circles removed onto the closure of Ω, which is holomorphic in the interior.

==Simultaneous uniformization==
{{See also|Simultaneous uniformization theorem}}
{{harvtxt|Bers|1961}} showed that two compact Riemannian 2-manifolds ''M''<sub>1</sub>, ''M''<sub>2</sub> of genus ''g'' > 1 can be simultaneously uniformized.

As topological spaces ''M''<sub>1</sub> and ''M''<sub>2</sub> are homeomorphic to a fixed quotient of the upper half plane '''H''' by a discrete cocompact subgroup Γ of PSL(2,'''R'''). Γ can be identified with the [[fundamental group]] of the manifolds and '''H''' is a [[universal covering space]]. The homeomorphisms can be chosen to be piecewise linear on corresponding  triangulations. A result of {{harvtxt|Munkres|1960}} implies that the homeomorphisms can be adjusted near the edges and the vertices of the triangulation to produce diffeomorphisms.  The metric on ''M''<sub>1</sub> induces a metric on '''H''' which is  Γ-invariant. Let ''μ'' be the corresponding Beltrami coefficient on '''H'''. It can be extended to '''C''' by reflection

:<math>\displaystyle \widehat{\mu}(z) = \mu(z)\,\, (z\in \mathbf{H}),\,\, \widehat{\mu}(z)=0 \,\, (z\in\mathbf{R}),\,\, \widehat{\mu}(z)
=\overline{\mu(\overline{z})}\,\, (\overline{z} \in \mathbf{H}). </math>

It satisfies the invariance property

:<math>\displaystyle \widehat{\mu}(g(z))=  {\overline{g_z}\over g_z}\, \widehat{\mu}(z), </math>

for ''g'' in Γ. The solution ''f'' of the corresponding Beltrami equation defines a homeomorphism of '''C''', preserving the real axis and the upper and lower half planes. Conjugation of the group elements by ''f''<sup>−1</sup> gives a new cocompact subgroup Γ<sub>1</sub> of PSL(2,'''R'''). Composing the original diffeomorphism with the inverse of ''f'' then yield zero as the Beltrami coefficient. Thus the metric induced on '''H''' is invariant under Γ<sub>1</sub> and conformal to the [[Poincaré metric]] on '''H'''. It must therefore be given by multiplying by a positive smooth function that is Γ<sub>1</sub>-invariant. Any such function corresponds to a smooth function on ''M''<sub>1</sub>. Dividing the metric on ''M''<sub>1</sub> by this function results in a conformally equivalent metric on ''M''<sub>1</sub> which agrees with the Poincaré metric on '''H''' /  Γ<sub>1</sub>. In this way ''M''<sub>1</sub> becomes a [[compact Riemann surface]], i.e. is uniformized and inherits a natural complex structure.

With this conformal change in metric ''M''<sub>1</sub> can be identified with '''H''' /  Γ<sub>1</sub>. The diffeomorphism between onto ''M''<sub>2</sub> induces another metric on '''H''' which is invariant under Γ<sub>1</sub>. It defines a Beltrami coefficient λ
omn '''H''' which this time is extended to '''C''' by defining λ to be 0 off '''H'''. The solution ''h'' of the Beltrami equation is a homeomorphism of '''C''' which is holomorphic on the lower half plane and smooth on the upper half plane. The image of the real axis is a [[Jordan curve]] dividing '''C''' into two components. Conjugation of Γ<sub>1</sub> by ''h''<sup>−1</sup> gives a [[Quasi-Fuchsian group|quasi-Fuchsian subgroup]] Γ<sub>2</sub> of PSL(2,'''C'''). It leaves invariant the Jordan curve and acts properly discontinuously on each of the two components. The quotients of the two components by Γ<sub>2</sub> are naturally identified with ''M''<sub>1</sub> and ''M''<sub>2</sub>. This identification is compatible with the natural complex structures on both ''M''<sub>1</sub> and ''M''<sub>2</sub>.

==Conformal welding==
{{See also|Conformal welding|Douady–Earle extension}}
An orientation-preserving homeomorphism ''f'' of the circle is said to be [[quasisymmetric map|quasisymmetric]] if there are positive constants ''a'' and ''b'' such that

:<math>\displaystyle{{{|f(z_1)-f(z_2)|\over |f(z_1)-f(z_3)|} \le a\cdot {|z_1-z_2|^b\over|z_1-z_3|^b}}.}</math>

If

:<math>\displaystyle{|z_1-z_2|=|z_1-z_3|,}</math>

then the condition becomes

:<math>\displaystyle{a^{-1}\le {|f(z_1)-f(z_2)|\over |f(z_1)-f(z_3)|} \le a.}</math>

Conversely if this condition is satisfied for all such triples of points, then ''f'' is quasisymmetric.<ref>{{harvnb|Tukia|Väisälä|1980}}</ref>

An apparently weaker condition on a homeomorphism ''f'' of the circle is that it be ''quasi-Möbius'', that is there are constants ''c'', ''d'' > 0 such that

:<math>\displaystyle{|(f(z_1),f(z_2);f(z_3),f(z_4))| \le c |(z_1,z_2;z_3,z_4)|^d,}</math>

where

:<math>\displaystyle{ (z_1,z_2;z_3,z_4)={(z_1-z_3)(z_2-z_4)\over(z_2-z_3)(z_1-z_4)}}</math>

denotes the [[cross-ratio]]. In fact if ''f'' is quasisymmetric then it is also quasi-Möbius, with ''c'' = ''a''<sup>2</sup> and ''d'' = ''b'': this follows by multiplying the first inequality above for (''z''<sub>1</sub>,''z''<sub>3</sub>,''z''<sub>4</sub>) and (''z''<sub>2</sub>,''z''<sub>4</sub>,''z''<sub>3</sub>).

Conversely if ''f'' is a quasi-Möbius homeomorphism then it is also quasisymmetric.<ref>{{harvnb|Väisälä|1984}}</ref> Indeed, it is immediate that if ''f'' is quasi-Möbius so is its inverse. It then follows that ''f''  (and hence ''f''<sup>−1</sup>) is [[Hölder continuous]]. To see this let ''S'' be the set of cube roots of unity, so that if ''a'' ≠ ''b'' in ''S'', then |''a'' − ''b''| = 2 sin {{pi}}/3 = {{radic|3}}. To prove a Hölder estimate, it can be assumed that ''x'' – ''y'' is uniformly small. Then both ''x'' and ''y'' are greater than a fixed distance away from ''a'', ''b'' in ''S'' with ''a'' ≠ ''b'', so the estimate follows by applying the quasi-Möbius inequality to ''x'', ''a'', ''y'', ''b''. To check that ''f'' is quasisymmetric, it suffices to find a uniform upper bound for |''f''(''x'') − ''f''(''y'')| / |''f''(''x'') − ''f''(''z'')| in the case of a triple with |''x'' − ''z''| = |''x'' − ''y''|, uniformly small. In this case there is a point ''w'' at a distance greater than 1 from ''x'', ''y'' and ''z''. Applying the quasi-Möbius inequality to ''x'', ''w'', ''y'' and ''z'' yields the required upper bound.

A homeomorphism ''f'' of the unit circle can be extended to a homeomorphism ''F'' of the closed unit disk which is diffeomorphism on its interior. {{harvtxt|Douady|Earle|1986}}, generalizing earlier results of Ahlfors and Beurling, produced such [[Douady–Earle extension|an extension]] with the additional properties that it commutes with the action of SU(1,1) by Möbius transformations and is quasiconformal if ''f'' is quasisymmetric. (A less elementary method was also found independently by {{harvtxt|Tukia|1985}}: Tukia's approach has the advantage of also applying in higher dimensions.) When ''f'' is a diffeomorphism of the circle, the [[Alexander trick|Alexander extension]] provides another way of extending ''f'':

:<math>\displaystyle{F(r,\theta)=r \exp [i\psi(r)g(\theta) + i(1-\psi(r))\theta] ,}</math>

where ψ is a smooth function with values in [0,1], equal to 0 near 0 and 1 near 1, and

:<math>\displaystyle{f(e^{i\theta})=e^{ig(\theta)},}</math>

with ''g''(''θ'' + 2{{pi}}) = ''g''(''θ'') + 2{{pi}}. {{harvtxt|Partyka|Sakan|Zając|1999}} give a survey of various methods of extension, including  variants of the Ahlfors-Beurling extension which are smooth or analytic in the open unit disk.

In the case of a diffeomorphism, the Alexander extension ''F'' can be continued to any larger disk |''z''| < ''R'' with ''R'' > 1. Accordingly, in the unit disc

:<math>\displaystyle{\|\mu\|_\infty < 1, \,\,\, \mu(z)=F_{\overline{z}}/F_z.}</math>

This is also true for the other extensions when ''f'' is only quasisymmetric.

Now extend ''μ'' to a Beltrami coefficient on the whole of '''C''' by setting it equal to 0 for |''z''| ≥ 1. Let ''G'' be the corresponding solution of the Beltrami equation. Let ''F''<sub>1</sub>(''z'') = ''G'' ∘ ''F''<sup>−1</sup>(''z'') for |''z''| ≤ 1 and
''F''<sub>2</sub>(''z'') = ''G'' (''z'') for |''z''| ≥ 1. Thus ''F''<sub>1</sub> and ''F''<sub>2</sub> are univalent holomorphic maps of |''z''| < 1 and |''z''| > 1 onto the inside and outside of a Jordan curve. They extend continuously to homeomorphisms ''f''<sub>''i''</sub> of the unit circle onto the Jordan curve on the boundary. By construction they satisfy the
'''conformal welding''' condition:

:<math>\displaystyle{f=f_1^{-1}\circ f_2.}</math>

==See also==
* [[Quasiconformal mapping]]
* [[Measurable Riemann mapping theorem]]
* [[Isothermal coordinates]]

==Notes==
{{Reflist|3}}

==References==
*{{citation|last=Ahlfors|first= Lars V.|title=Conformality with respect to Riemannian metrics|
series=Ann. Acad. Sci. Fenn. Ser. A. I.|year= 1955|volume= 206}}
*{{citation|last=Ahlfors|first=Lars V.|authorlink=Lars Ahlfors|title=Lectures on quasiconformal mappings|publisher=Van Nostrand|year=1966}}
*{{citation|last1=Astala|first1=Kari|last2= Iwaniec|first2= Tadeusz|author2-link=Tadeusz Iwaniec|last3=Martin|first3= Gaven|author3-link=Gaven Martin|title=Elliptic partial differential equations and quasiconformal mappings in the plane|series= Princeton Mathematical Series|volume= 48|publisher= Princeton University Press|year= 2009|isbn= 978-0-691-13777-3}}
*{{Citation | last1=Beltrami | first1=Eugenio | title=Saggio di interpretazione della geometria non euclidea (Essay on the interpretation of noneuclidean geometry) | url=http://www.caressa.it/pdf/beltrami01.pdf | language=Italian | jfm=01.0275.02 | year=1867 | journal=Giornale di Mathematica | volume=6}} English translation in {{harvtxt|Stillwell|1996}}
*{{citation|last=Bers|first=Lipman|title= Riemann surfaces|year=1958|publisher=Courant Institute}}
*{{citation|last1=Bers|first1=Lipman|last2=John|first2=Fritz|last3= Schechter|first3= Martin|title=Partial differential equations, with supplements by Lars Gårding and A. N. Milgram|series= Lectures in Applied Mathematics|volume= 3A|publisher= American Mathematical Society|year=1979|isbn=0-8218-0049-3}}, Chapter VI.
*{{citation|last=Bers|first= Lipman|title= Uniformization by Beltrami equations|journal= Comm. Pure Appl. Math. |volume=14|year=1961|issue= 3|pages= 215–228|doi=10.1002/cpa.3160140304}}
*{{citation|last1=Douady|first1= Adrien|last2= Earle|first2 =Clifford J.|title=Conformally natural extension of homeomorphisms of the circle|journal=Acta Math.|volume= 157|year=1986|pages=23–48|doi=10.1007/bf02392590|doi-access=free}}
*{{citation|last1=Douady|first1= Adrien|authorlink=Adrien Douady|last2= Buff|first2= X.|title=Le théorème d'intégrabilité des structures presque complexes. [Integrability theorem for almost complex structures]|pages= 307–324|
series=London Math. Soc. Lecture Note Ser.|volume= 274|year= 2000|publisher =Cambridge Univ. Press}} 
*{{citation|last=Glutsyuk|first=Alexey A.|title=Simple proofs of uniformization theorems|pages= 125–143|journal=Fields Inst. Commun.|volume= 53|year=2008}}
*{{Citation | last1=Hubbard | first1=John Hamal | title=Teichmüller theory and applications to geometry, topology, and dynamics. Vol. 1 | url=http://matrixeditions.com/TeichmullerVol1.html | publisher=Matrix Editions, Ithaca, NY | isbn=978-0-9715766-2-9 | mr=2245223 | year=2006}}
*{{citation|first1=Y. |last1=Imayoshi|first2=M.|last2=Taniguchi|title=An Introduction to Teichmüller spaces|publisher=Springer-Verlag|year=1992|isbn=0-387-70088-9}}
*{{Citation | last1=Iwaniec | first1=Tadeusz | last2=Martin | first2=Gaven | title=The Beltrami equation | url=https://books.google.com/books?id=8iZczAQc59cC | series=Memoirs of the American Mathematical Society | isbn=978-0-8218-4045-0 | mr=2377904 | year=2008 | volume=191 | doi=10.1090/memo/0893| doi-access=free }}
*{{citation|first=Erwin|last=Kreyszig|authorlink=Erwin Kreyszig|year=1991|title=Differential Geometry|isbn=0-486-66721-9|publisher=Dover}}
*{{citation|last1=Lehto|first1= Olli|last2= Virtanen|first2=  K.I. |title=Quasiconformal mappings in the plane|edition=2nd|series= Die Grundlehren der mathematischen Wissenschaften|volume= 126|publisher= Springer-Verlag|year= 1973}}
*{{citation|last=Lehto|first= Olli|title=Univalent functions and Teichmüller spaces|series=Graduate Texts in Mathematics|volume= 109|publisher= Springer-Verlag|year= 1987|isbn=0-387-96310-3}}
*{{Citation | last1=Morrey | first1=Charles B. | title=On the solutions of quasi-linear elliptic partial differential equations. | doi=10.1090/S0002-9904-1936-06297-X | year=1936 | journal=[[Bulletin of the American Mathematical Society]] | issn=0002-9904 | volume=42 | issue=5 | pages=316 | jfm=62.0565.02| doi-access=free }}
*{{citation|title=On the Solutions of Quasi-Linear Elliptic Partial Differential Equations|first=Charles B. Jr.|last=Morrey|authorlink=Charles B. Morrey Jr.|journal=Transactions of the American Mathematical Society|volume=43|year=1938|pages=126–166|doi=10.2307/1989904|issue=1|jstor=1989904|publisher=American Mathematical Society|zbl=0018.40501|doi-access=free}}
*{{citation|last=Munkres|first= James|title=Obstructions to the smoothing of piecewise-differentiable homeomorphisms|
journal=Ann. of Math. |volume= 72|year= 1960 |issue= 3|pages=521–554|doi=10.2307/1970228|jstor= 1970228|url= http://projecteuclid.org/euclid.bams/1183523321}}
*   Papadopoulos, Athanase, ed. (2007), Handbook of Teichmüller theory. Vol. I, IRMA Lectures in Mathematics and Theoretical Physics, 11, European Mathematical Society (EMS), Zürich, {{doi|10.4171/029}}, {{ISBN|978-3-03719-029-6}}, {{MR|2284826}}
*   Papadopoulos, Athanase, ed. (2009), Handbook of Teichmüller theory. Vol. II, IRMA Lectures in Mathematics and Theoretical Physics, 13, European Mathematical Society (EMS), Zürich, {{doi|10.4171/055}}, {{ISBN|978-3-03719-055-5}}, {{MR|2524085}}
*   Papadopoulos, Athanase, ed. (2012), Handbook of Teichmüller theory. Vol. III, IRMA Lectures in Mathematics and Theoretical Physics, 19, European Mathematical Society (EMS), Zürich, {{doi|10.4171/103}}, {{ISBN|978-3-03719-103-3}}
*{{citation|last1=Partyka|first1=Dariusz|last2= Sakan|first2= Ken-Ichi|last3= Zając|first3= Józef|title=The harmonic and quasiconformal extension operators|pages= 141–177|journal=Banach Center Publ.|volume= 48|year= 1999|doi=10.4064/-48-1-141-177 |doi-access= free}}
*{{citation|last=Sibner|first= Robert J.|title=Uniformization of symmetric Riemann surfaces by Schottky groups|journal=Trans. Amer. Math. Soc.|volume= 116|year= 1965|pages= 79–85|doi=10.1090/s0002-9947-1965-0188431-2|doi-access=free}}
*{{citation|last=Spivak|first= Michael|title= A comprehensive introduction to differential geometry. Vol. IV|edition=3nd|publisher= Publish or Perish|year= 1999|isbn=0-914098-70-5}}
*{{Citation | last1=Stillwell | first1=John | author1-link=John Stillwell | title=Sources of hyperbolic geometry | url=https://books.google.com/books?id=ZQjBXxxQsucC | publisher=[[American Mathematical Society]] | location=Providence, R.I. | series=History of Mathematics | isbn=978-0-8218-0529-9 | mr=1402697 | year=1996 | volume=10}}
*{{citation|last1=Tukia|first1= P.|last2= Väisälä|first2= J.|title=Quasisymmetric embeddings of metric spaces|
journal=Ann. Acad. Sci. Fenn. Ser. A I Math.|volume= 5 |year=1980|pages= 97–114|doi= 10.5186/aasfm.1980.0531|doi-access= free}}
*{{citation|last=Tukia|first= Pekka|title=Quasiconformal extension of quasisymmetric mappings compatible with a Möbius group|journal=Acta Math. |volume=154|year=1985|issue= 3–4|pages= 153–193|doi=10.1007/bf02392471|doi-access=free}}
*{{citation|last=Väisälä|first= Jussi|title=Quasi-Möbius maps|journal=[[Journal d'Analyse Mathématique]] |volume=44|year=1984|pages= 218–234|doi=10.1007/bf02790198|doi-access=free|s2cid= 189767039}}
*{{citation|last=Vekua|first= I. N.|title= Generalized analytic functions|publisher= Pergamon Press|year= 1962}}

{{Authority control}}

[[Category:Partial differential equations]]
[[Category:Complex analysis]]
[[Category:Operator theory]]
[[Category:Moduli theory]]