Parabolic implosion in dimension 2

Parabolic implosion is the study of the bifurcation phenomena which occur when a multiple (i.e., parabolic) fixed point is perturbed and splits into several fixed points or periodic cycles. It was first developed by Lavaurs in his PhD thesis ([Lav89]). A first consequence of this theory is a precise description of the discontinuity (enrichment) of Julia sets with respect to the parameter, in the presence of a non-persistent parabolic cycle. This was used by Shishikura ([Shi98]) to prove that the boundary of the Mandelbrot set has Hausdorff dimension 2. A refinement of parabolic implosion (near-parabolic renormalization), developed by Inou and Shishikura ([IS06]), has led to remarkable results, such as the construction of quadratic Julia sets with positive area by Buff and Chéritat ([BC12]) or progress towards the hyperbolicity conjecture ([CS15]).

More recently, the theory of parabolic implosion has started to develop and to find successful applications in higher dimension. In [BSU17], Bedford, Smillie and Ueda develop a parabolic implosion theory in the setting of semi-parabolic diffeomorphisms in dimension 2, i.e., in the case of a fixed point with one attracting direction and the other one with multiplier equal to 1. In particular, in the important case of a dissipative Hénon map, they were able to deduce the discontinuity of several dynamically defined sets (including the forward Julia set $J^{+}$ and the closure $J^{*}$ of the saddle periodic points) with respect to the parameter. Building on their result, Bianchi and the first named author proved in [AB25] the existence of perturbations of such Hénon maps whose forward Julia set $J^{+}$ has large Hausdorff dimension.

In [DL15], Dujardin and Lyubich adapted the results of [BSU17] to construct homoclinic tangencies for perturbations of dissipative Hénon maps with a semi-parabolic periodic cycle, with applications to bifurcation theory. In [ABD⁺16], the authors used parabolic implosion techniques to construct the first examples of polynomial maps (in dimension 2) with a wandering Fatou component (see also [ABTP23], [ABT26]). The dynamical systems under consideration are polynomial skew-products, hence the techniques employed can be seen as a non-autonomous version of one-dimensional parabolic implosion. Finally, in [Bia19b], Bianchi obtained results analogous to those of [BSU17] but for the more difficult case of maps with a fully parabolic fixed point, i.e., in the case where the differential at the fixed point is the identity. The purpose of this article is to extend the results from [Bia19b].

Let us now provide a quick overview of classical parabolic implosion in dimension one, in the simplest case of a parabolic fixed point with just one attracting and one repelling petal. Let $f:\mathbb{C}\to\mathbb{C}$ be a holomorphic map of the form $f(z)=z+z^{2}+\mathcal{O}(z^{3})$. By the classical Leau-Fatou theorem, there exists $r>0$ and univalent maps $\phi^{\iota}:\mathbb{D}(-r,r)\to\mathbb{C}$ and $\phi^{o}:\mathbb{D}(r,r)\to\mathbb{C}$ such that

1. (1)

   $P^{\iota}:=\mathbb{D}(-r,r)$ is forward invariant under $f$, and $P^{o}:=\mathbb{D}(r,r)$ is invariant under the branch of $f^{-1}$ fixing the origin;

2. (2)

   $\phi^{\iota}\circ f=\phi^{\iota}+1$, and $\phi^{o}\circ f=\phi^{o}+1$.

The domains $P^{\iota}$ and $P^{o}$ are respectively called incoming and outgoing petals and the maps $\phi^{\iota}$ and $\phi^{o}$ are called incoming and outgoing Fatou coordinates. Let

be the parabolic basin. Then $P^{\iota}\subset\mathcal{B}$, and moreover $\mathcal{B}=\bigcup_{n\geq 0}f^{-n}(P^{\iota})$. The incoming Fatou coordinate extends to a holomorphic map $\phi^{\iota}:\mathcal{B}\to\mathbb{C}$ and the inverse $(\phi^{o})^{-1}$ of the outgoing Fatou coordinate extends to a holomorphic map $\psi^{o}:\mathbb{C}\to\mathbb{C}$, called the outgoing Fatou parametrization. We refer the reader to [Mil11] for details. In particular, for any $\sigma\in\mathbb{C}$, the change of coordinate $\mathcal{L}_{\sigma}:=\psi^{o}\circ T_{\sigma}\circ\phi^{\iota}$ is well-defined on $\mathcal{B}$, where $T_{\sigma}(z):=z+\sigma$ is the translation of vector $\sigma$. It is called the Lavaurs map of phase $\sigma$.¹¹1Fatou coordinates are in fact only unique up to addition by a constant; thus $\mathcal{L}_{\sigma}$ may be thought of as $(\phi^{o})^{-1}\circ\phi^{\iota}$, for a different choice of $\phi^{\iota}$. In practice, we will work with some unique normalizations of Fatou coordinates based on their formal expansion at $0$.

Consider now the family of perturbations $f_{\varepsilon}(z)=f(z)+\varepsilon^{2}$, $\varepsilon\in\mathbb{C}$. For $\varepsilon$ small but non-zero, the double fixed point at the origin for $f$ splits into 2 simple fixed points for $f_{\varepsilon}$, of the form $z^{\pm}(\varepsilon)=\pm i\varepsilon+\mathcal{O}(\varepsilon^{2})$. If we take say $\varepsilon>0$, then Lavaurs proved that orbits under $f_{\varepsilon}$ starting from a point in $P^{\iota}$ will approach the origin, then cross the "gate" given by the vertical segment $[z^{-}(\varepsilon),z^{+}(\varepsilon)]$ between $z^{-}(\varepsilon)$ and $z^{+}(\varepsilon)$ and then move away from the origin inside $P^{o}$. Since $f_{\varepsilon}$ is close to the identity near $0$, as $\varepsilon\to 0$ it takes more and more iterations to do this. At the limit, we obtain in this way a "transit map" from $P^{\iota}$ to $P^{o}$, which is useful for studying the dynamics of $f_{\varepsilon}$. It turns out that this transit map is exactly the Lavaurs map defined above. More precisely:

We now move on to the setting of complex dimension 2. Let $f:U{\subset\mathbb{C}^{2}}\to\mathbb{C}^{2}$ be a holomorphic map on a neighborhood $U$ of the origin, with a power series expansion of the form

where the $P_{j}$ are homogeneous degree $j$ polynomial maps from $\mathbb{C}^{2}$ to $\mathbb{C}^{2}$, and $P_{2}\not\equiv 0$. Such a map is called tangent to the identity at order 2. Following Hakim [Hak98] and Écalle [É85], we say that $v\in\mathbb{C}^{2}\setminus\{(0,0)\}$ is a characteristic direction for $f$ if there exists $\lambda\in\mathbb{C}$ so that $P_{k}(v)=\lambda v$. If $\lambda\neq 0$ then $v$ is said to be non-degenerate. We shall denote by $v\mapsto[v]$ the canonical projection of $\mathbb{C}^{2}\backslash\{(0,0)\}$ onto $\mathbb{P}^{1}$. The director of a non-degenerate characteristic direction $v$ is the eigenvalue of the linear operator

If the real part of the director of a non-degenerate characteristic direction $v$ is strictly positive, Hakim proved in [Hak97] that for any $C>0$ there exist incoming and outgoing petals $P_{C}^{\iota}$ and $P_{C}^{o}$ and incoming/outgoing Fatou coordinates $\Phi^{\iota/o}:P_{C}^{\iota/o}\to\mathbb{C}^{2}$ conjugating $f$ to a translation of vector $(1,0)$ (see Propositions 3.2 and 3.3 for details). Let

denote the parabolic basin associated to $v$. Similarly to the one-dimensional case, we have that $\mathcal{B}_{U,v}=\bigcup_{C>0}\bigcup_{n\geq 0}f^{-n}(P_{C}^{\iota})$ (see [LHR25]).

We say that a formal non-singular curve $\mathcal{C}$ is invariant for $f$ if given a parametrization $\gamma(t)$ of $\mathcal{C}$ (i.e., $\gamma(t)\in\mathbb{C}[[t]]^{2}$ with $\gamma(0)=(0,0)$ and $\gamma^{\prime}(0)\neq(0,0)$) there exists $h\in t\mathbb{C}[[t]]$ with $h^{\prime}(0)\neq 0$ such that

The tangent of $\mathcal{C}$ is, by definition, $\mathbb{C}\cdot\gamma^{\prime}(0)$.

From now on, we will assume that $f:U{\subset\mathbb{C}^{2}}\to\mathbb{C}^{2}$ is a holomorphic map defined on a neighborhood $U$ of the origin and which satisfies the following assumption:

1. $(H_{1})$

   The map $f$ is tangent to the identity at order $2$ and has a non-degenerate characteristic direction $v$, with a formal non-singular invariant curve $\mathcal{C}$ tangent to $v$, and a director $\alpha$ such that $\mathrm{Re}\,\alpha>2$.

It is worth mentioning that the existence of a non-singular formal invariant curve tangent to a non-degenerate characteristic direction $v$ is a generic hypothesis. It is equivalent to the existence of an analytic curve $C$ tangent to $v$ which is preserved up to order $k=\lfloor\mathrm{Re}\,\alpha\rfloor+2$ in the following sense: if $\gamma(t)$ is a parametrization of $C$ then there exists $h\in t\mathbb{C}\{t\}$ such that

and this condition is always satisfied when $\alpha\not\in\mathbb{N}$ (see [[Hak98], Section 3]).

We now state a first, non-technical version of our main result:

Since the maps $f$ and $f_{\varepsilon}$ are only defined on $U$, in general the iterates $f_{\varepsilon_{n}}^{n}$ may not be well-defined; however Theorem 1 implies that there is a constant integer $N$ such that $f_{\varepsilon_{n}}^{n-N}$ is well-defined on $K$ for all $n$ large enough. Moreover, since the map

satisfies

it is not difficult to see that we also have $f_{\varepsilon_{n}}^{n-N^{\prime}}\to\mathcal{L}_{\sigma-N^{\prime},q}$ for any $N^{\prime}\in\mathbb{Z}$ such that $\mathcal{L}_{\sigma-N^{\prime},q}(K)\subset U$. In particular, in the case where the maps $f_{\varepsilon}$ are global endomorphisms of a complex manifold we may simply take $N=0$.

We can be more precise on the requirements for the family of perturbations $(f_{\varepsilon})_{\varepsilon\in\mathbb{D}}$, and interpret the constants $\sigma$ and $q$ in terms of $f_{\varepsilon}$. We will make the following assumptions:

1. $(H_{2})$

   The family $(f_{\varepsilon})_{\varepsilon\in\mathbb{D}}$ preserves the formal curve $\mathcal{C}$ to order $m+1$, where $m:=\lfloor\mathrm{Re}\,\alpha\rfloor+1$, in the following sense: if $\gamma(t)$ is a parametrization of $\mathcal{C}$, there exists $h_{\varepsilon}\in\mathcal{O}(\mathbb{D})[[t]]$ such that f_ε∘γ- γ∘h_ε∈⟨t^m+2, εt^m+1, …, ε^m+2 ⟩.

2. $(H_{3})$

   For all $\varepsilon$ in a neighborhood of $0$, $f_{\varepsilon}$ has exactly 4 fixed points $z_{i}(\varepsilon)$ near the origin ($1\leq i\leq 4$) and counted with multiplicity, which depend holomorphically on $\varepsilon$ and such that $z_{i}^{\prime}(0):=\frac{d}{d\varepsilon}_{|\varepsilon=0}z_{i}(\varepsilon)$ are non-zero and pairwise distinct.

Note that the existence of 4 fixed points depending holomorphically on $\varepsilon$ is always satisfied up to passing to a branched cover in parameter space (although doing so affects the derivatives $z_{i}^{\prime}(0)$).

With a slight abuse of notation we will say that $(f_{\varepsilon})_{\varepsilon\in\mathbb{D}}$ satisfies $(H_{1})-(H_{3})$ when $f_{0}$ satisfies $(H_{1})$ and $(f_{\varepsilon})_{\varepsilon\in\mathbb{D}^{*}}$ satisfies $(H_{2})$ and $(H_{3})$.

We can interpret the constants $\sigma$ and $q$ in terms of the multipliers of the fixed points of $f_{\varepsilon}$. To make this precise, we first need the following Proposition, whose proof is deferred until the next section:

We can now state a second version of our results, which is more explicit on the requirements on the family of perturbations $(f_{\varepsilon})_{\varepsilon\in\mathbb{D}}$. To do so, we will use the following convention: let $z_{1}(\varepsilon)$ and $z_{2}(\varepsilon)$ be the fixed points from Proposition 1 and let $\lambda_{1}(\varepsilon),\mu_{1}(\varepsilon)$ and $\lambda_{2}(\varepsilon),\mu_{2}(\varepsilon)$ be their eigenvalues. If $\lambda_{i}^{\prime}(0)\neq\mu_{i}^{\prime}(0)$, we denote by $\rho_{T}^{i}(\varepsilon)$ the eigenvalue whose eigenspace tends to $\mathbb{C}v$ as $\varepsilon\to 0$ and by $\rho_{N}^{i}(\varepsilon)$ the other one; if $\lambda_{i}^{\prime}(0)=\mu_{i}^{\prime}(0)$, we assign the names $\rho_{T}^{i}(\varepsilon)$ and $\rho_{N}^{i}(\varepsilon)$ indifferently to the two eigenvalues $\lambda_{i}(\varepsilon)$ and $\mu_{i}(\varepsilon)$.

We finally state a third version of our main result, explicitly expressed in coordinates:

Theorem 3 can be interpreted as a generalization of Bianchi’s main result in [Bia19b], Theorem 1.4. Let us comment here on the differences between Theorem 3 and [Bia19b, Theorem 1.4]. First, [Bia19b, Theorem 1.4] applies only to maps satisfying a strong assumption, namely that they leave invariant the 3 complex lines $x=\pm{i\varepsilon}$ and $y=0$ (which amounts to taking $b_{\varepsilon}(x,y)=d_{\varepsilon}(x)=0$ with our notations). If that is the case, then the formal curve corresponding to our assumptions $(H_{1})$ and $(H_{2})$ is the curve $y=0$, and it is invariant for the whole family of perturbations. Secondly, [Bia19b, Theorem 1.4] only proves convergence near the line $y=0$ (instead of the whole parabolic basin associated to $v=(1,0)$), and only up to extraction. In particular, it does not rule out the possibility of the sequence $(g_{\varepsilon_{n}}^{n})_{n\in\mathbb{N}}$ having more than one limit value. Finally, in [Bia19b, Theorem 1.4] the possible limits of $(g_{\varepsilon_{n}}^{n})_{n\in\mathbb{N}}$ are not described explicitly as maps of the form $(\Phi^{o})^{-1}\circ A_{\sigma,q}\circ\Phi^{\iota}$, and only depend on the parameter $\sigma$ since in his case $q=0$. On the other hand, we must note that [Bia19b, Theorem 1.4] only assumes $\mathrm{Re}\,\eta>1$, compared to our assumption that $\mathrm{Re}\,\eta>3$, so our results do not strictly imply his.

Let us now give an application case of our main results. A holomorphic endomorphism of $\mathbb{P}^{2}$ is a map which may be written in homogeneous coordinates as

where $P_{0},P_{1},P_{2}:\mathbb{C}^{3}\to\mathbb{C}$ are homogeneous polynomials of degree $d$ with no common factors. The integer $d$ is called the algebraic degree of $f$. Given such an endomorphism $f$, one can define two distinct notions of Julia sets: the set $J_{1}(f)$, which may be defined as the non-normality locus; and the set $J_{2}(f)$, also sometimes called the small Julia set, which may be defined as the support of the unique measure of maximal entropy. Alternatively, $J_{m}(f)$ ($1\leq m\leq 2$) may be defined as the support of $T_{f}^{\wedge m}$, where $T_{f}$ is the so-called Green current of $f$ (this construction is not specific to the case of dimension 2). In general, $J_{2}(f)\subsetneq J_{1}(f)$, and even for very simple maps (such as product maps) there is no equality. We refer the reader to the survey [DS10] for more details.

In the following, for any $(\sigma,q)\in\mathbb{C}^{2}$ we will let $\mathcal{L}_{\sigma,q}:=\Psi^{o}\circ A_{\sigma,q}\circ\Phi^{\iota}$, where $\Psi^{o}$ is the extension of $(\Phi^{o})^{-1}$ to $\mathbb{C}^{2}$. For any endomorphism $f$ of $\mathbb{P}^{2}$ satisfying $(H_{1})$, we will set

Even if we drop the assumption that $\Psi^{o}(\mathbb{C}^{2})\cap J_{2}(f)\neq\emptyset$, our arguments still prove the discontinuity of the closure of the set of repelling periodic points; however, as mentioned above, the Julia set $J_{2}$ may be smaller than this closure. This hypothesis is not easy to check in practice; let us however give a concrete example. Let $\eta\in\mathbb{C}$ with $\mathrm{Re}\,\eta>3$, and let $d>\mathrm{Re}\,\eta$. Let

The polynomial map $f:\mathbb{C}^{2}\to\mathbb{C}^{2}$ extends to an endomorphism of $\mathbb{P}^{2}$. If $(a,b)=(0,0)$, then the map $f$ is a polynomial skew-product. In that case, $\Psi^{o}$ is of the form $\Psi^{o}(x,y)=(\Psi_{p}^{o}(x),\Psi_{2}^{o}(x,y))$, where $\Psi_{p}^{o}$ is an outgoing Fatou parametrization of the base polynomial map $p(z):=z+z^{2}+z^{d}$. By [Jon99], $J_{2}(f):=\overline{\bigcup_{z\in J(p)}J_{z}}$, where $J(p)$ denotes the Julia set of $p$ and $J_{z}$ is the non-normality locus of $\{f^{n}:n\in\mathbb{N}\}$ restricted to the vertical line $x=z$. Since $\Psi_{p}^{o}$ is non-constant and entire, it omits at most one value, so there exists $x_{0}\in J(p)$ and $X_{0}\in\mathbb{C}$ such that $\Psi_{p}^{o}(X_{0})=x_{0}$. Similarly, the map $Y\mapsto{\Psi_{2}^{o}}(X_{0},Y)$ is entire and non-constant and $J_{x_{0}}$ is uncountable; so $\Psi^{o}(\mathbb{C}^{2})\cap(J_{2}(f))\neq\emptyset$. Now, the set $J_{2}(f)$ varies lower semi-continuously with respect to the parameters $(a,b,c)$; and the map $\Psi^{o}$ depends holomorphically (hence continuously) on $(a,b,c)$. (For a proof of this fact in dimension one, see the Appendix in [ABD⁺16]; the argument remains valid in higher dimension). Therefore, there exists some open set $W\subset\mathbb{C}^{3}$ such that for all $(a,b,c)\in\mathbb{C}^{3}$, the map $f$ satisfies $\Psi_{f}^{o}(\mathbb{C}^{2})\cap J_{2}(f)\neq\emptyset$.

Let us first show how Theorem 3 implies Theorem 1. Let $f:U\subset\mathbb{C}^{2}\to\mathbb{C}^{2}$ be a holomorphic map satisfying $(H_{1})$. If we choose coordinates $(x,y)$ such that $v=(1,0)$ we have that

with $\lambda\neq 0$ and $\eta/\lambda=\alpha+1$, where $\alpha$ is the director of $v$. Up to conjugating by the linear map $(x,y)\mapsto(\lambda x,y)$, we can assume that $\lambda=1$, so $\eta=\alpha+1$. In those coordinates, the formal invariant curve $\mathcal{C}$ has a parametrization $\gamma(t)=(t,\zeta(t))$, with $\zeta(t)\in t^{2}\mathbb{C}[[t]]$. If we take $m:=\lfloor\mathrm{Re}\,\eta\rfloor$ and $\Psi(x,y):=(x,y-J_{m+2}\zeta(x))$, where $J_{m+2}\zeta$ is the jet of order $m+2$ of $\zeta$, we have that $g_{0}(x,y)=\Psi\circ f\circ\Psi^{-1}$ has the form

with $a_{0}(x)=1+\mathcal{O}(x)$, $b_{0}(0,0)=0$ and $c_{0}(x,y)=\eta x+\mathcal{O}(x^{2},y)$. Then Theorem 1 follows immediately considering a family $(g_{\varepsilon})$ and a sequence $(\varepsilon_{n})$ as in Theorem 3 and taking $f_{\varepsilon}:=\Psi^{-1}\circ g_{\varepsilon}\circ\Psi$.

Let us now obtain Theorem 2 from Theorem 3. Consider a family $(f_{\varepsilon}:U\subset\mathbb{C}^{2}\to\mathbb{C}^{2})_{\varepsilon\in\mathbb{D}}$ of holomorphic maps satisfying $(H_{1})-(H_{3})$. We will start making several successive changes of coordinates until we obtain the form of Theorem 3.

As above, we first choose coordinates in which $v=(1,0)$ and the curve $\mathcal{C}$ has a parametrization $\gamma(t):=(t,\zeta(t))$, with $\zeta(t)\in t^{2}\mathbb{C}[[t]]$. Let $m:=\lfloor\mathrm{Re}\,\alpha\rfloor+1$ and $\Psi(x,y):=(x,y-J_{m+1}\zeta(x))$, where $J_{m+1}\zeta$ is the jet of order $m+1$ of $\zeta$, and set $\widetilde{g}_{\varepsilon}:=\Psi\circ f_{\varepsilon}\circ\Psi^{-1}$. Since the family $(f_{\varepsilon})$ preserves $\mathcal{C}$ up to order $m+1$ by hypothesis $(H_{2})$, we have

where $p_{\varepsilon}$ depends holomorphically on $\varepsilon$, and $p_{0}(x)=x+\lambda x^{2}+\mathcal{O}(x^{3})$ for some constant $\lambda\neq 0$ (as above, the fact that $\lambda\neq 0$ follows from the assumptions of $f_{0}$ having order 2 and $v$ being non-degenerate). Up to conjugating by the linear map $(x,y)\mapsto(\lambda x,y)$, we may assume without loss of generality that $\lambda=1$.

In particular, there exist holomorphic functions $\widetilde{b}_{\varepsilon}$ and $\widetilde{c}_{\varepsilon}$ such that

Since $\widetilde{g}_{0}$ is tangent to the identity, we also have $\widetilde{b}_{0}(0,0)=\widetilde{c}_{0}(0,0)=0$. Moreover, by hypotheses $(H_{3})$, $\widetilde{g}_{\varepsilon}$ has 4 fixed points $\widetilde{z}_{i}(\varepsilon)=\Psi(z_{i}(\varepsilon))$ ($1\leq i\leq 4$) near the origin, which depend holomorphically on $\varepsilon$ and satisfy that $\widetilde{z}_{i}^{\prime}(0)$ are non-zero and pairwise different.