The quadratic family f c ( z ) = z 2 + c ,   c ∈ ℂ introduced formally in Definition 3.1, occupies a central position in one-dimensional holomorphic dynamics. For a fixed parameter c , its dynamical-plane structure is encoded by the filled Julia set K c , the Julia set J c = ∂ K c , and the basin of infinity A c ( ∞ ) = ℂ \ K c . Its parameter-space organization is governed by the Mandelbrot set

The formal definition is given in Definition 3.3. The critical orbit 0 , f c ( 0 ) = c , f c ∘ 2 ( 0 ) , ⋯ therefore determines whether the parameter belongs to ℳ or to its complement.

In the classical theory developed by Douady and Hubbard, the exterior ℂ \ ℳ admits a normalized conformal uniformization

Under the uniformization (2), the associated parameter Green function is represented by

The representation (3), together with the classical properties of Φ ℳ , yields the harmonicity and strict positivity of G ℳ on ℂ \ ℳ , its continuous extension by zero on ℳ , and its normalized logarithmic growth

The boundary condition and the growth condition (4) enter the classical uniqueness characterization of G ℳ as the normalized exterior Green function with pole at infinity. The same conformal coordinate determines the classical parameter equipotentials and external rays. These facts are standard consequences of classical exterior uniformization and logarithmic potential theory. They are recalled here to establish the analytic framework and normalization required for the subsequent quantitative analysis; see [1]-[5]. The relationship between the dynamical-plane and parameter-plane normalizations is likewise classical. Let g c denote the dynamical Green function of f c . Its functional equation is

Applying (5) at z = 0 , together with f c ( 0 ) = c , gives

Combining (3) and (6) gives the complete normalization identity

The classical identity (7) provides the bridge between critical orbit escape and the parameter Böttcher potential. The representation of G ℳ by the parameter Böttcher map, its standard potential-theoretic properties, its uniqueness characterization, and the preceding dynamical–parameter identity are not claimed as new results. They are included to provide a self-contained structural setting and to establish the precise normalization used throughout the paper.

For the critical orbit

consider the finite escape-rate approximants

The quantitative escape-rate estimate governing the approximants in (9) was established in the author’s earlier work [6]. Accordingly, the estimate considered in isolation is not claimed as a new result of the present paper.

For completeness, the present manuscript gives a self-contained telescoping derivation directly from the critical orbit recurrence in (8). More importantly, the present manuscript identifies the parameter-plane conformal meaning of the estimate through the classical normalization identity (7). Through the identification (7), the normalized finite critical orbit quantities

become quantitative approximations to the parameter Böttcher potential. Once the critical orbit has entered a prescribed escape region, Theorem 8.1 gives an explicit estimate of the form

In (11), the constant C ( c , R ) is independent of n . Thus, through the conformal interpretation developed here, the previously established escape-rate estimate becomes an explicit approximation result for the exterior Green function of the Mandelbrot set. The pointwise estimate (11) admits a uniform compact-subset consequence. The present paper also obtains uniform exponential convergence on every compact subset of ℂ \ ℳ . More precisely, if K ⊂ ℂ \ ℳ is compact, then Corollary 8.4 provides common constants n K ≥ 1 and C K > 0 such that

The compact-uniform estimate (12) clarifies why a common finite iteration depth can approximate G ℳ reliably on parameter regions lying a positive distance from ∂ ℳ . It also identifies the boundary as the region in which no global common escape time should be expected.

The numerical methodology is formulated to preserve the distinction between fixed-index approximation and parameter-dependent escape detection. For a prescribed index N , the fixed-index field is

This quantity uses the same iteration index for every parameter. By contrast, the variable-time first-escape estimator is

where τ N max ( c ) is the first index, not exceeding the iteration budget N max , at which the critical orbit crosses the chosen escape radius. The estimators in (13) and (14) coincide only when N = τ N max ( c ) . They are therefore kept notationally and computationally distinct. For every parameter whose escape is detected, the first-escape estimator is accompanied by the explicit pointwise certificate

The certificate (15) applies only when escape is detected and the stated admissibility condition on the escape radius is satisfied. If escape is not detected within the prescribed iteration budget, the parameter is placed in a finite-time unresolved mask. This numerical outcome records only the failure to detect escape and does not certify membership in ℳ . The corresponding dark regions in the figures must therefore not be interpreted as rigorous numerical representations of the Mandelbrot set.

The convergence visualization similarly distinguishes the exact Green function from a certified numerical reference. If G ref ( c ) is a first-escape reference computed using a large reference radius and

then the triangle inequality gives

Here B n ( c ; R ) is the bound supplied by Theorem 8.1. The right-hand side of (17) is the certified comparison bound represented by the dotted curves in the convergence Figure 5.

The exterior conformal representation also gives the identity for every positive equipotential

The present paper proves that the family of equipotentials varies continuously in the Hausdorff metric over every compact positive level range. More precisely, Theorem 9.1 shows that if 0 < a < b , t n , t ∈ [ a , b ] , and t n → t , then

The convergence in (19) concerns the exact equipotentials of G ℳ . The numerical level curves displayed later are finite-resolution illustrations of the theorem and are not used in its proof.

The computed exterior-geometry figure also uses a truncated approximation to the Böttcher logarithm in order to visualize numerical external-angle fields and parameter-ray segments. These objects are included only as finite-resolution approximations. No global branch-certified error estimate for arg Φ ℳ is claimed, and the extracted numerical ray components are not identified with exact external rays.

The contributions developed in the present paper may therefore be summarized as follows:

(i) A self-contained telescoping derivation of the previously established quantitative escape-rate estimate;

(ii) An explicit interpretation of that estimate as an approximation to the parameter Böttcher potential log | Φ ℳ ( c ) | ;

(iii) Uniform exponential convergence of the normalized finite critical orbit approximants on compact subsets of ℂ \ ℳ ;

(iv) A precise distinction between fixed-index approximation and variable-time first-escape estimation, including a pointwise error certificate for detected escaping parameters;

(v) A certified comparison bound for the discrepancy between fixed-index approximants and the numerical first-escape reference;

(vi) Hausdorff continuity of the exact equipotentials over compact positive level ranges;

(vii) Finite-resolution numerical procedures whose notation and normalization are consistent with the analytic quantities being approximated.

The emphasis of the present article therefore differs from that of the author’s earlier work [6]. The earlier paper developed the effective quantitative behaviour of regularized critical orbit escape. The present paper places that estimate within the classical exterior conformal geometry of the Mandelbrot set, supplies the telescoping derivation for self-containment, derives the compact-uniform consequence, formulates the certified first-escape comparison, and establishes the positive-level equipotential-stability result. Accordingly, the contributions of the present article lie in the quantitative–conformal interpretation, its compact-uniform consequence, the certified formulation of the numerical approximations, and the Hausdorff stability analysis. The classical construction of the parameter Green function, its standard potential-theoretic properties, its uniqueness characterization, the dynamical–parameter normalization, and the escape-rate estimate considered independently are not claimed as new results.

The paper is organized as follows. Section 2 reviews the relevant literature and places the present work in its mathematical context. Section 3 introduces the definitions and normalizations used throughout the paper. Section 4 recalls the classical dynamical Böttcher coordinate, its holomorphic dependence near infinity, and the Douady–Hubbard parameter uniformization, summarized schematically in Figure 1. Section 5 records the classical exterior conformal representation of the parameter Green function and derives its standard potential-theoretic properties. Section 6 gives, for completeness, a direct maximum-principle proof of the classical uniqueness characterization. Section 7 recalls the exact conformal description of equipotentials and external rays. Section 8 gives the self-contained telescoping derivation, establishes the quantitative–conformal interpretation, derives uniform exponential convergence on compact exterior parameter sets, and introduces the numerical procedures used in Figures 2-5. Section 9 proves Hausdorff continuity of equipotentials over compact positive level ranges, illustrated in Figure 6. Section 10 explains the relation between the present conformal interpretation and the author’s earlier quantitative study. Section 11 records the limitations of the present analysis, and Section 12 formulates the proposed directions for future research. The paper concludes in Section 13.

The modern study of quadratic complex dynamics is rooted in the work of Douady and Hubbard, who developed the foundational theory of polynomial-like mappings, external uniformization, and the parameter-space description of the Mandelbrot set [2][3]. In their framework, the Mandelbrot set appears as the connectedness locus for the quadratic family, while its complement is naturally described by a Green function with pole at infinity. Earlier potential-theoretic foundations of polynomial iteration, including invariant sets and capacity, go back to Brolin’s work [7].

A standard and geometrically transparent account of these ideas is given by Milnor [5], who emphasizes the role of critical orbits, Julia sets, Böttcher coordinates, and external rays in one-dimensional holomorphic dynamics. Closely related presentations appear in Carleson and Gamelin [1] and in Beardon’s classical treatment of iteration theory [8]. For broader analytic and geometric background, including Teichmüller-theoretic aspects and the global organization of parameter space, one may also consult Hubbard [4].

From the standpoint of analysis, the Green function of ℂ \ ℳ belongs to classical logarithmic potential theory. Its basic properties—harmonicity on the complement, vanishing on the compact set, and logarithmic growth at infinity—fit naturally into the theory of harmonic and subharmonic functions, planar Green functions, Dirichlet problems, and conformal mapping. Useful background for these topics may be found in Ransford [9] and in Axler, Bourdon, and Ramey [10]. In parameter spaces of rational maps, potential-theoretic and Green-function ideas also underlie bifurcation currents and Lyapunov-exponent methods, as developed for example by DeMarco [11] and Dujardin–Favre [12]. For Green functions in varying polynomial families, see also Favre–Gauthier [13].

Recent developments show that the geometry of the Mandelbrot set and of related parameter spaces continues to evolve in several active directions. These include renormalization and local-connectivity phenomena associated with the Mandelbrot set and the MLC programme [14], self-similarity near Siegel parameters and pacman renormalization [15], positive-area phenomena for quadratic Julia sets [16], continuity and limiting behaviour for Julia sets and invariant external rays [17], richer parameter-space structures in cubic and parabolic settings [18], and Mandelbrot-type analogues in broader polynomial and transcendental families [19][20]. In this broader context, the present paper remains within the classical quadratic setting, but its emphasis lies at the intersection of exterior conformal uniformization, logarithmic potential theory, and the structural interpretation of escape-rate regularization.

The present article should be read together with [6]. In that work, the principal object was the regularized critical orbit escape-rate function

The function in (20) was studied as a quantitatively tractable observable of escape. The present article is a geometric companion: the earlier paper studies the effective quantitative behaviour of regularized escape, whereas the present paper identifies the structural analytic object underlying that theory and proves a bridge theorem between the two viewpoints.

Definition 3.1 (Quadratic Family).For c ∈ ℂ , define

Its n -fold iterate is denoted by f c ∘ n .

Definition 3.2 (Filled Julia set, Julia set, and Basin of Infinity [1], Chapter VI and[5]Chapter 8).For fixed c ∈ ℂ , the filled Julia set of f c is

the Julia set is

and the basin of infinity is

Definition 3.3 (Mandelbrot Set,[2][3]and[5], Chapter 14).The Mandelbrot set is

Equivalently, c ∈ ℳ if and only if the critical orbit of f c is bounded.

Definition 3.4 (Positive Logarithm).For r ≥ 0 , define

For r > 0 , the definition is equivalently written as

Definition 3.5 (Dynamical Green Function) ([1] Chapter VI and [5] Chapter 8).For each c ∈ ℂ , the dynamical Green function of f c ( z ) = z 2 + c is

Here log + is defined in (26). It is a classical theorem that the limit in (28) exists for every z ∈ ℂ .

Definition 3.6 (FiniteCritical OrbitEscapeRateApproximants).For c ∈ ℂ , let

For n ≥ 0 , define

The quantities defined in (30) are the finite critical orbit regularizations used later in Theorem 8.1 and in the numerical schemes of Section 8.1. The recurrence (29) is used in the telescoping argument of Theorem 8.1.

By the defining limit (28), evaluated at z = 0 , and the definition in (30), one has

Remark3.7 (Classical Properties ofg_c).It is classical that for every c ∈ ℂ , the function g c defined in (28) exists on all of ℂ , vanishes precisely on K c , is harmonic on A c ( ∞ ) , and satisfies the functional equation

See [1], Chapter VI, or [5], Chapter 8.

Definition 3.8 (Logarithmic Pole at Infinity).Let Ω ⊂ ℂ be an unbounded domain whose complement is compact. A real-valued function u on Ω has a logarithmic pole at infinity with the normalization used in this paper if

Equivalently, u ( c ) = log | c | + o ( 1 ) as c → ∞ . This is the normalization used for the Green function of ℂ \ ℳ .

Definition 3.9 (Parameter Green Function, Dynamical Normalization).For c ∈ ℂ \ ℳ , define initially

The dynamical normalization (34) will be identified in Proposition 5.3 with the exterior conformal formula introduced in Definition 5.1.

Remark3.10 (Normalization Bridge).The normalization in Definition 3.9 is naturally compatible with the parameter-plane Böttcher map introduced in Notation 3.11. It differs by a factor 2 from the critical-point normalization G ( c ) : = g c ( 0 ) used in [6]. Since f c ( 0 ) = c , the functional Equation (32) gives

The identity (35) shows that the two normalizations encode the same escape geometry, but they are adapted to different viewpoints: G is suited to critical orbit regularization, whereas G ℳ is suited to parameter-plane exterior uniformization.

Notation 3.11 (Parameter Böttcher Map).Throughout the paper, Φ ℳ denotes the conformal map Φ ℳ : ℂ \ ℳ → ℂ \ D ¯ normalized by

This map is recalled more formally in Theorem 4.8.

Definition 3.12 (Hausdorff Distance).Let A , B ⊂ ℂ be nonempty compact sets. Their Hausdorff distance is

The metric (37) is used in Theorem 9.1 to measure stability of equipotentials.

Notation 3.13 (Exterior Circles and Compact Annuli). For t > 0 , write

For 0 < a < b , write

Remark3.14 (Compact Subsets of the Exterior). If K ⊂ ℂ \ ℳ is a nonempty compact set, then K lies at a positive Euclidean distance from ℳ . Equivalently,

This positive separation does not by itself provide a common escape time for all parameters in K . In the proof of Corollary 8.4, continuity of the critical orbit polynomials from Notation 4.4, together with a finite-subcover argument, is used to construct the common entry index n K defined in (129).

We recall the classical Böttcher coordinate near infinity and the corresponding parameter-space uniformization. These results provide the conformal input used in Definition 5.1 and Proposition 5.3.

Theorem 4.1 (Böttcher CoordinatenearInfinity)Fix c ∈ ℂ . There exist a neighbourhood U c of infinity, a number R c > 1 , and a unique conformal map ϕ c : U c → { w ∈ ℂ : | w | > R c } such that

whenever z , f c ( z ) ∈ U c , and ϕ c ( z ) z → 1 as z → ∞ . Moreover,

Proof. This is the classical Böttcher theorem applied to the super-attracting fixed point at infinity of f c ( z ) = z 2 + c ; see [8], Chapter 3, [1], Chapter VI, and [5], Chapter 8. For z sufficiently large,

The factorization (43) shows that the correction term c / z 2 is small near infinity. The standard Böttcher iteration construction produces a holomorphic function ϕ c , defined in an exterior neighbourhood of infinity, satisfying (41) and ϕ c ( z ) = z + O ( 1 ) as ( z → ∞ ). The stronger normalization ϕ c ( z ) z → 1 determines ϕ c uniquely.

For completeness, the construction can be written locally in the form

In the representation (44), the neighbourhood of infinity is chosen so that

Here, Log denotes the branch satisfying Log   1 = 0 . Condition (45) ensures that the logarithmic terms are evaluated in a common neighbourhood of 1. The series in (44) converges locally uniformly near infinity and therefore defines a holomorphic function there. Finally, by the defining limit (28), the dynamical Green function is obtained from the normalized asymptotic growth of f c ∘ n ( z ) . For z ∈ U c , the orbit tends to infinity, and the Böttcher functional equation gives

Taking absolute values and logarithms in (46) gives

Since ϕ c ( w ) w → 1 as ( w → ∞ ), the difference between log | ϕ c ( f c ∘ n ( z ) ) | and log | f c ∘ n ( z ) | tends to zero. Combining the defining limit (28) with (47) therefore proves the identity (42). □

The next result explains why the critical value c = f c ( 0 ) is the correct parameter-plane point at which to evaluate the dynamical Böttcher coordinate.

Proposition 4.2.If c ∈ ℂ \ ℳ , then both 0 and c = f c ( 0 ) belong to the basin of infinity A c ( ∞ ) from Definition3.2.Moreover, thecritical orbiteventually enters every sufficiently small neighbourhood of infinity. In particular, there exists an integer N ≥ 1 such that f c ∘ N ( 0 ) ∈ U c , where U c is the exterior neighbourhood on which the local Böttcher coordinate in Theorem 4.1 is defined.

Proof. Assume that c ∉ ℳ . By Definition 3.3, the critical orbit 0 , f c ( 0 ) = c , f c ∘ 2 ( 0 ) , ⋯ is unbounded.

For a polynomial, an unbounded forward orbit necessarily tends to infinity. Indeed, the reverse triangle inequality gives

The estimate (48) shows that, once an orbit enters a sufficiently large exterior region, its modulus continues to increase. Consequently, f c ∘ n ( 0 ) → ∞ as n → ∞ . It follows that the critical point 0 does not belong to the filled Julia set K c . Hence

Since c = f c ( 0 ) , the forward orbit of c is the tail of the critical orbit:

Therefore,

The limit (51) proves that the critical value c = f c ( 0 ) also belongs to A c ( ∞ ) . Finally, because f c ∘ n ( 0 ) → ∞ , the critical orbit eventually enters the exterior neighbourhood U c supplied by Theorem 4.1. Hence there exists N ≥ 1 such that f c ∘ N ( 0 ) ∈ U c . □

Remark4.3. Proposition 4.2 shows that the critical value c = f c ( 0 ) belongs to the basin of infinity and that its forward orbit eventually enters the exterior neighbourhood U c . It does not, however, imply that the critical value c itself lies in the particular neighbourhood U c on which the local coordinate of Theorem 4.1 is initially defined. Moreover, Theorem 4.1 is a fixed-parameter statement and does not by itself establish holomorphic dependence of the Böttcher coordinate on c . These two issues are addressed by the family construction near infinity and by the classical Douady–Hubbard parameter uniformization stated below.

Notation 4.4 (Critical orbitPolynomials). For c ∈ ℂ , define

Equivalently, the recurrence (52) gives

For every n ≥ 1 , the function c ↦ z n ( c ) is a polynomial in c . In particular,

If c ∈ ℂ \ ℳ , then z n ( c ) → ∞ as n → ∞ by Proposition 4.2.

Definition 4.5 (Local Parameter Böttcher Branch). Let V ⊂ ℂ \ ℳ be a simply connected parameter domain. Suppose that N ≥ 1 and R > 1 are such that | z N ( c ) | > R for c ∈ V , where z N ( c ) = f c ∘ N ( 0 ) .

Let ϕ : V × { z ∈ ℂ : | z | > R } → ℂ be the jointly holomorphic Böttcher coordinate supplied by Proposition 4.7. Alocal parameter Böttcher branch on V is a holomorphic function Φ V : V → ℂ \ D ¯ , satisfying

Equation (55) is understood with the root branch selected compatibly with the normalization (36).

Notation 4.6 (NormalizedCritical orbitRoot Approximants). On a simply connected parameter domain V ⊂ ℂ \ ℳ where compatible holomorphic roots have been selected, write

The root expressions in (56) are understood locally through compatible holomorphic branches. This notation is local until the global parameter Böttcher map has been established by Theorem 4.8.

The following standard family version of the local Böttcher construction is included with proof in the notation required below; compare [1], Chapter VI and [5], Chapter 8.

Proposition 4.7 (Joint Holomorphic Dependence near Infinity).Let V ⊂ ℂ be a relatively compact parameter domain. Then there exists R > 1 and a jointly holomorphic function ϕ : V × { z ∈ ℂ : | z | > R } → ℂ such that, for every c ∈ V , the map z ↦ ϕ c ( z ) : = ϕ ( c , z ) is the normalized Böttcher coordinate of f c on { | z | > R } . In particular, ϕ ( c , f c ( z ) ) = ϕ ( c , z ) 2 whenever | z | > R , and ϕ ( c , z ) z → 1 as z → ∞ , locally uniformly with respect to c ∈ V .

Proof. Since V is relatively compact, its closure V ¯ is compact. Hence there exists B > 0 such that | c | ≤ B for every c ∈ V ¯ , and therefore for every c ∈ V . Choose R > 1 sufficiently large that R 2 − B > R and B R 2 < 1 2 . Then, for c ∈ V and | z | > R ,

Estimate (57) shows that the common exterior region { z ∈ ℂ : | z | > R } is forward invariant under f c for every c ∈ V . Consequently, for c ∈ V and | z | > R , define z j ( c , z ) : = f c ∘ j ( z ) . By (57),

For c ∈ V and | z | > R , define

where Log denotes the branch of the logarithm defined near 1 and normalized by Log   1 = 0 .

Since | c | ≤ B , the orbit estimate (58) and the choice B / R 2 < 1 / 2 give

Condition (60) ensures that every logarithmic term in (59) is evaluated in a common neighbourhood of 1. Moreover, the forward-invariance estimate (57) implies rapid growth of the iterates in the common exterior region. Consequently, the series in (59) converges locally uniformly on V × { z ∈ ℂ : | z | > R } . Each summand is jointly holomorphic in ( c , z ) , and therefore H ( c , z ) is jointly holomorphic.

Define

Since H is jointly holomorphic, the definition (61) shows that ϕ is jointly holomorphic.

Separating the first term of the series (59) and using the recurrence z j + 1 ( c , z ) = z j ( c , f c ( z ) ) gives

By (61), ϕ ( c , z ) z = exp ( H ( c , z ) ) . Since the series (59) tends locally uniformly to zero as z → ∞ , one obtains

locally uniformly with respect to c ∈ V .

For each fixed c ∈ V , the functional Equation (62) and the normalization (63), together with uniqueness in Theorem 4.1, show that ϕ ( c , ⋅ ) is the normalized local Böttcher coordinate of f c . □

Theorem 4.8 (Douady--Hubbard Parameter Uniformization; classical, see[2][3],[4]Chapter 13, and[5]Chapter 14).There exists a unique conformal isomorphism Φ ℳ : ℂ \ ℳ → ℂ \ D ¯ satisfying the normalization Φ ℳ ( c ) c → 1 as c → ∞ . Let c 0 ∈ ℂ \ ℳ . Then there exist a simply connected neighbourhood V ⊂ ℂ \ ℳ of c 0 , an integer N ≥ 1 , and a number R > 1 such that | z N ( c ) | > R , c ∈ V where z N ( c ) = f c ∘ N ( 0 ) is thecritical orbitpolynomial defined in Notation 4.4. Moreover, if ϕ : V × { z ∈ ℂ : | z | > R } → ℂ is the jointly holomorphic exterior Böttcher coordinate supplied by Proposition 4.7, then

The relation(64) identifies the restriction Φ ℳ | V with the local parameter Böttcher branch of Definition4.5that is compatible with the normalization at infinity. In particular, Φ ℳ is holomorphic and nonvanishing on ℂ \ ℳ , and

Proof.The existence of the conformal isomorphism Φ ℳ : ℂ \ ℳ → ℂ \ D ¯ , with the normalization (36), is the classical Douady–Hubbard external uniformization theorem for the quadratic connectedness locus.

We verify the stated local representation in the notation used in this paper. Fix c 0 ∈ ℂ \ ℳ . By Proposition 4.2, and the critical orbit representation (53),

Choose a relatively compact parameter neighbourhood W ⊂ ℂ \ ℳ of c 0 . By Proposition 4.7, there exist R > 1 and a jointly holomorphic function ϕ : W × { z ∈ ℂ : | z | > R } → ℂ such that, for each c ∈ W , the function ϕ c ( z ) : = ϕ ( c , z ) is the normalized dynamical Böttcher coordinate of f c on the common exterior region { | z | > R } . By Proposition 4.7, the map ϕ satisfies the functional Equation (62) and the normalization (63). Since z n ( c 0 ) → ∞ , there exists N ≥ 1 such that | z N ( c 0 ) | > R . By Notation 4.4, z N ( c ) is a polynomial in c , and hence is holomorphic and continuous. After shrinking W to a sufficiently small simply connected neighbourhood V of c 0 , we therefore have

Define

By (67), the composition in (68) is holomorphic on V . Moreover, because each ϕ c is a Böttcher coordinate with values in ℂ \ D ¯ , the function F N is nonvanishing on V .

Because the function F N defined in (68) is nonvanishing on the simply connected domain V , it admits a holomorphic logarithm and, consequently, a holomorphic 2 N − 1 -st root. By Definition 4.5, the root branch compatible with the normalization at infinity defines Φ V and satisfies (55).

The definition is independent of the chosen entry time. Indeed, if m ≥ N , then

Repeated application of the functional Equation (62) along the orbit identity (69) gives

Therefore, for the compatible local root branches,

Identity (71) proves that the normalized pullback is unchanged when a later entry time is used.

The classical Douady–Hubbard theorem guarantees that the normalized local branches obtained in this manner are compatible on overlaps and patch together to form a single-valued holomorphic map on ℂ \ ℳ . This global map is precisely Φ ℳ . Consequently, the globally defined parameter Böttcher map satisfies the local representation (64) on V .

The same classical theorem establishes that Φ ℳ is one-to-one and maps ℂ \ ℳ onto ℂ \ D ¯ . Therefore Φ ℳ is a conformal isomorphism. Since Φ ℳ maps ℂ \ ℳ onto ℂ \ D ¯ , the modulus estimate (65) follows.

It remains to prove uniqueness under the stated normalization. Suppose that Ψ : ℂ \ ℳ → ℂ \ D ¯ is another conformal isomorphism satisfying the normalization Ψ ( c ) / c → 1 as c → ∞ , analogous to (36). Then A : = Ψ ∘ Φ ℳ − 1 is an automorphism of ℂ \ D ¯ that fixes infinity. Every such automorphism is a rotation, A ( w ) = e i θ w . The normalization at infinity gives e i θ = 1 . Hence A is the identity map, and therefore Ψ = Φ ℳ . □

Corollary 4.9 (Critical OrbitRepresentation).Let V ⊂ ℂ \ ℳ be a simply connected parameter neighbourhood on which the compatible holomorphic root branches are chosen according to Definition4.5. Then

locally uniformly for c ∈ V . Equivalently,

where the roots denote the local holomorphic branches compatible with the normalization (36).

Proof. Let N and R be as in Theorem 4.8. For n ≥ N , the orbit point z n ( c ) lies in the common exterior Böttcher region. Since

locally uniformly with respect to c , applying (74) along the escaping critical orbit gives

locally uniformly on V . By the functional equation and Theorem 4.8,

Consequently,

By (75), the second factor in (77) converges locally uniformly to 1 . Hence the convergence asserted in (72) follows. □

Remark4.10. The root-limit formula in Corollary 4.9 is not used as a branch-free global definition of Φ ℳ . For arbitrary root choices, the expressions z n ( c ) 1 / 2 n − 1 are multivalued and do not by themselves establish either single-valuedness or holomorphic dependence on c . The global single-valued holomorphic map Φ ℳ is supplied by the classical Douady–Hubbard uniformization theorem. The corollary expresses this already-established map locally through the uniquely compatible holomorphic root branches determined by its normalization at infinity.

Remark4.11. The notation Φ ℳ ( c ) = ϕ c ( c ) is commonly used as shorthand for the classical identification between the parameter Böttcher map and the normalized dynamical exterior coordinate evaluated along the escaping critical orbit. It must not be interpreted as a direct evaluation of the local map supplied by Theorem 4.1, because the critical value c need not belong to the particular neighbourhood U c on which that local coordinate is initially defined. More precisely, if N ≥ 1 is chosen so that z N ( c ) = f c ∘ N ( 0 ) lies in the exterior Böttcher region, then

with the root selected by the normalization at infinity. Consequently,

Taking absolute values and logarithms in (78) gives (79). By the local Green-function identity (42), evaluated at z N ( c ) ,

Combining (79), (80), and the dynamical functional Equation (32) gives

The geometric content of Theorem 4.8 is summarized in Figure 1. The exterior coordinate sends equipotentials to concentric circles and parameter external rays to radial half-lines.

Figure 1. Schematic representation of the classical Douady–Hubbard exterior uniformization. By Theorem 4.8, the parameter Böttcher map Φ ℳ conformally identifies ℂ \ ℳ with ℂ \ D ¯ . Under this map, the equipotential G ℳ ( c ) = t , defined using Definition 5.1, corresponds to the circle | w | = e t , as stated in Proposition 7.2. The parameter external ray of angle θ , introduced in Definition 7.1, corresponds to the radial half-line w = r e 2 π i θ , r > 1 . The parameter-plane ray is schematic and is included to illustrate the conformal correspondence rather than to represent a numerically computed inverse Böttcher ray.

We now construct the Green function of ℂ \ ℳ explicitly from the parameter Böttcher map of Notation 3.11 and Theorem 4.8. The results of this section are classical consequences of the Douady-Hubbard parameter uniformization and standard logarithmic potential theory; see [1], Chapter VI, [2][3], [5], Chapter 14, and [9]. They are included with proofs to fix the precise normalization used in the quantitative analysis.

Definition 5.1 (Exterior Green function). Let Φ ℳ : ℂ \ ℳ → ℂ \ D ¯ be the normalized parameter Böttcher map from Theorem 4.8. Define

Equation (82) is the exact conformal definition used throughout the remainder of the paper.

Remark5.2. Definition 5.1 is the classical conformal representation of the exterior Green function:

The composition formula (83) shows that G ℳ is defined exactly on ℂ \ ℳ , rather than only through its asymptotic behaviour at infinity.

We next identify this parameter-plane function with the dynamical Green function evaluated at the critical value.

Proposition 5.3.For every c ∈ ℂ \ ℳ , G ℳ ( c ) = g c ( c ) .

Proof. Fix c ∈ ℂ \ ℳ , and write z n ( c ) = f c ∘ n ( 0 ) . Choose N ≥ 1 so that z N ( c ) belongs to an exterior neighbourhood on which the dynamical Böttcher coordinate ϕ c is defined. By Theorem 4.8, the local representation (64),

Taking absolute values and logarithms in (84) gives

By Theorem 4.1, and the local identity (42),

Since the critical value is c = f c ( 0 ) , one has

Repeated application of the functional Equation (32) along (87) gives

Combining (85), (86), and (88) yields

Using Definition 5.1 and (82), we obtain G ℳ ( c ) = g c ( c ) . □

Corollary 5.4.For every c ∈ ℂ \ ℳ , G ℳ ( c ) = 2 g c ( 0 ) . Consequently,

Proof. Since c = f c ( 0 ) , the functional Equation (32) from Remark 3.7 gives

Using Proposition 5.3, we conclude that G ℳ ( c ) = 2 g c ( 0 ) . The limit formula (90) then follows from Definition 3.5 and (28), evaluated at z = 0 . □

Remark5.5. Corollary 5.4 is the precise normalization bridge between the present paper and [6]. The earlier quantitative paper studied the critical-point regularization g c ( 0 ) , whereas the present paper is organized around the parameter-plane Green function G ℳ ( c ) = g c ( c ) . The factor 2 reflects the single-step passage from the critical point 0 to the critical value c = f c ( 0 ) .

Proposition 5.6.The function G ℳ is harmonic on ℂ \ ℳ .

Proof. By Theorem 4.8, the map Φ ℳ is holomorphic and maps ℂ \ ℳ into ℂ \ D ¯ . In particular, Φ ℳ is nowhere zero on ℂ \ ℳ . Hence, near every point of ℂ \ ℳ , one may choose a holomorphic branch of Log   Φ ℳ . On such a neighbourhood,

Identity (92) shows that log | Φ ℳ | is harmonic locally. Since harmonicity is a local property and G ℳ = log | Φ ℳ | by Definition 5.1 and (82), G ℳ is harmonic on ℂ \ ℳ . □

Proposition 5.7.For every c ∈ ℂ \ ℳ , one has G ℳ ( c ) > 0.

Proof. If c ∈ ℂ \ ℳ , then Theorem 4.8 gives Φ ℳ ( c ) ∈ ℂ \ D ¯ . The modulus estimate (65) gives | Φ ℳ ( c ) | > 1 . Taking logarithms and using Definition 5.1 together with (82) gives G ℳ ( c ) > 0 . □

Theorem 5.8 (Logarithmic growth at infinity).As c → ∞ ,

Proof. By Theorem 4.8 and the normalization (36), Φ ℳ ( c ) c → 1 as c → ∞ . Taking moduli gives

Since the real logarithm is continuous at 1, it follows that

Using (82) in (95) gives (93). □

Corollary 5.9.The function G ℳ has a logarithmic pole at infinity in the sense of Definition 3.8.

Proof. By Theorem 5.8 and (93), G ℳ ( c ) = log | c | + o ( 1 ) as c → ∞ , which is exactly the normalization in Definition 3.8 and (33). □

Proposition 5.10.The function G ℳ extends continuously to ℂ by setting G ℳ ( c ) = 0 for all c ∈ ℳ . In particular, G ℳ = 0 on ∂ ℳ .

Proof. By Corollary 5.4, G ℳ ( c ) = 2 g c ( 0 ) for all c ∈ ℂ \ ℳ . The escape-rate function ( c , z ) ↦ g c ( z ) is continuous for the quadratic family; see, for example, [1], Chapter VI, or [5], Chapter 8. In particular, the function c ↦ g c ( 0 ) is continuous on ℂ . Moreover,

Let c j ∈ ℂ \ ℳ and suppose that c j → c * ∈ ∂ ℳ . By continuity,

The final equality in (97) follows from (96), since c * ∈ ∂ ℳ ⊂ ℳ . Define

The convergence (97) shows that the extension (98) is continuous on all of ℂ . In particular, G ℳ = 0 on ∂ ℳ . □

Theorem 5.11 (Structural Characterization of G ℳ ).The function G ℳ constructed in Definition5.1 has the following properties:

(i) G ℳ is harmonic on ℂ \ ℳ ;

(ii) G ℳ ( c ) > 0 for every c ∈ ℂ \ ℳ ;

(iii) G ℳ extends continuously to ℂ and vanishes on ℳ ;

(iv) G ℳ ( c ) − log | c | → 0 as c → ∞ .

Proof. Part (i) is Proposition 5.6. Part (ii) is Proposition 5.7. Part (iii) is Proposition 5.10. Part (iv) is Theorem 5.8, equivalently the logarithmic-pole normalization of Definition 3.8. □

We now prove the classical uniqueness characterization of G ℳ as the Green function of ℂ \ ℳ with pole at infinity in the normalization of Definition 3.8; compare [9] and [10]. The proof is included for completeness.

Theorem 6.1 (Uniqueness of the Normalized Exterior Green Function)Let u : ℂ \ ℳ → ℝ be a function satisfying:

(a) u is harmonic on ℂ \ ℳ ;

(b) u extends continuously to ∂ ℳ and u = 0 on ∂ ℳ ;

(c) u has a logarithmic pole at infinity with the normalization of Definition 3.8, namely the condition (33).

Then u ( c ) = G ℳ ( c ) for every c ∈ ℂ \ ℳ .

Proof. Define

By the harmonicity of u and Proposition 5.6, the function h defined in (99) is harmonic on ℂ \ ℳ .

By the boundary condition on u and Proposition 5.10, the function h extends continuously to ∂ ℳ and satisfies h = 0 on ∂ ℳ . Furthermore, by the growth assumption on u , Theorem 5.8, and (93),

Fix c 0 ∈ ℂ \ ℳ . Since ℳ is compact, choose R 0 > 0 such that ℳ ⊂ { | c | < R 0 } . For R > R 0 with | c 0 | < R , set

The function h is harmonic on Ω R , continuous up to the boundary pieces ∂ ℳ and { | c | = R } , and the boundary of Ω R is contained in

Set

By (100), ε ( R ) → 0 as R → ∞ . The maximum principle applied to both h and − h on Ω R gives

Applying (104) at c = c 0 gives | h ( c 0 ) | ≤ ε ( R ) . Letting R → ∞ and using ε ( R ) → 0 , which follows from (100), gives h ( c 0 ) = 0 . Since c 0 ∈ ℂ \ ℳ was arbitrary, it follows that h ≡ 0 on ℂ \ ℳ . Therefore, u = G ℳ on ℂ \ ℳ . □