An extremely large portion of mathematics is concerned with locating solutions to equations such as

\displaystyle  f(x) = 0


\displaystyle  \Phi(x) = x \ \ \ \ \ (1)

for {x} in some suitable domain space (either finite-dimensional or infinite-dimensional), and various maps {f} or {\Phi}. To solve the fixed point iteration equation (1), the simplest general method available is the fixed point iteration method: one starts with an initial approximate solution {x_0} to (1), so that {\Phi(x_0) \approx x_0}, and then recursively constructs the sequence {x_1, x_2, x_3, \dots} by {x_n := \Phi(x_{n-1})}. If {\Phi} behaves enough like a “contraction”, and the domain is complete, then one can expect the {x_n} to converge to a limit {x}, which should then be a solution to (1). For instance, if {\Phi: X \rightarrow X} is a map from a metric space {X = (X,d)} to itself, which is a contraction in the sense that

\displaystyle  d( \Phi(x), \Phi(y) ) \leq (1-\eta) d(x,y)

for all {x,y \in X} and some {\eta>0}, then with {x_n} as above we have

\displaystyle  d( x_{n+1}, x_n ) \leq (1-\eta) d(x_n, x_{n-1} )

for any {n}, and so the distances {d(x_n, x_{n-1} )} between successive elements of the sequence decay at at least a geometric rate. This leads to the contraction mapping theorem, which has many important consequences, such as the inverse function theorem and the Picard existence theorem.

A slightly more complicated instance of this strategy arises when trying to linearise a complex map {f: U \rightarrow {\bf C}} defined in a neighbourhood {U} of a fixed point. For simplicity we normalise the fixed point to be the origin, thus {0 \in U} and {f(0)=0}. When studying the complex dynamics {f^2 = f \circ f}, {f^3 = f \circ f \circ f}, {\dots} of such a map, it can be useful to try to conjugate {f} to another function {g = \psi^{-1} \circ f \circ \psi}, where {\psi} is a holomorphic function defined and invertible near {0} with {\psi(0)=0}, since the dynamics of {g} will be conjguate to that of {f}. Note that if {f(0)=0} and {f'(0)=\lambda}, then from the chain rule any conjugate {g} of {f} will also have {g(0)=0} and {g'(0)=\lambda}. Thus, the “simplest” function one can hope to conjugate {f} to is the linear function {z \mapsto \lambda z}. Let us say that {f} is linearisable (around {0}) if it is conjugate to {z \mapsto \lambda z} in some neighbourhood of {0}. Equivalently, {f} is linearisable if there is a solution to the Schröder equation

\displaystyle  f( \psi(z) ) = \psi(\lambda z) \ \ \ \ \ (2)

for some {\psi: U' \rightarrow {\bf C}} defined and invertible in a neighbourhood {U'} of {0} with {\psi(0)=0}, and all {z} sufficiently close to {0}. (The Schröder equation is normalised somewhat differently in the literature, but this form is equivalent to the usual form, at least when {\lambda} is non-zero.) Note that if {\psi} solves the above equation, then so does {z \mapsto \psi(cz)} for any non-zero {c}, so we may normalise {\psi'(0)=1} in addition to {\psi(0)=0}, which also ensures local invertibility from the inverse function theorem. (Note from winding number considerations that {\psi} cannot be invertible near zero if {\psi'(0)} vanishes.)

We have the following basic result of Koenigs:

Theorem 1 (Koenig’s linearisation theorem) Let {f: U \rightarrow {\bf C}} be a holomorphic function defined near {0} with {f(0)=0} and {f'(0)=\lambda}. If {0 < |\lambda| < 1} (attracting case) or {1 < |\lambda| < \infty} (repelling case), then {f} is linearisable near zero.

Proof: Observe that if {f, \psi, \lambda} solve (2), then {f^{-1}, \psi^{-1}, \lambda^{-1}} solve (2) also (in a sufficiently small neighbourhood of zero). Thus we may reduce to the attractive case {0 < |\lambda| < 1}.

Let {r>0} be a sufficiently small radius, and let {X} denote the space of holomorphic functions {\psi: B(0,r) \rightarrow {\bf C}} on the complex disk {B(0,r) := \{z \in {\bf C}: |z| < r \}} with {\psi(0)=0} and {\psi'(0)=1}. We can view the Schröder equation (2) as a fixed point equation

\displaystyle  \psi = \Phi(\psi)

where {\Phi: X' \rightarrow X} is the partially defined function on {X} that maps a function {\psi: B(0,r) \rightarrow {\bf C}} to the function {\Phi(\psi): B(0,r) \rightarrow {\bf C}} defined by

\displaystyle  \Phi(\psi)(z) := f^{-1}( \psi( \lambda z ) ),

assuming that {f^{-1}} is well-defined on the range of {\psi(B(0,\lambda r))} (this is why {\Phi} is only partially defined).

We can solve this equation by the fixed point iteration method, if {r} is small enough. Namely, we start with {\psi_0: B(0,r) \rightarrow {\bf C}} being the identity map, and set {\psi_1 := \Phi(\psi_0), \psi_2 := \Phi(\psi_1)}, etc. We equip {X} with the uniform metric {d( \psi, \tilde \psi ) := \sup_{z \in B(0,r)} |\psi(z) - \tilde \psi(z)|}. Observe that if {d( \psi, \psi_0 ), d(\tilde \psi, \psi_0) \leq r}, and {r} is small enough, then {\psi, \tilde \psi} takes values in {B(0,2r)}, and {\Phi(\psi), \Phi(\tilde \psi)} are well-defined and lie in {X}. Also, since {f^{-1}} is smooth and has derivative {\lambda^{-1}} at {0}, we have

\displaystyle  |f^{-1}(z) - f^{-1}(w)| \leq (1+\varepsilon) |\lambda|^{-1} |z-w|

if {z, w \in B(0,r)}, {\varepsilon>0} and {r} is sufficiently small depending on {\varepsilon}. This is not yet enough to establish the required contraction (thanks to Mario Bonk for pointing this out); but observe that the function {\frac{\psi(z)-\tilde \psi(z)}{z^2}} is holomorphic on {B(0,r)} and bounded by {d(\psi,\tilde \psi)/r^2} on the boundary of this ball (or slightly within this boundary), so by the maximum principle we see that

\displaystyle  |\frac{\psi(z)-\tilde \psi(z)}{z^2}| \leq \frac{1}{r^2} d(\psi,\tilde \psi)

on all of {B(0,r)}, and in particular

\displaystyle  |\psi(z)-\tilde \psi(z)| \leq |\lambda|^2 d(\psi,\tilde \psi)

on {B(0,\lambda r)}. Putting all this together, we see that

\displaystyle  d( \Phi(\psi), \Phi(\tilde \psi)) \leq (1+\varepsilon) |\lambda| d(\psi, \tilde \psi);

since {|\lambda|<1}, we thus obtain a contraction on the ball {\{ \psi \in X: d(\psi,\psi_0) \leq r \}} if {\varepsilon} is small enough (and {r} sufficiently small depending on {\varepsilon}). From this (and the completeness of {X}, which follows from Morera’s theorem) we see that the iteration {\psi_n} converges (exponentially fast) to a limit {\psi \in X} which is a fixed point of {\Phi}, and thus solves Schröder’s equation, as required. \Box

Koenig’s linearisation theorem leaves open the indifferent case when {|\lambda|=1}. In the rationally indifferent case when {\lambda^n=1} for some natural number {n}, there is an obvious obstruction to linearisability, namely that {f^n = 1} (in particular, linearisation is not possible in this case when {f} is a non-trivial rational function). An obstruction is also present in some irrationally indifferent cases (where {|\lambda|=1} but {\lambda^n \neq 1} for any natural number {n}), if {\lambda} is sufficiently close to various roots of unity; the first result of this form is due to Cremer, and the optimal result of this type for quadratic maps was established by Yoccoz. In the other direction, we have the following result of Siegel:

Theorem 2 (Siegel’s linearisation theorem) Let {f: U \rightarrow {\bf C}} be a holomorphic function defined near {0} with {f(0)=0} and {f'(0)=\lambda}. If {|\lambda|=1} and one has the Diophantine condition {\frac{1}{|\lambda^n-1|} \leq C n^C} for all natural numbers {n} and some constant {C>0}, then {f} is linearisable at {0}.

The Diophantine condition can be relaxed to a more general condition involving the rational exponents of the phase {\theta} of {\lambda = e^{2\pi i \theta}}; this was worked out by Brjuno, with the condition matching the one later obtained by Yoccoz. Amusingly, while the set of Diophantine numbers (and hence the set of linearisable {\lambda}) has full measure on the unit circle, the set of non-linearisable {\lambda} is generic (the complement of countably many nowhere dense sets) due to the above-mentioned work of Cremer, leading to a striking disparity between the measure-theoretic and category notions of “largeness”.

Siegel’s theorem does not seem to be provable using a fixed point iteration method. However, it can be established by modifying another basic method to solve equations, namely Newton’s method. Let us first review how this method works to solve the equation {f(x)=0} for some smooth function {f: I \rightarrow {\bf R}} defined on an interval {I}. We suppose we have some initial approximant {x_0 \in I} to this equation, with {f(x_0)} small but not necessarily zero. To make the analysis more quantitative, let us suppose that the interval {[x_0-r_0,x_0+r_0]} lies in {I} for some {r_0>0}, and we have the estimates

\displaystyle  |f(x_0)| \leq \delta_0 r_0

\displaystyle  |f'(x)| \geq \eta_0

\displaystyle  |f''(x)| \leq \frac{1}{\eta_0 r_0}

for some {\delta_0 > 0} and {0 < \eta_0 < 1/2} and all {x \in [x_0-r_0,x_0+r_0]} (the factors of {r_0} are present to make {\delta_0,\eta_0} “dimensionless”).

Lemma 3 Under the above hypotheses, we can find {x_1} with {|x_1 - x_0| \leq \eta_0 r_0} such that

\displaystyle  |f(x_1)| \ll \delta_0^2 \eta_0^{-O(1)} r_0.

In particular, setting {r_1 := (1-\eta_0) r_0}, {\eta_1 := \eta_0/2}, and {\delta_1 = O(\delta_0^2 \eta_0^{-O(1)})}, we have {[x_1-r_1,x_1+r_1] \subset [x_0-r_0,x_0+r_0] \subset I}, and

\displaystyle  |f(x_1)| \leq \delta_1 r_1

\displaystyle  |f'(x)| \geq \eta_1

\displaystyle  |f''(x)| \leq \frac{1}{\eta_1 r_1}

for all {x \in [x_1-r_1,x_1+r_1]}.

The crucial point here is that the new error {\delta_1} is roughly the square of the previous error {\delta_0}. This leads to extremely fast (double-exponential) improvement in the error upon iteration, which is more than enough to absorb the exponential losses coming from the {\eta_0^{-O(1)}} factor.

Proof: If {\delta_0 > c \eta_0^{C}} for some absolute constants {C,c>0} then we may simply take {x_0=x_1}, so we may assume that {\delta_0 \leq c \eta_0^{C}} for some small {c>0} and large {C>0}. Using the Newton approximation {f(x_0+h) \approx f(x_0) + h f'(x_0)} we are led to the choice

\displaystyle  x_1 := x_0 - \frac{f(x_0)}{f'(x_0)}

for {x_1}. From the hypotheses on {f} and the smallness hypothesis on {\delta} we certainly have {|x_1-x_0| \leq \eta_0 r_0}. From Taylor’s theorem with remainder we have

\displaystyle  f(x_1) = f(x_0) - \frac{f(x_0)}{f'(x_0)} f'(x_0) + O( \frac{1}{\eta_0 r_0} |\frac{f(x_0)}{f'(x_0)}|^2 )

\displaystyle  = O( \frac{1}{\eta_0 r_0} (\frac{\delta_0 r_0}{\eta_0})^2 )

and the claim follows. \Box

We can iterate this procedure; starting with {x_0,\eta_0,r_0,\delta_0} as above, we obtain a sequence of nested intervals {[x_n-r_n,x_n+r_n]} with {f(x_n)| \leq \delta_n}, and with {\eta_n,r_n,\delta_n,x_n} evolving by the recursive equations and estimates

\displaystyle  \eta_n = \eta_{n-1} / 2

\displaystyle  r_n = (1 - \eta_{n-1}) r_{n-1}

\displaystyle  \delta_n = O( \delta_{n-1}^2 \eta_{n-1}^{-O(1)} )

\displaystyle  |x_n - x_{n-1}| \leq \eta_{n-1} r_{n-1}.

If {\delta_0} is sufficiently small depending on {\eta_0}, we see that {\delta_n} converges rapidly to zero (indeed, we can inductively obtain a bound of the form {\delta_n \leq \eta_0^{C (2^n + n)}} for some large absolute constant {C} if {\delta_0} is small enough), and {x_n} converges to a limit {x \in I} which then solves the equation {f(x)=0} by the continuity of {f}.

As I recently learned from Zhiqiang Li, a similar scheme works to prove Siegel’s theorem, as can be found for instance in this text of Carleson and Gamelin. The key is the following analogue of Lemma 3.

Lemma 4 Let {\lambda} be a complex number with {|\lambda|=1} and {\frac{1}{|\lambda^n-1|} \ll n^{O(1)}} for all natural numbers {n}. Let {r_0>0}, and let {f_0: B(0,r_0) \rightarrow {\bf C}} be a holomorphic function with {f_0(0)=0}, {f'_0(0)=\lambda}, and

\displaystyle  |f_0(z) - \lambda z| \leq \delta_0 r_0 \ \ \ \ \ (3)

for all {z \in B(0,r_0)} and some {\delta_0>0}. Let {0 < \eta_0 \leq 1/2}, and set {r_1 := (1-\eta_0) r_0}. Then there exists an injective holomorphic function {\psi_0: B(0, r_1) \rightarrow B(0, r_0)} and a holomorphic function {f_1: B(0,r_1) \rightarrow {\bf C}} such that

\displaystyle  f_0( \psi_1(z) ) = \psi_1(f_1(z)) \ \ \ \ \ (4)

for all {z \in B(0,r_1)}, and such that

\displaystyle  |\psi_1(z) - z| \ll \delta_0 \eta_0^{-O(1)} r_1


\displaystyle  |f_1(z) - \lambda z| \leq \delta_1 r_1

for all {z \in B(0,r_1)} and some {\delta_1 = O(\delta_0^2 \eta_0^{-O(1)})}.

Proof: By scaling we may normalise {r_0=1}. If {\delta_0 > c \eta_0^C} for some constants {c,C>0}, then we can simply take {\psi_1} to be the identity and {f_1=f_0}, so we may assume that {\delta_0 \leq c \eta_0^C} for some small {c>0} and large {C>0}.

To motivate the choice of {\psi_1}, we write {f_0(z) = \lambda z + \hat f_0(z)} and {\psi_1(z) = z + \hat \psi(z)}, with {\hat f_0} and {\hat \psi_1} viewed as small. We would like to have {f_0(\psi_1(z)) \approx \psi_1(\lambda z)}, which expands as

\displaystyle  \lambda z + \lambda \hat \psi_1(z) + \hat f_0( z + \hat \psi_1(z) ) \approx \lambda z + \hat \psi_1(\lambda z).

As {\hat f_0} and {\hat \psi} are both small, we can heuristically approximate {\hat f_0(z + \hat \psi_1(z) ) \approx \hat f_0(z)} up to quadratic errors (compare with the Newton approximation {f(x_0+h) \approx f(x_0) + h f'(x_0)}), and arrive at the equation

\displaystyle  \hat \psi_1(\lambda z) - \lambda \hat \psi_1(z) = \hat f_0(z). \ \ \ \ \ (5)

This equation can be solved by Taylor series; the function {\hat f_0} vanishes to second order at the origin and thus has a Taylor expansion

\displaystyle  \hat f_0(z) = \sum_{n=2}^\infty a_n z^n

and then {\hat \psi_1} has a Taylor expansion

\displaystyle  \hat \psi_1(z) = \sum_{n=2}^\infty \frac{a_n}{\lambda^n - \lambda} z^n.

We take this as our definition of {\hat \psi_1}, define {\psi_1(z) := z + \hat \psi_1(z)}, and then define {f_1} implicitly via (4).

Let us now justify that this choice works. By (3) and the generalised Cauchy integral formula, we have {|a_n| \leq \delta_0} for all {n}; by the Diophantine assumption on {\lambda}, we thus have {|\frac{a_n}{\lambda^n - \lambda}| \ll \delta_0 n^{O(1)}}. In particular, {\hat \psi_1} converges on {B(0,1)}, and on the disk {B(0, (1-\eta_0/4))} (say) we have the bounds

\displaystyle  |\hat \psi_1(z)|, |\hat \psi'_1(z)| \ll \delta_0 \sum_{n=2}^\infty n^{O(1)} (1-\eta_0/4)^n \ll \eta_0^{-O(1)} \delta_0. \ \ \ \ \ (6)

In particular, as {\delta_0} is so small, we see that {\psi_1} maps {B(0, (1-\eta_0/4))} injectively to {B(0,1)} and {B(0,1-\eta_0)} to {B(0,1-3\eta_0/4)}, and the inverse {\psi_1^{-1}} maps {B(0, (1-\eta_0/2))} to {B(0, (1-\eta_0/4))}. From (3) we see that {f_0} maps {B(0,1-3\eta_0/4)} to {B(0,1-\eta_0/2)}, and so if we set {f_1: B(0,1-\eta_0) \rightarrow B(0,1-\eta_0/4)} to be the function {f_1 := \psi_1^{-1} \circ f_0 \circ \psi_1}, then {f_1} is a holomorphic function obeying (4). Expanding (4) in terms of {\hat f_0} and {\hat \psi_1} as before, and also writing {f_1(z) = \lambda z + \hat f_1(z)}, we have

\displaystyle  \lambda z + \lambda \hat \psi_1(z) + \hat f_0( z + \hat \psi_1(z) ) = \lambda z + \hat f_1(z) + \hat \psi_1(\lambda z + \hat f_1(z))

for {z \in B(0, 1-\eta_0)}, which by (5) simplifies to

\displaystyle  \hat f_1(z) = \hat f_0( z + \hat \psi_1(z) ) - \hat f_0(z) + \hat \psi_1(\lambda z) - \hat \psi_1(\lambda z + \hat f_1(z)).

From (6), the fundamental theorem of calculus, and the smallness of {\delta_0} we have

\displaystyle  |\hat \psi_1(\lambda z) - \hat \psi_1(\lambda z + \hat f_1(z))| \leq \frac{1}{2} |\hat f_1(z)|

and thus

\displaystyle  |\hat f_1(z)| \leq 2 |\hat f_0( z + \hat \psi_1(z) ) - \hat f_0(z)|.

From (3) and the Cauchy integral formula we have {\hat f'_0(z) = O( \delta_0 \eta_0^{-O(1)})} on (say) {B(0,1-\eta_0/4)}, and so from (6) and the fundamental theorem of calculus we conclude that

\displaystyle  |\hat f_1(z)| \ll \delta_0^2 \eta_0^{-O(1)}

on {B(0,1-\eta_0)}, and the claim follows. \Box

If we set {\eta_0 := 1/2}, {f_0 := f}, and {\delta_0>0} to be sufficiently small, then (since {f(z)-\lambda z} vanishes to second order at the origin), the hypotheses of this lemma will be obeyed for some sufficiently small {r_0}. Iterating the lemma (and halving {\eta_0} repeatedly), we can then find sequences {\eta_n, \delta_n, r_n > 0}, injective holomorphic functions {\psi_n: B(0,r_n) \rightarrow B(0,r_{n-1})} and holomorphic functions {f_n: B(0,r_n) \rightarrow {\bf C}} such that one has the recursive identities and estimates

\displaystyle  \eta_n = \eta_{n-1} / 2

\displaystyle  r_n = (1 - \eta_{n-1}) r_{n-1}

\displaystyle  \delta_n = O( \delta_{n-1}^2 \eta_{n-1}^{-O(1)} )

\displaystyle  |\psi_n(z) - z| \ll \delta_{n-1} \eta_{n-1}^{-O(1)} r_n

\displaystyle  |f_n(z) - \lambda z| \leq \delta_n r_n

\displaystyle  f_{n-1}( \psi_n(z) ) = \psi_n(f_n(z))

for all {n \geq 1} and {z \in B(0,r_n)}. By construction, {r_n} decreases to a positive radius {r_\infty} that is a constant multiple of {r_0}, while (for {\delta_0} small enough) {\delta_n} converges double-exponentially to zero, so in particular {f_n(z)} converges uniformly to {\lambda z} on {B(0,r_\infty)}. Also, {\psi_n} is close enough to the identity, the compositions {\Psi_n := \psi_1 \circ \dots \circ \psi_n} are uniformly convergent on {B(0,r_\infty/2)} with {\Psi_n(0)=0} and {\Psi'_n(0)=1}. From this we have

\displaystyle  f( \Psi_n(z) ) = \Psi_n(f_n(z))

on {B(0,r_\infty/4)}, and on taking limits using Morera’s theorem we obtain a holomorphic function {\Psi} defined near {0} with {\Psi(0)=0}, {\Psi'(0)=1}, and

\displaystyle  f( \Psi(z) ) = \Psi(\lambda z),

obtaining the required linearisation.

Remark 5 The idea of using a Newton-type method to obtain error terms that decay double-exponentially, and can therefore absorb exponential losses in the iteration, also occurs in KAM theory and in Nash-Moser iteration, presumably due to Siegel’s influence on Moser. (I discuss Nash-Moser iteration in this note that I wrote back in 2006.)