Nicolas Marque
email: nicolas.marque@imj-prg.fr

Minimal bubbling for Willmore surfaces

In this paper we build an explicit example of a minimal bubble on a Willmore surface, showing there cannot be compactness for Willmore immersions of Willmore energy above 16π. Additionally we prove an inequality on the second residue for limits sequences of Willmore immersions with simple minimal bubbles. Doing so, we exclude some gluing configurations and prove compactness for immersed Willmore tori of energy below 12π.

Introduction

The following is primarily concerned with the study of Willmore immersions in R 3 . Let Φ be an immersion from a closed Riemann surface Σ into R 3 . We denote by g := Φ * ξ the pullback by Φ of the Euclidean metric ξ of R 3 , also called the first fundamental form of Φ or the induced metric. Let dvol g be the volume form associated with g. The Gauss map n of Φ is the normal to the surface. In local coordinates (x, y):

n := Φ x × Φ y Φ x × Φ y , 2 N. Marque
where Φ x = ∂ x Φ, Φ y = ∂ y Φ and × is the usual cross product in R 3 . Denoting π n the orthonormal projection on the normal (meaning π n (v) = n, v n), the second fundamental form of Φ at the point p ∈ Σ is defined as follows.

A p (X, Y ) := A p (X, Y ) n := π n d 2 Φ (X, Y ) for all X, Y ∈ T p Σ.

The mean curvature of the immersion at p is then

H(p) = H(p) n = 1 2 T r g (A) n,
while its tracefree second fundamental form is

Åp (X, Y ) = A p (X, Y ) - 1 2 H(p)g p (X, Y ).
The Willmore energy is defined as

W ( Φ) := Σ H 2 dvol g .
Willmore immersions are critical points of this Willmore energy, and satisfy the Willmore equation:

∆ g H + Å 2 H = 0. (1)
The Willmore energy was already under scrutiny in the XIXth century in the study of elastic plates, but to our knowledge W. Blaschke was the first to state (see [START_REF] Blaschke | Vorlesungen über Integralgeometrie[END_REF]) its invariance by conformal diffeomorphisms of R 3 (which was later rediscovered by T. Willmore, see [START_REF] Willmore | Riemannian geometry[END_REF]) and to study it in the context of conformal geometry. While the Willmore energy is the canonically studied Lagrangian, and serves as a natural measure of the complexity of a given immersion, its invariance is contextual. Indeed W is not invariant by inversions whose center is on the surface, with the simplest example being the Euclidean sphere which is sent to a plane once inverted at one of its points. The true pointwise conformal invariant (as shown by T. Willmore, [START_REF] Willmore | Riemannian geometry[END_REF]) is in fact Åp dvol gp . The total curvature and tracefree curvature are then two relevant energies, respectively defined as follows:

E( Φ) := Σ A 2 g dvol g = Σ |∇ g n| 2 dvol g , E( Φ) := Σ Å 2 g dvol g .
Quick and straightforward computations (done for instance in appendix A.1 of [START_REF] Marque | An ε-regularity result with mean curvature control for Willmore immersions and application to minimal bubbling[END_REF] in a conformal chart) ensure that both

E( Φ) = 4W ( Φ) -4πχ(Σ) (2)
with χ(Σ) the Euler characteristic of Σ, and

E( Φ) = 2W ( Φ) -4πχ(Σ). (3)
The conformal invariance of W when the topology of the surface is not changed then follows from [START_REF] Bernard | Singularity removability at branch points for Willmore surfaces[END_REF]. A Willmore surface is thus a critical point of W , E and E.

In the study of the moduli spaces of Willmore immersions, the compactness question has proven pivotal. E. Kuwert and R. Schätzle (see [START_REF] Kuwert | The Willmore flow with small initial energy[END_REF]) and later T. Rivière (in arbitrary codimension see for instance theorem I.5 in [START_REF] Rivière | Analysis aspects of Willmore surfaces[END_REF]) showed that Willmore immersions follow an ε-regularity result. These induce a now classical concentration of compactness dialectic, as originally developed by J. Sacks and K. Uhlenbeck, for Willmore surfaces with bounded total curvature (or alternatively, given [START_REF] Bernard | Noether's theorem and the Willmore functional[END_REF], bounded Willmore energy and topology). In essence, sequences of Willmore surfaces converge smoothly away from concentration points, on which trees of Willmore spheres are blown (see [START_REF] Brezis | Convergence of solutions of H-systems or how to blow bubbles[END_REF] for an exploration of the bubble tree phenomenon in another simpler case). Y. Bernard and T. Rivière developed an energy quantization result for such sequences of Willmore immersions assuming their conformal class is in a compact of the Teichmuller space (see theorem I.2 in [START_REF] Bernard | Energy quantization for Willmore surfaces and applications[END_REF]). P. Laurain and T. Rivière then showed one could replace the bounded conformal class hypothesis by a weaker convergence of residues linked with the conservation laws. Since we will work with bounded conformal classes we here give abridged versions of theorems I.2 and I.3 of [START_REF] Bernard | Energy quantization for Willmore surfaces and applications[END_REF].

Minimal bubbling for Willmore surfaces 3 Theorem 1.1. Let Φ k be a sequence of Willmore immersions of a closed surface Σ. Assume that lim sup k→∞ W ( Φ k ) < ∞, and that the conformal class of Φ * k ξ remains within a compact subdomain of the moduli space of Σ. Then modulo extraction of a subsequence, the following energy identity holds

lim k→∞ W ( Φ k ) = W ( Φ ∞ ) + p s=1 W ( η s ) + q t=1 W ( ζ t ) -4πθ t ,
where Φ ∞ (respectively η s , ζ t ) is a possibly branched smooth immersion of Σ (respectively S 2 ) and θ t ∈ N.

Further there exists a 1 . . . a n ∈ Σ such that

Φ k → Φ ∞ in C ∞ loc Σ\{a 1 , . . . , a n }
up to conformal diffeomorphisms of R 3 ∪ {∞}. Moreover there exists a sequence of radii ρ s k , points x s k ∈ C converging to one of the a i such that up to conformal diffeomorphisms of R 3

Φ k (ρ s k y + x s k ) → η s • π -1 (y) in C ∞ loc (C\{finite set}) .
Finally there exists a sequence of radii ρ t k , points x t k ∈ C converging to one of the a i such that up to conformal diffeomorphisms of R 3

Φ k ρ t k y + x t k → ι pt • ζ t • π -1 (y) in C ∞ loc (C\{finite set}) .
Here ι pt is an inversion at p ∈ ζ t (S 2 ). The integer θ t is the density of ζ t at p t .

While theorem 1.1 states an energy quantization for W , equality VIII.8 in [START_REF] Bernard | Energy quantization for Willmore surfaces and applications[END_REF] offers in fact a stronger energy quantization for E (and one for E follows). The a i are the aforementioned concentration points and the η s and ι pt • ζ t are the bubbles blown on those concentration points. More precisely, the η s are the compact bubbles, while the ι pt • ζ t are the non compact ones. Non-compact bubbles stand out as a consequence of the conformal invariance of the problem (see [START_REF] Laurain | Asymptotic analysis for surfaces with large constant mean curvature and free boundaries[END_REF] to compare with the bubble tree extraction in the constant mean curvature framework). One might notice that W (ι pt • ζ t ) = W ( ζ t ) -4πθ t , and deduce that if W ( ζ t ) = 4πθ t , then the bubble ι pt • ζ t is minimal. This case, which we will refer to as minimal bubbling will be of special interest to us in this article. Furthermore, if there is only one bubble at a given concentration point we will call the bubbling simple.

Works from Y. Li in [START_REF] Li | Some remarks on Willmore surfaces embedded in R 3[END_REF] (see also [START_REF] Laurain | Energy quantization of Willmore surfaces at the boundary of the moduli space[END_REF]) ensure that compact simple bubbles cannot appear. These studies, furthered by P. Laurain and T. Rivière (see theorem 0.2 of [START_REF] Laurain | Energy quantization of Willmore surfaces at the boundary of the moduli space[END_REF], written just below) have yielded a compactness result for Willmore immersions of energy strictly below 12π. Theorem 1.2. Let Σ be a closed surface of genus g ≥ 1 and Φ k : Σ → R 3 a sequence of Willmore immersions such that the induced metric remains in a compact set of the moduli space and

lim sup k→∞ W Φ k < 12π.
Then there exists a diffeomorphism ψ k of Σ and a conformal transformation

Θ k of R 3 ∪ {∞}, such that Θ k • Φ k • ψ k converges up to a subsequence toward a smooth Willmore immersion Φ ∞ : Σ → R 3 in C ∞ (Σ).
In the aforementioned paper P. Laurain and T. Rivière put forth a potential candidate for Willmore bubbling, consisting of an Enneper bubble glued on the branch point of an inverted Chen-Gackstatter torus, with an energy of exactly 12π (see [START_REF] Chen | Elliptische und hyperelliptische Funktionen und vollständige Minimalflächen vom Enneperschen Typ[END_REF] for the definition of the Chen-Gackstatter torus).

However before considering genus one sequences, a study of the spherical case offers interesting perspectives. Indeed in his seminal work [START_REF] Bryant | A duality theorem for Willmore surfaces[END_REF], R. Bryant offered a classification of Willmore immersions of a sphere in R 3 , showed they were conformal transforms of minimal immersions (see theorem F), and thus that their Willmore energy was 4π-quantized. Moreover while giving a complete description of the Willmore immersions of energy 16π (part 5), R. Bryant remarked: "Surprisingly, this space [of Willmore immersions of energy 16π] is not compact."

It is then interesting to consider whether one can degenerate a sequence of 16π immersions into a bubble blown on a Willmore sphere. A quick study direct our search toward the most likely case: a sequence degenerating into an Enneper immersion glued on the branch point of the inverse of a López minimal surface (that is a minimal sphere with one branched end of multiplicity 3 and one end of multiplicity 1, see [START_REF] Lamm | Branched Willmore spheres[END_REF] below for a parametrization). This will be our first result: Theorem 1.3. There exists Φ k : S 2 → R 3 a sequence of Willmore immersions such that

W ( Φ k ) = 16π, and Φ k → Φ ∞ , smoothly on S 2 \{0}
, where Φ ∞ is the inversion of a López surface. Moreover,

lim k→∞ E( Φ k ) = E( Φ ∞ ) + E(Ψ ∞ ),
where Ψ ∞ : C → R 3 is the immersion of an Enneper surface.

Theorem 1.3 proves that minimal bubbles can appear and thus that Willmore immersions are not compact. It might also indicate the possibility of gluing an Enneper bubble on an inverted Chen-Gackstatter torus. However R. Bryant's classification result proves that one cannot glue an Enneper bubble on an inverted Enneper surface (the resulting surface would be of energy 12π, and thus limit of Willmore immersions of equal energy, which R. Bryant showed did not exist). The local behavior of the limit surface around its branch point needs then to be constrained in order to forbid this case. Since the Chen-Gackstatter torus and the Enneper surface are asymptotic near their branched end there is hope yet to eliminate this configuration.

The behavior of a Willmore surface around a branch point can be fully described by an expansion, proven by Y. Bernard and T. Rivière in theorem 1.8 of [START_REF] Bernard | Singularity removability at branch points for Willmore surfaces[END_REF].

Theorem 1.4. Let Φ ∈ C ∞ (D\{0}) ∩ W 2,2 ∩ W 1,∞ (
D) be a Willmore conformal branched immersion whose Gauss map n lies in W 1,2 (D) and with a branch point at 0 of multiplicity θ + 1. Let λ be its conformal factor, γ 0 the first residue defined as

γ 0 := 1 4π ∂D ν. ∇ H -3π n ∇ H + ∇ ⊥ n × H .
Then there exists α ∈ Z such that α ≤ θ and locally around the origin, Φ has the following asymptotic expansion:

Φ(z) = Az θ+1 + θ+1-α j=1 B j z θ+1+j + C α |z| 2(θ+1) z -α -C γ 0 ln |z| 2(m+1) -4 + ξ (z) ,
where B j , C α ∈ C 3 are constant vectors, A ∈ C\{0}, and C ∈ R. Furthermore ξ satisfies the estimates

∇ j ξ(z) = O |z| 2(θ+1)-α-j+1-κ for all κ > 0 and j ≤ θ + 2 -α, |z| -θ ∇ θ-α+3 ξ ∈ L p for all p < ∞.
In particular:

H(z) = E α z-α -γ 0 ln |z| + η(z),
where E α ∈ C 3 \{0}. The function η satisfies

∇ j η(z) = O |z| 1-j-α-κ for all κ > 0 and j ≤ θ -α, |z| θ ∇ θ+1-α η ∈ L p for all p < ∞.
Minimal bubbling for Willmore surfaces 5

In the specific case of limits of Willmore immersions, the punctured disk described in theorem 1.4 is in fact the limit of simply connected disks, on which the first residue is null (see remark 1.1 in [START_REF] Laurain | Energy quantization of Willmore surfaces at the boundary of the moduli space[END_REF]). Since away from

the concentration point ∇ H k -3π n k ∇ H k + ∇ ⊥ n k × H k converges
, the first residue around branch points of limit Willmore surfaces is always null. Such surfaces are called true Willmore surfaces. The quantity α, although called the second residue (see definition 1.7 in [START_REF] Bernard | Singularity removability at branch points for Willmore surfaces[END_REF]), is not actually a residue and is thus not necessarily null. It will then take center stage in the study of limit Willmore surfaces. With this tool we can refine our understanding of the behavior around minimal concentration points with the following theorem: Theorem 1.5. Let Φ k be a sequence of Willmore immersions of a closed surface Σ satisfying the hypotheses of theorem 1.1. Then at each concentration point p ∈ Σ of multiplicity θ p + 1 on which a simple minimal bubble is blown, the second residue α p of the limit immersion Φ ∞ satisfies

α p ≤ θ p -1.
One should be aware that minimal bubbling is not necessarily simple. Indeed one could for instance imagine an Enneper surface bubbling on the branch point of a minimal surface of Enneper-Weierstrass data (f, g) = (z 2 , z), itself glued on a branch point of multiplicity 5. However piling minimal spheres that way increases the total multiplicity. Minimal bubbling on branch points of multiplicity 3 is thus simple. Consequently theorem 1.5 allows us to eliminate some surfaces as a support for minimal bubbling.

Corollary 1.6. The convergence of Willmore immersions cannot lead to a minimal bubble and an inverted Chen-Gackstatter torus.

Let us recall that a Chen-Gackstatter torus is a minimal torus with one branched end, which is asymptotic to an Enneper end of multiplicity 3. Its inversion is thus a branched Willmore torus with one branch point of multiplicity 3 and Willmore energy 12π. We refer the reader to [START_REF] Chen | Elliptische und hyperelliptische Funktionen und vollständige Minimalflächen vom Enneperschen Typ[END_REF] or [START_REF] López | Complete minimal surfaces in R 3[END_REF] for a complete study.

We can now extend theorem 1.2:

Theorem 1.7. Let Σ be a closed surface of genus 1 and Φ k : Σ → R 3 a sequence of Willmore immersions such that the induced metric remains in a compact set of the moduli space and lim sup k→∞

W Φ k ≤ 12π.
Then there exists a diffeomorphism ψ k of Σ and a conformal transformation

Θ k of R 3 ∪ {∞}, such that Θ k • Φ k • ψ k converges up to a subsequence toward a smooth Willmore immersion Φ ∞ : Σ → R 3 in C ∞ (Σ).
Proof . We only have to exclude the case lim sup

k→∞ W Φ k = 12π. ( 4
)
Consider then Φ k satisfying (4) and converging toward Φ ∞ away from a finite number of concentration points. We consider a concentration point and reason on its multiplicity θ 0 + 1. If θ 0 ≥ 1, using corollary 3.3 (see below), the bubble glued on its concentration point is branched, with the same multiplicity. Using proposition C.1 in [START_REF] Laurain | Energy quantization of Willmore surfaces at the boundary of the moduli space[END_REF] ensures that the multiplicity is odd, and then θ 0 ≥ 2. Given P. Li and S. Yau's inequality (see [START_REF] Li | A new conformal invariant and its applications to the Willmore conjecture and the first eigenvalue of compact surfaces[END_REF]) and (4), Φ ∞ has a Willmore energy of exactly 12π meaning that the branch point is of multiplicity exactly 3, that the bubbles have no Willmore energy (i.e. they are minimal and more accurately Enneper). Using formulas from [START_REF] Lamm | Branched Willmore spheres[END_REF], detailed in appendix (see propositions A.9 and A.10), Φ ∞ is the inverse of a minimal torus of total curvature -8π. The main result of [START_REF] López | The classification of complete minimal surfaces with total curvature greater than -12π[END_REF] ensures that this minimal torus is a Chen-Gackstatter immersion. We are then in the case excluded by corollary 1.6.

If the concentration point is not branched, we refer the reader to the concluding remark of P. Laurain and T. Rivière's [START_REF] Laurain | Energy quantization of Willmore surfaces at the boundary of the moduli space[END_REF] (found just before the appendix) which states that the energy is then at least β 1 + 12π, where β 1 is the infimum of the Willmore energy of Willmore tori. We would then be above our 12π ceiling, which concludes the proof.

The compactness of Willmore tori could fail with an energy strictly above 12π. One could imagine a sequence of non conformally minimal tori (similar to the ones described by U. Pinkall in [START_REF] Pinkall | Hopf tori in S 3[END_REF]) of Willmore energy 12π + δ which degenerates into a branched torus of same energy, with an Enneper bubble. To avoid contradicting theorem 1.5, the branch point would need a second residue α ≤ 1. The existence of such a true Willmore torus is key in understanding compactness above 12π. In addition, one must notice that under the conclusion of theorem 1.5, A. Michelat and T. Rivière, in [START_REF] Michelat | The classification of branched Willmore spheres in the 3-sphere and the 4-sphere[END_REF], have proven that the Bryant's quartic, denoted Q, of the limit surface is then holomorphic, meaning constant in the torus case. Since, according to R. Bryant's [START_REF] Bryant | A duality theorem for Willmore surfaces[END_REF], Q = 0 implies that the surface is the inversion of a minimal surface, one could hope to push A. Michelat and T. Rivière's reasoning to the next order in the special case of a minimal bubble and conclude that the limit surface is an inversion of a minimal surface. Its energy would then be at least 16π, which would get us closer to the compactness strictly below 16π to surfaces of genus lower than 1. The possible counter example would be a 12π branched minimal torus on which a non conformally minimal, non compact Willmore sphere is blown.

To extend the compactness below 12π to immersions of a higher genus, it would be enough to show that the only minimal immersions of critical curvature are asymptotic to the Enneper surface near their branched end. This would be in agreement with the conjecture that the Chen-Gackstatter immersions are the only one with critical curvature for a given genus.

Section 2 will prove theorem 1.3 and build an exemple of minimal bubbling. Then section 3 will be devoted to translating theorem 1.5 in local conformal charts and to the possible adjustments, in notation or with the conformal group, that can be done to simplify the problem. Section 4 will give the first expansions on the conformal factor. Section 5 will prove theorem 1.5 and its corollaries.

Acknowledgments: The author would like to thank his advisor Paul Laurain for his support and precious advices. This work was partially supported by the ANR BLADE-JC ANR-18-CE40-0023.

2 Proof of theorem 1.3

We will build a sequence of Willmore immersions whose energy E concentrates on a point where an Enneper bubble blows up. Working from section 5 of R. Bryant's [START_REF] Bryant | A duality theorem for Willmore surfaces[END_REF], we study a family of four ended minimal immersions

Ψ µ : C\{a 1 , a 2 , a 2 } → R 3 : Ψ µ = 2 (f µ ) f µ = a 1 z -µ + a 2 z -µj + a 3 z -µj 2 + a 4 z, (5)
with a 1 , a 2 , a 3 , a 4 ∈ C 3 , j 3 = 1, and µ a real parameter that will go toward 0. As explained in [START_REF] Bryant | A duality theorem for Willmore surfaces[END_REF] the ( a i ) must be constrained for Ψ µ to be a conformal immersion. Indeed:

(Ψ µ ) z , (Ψ µ ) z = (f µ ) z , (f µ ) z = a 1 , a 1 (z -µ) 4 + a 2 , a 2 (z -µj) 4 + a 3 , a 3 (z -µj 2 ) 4 + a 4 , a 4 + 2 a 1 , a 2 (z -µ) 2 (z -µj) 2 + 2 a 1 , a 3 (z -µ) 2 (z -µj 2 ) 2 + 2 a 2 , a 3 (z -µj) 2 (z -µj 2 ) 2 - 2 a 1 , a 4 (z -µ) 2 - 2 a 2 , a 4 (z -µj) 2 - 2 a 3 , a 4 (z -µj 2 ) 2 .
Further since given u, v ∈ C:

1 (z -u) 2 (z -v) 2 = 1 (u -v) 2 1 (z -u) 2 + 1 (u -v) 2 1 (z -v) 2 - 2 (u -v) 3 1 z -u + 2 (u -v) 3 1 z -v ,
we deduce that (Ψ µ ) z , (Ψ µ ) z = 0 if and only if

a 1 , a 1 = a 2 , a 2 = a 3 , a 3 = a 4 , a 4 = 0, a 1 , a 2 = a 1 , a 3 = a 2 , a 3 , a 4 = - 1 3µ 2 a 1 + j a 2 + j 2 a 3 . (6)
One can check that under the conditions ( 6), ( a 1 , a 2 , a 3 ) is a linearly independant family of C 3 and thus that Ψ µ is an immersion.

Minimal bubbling for Willmore surfaces 7

Here we take, with b ∈ C a parameter to be adjusted later,

a 1 = 1 2µ 2 e 0 , a 2 = j 2µ 2 e 0 - µ 2 b 2 2 e 0 + bj 2 e 3 , a 3 = j 2 2µ 2 e 0 - µ 2 b 2 2 e 0 -bj e 3 ,
where

e 0 =   1 i 0   , e 3 =   0 0 1   .
One can check that these (a i ) satisfy [START_REF] Brezis | Convergence of solutions of H-systems or how to blow bubbles[END_REF]. Computing, we find:

f µ = 3 z 3 -µ 3 1 2 e 0 -b 2 µ 2 2z + µ z 2 + µz + µ 2 + z 3 1 2 e 0 + bj(j -1) z + µ z 2 + µz + µ 2 e 3 .
To simplify this expression we set b = 3a 2j(j-1) with a ∈ C to be fixed at the end of the reasoning, and reach:

f µ = 3 z 3 -µ 3 1 2 e 0 + a 2 4 3µ 2 2z + µ z 2 + µz + µ 2 + z 1 2 e 0 + 3a 2
z + µ z 2 + µz + µ 2 e 3 . (7)
Then Ψ µ : S 2 → R 3 is a sequence of minimal immersions with four simple planar ends. Applying propositions A.9 and A.10 (see in the appendix) we find:

S 2 K Ψµ dvol g Ψµ = -12π, (8)
S 2 Å 2 Ψµ dvol g Ψµ = 24π. (9)
Letting µ → 0 in [START_REF] Bryant | A duality theorem for Willmore surfaces[END_REF] we find that, away from 0,

f µ → f 0 = 3 2z 3 e 0 + a 2 z 8 e 0 + 3a 2z e 3 , (10)
and deduce that Ψ µ → Ψ 0 := 2 ( f 0 ) smoothly away from 0, where Ψ 0 is a branched minimal immersion of the sphere with one simple planar end and one planar end of multiplicity 3. This immersion is in fact the López minimal surface mentioned in theorem 1.3. Then

S 2 K Ψ0 dvol g Ψ 0 = -8π, (11)
S 2 Å 2 Ψ0 dvol g Ψ 0 = 16π. ( 12
)
Let p be a point in R 3 such that d(p, Ψ µ ) > 1. We now introduce Φ µ := ι p • Ψ µ and Φ 0 := ι p • Ψ 0 , with ι(x) = x-p |x-p| 2 the inversion in R 3 centered at p. Then Φ µ is a sequence of closed Willmore conformal immersions of the sphere converging toward Φ 0 smoothly away from 0, and Φ 0 is a closed Willmore conformal branched immersion of the sphere with a single branch point of multiplicity 3 at 0. Thus

S 2 K Φµ dvol g Φµ = 4π, ( 13
)
S 2 K Φ0 dvol g Φ 0 = 8π. ( 14
)
Since Å 2 dvol g is a conformal invariant, we deduce from ( 9) and ( 12):

S 2 Å 2 Φµ dvol g Φµ = 24π, (15)
S 2 Å 2 Φ0 dvol g Φ 0 = 16π. ( 16
)
With proposition A.10 (see in the appendix) we conclude with ( 13) and ( 15):

S 2 H 2 Φµ dvol g Φµ = 1 2 S 2 Å 2 Φµ dvol g Φµ + S 2 K Φµ dvol g Φµ = 16π, (17)
and with ( 14) and ( 16):

S 2 H 2 Φ0 dvol g Φ 0 = 1 2 S 2 Å 2 Φ0 dvol g Φ 0 + S 2 K Φ0 dvol g Φ 0 = 16π. (18)
Comparing ( 15)-( 18) reveals that while: W ( Φ µ ) → W ( Φ 0 ), there is an energy gap of 8π in E (or equivalently in E). From this, and the energy quantization theorem (theorem 1.1, written above), we deduce that a simple minimal bubble of energy E = 8π is blown. The only possible bubble is then an Enneper surface (see for instance [START_REF] Osserman | A survey of minimal surfaces[END_REF]), given by:

E(z) = 2 z 2 e 0 + z 2 2 e 3 - z 3 6 e 0 . (19)
This is enough to ensure that the immersions Φ µ offer an exemple of an Enneper bubble appearing on a sequence of Willmore immersions, which proves theorem 1.3. We however wish to make the appearance of the Enneper bubble explicit in the computations. To do that we will perform a blow-up at the origin at scale µ 3 . This concentration scale has been determined the classical way (see the bubble tree extraction procedure in [START_REF] Laurain | Energy quantization of Willmore surfaces at the boundary of the moduli space[END_REF] or [START_REF] Bernard | Energy quantization for Willmore surfaces and applications[END_REF]) by computing

∇ n Φµ L ∞ (S 2 )
. Since these computations do not by themselves further the understanding of the bubbling phenomena, they are omitted. Considering [START_REF] Bryant | A duality theorem for Willmore surfaces[END_REF] we find

f µ µ 3 z = 3 2µ 3 (µ 6 z 3 -1) e 0 + a 2 µ 8 3 2zµ 2 + 1 1 + µ 2 z + µ 4 z 2 + µ 2 z e 0 + 3a 2µ 1 + µ 2 z 1 + µ 2 z + µ 4 z 2 e 3 = - 3 2µ 3 - 3µ 3 z 3 2 + O µ 9 e 0 + a 2 2 3µ 4 + zµ 3 -+O(µ 5 ) e 0 + a 3 2µ - 3 2 µ 3 z 2 + O(µ 5 ) e 3 .
Which means that

Ψ µ (µ 3 z) = 1 2 3a 2 µ 4 - 3 µ 3 + µ 3 a 2 z -3z 3 + O(µ 5 ) e 0 + 1 2 3a 2 µ 4 - 3 µ 3 + µ 3 a 2 z -3z 3 + O(µ 5 ) e 0 + 3 (a + a) 2µ - 3µ 3 2 az 2 + az 2 + O(µ 5 ) e 3 .
With p = p1 2 e 0 + p1 2 e 0 + p 3 e 3 defined previously we conclude:

Ψ µ (µ 3 z) -p = 1 µ 3   1 2 -3 -µ 3 p 1 + 3a 2 µ 4 4 + µ 6 a 2 z -3z 3 + O(µ 8 )   1 i 0   + 1 2 -3 -µ 3 p 1 + 3a 2 µ 4 4 + µ 6 a 2 z -3z 3 + O(µ 8 )   1 -i 0   + 3(a + a)µ 2 2 -µ 3 p 3 - 3µ 6 2 (az 2 + az 2 ) + O(µ 8 ) e 3 . (20)
Here the only relevant terms are the first non constant ones i.e. those in µ 6 . This yields:

Ψ µ (µ 3 z) -p 2 = 1 µ 6 9 + 3µ 3 (p 1 + p 1 ) + µ 4 9(a 2 + a 2 ) 4 + µ 4 9(a + a) 2 4 - 3µ 5 (a + a) 2 (p 3 + p 3 ) -3µ 6 (a 2 z + a 2 z-3z 3 -3z 3 + |p 1 | 2 + |p 3 | 2 ) +O(µ 7 ) . (21)
Minimal bubbling for Willmore surfaces 9

We can combine ( 20) and ( 21):

Φ µ µ 3 z = Ψ µ -p Ψ µ -p 2 = µ 3 1 2 - 1 3 - µ 3 9 p 1 + a 2 µ 4 12 + µ 6 9 a 2 z -3z 3 + O(µ 7 ) e 0 + 1 2 - 1 3 - µ 3 9 p 1 + a 2 µ 4 12 + µ 6 9 a 2 z -3z 3 + O(µ 7 ) e 0 + a + aµ 2 6 - µ 3 9 p 3 - µ 6 6 (az 2 + az 2 ) + O(µ 7 ) e 3 1 - 1 3 µ 3 (p 1 + p 1 ) -µ 4 (a 2 + a 2 ) 4 -µ 4 (a + a) 2 4 + µ 5 (a + a) 6 (p 3 + p 3 ) + 1 3 µ 6 (a 2 z + a 2 z -3z 3 -3z 3 + |p 1 | 2 + |p 3 | 2 + 1 3 (p 1 + p 1 ) 2 ) + O(µ 7 ) = Φ µ (0) + µ 9 1 9 a 2 z -3z 3 -a 2 z -a 2 z + 3z 3 + 3z 3 1 2 e 0 + 1 9 a 2 z -3z 3 -a 2 z -a 2 z + 3z 3 + 3z 3 1 2 e 0 - az 2 + az 2 6 e 3 + O(µ 10 ) = Φ µ (0) + µ 9 z3 3 - a 2 9 z 1 2 e 0 + z 3 3 - a 2 9 z 1 2 e 0 - a 3
z 2 2 + a 3 z2 2 e 3 + O(µ 10 ).
Taking a = 3 we find exactly

Φ µ µ 3 z = Φ µ (0) -µ 9 E(z) + O µ 10 . ( 22
)
Hence we do have:

Φ µ (µ 3 z) -Φ µ (0) -µ 9 → E(z)
smoothly on every compact of C, which does illustrate theorem 1.3.

Remark 2.0.1. Chosing another value for a would have led to another Enneper surface, with Enneper-

Weierstrass data (f, g) = (1, 3 a z) instead of simply (f, g) = (1, z).
Remark 2.0.2. One must notice the fundamentally asymetric role of Φ 0 (the surface) and E (the bubble). Indeed while we have compactly glued E on Φ 0 we cannot compactly glue Ψ 0 = ι • Φ 0 on an inverted Enneper using the same construction, since Ψ 0 has an end which is not on the concentration point. Doing so would require to glue a closed bubble tree on said planar end (and would necessarily add Willmore energy to the concentration point). Furthermore, theorem 1.5 ensures that no construction will ever enable us to do so, given that the second residue of the inverted Enneper surface is α = 2 (see [START_REF] Bernard | Singularity removability at branch points for Willmore surfaces[END_REF]). While theorem 1.5 is stated globally, its conclusion is localized on the neighborhoods of concentration points on which a simple minimal bubble is blown. We will then work in conformal charts around such points. The aim of this section is to draw a set of hypotheses that these maps may satisfy, sometimes up to slight adjustments that can be done without loss of generality. We will also delve into the first consequences of these hypotheses to clarify our framework.

Convergence in local conformal charts:

We here wish to show the following:

Lemma 3.1. Let ξ k be a sequence of Willmore immersions of a closed surface Σ satisfying the hypotheses of theorem 1.1. Then, in proper conformal charts around a concentration point on which a simple minimal bubble is blown, ξ k yields a sequence of Willmore conformal immersions Φ ε : D → R 3 , of conformal factor

λ ε = 1 2 ln |∇ Φ ε | 2 2
, Gauss map n ε , mean curvature H ε and tracefree curvature Ω ε := 2 Φ ε zz , n ε , satisfying the following set of hypotheses:

1. There exists C 0 > 0 such that Φ ε L ∞ (D) + ∇ Φ ε L 2 (D) + ∇λ ε L 2,∞ (D) + ∇ n ε L 2 (D) ≤ C 0 . 2. Φ ε → Φ 0 C ∞ loc (D\{0})
, where Φ 0 is a true branched Willmore conformal immersion, with a unique branch point of multiplicity θ 0 + 1 at 0, meaning that

Φ 0 z ∼ 0 Az θ0 . ( 23
)
We denote λ 0 its conformal factor, n 0 its Gauss map, H 0 its mean curvature and Ω 0 its tracefree curvature.

Minimal bubbling for Willmore surfaces 11

3. There exists a sequence of real numbers C ε > 0 such that

Φ ε := Φ ε (ε.) -Φ ε (0) C ε → Φ 1 z C ∞ loc (C)
, where Φ 1 is assumed to be a minimal conformal immersion of C with a branched end of multiplicity θ 1 + 1, meaning that: Φ 1 z ∼ ∞ Az θ1 . We denote λ 1 its conformal factor, n 1 its Gauss map, H 1 its mean curvature and Ω 1 its tracefree curvature.

lim

R→∞   lim ε→0 D 1 R \D εR |∇ n ε | 2 dz   = 0.
5. Ω ε e -λ ε reaches its maximum at 0 and

Ω ε e -λ ε (0) = 2 ε .
Proof . Such assumptions are natural if we consider ξ k satisfying the hypotheses of theorem 1.1. Thanks to theorems I.2 and I.3 of [START_REF] Bernard | Energy quantization for Willmore surfaces and applications[END_REF] such a sequence ξ k converges smoothly away from concentration points. In a conformal chart centered on a concentration point, ξ k yields a sequence of conformal, weak Willmore immersions Φ k : D → R 3 converging smoothly away from the origin toward a true Willmore surface (i.e. hypothesis 2). Hypothesis 1 stands if we choose proper conformal charts (see theorem 3.1 of P. Laurain and T. Rivière's [START_REF] Laurain | Optimal estimate for the gradient of Green's function on degenerating surfaces and applications[END_REF] for a detailed explanation). Hypothesis 3 then specifies that we consider the case where there is only one simple minimal bubble which concentrates on 0 in the aforementioned chart. Hypothesis 4 is just the energy quantization once the whole bubble tree is extracted and corresponds to inequality VIII.8 in [START_REF] Bernard | Energy quantization for Willmore surfaces and applications[END_REF]. Moreover, by definition of a concentration point

∇ n k L ∞ (D) → ∞.
On the other hand, the main result of [START_REF] Marque | An ε-regularity result with mean curvature control for Willmore immersions and application to minimal bubbling[END_REF] (namely inequality (96)) states that

H k ∇ Φ k L ∞ (D) ≤ C(C 0 ). Since ∇ n k 2 = H k ∇ Φ k 2 + Ω k e -λ k 2 , necessarily Ω k e -λ k L ∞ (D) → ∞.
We then define the concentration speed as

ε k = 2 Ω k e -λ k L ∞ (D)
, and we assume it is reached at the origin. For simplicity's sake we reparametrize this sequence by the concentration speed which we denote ε. Hypothesis 5 is then a consequence of this slight adjustment.

An immediate consequence of hypotheses 2-4 is the following energy quantization result

D |∇ n ε | 2 dz → D ∇ n 0 2 dz + C ∇ n 1 2 dz. ( 24
)
while hypothesis 1 ensures

∇ n 0 L 2 (D) + ∇ n 1 L 2 (C) ≤ C(C 0 ). ( 25
)
Further, if we denote λ ε the conformal factor of Φ ε , H ε its mean curvature, Ω ε its tracefree curvature and n ε its Gauss map, we have

λ ε = λ ε (ε.) -ln C ε ε , (26)
and

Ω ε e -λ ε = ε Ω ε e -λ ε (ε.) . ( 27
)
Using hypothesis 5 one may conclude that

Ω ε e -λ ε (0) = 2. ( 28
)
Hypothesis 3 then yields

Ω 1 e -λ 1 (0) = 2. ( 29
)
Similarly we know, thanks to hypotheses 3 and 5,

Ω ε e -λ ε z (0) = 0, (30)
and:

Ω 1 e -λ 1 z (0) = 0. ( 31
)
Finally hypothesis 4 allows us to apply theorem 1.4 to Φ 0 : there exists

α ≤ θ 0 , A ∈ C 3 \{0}, B j j=1..θ0+1-α ∈ C 3 , C α ∈ C 3 \{0} and ξ : D → R 3 such that Φ 0 z = Az θ0 + θ0+1-α j=1 B j z θ0+j + C α θ 0 + 1 z θ0-α z θ0+1 + C α θ 0 + 1 -α z θ0 z θ0+1-α + ξ z , (32)
where ξ satisfies:

∇ j ξ = O r 2θ0+3-α-j-κ ,
for all κ > 0 and j ≤ θ 0 + 2 -α. The second residue α is in fact defined as follows (see theorem I.8 in [START_REF] Bernard | Singularity removability at branch points for Willmore surfaces[END_REF]):

H 0 ∼ C α |z| -α . ( 33
)
Our proofs will use the quantities L, S and R, stemming from the Willmore conservation laws (see for instance theorem I.4 in [START_REF] Rivière | Analysis aspects of Willmore surfaces[END_REF]), which at the core, are a consequence of the conformal invariance of W (see [START_REF] Bernard | Noether's theorem and the Willmore functional[END_REF]). More precisely L, S and R are defined as follows:

∇ ⊥ L = ∇ H -3π n ∇ H + ∇ ⊥ n × H, ∇ ⊥ S = L, ∇ ⊥ Φ , ∇ ⊥ R = L × ∇ ⊥ Φ + 2H∇ ⊥ Φ. ( 34
)
Exploiting these was key in T. Rivière's proof of the ε-regularity for Willmore surfaces.

Under hypotheses 1-5, the conclusion of [START_REF] Marque | An ε-regularity result with mean curvature control for Willmore immersions and application to minimal bubbling[END_REF] stands and yields (see ( 96)-(98) in the aforementioned paper):

H ε ∇ Φ ε L ∞ (D) ≤ C(C 0 ), (35)
∇ Φ ε W 3,p (D) ≤ C(C 0 ), (36)
while the second and third Willmore quantities satisfy

∇S ε W 1,p (D) + ∇ R ε W 1,p (D) ≤ C(C 0 ) (37)
for all p < ∞. Up to an inconsequential translation one can further assume Φ ε (0) = 0.

Minimal bubbling for Willmore surfaces 13

3.2 Branch point-branched end correspondence:

The goal of this subsection is to show the equality of the multiplicity of the end of the bubble and the multiplicity of the branch point of the surface:

Theorem 3.2. Let Φ ε : D → R 3 satisfying 1-5. Then θ 0 = θ 1 = θ.
Proof . Since Φ ε is conformal, the Liouville equation states

∆λ ε = K ε e 2λ ε , (38)
where

K ε is the Gauss curvature of Φ ε . Then given R ∈ R + D 1 R \D εR K ε e 2λ ε dz = D 1 R \D εR ∆λ ε dz = ∂D 1 R ∂ r λ ε dσ - ∂D εR ∂ r λ ε dσ = ∂D 1 R ∂ r λ ε dσ - ∂D R ε∂ r λ ε (ε.)dσ = ∂D 1 R ∂ r λ ε dσ - ∂D R ∂ r λ ε dσ. (39)
Besides, hypotheses 2 and 3 ensure that λ ε → λ 0 on ∂D 1 R and λ ε → λ 1 on ∂D R . In addition, since Φ 0 has a branch point of multiplicity θ 0 + 1 at 0, lim

R→∞ D 1 R ∂ r λ 0 dσ → 2πθ 0 . (40)
Similarly Φ 1 has an end of multiplicity θ 1 + 1 at ∞, which implies:

lim R→∞ D R ∂ r λ 1 dσ → 2πθ 1 . (41)
Injecting ( 40) and (41) in (39) yields

2π |θ 0 -θ 1 | ≤ lim R→∞   lim ε→0 ∂D 1 R ∂ r λ ε dσ - ∂D R ∂ r λ ε dσ   ≤ lim R→∞   lim ε→0 D 1 R \D εR K ε e 2λ ε dz   ≤ lim R→∞   lim ε→0 D 1 R \D εR K ε e 2λ ε dz   ≤ lim R→∞   lim ε→0 D 1 R \D εR |∇ n ε | 2 dz   = 0, using hypothesis 1. As a conclusion θ 0 = θ 1 = θ.
While we wrote the proof in the case specific to our paper, it remains valid for any behavior of Φ 0 and Φ 1 (branched points or ends) and relies solely on the energy quantization result. In a broader frame this corresponds to the following: Corollary 3.3. A Willmore bubble with a branched end of multiplicity θ + 1 at infinity can only appear on a branch point of multiplicity θ + 1. A Willmore bubble with a branch point of multiplicity θ -1 at infinity can only appear on a branched end of multiplicity θ -1.

Here, it must be emphasized that the true invariant is the behavior of the conformal factor, transmitted from the bubble to the surface. It then translates as information on the multiplicity of the respective singular points, depending on their nature (branch point or end) and their location (at ∞ or in a compact).

Macroscopic adjustments

This section is dedicated to the proof of the following result: Lemma 3.4. Let Φ ε : D → R 3 be a sequence of Willmore conformal immersions satisfying hypotheses 1-5. Then θ is even, and up to macroscopic adjustments we can assume that: 6.

Φ 1 z = 1 2   Q 2 -P 2 i Q 2 + P 2 2P Q   = Q 2 2 e 0 - P 2 2 e 0 + P Q e 3 ,
where P, Q ∈ C θ 2

[X], P ∧ Q = 1, and

P (0) = P (0) -2Q (0) = 0, Q(0) = P (0) = 1 (42)
The end of multiplicity θ + 1 of Φ 1 can be highlighted as follows: there exists

A ∈ C 3 \{0} and V ∈ C θ-1 [X] such that Φ 1 z = Az θ + V .
Here we used C[X] to denote polynomials with complex coefficients, while for any m ∈ N, C m [X] denotes the space of polynomials of degree lower or equal to m.

Proof . Adjusting with homothetic transformations of R 3 : Since Φ 1 has no branch point on C, up to a fixed rotation and a dilation in R 3 one can assume:

Φ 1 z (0) = 1 2   1 i 0   . (43)
Adjusting the parametrization:

Taking M -ω =   cos ω -sin ω 0 sin ω cos ω 0 0 0 1
  with ω ∈ R to be determined at a later point, we set Ψ := M -ω Φ 1 e iω . .

We denote respectively λ Ψ , H Ψ , Ω Ψ and n Ψ its conformal factor, its mean curvature, its tracefree curvature and its Gauss map. Then Ψ z = e iω M -ω Φ 1 z e iω . which implies n Ψ = M -ω n and e λ Ψ = e λ 1 . Consequently using (43)

Ψ z (0) = e iω 2 M ω   1 i 0   = e iω 2   e -iω ie -iω 0   = 1 2 e 0 . (44)
Thus, one still has e λΨ (0) = 1

. Moreover, we can compute Ψ zz = e 2iω M -ω Φ 1 zz e iω and deduce

Ω Ψ = e 2iω Ω 1 , which in turn implies Ω Ψ e -λΨ (0) = e 2iω Ω 1 e -λ 1 (0).
Thanks to (29), one can choose ω such that e 2iω = -2

[Ω 1 e -λ 1 ](0) . In that case Ω Ψ (0)e -λΨ(0) = -2, which yields thanks to (45)

Ω Ψ (0) = -2. ( 46
)
Comparing (29) to (46), one can see that this adjustment done with rotations ensured a jump from an equality in moduli (namely (29)) to an equality for complex numbers (that is (46)).

The sequence Ψ ε = M ω Φ ε e iω . satisfies hypotheses 1, 4 and 5, while

Ψ ε → M ω Φ 0 e iω . C ∞ loc (D\{0}) , Ψ ε → Ψ C ∞ loc (C) .
We will not change notations for simplicity's sake, and will merely assume, without loss of generality, that

Ω 1 (0) = -2. ( 47
)
Summing up:

Φ 1 z (0) = 1 2 e 0 , Ω 1 (0) = -2, Ω 1 2 e -2λ 1 z (0) = 0. ( 48
)
Consequences on the Enneper-Weierstrass representation: Since Φ 1 is a minimal immersion we can use the Enneper-Weierstrass representation:

Φ 1 z = f 2   1 -g 2 i 1 + g 2 2g   ( 49
)
where f is a holomorphic function on C and g a meromorphic one. Since (according to ( 25)) Φ 1 has finite total curvature, g is a meromorphic function of finite degree on C. Thus there exists two polynomials P, Q ∈ C [X] such that P ∧ Q = 1 and g = P Q . Since Φ 1 has no end on C, f has a zero of order at least 2k at each pole of order k of g (if not, (49) shows that at such a pole Φ z → ∞, that is: this point is an end of Φ). Consequently there exists a holomorphic function f such that f = Q 2 f . In addition, Φ 1 has no branch point on C and one finite end at ∞, thus f is a holomorphic function without zeros and of finite order at infinity, i.e. a constant. We can then write

Φ 1 z = 1 2   Q 2 -P 2 i Q 2 + P 2 2P Q   = Q 2 2   1 i 0   - P 2 2   1 -i 0   + P Q   0 0 1   . (50)
Further, since Φ 1 is assumed to have an end of multiplicity θ + 1 one can expand (50) as

Φ 1 z = Az θ + O z θ-1 , (51)
where A ∈ C 3 \{0}. Comparing (50) and (51) notably implies that θ is even and

P, Q ∈ C θ 2 [X]
. From (50) we then successively deduce

n 1 = 1 |P | 2 + |Q| 2   P Q + P Q i P Q -P Q |P | 2 -|Q| 2   , (52)
e 2λ 1 = |P | 2 + |Q| 2 2 , (53)
Φ 1 zz = Q   Q Qi P   -P   P -iP -Q   , (54)
which implies

Ω 1 = 2 (P Q -P Q) , (55)
and in turn

Ω 1 e -λ 1 = 2 P Q -P Q |P | 2 + |Q| 2 , (56)
Ω 1 e -λ 1 2 = 4 |P | 2 |Q | 2 + |P | 2 |Q| 2 -P P Q Q -P QP Q (|P | 2 + |Q| 2 ) 2 . ( 57
)
Comparing the first equality of (48) to (50) with z = 0 ensures that P (0) = 0 and Q(0) = 1, while taking (55) at z = 0 and comparing it to the second equality of (48) yields P (0) = 1. Differentiating (57) at z = 0 and comparing it to the third equality of (48) yields the last equality of (42). This concludes the proof.

16 N. Marque

Infinitesimal adjustments:

This section will prove the infinitesimal counterpart of theorem 3.4.

Lemma 3.5. Let Φ ε : D → R 3 be a sequence of Willmore conformal immersions satisfying hypotheses 1-6.

Up to infinitesimal adjustments we can assume that:

7. Φ ε z (0) = C ε ε Φ 1 z (0) = C ε 2ε e 0 , Ω ε e -λ ε (0) = 1 ε Ω 1 e -λ 1 (0) = - 2 ε .
Proof . Using homothetic transformations of R 3 : By hypothesis 3, Φ ε z (0) → Φ 1 z (0), thus there exists a sequence of homothetic transformations

σ ε → Id such that σ ε Φ ε z (0) = Φ 1 z (0)
. Since σ ε tends toward the identity, hypotheses 1-3 are still satisfied, and hypothesis 5 still stands due to the conformal invariance properties of the tracefree curvature. We will then apply this sequence of transformations without changing the notations for simplicity's sake and assume Φ ε z (0) = Φ 1 z (0). Considering

Φ ε z = ε C ε Φ ε z ,
we deduce with (48),

Φ ε z (0) = C ε 2ε e 0 . ( 58
)
Adjusting the parametrization: Using (29) and ( 27), we can set

e 2iω ε = Ω ε e -λ ε (0) Ω 1 e -λ 1 (0) = ε Ω ε e -λ ε (0) Ω 1 e -λ 1 (0) → 1.
We consider Ψ ε = M ω ε Φ ε e iω ε . . Since e iω ε → 1, Ψ ε satisfies 1-6 (5 is still satisfied due to the invariance properties of the tracefree curvature). As detailed in the previous section we have

Ω ε Ψ e -λ ε Ψ = -2, which implies Ω ε Ψ e -λ ε Ψ (0) = - 2 ε .
For simplicity's sake we will not change the notations and assume that Φ ε satisfies

Ω ε e -λ ε (0) = - 2 ε . (59)
This gives us the desired result.

Expanding the conformal factor

This section will prove the following expansion on the conformal factor, which will serve as a stepping-stone in the proof of theorem 1.5.

Theorem 4.1. Let Φ ε be a sequence of Willmore conformal immersions satisfying 1-7. Then there exists l ε ∈ L ∞ (D) such that:

λ ε = ln ε θ + r θ + l ε , l ε L ∞ (D) ≤ C(C 0 ).
As a result if we denote χ = √ ε 2 + r 2 , the immersion satisfies the following Harnack inequality:

χ θ C(C 0 ) ≤ e λ ε ≤ C(C 0 )χ θ .
Minimal bubbling for Willmore surfaces 17

Proof .

Step 1: Controls in the neck area Given hypothesis 4, for any ε 0 > 0 arbitrarily small there exists R big enough such that

lim ε→0   D 1 R \D εR |∇ n ε | 2 dz   ≤ ε 0 . (60)
We first recall lemma V.3 of [START_REF] Bernard | Energy quantization for Willmore surfaces and applications[END_REF].

Lemma 4.2. There exists a constant η > 0 with the following property. Let 0 < 4r < R < ∞. If Φ is any (weak) conformal immersion of Ω := D R \D r into R 3 with L 2 -bounded second fundamental form and satisfying

∇ n L 2,∞ (Ω) < √ η,
then there exist 1 2 < α < 1 and A ∈ R depending on R, r, m and Φ such that

λ(x) -d ln |x| -A L ∞ D αR \D r α ≤ C ∇λ L 2,∞ (Ω) + Ω |∇ n| 2 , (61)
where

d satisfies 2πd - ∂Dr ∂ r λdl ∂Dr ≤C D2r\Dr |∇ n| 2 dz + 1 ln R r ∇λ L 2,∞ (Ω) + Ω |∇ n| 2 .
(62) Thus, according to (60), there exists R 0 such that for all R ≥ R 0 and ε small enough, we can apply lemma 4.2 on D 1 R \D εR and conclude that there exists

d ε R and A ε R ∈ R such that λ ε (x) -d ε R ln r -A ε R L ∞ D 1 2R \D 2εR ≤ C 0 , (63)
d ε R - 1 2π ∂D εR ∂ r λ ε dl ∂D εR ≤C D 2εR \D εR |∇ n| 2 dz + C 0 -ln (εR 2 ) . ( 64
)
Here C 0 is the uniform bound given by hypothesis 1 (up to a multiplicative uniform constant). We saw in (41) that

lim R→∞ lim ε→0 1 2π ∂D εR ∂ r λ ε dl ∂D εR = θ, while hypothesis 4 ensures that lim R→∞ lim ε→0 D 2εR \D εR |∇ n| 2 dz = 0.
Hence we can fix R 1 > 0 such that for ε small enough:

|d ε -θ| ≤ 1 10 3 . ( 65
)
Since R 1 is fixed, we will get rid of the subscript on d ε and A ε . Then for any ε small enough:

λ ε (x) -d ε ln r -A ε L ∞ D 1 2R 1 \D 2εR 1 ≤ C(C 0 , R 1 ). ( 66
)
Step 2: Estimates on the exterior boundary: Hypothesis 2 ensures that on ∂D 1 2R 1

, λ ε → λ 0 smoothly, and that λ 0 is a bounded function away from 0, which implies

λ ε L ∞ ∂D 1 2R 1 ≤ C(C 0 , R 1 ). ( 67
)
On the other hand (65) ensures

|d ε ln R 1 | ≤ C(C 0 , R 1 ). ( 68
)
As a result, combining (66) on ∂D 1 2R 1

, (67), and (68) yields

|A ε | ≤ C(C 0 , R 1 ),
which we can inject in (66) to obtain

λ ε (x) -d ε ln r L ∞ D 1 2R 1 \D 2εR 1 ≤ C(C 0 , R 1 ). ( 69
)
Step 3: Estimates on the interior boundary:

Estimate (69) implies λ ε (x) -d ε ln r L ∞ (∂D2εR 1 ) ≤ C(C 0 , R 1 ). ( 70
)
In addition, (26) yields

λ ε (x) -d ε ln r L ∞ (∂D2εR 1 ) = λ ε (x) -d ε ln r + ln C ε ε -d ε ln ε L ∞ (∂D2R 1 )
.

Hypothesis 3 then ensures that λ ε L ∞ (D2R 1 ) ≤ C(C 0 , R 1 ), (71)
and (65) that

d ε ln r L ∞ (D2R 1 ) ≤ C(C 0 , R 1 ). ( 72
)
Together ( 70), ( 71) and (72) yield

ln C ε ε d ε +1 ≤ C(C 0 , R 1 ). (73)
A direct consequence of ( 65) and ( 73) is the following estimate:

ε θ+1-10 -3 C(R 1 , C 0 ) ≤ ε d ε +1 C(C 0 , R 1 ) ≤ C ε ≤ C(C 0 , R 1 )ε d ε +1 ≤ C(C 0 , R 1 )ε θ+1+10 -3 . ( 74
)
Step 4: Expanding the conformal factor on the whole disk:

We forcefully write λ ε = ln ε d ε + r d ε + l ε . We aim to show that

l ε L ∞ (D) ≤ C(C 0 , R 1 ). On D\D 1 4R 1 : Using hypothesis 2, λ ε L ∞ D\D 1 4R 1 ≤ C(C 0 , R 1 ). ( 75
)
One might also notice that, thanks to (65), on D\D 1 4R 1

:

ln ε d ε + r d ε ≤ C(C 0 , R 1 ) ln 1 2R 1
m+ 1 100 . ( 76
)
Then using (75) and (76):

l ε L ∞ D\D 1 4R 1 ≤ λ ε L ∞ D\D 1 4R 1 + ln ε d ε + r d ε L ∞ D\D 1 4R 1 ≤ C(C 0 , R 1 ). ( 77
)
On D 4εR1 : Using hypothesis 3

λ ε L ∞ (D4R 1 ) ≤ C(C 0 , R 1 ), (78)
while thanks to (73)

ln ε d ε + r d ε -ln C ε ε L ∞ (D4εR 1 ) ≤ ln ε d ε + (εr) d ε -ln C ε ε L ∞ (D4R 1 ) ≤ ln 1 + r d ε L ∞ (D4R 1 ) + ln C ε ε d ε +1 L ∞ (D4R 1 ) ≤ C(C 0 , R 1 ). (79)
To conclude we deduce thanks to (78) and (79):

l ε L ∞ (D4εR 1 ) ≤ λ ε -ln ε d ε + r d ε + ln C ε ε -ln C ε ε L ∞ (D4εR 1 ) ≤ λ ε -ln C ε ε L ∞ (D4εR 1 ) + ln ε d ε + r d ε -ln C ε ε L ∞ (D4εR 1 ) ≤ λ ε L ∞ (D4R 1 ) + ln ε d ε + r d ε -ln C ε ε L ∞ (D4εR 1 ) ≤ C(C 0 , R 1 ). (80) On D 1 2R 1 \D 2εR1 :
Thanks to (69):

l ε L ∞ D 1 2R 1 \D 2εR 1 ≤ λ ε -d ε ln r L ∞ D 1 2R 1 \D 2εR 1 + ln r d ε ε d ε + r d ε L ∞ D 1 2R 1 \D 2εR 1 ≤ C(C 0 , R 1 ). (81)
Combining (77), ( 80) and (81) yields

l ε L ∞ (D) ≤ C(C 0 , R 1 ), (82)
which is as desired. We now wish to refine this first expansion by showing that d ε converges toward θ fast enough to be updated in (69).

Step 5: Refinement: A consequence of estimate (82) is the following Harnack inequality on the conformal factor:

ε d ε + r d ε C(C 0 , R 1 ) ≤ e λ ε ≤ C(C 0 , R 1 ) ε d ε + r d ε
which, using the notation χ = √ ε 2 + r 2 , we will rewrite in the more convenient form

χ d ε C(C 0 , R 1 ) ≤ e λ ε ≤ C(C 0 , R 1 )χ d ε . ( 83
)
Injecting (65) into (83) yields

e λ ε ≤ C(C 0 , R 1 )χ θ-10 -3 . ( 84
) Since Φ ε is conformal, ∆ Φ ε = 2H ε e 2λ ε n ε = χ θ-10 -3 2H ε e λ ε e λ ε χ θ-10 -3 n ε . ( 85
)
Noticing that (35) and(84) imply

2H ε e λ ε e λ ε χ θ-10 -3 n ε L ∞ (D) ≤ C(C 0 , R 1 ), ( 86
)
we can apply theorem A.3 (see in the appendix) to equation ( 85) and find: where

Φ ε z = P ε (z) + ϕ ε 0 , (87)
P ε = θ q=0 p ε q z q ∈ C θ [X] with p ε q ≤ C(C 0 , R 1 ) for all q ≤ θ and ϕ ε 0 : D → R 3 satisfies ∀κ > 0 ϕ ε 0 χ θ+1-10 -3 -κ L ∞ (D) ≤ C κ (C 0 , R 1 ) , ∀p < ∞ ∇ ϕ ε 0 χ θ-10 -3 L p (D) ≤ C p (C 0 , R 1 ) . ( 88
)
By convergence of Φ ε away from zero (hypothesis 2), p ε q → p q ∈ C as ε goes to 0. Furthermore, (88) yields

ϕ ε 0 → ϕ 0 W 1,p (D), with ϕ 0 satisfying ∀κ > 0 ϕ 0 r θ+1-10 -3 -κ L ∞ (D) ≤ C κ (C 0 , R 1 ) , ∀p < ∞ ∇ ϕ 0 r θ-10 -3 L p (D) ≤ C p (C 0 , R 1 ) , (89) since χ → r as ε → 0. Then (87) ensures that Φ ε z → m q=0 p q z q + ϕ 0 . ( 90
)
Since we assumed that Φ ε → Φ 0 away from 0, comparing (32) and (90) yields

∀q < θ p ε q → 0 p ε θ → A = 0. ( 91
)
Besides, (87) gives the following

Φ ε z = θ q=0 p ε q ε q+1 C ε z q + ε ϕ ε 0 (ε.) C ε . ( 92
)
One might also notice using (88)

| ϕ ε 0 | (εz) ≤ C(C 0 , R 1 )χ θ+ 1 2 (εz) ≤ ε θ+ 1 2 C(C 0 , R 1 ) 1 + r 2 θ+ 1 2 . ( 93
)
This, along with (73), implies

ε ϕ ε 0 (ε.) C ε ≤ ε θ+1+ 1 2 C ε C(C 0 , R 1 ) 1 + r 2 θ+ 1 2 ≤ C(C 0 , R 1 ) 1 + r 2 θ+ 1 2 ε 1 2 . Consequently ε ϕ ε 0 (εz) C ε → 0 L ∞ loc (C) . ( 94
)
Since Φ ε is assumed to converge smoothly towards Φ 1 on compacts of C, we deduce from (51), ( 92) and ( 94)

θ q=0 p ε q ε q+1 C ε z q → Φ 1 z = Az θ + O(z θ-1 ). ( 95
) Hence ε θ+1 C ε p ε θ → A = 0.
Further, given that p ε θ → A = 0, there exists

C(C 0 , R 1 ) > 0 such that ln C ε ε θ+1 ≤ C(C 0 , R 1 ). ( 96
)
Combining ( 73) and (96) yields

ln ε d ε ε θ ≤ C(C 0 , R 1 ), which ensures |(d ε -θ) ln ε| ≤ C(C 0 , R 1 ). ( 97
)
Then, ( 69) and (97) combine and yield

λ ε -θ ln r L ∞ D 1 2R 1 \D 2εR 1 ≤ λ ε -θ ln r L ∞ D 1 2R 1 \D 2εR 1 + (θ -d ε ) ln r L ∞ D 1 2R 1 \D 2εR 1 ≤ C(C 0 , R 1 ). ( 98
)
Since inequality ( 98) is analogous to (69), we can do all the reasonings from ( 69) to (95) with d ε = θ. Then the conformal factor satisfies:

λ ε = ln ε θ + r θ + l ε , (99)
with

l ε such that l ε L ∞ (D) ≤ C(C 0 , R 1 ), and χ θ C(C 0 , R 1 ) ≤ e λ ε ≤ C(C 0 , R 1 )χ θ . ( 100
)
This concludes the proof of the desired result since R 1 is fixed. Moreover, for simplicity's sake we can, up to an inconsequential (thanks to (96)) adjustment, assume C ε = ε θ+1 . Then, exploiting (95) yields:

A = A. ( 101
) θ q=0 p ε q ε q-θ z q → Φ 1 z .
We can then decompose

P ε = ε θ Φ 1 z z ε + ε θ Q ε z ε , with Q ε ∈ C θ [X] such that Q ε → 0.
Φ ε then satisfies the following decomposition:

Φ ε z = ε θ Φ 1 z z ε + ε θ Q z ε + ϕ ε 0 , (102)
where

∀κ > 0 ϕ ε 0 χ θ+1-κ L ∞ (D) ≤ C κ (C 0 , R 1 ) , ∀p < ∞ ∇ ϕ ε 0 χ θ L p (D) ≤ C p (C 0 , R 1 ) . ( 103
)
Remark 4.2.1. Equality (101) can be seen as prolonging theorem 3.2: not only must the multiplicity of the end and the multiplicity of the branch point correspond, but so must the parametrization of the limit planes in both cases.

Remark 4.2.2. A. Michelat and T. Rivière have presented the author with another proof of the expansion which works in the more general framework of any simple bubble (in [START_REF] Michelat | [END_REF]).

5 Conditions on the limit surface:

The aim of this section is to prove both theorem 1.5 and corollary 1.6.

5.1 Proof of theorem 1.5:

As explained in section 3 we can equivalently work in conformal parametrizations under hypotheses 1-7. We will then instead prove:

Theorem 5.1. Let Φ ε : D → R 3 be a sequence of Willmore conformal immersions satisfying hypotheses 1-7.

Then the second residue of Φ 0 at 0 satisfies α ≤ θ -1.

Proof .

Step 1: Expansion of Φ ε zz : We consider Φ ε satisfying hypotheses 1-7, and thus (99)-(103). The system [START_REF] Bryant | A duality theorem for Willmore surfaces[END_REF] of [START_REF] Marque | An ε-regularity result with mean curvature control for Willmore immersions and application to minimal bubbling[END_REF] states

         ∆S ε = H ε ∇ Φ ε , ∇ ⊥ R ε ∆ R ε = -H ε ∇ Φ ε × ∇ ⊥ R ε -∇ ⊥ S ε H ε ∇ Φ ε ∆ Φ ε = 1 2 ∇ ⊥ S ε .∇ Φ ε + ∇ ⊥ R ε × ∇ Φ ε . ( 104
)
Then ( 35), ( 37) and (104) yield:

∆S ε L ∞ (D) + ∆ R ε L ∞ (D) ≤ C(C 0 ). ( 105
)
Applying theorem A.3 (in the appendix) gives the following decomposition on S ε and R ε :

S ε z = S ε z (0) + σ ε 0 R ε z = R ε z (0) + ρ ε 0 , (106)
with σ ε 0 and

ρ ε 0 satisfying ∀κ > 0 |σ ε 0 | (z) + | ρ ε 0 | (z) ≤ C κ (C 0 )χ 1-κ , ∀p < ∞ ∇σ ε 0 L p (D) + ∇ ρ ε 0 L p (D) ≤ C p (C 0 ). (107)
Injecting ( 102) and (106) into the third equation of (104) yields

∆ Φ ε = 2 R ε z (0) × ε θ Φ 1 z z ε + ε θ Q ε z ε + S ε z (0) ε θ Φ 1 z z ε + ε θ Q ε z ε + Ψ ε 0 , (108)
where

Ψ ε 0 := 2 ρ ε 0 × Φ ε z + σ ε 0 Φ ε z + R ε z × ϕ ε 0 + S ε z ϕ ε 0 satisfies ∀κ > 0 Ψ ε 0 (z) ≤ C κ (C 0 )χ θ+1-κ , ∀p < ∞ ∇ Ψ ε 0 χ θ L p (D) ≤ C p (C 0 ). ( 109
)
One may notice that

h ε := 2 R ε z (0) × ε θ Φ 1 z z ε + ε θ Q ε z ε + S ε z (0) ε θ Φ 1 z z ε + ε θ Q ε z ε
is the sum of a polynomial of degree θ in z and a polynomial of degree θ in z, whose coefficients are uniformly bounded by a constant depending on C 0 . Additionally it is of order O(χ θ ) thanks to theorem 4.1. We can then find a polynomial h ε in z and z of total degree θ + 2 such that

h ε (0) = h ε z (0) = h ε z (0) = 0, h ε = O χ θ+2 ∆ h ε = h ε . Then ∆ Φ ε -h ε = Ψ ε 0 . ( 110
)
Minimal bubbling for Willmore surfaces 23

Applying theorem A.8 (in the appendix) to (110), with a = θ + 1 -κ for κ arbitrarily small yields

Φ ε z = P ε (z) + h ε z + ϕ ε 1 , (111)
where P ε is a polynomial of degree θ + 1 that we can split P ε = P ε θ + p ε z θ+1 with P ε θ ∈ C θ [X], and ϕ ε 1 satisfies:

∀κ > 0 | ϕ ε 1 | χ θ+2-κ + |∇ ϕ ε 1 | χ θ+1-κ ≤ C κ (C 0 ), ∀p < ∞ ∇ 2 ϕ ε 1 χ θ L p (D) ≤ C p (C 0 ). ( 112
)
Comparing ( 102) and (111) as in the proof of lemma A.5 (in the appendix) yields:

P ε θ = ε θ Φ 1 z ε + Q ε z ε , ϕ ε 0 = p ε z θ+1 + h ε z + ϕ ε 1 .
Consequently ϕ ε 0 satisfies:

| ϕ ε 0 | χ θ+1 + |∇ ϕ ε 0 | χ θ ≤ C(C 0 ) ∀p < ∞ ∇ 2 ϕ ε 0 χ θ-1 L p (D) ≤ C p (C 0 ). (113)
Estimates (113) applied to (102) allow for a pointwise expansion of Φ ε zz :

Φ ε z = ε θ Φ 1 z z ε + Q ε z ε + ϕ ε 0 , Φ ε zz = ε θ-1 Φ 1 zz z ε + Q ε z z ε + (ϕ ε 0 ) z . (114)
Step 2: Initial conditions The relations (114) yield when evaluated at 0

Φ ε z (0) = ε θ Φ 1 z (0) + ε m Q ε (0) + O ε m+1 , Φ ε zz (0) = ε θ-1 Φ 1 zz (0) + ε θ-1 Q ε z (0) + O (ε m ) . (115)
Per hypothesis 7, we have:

Φ ε (0) = 0, Φ ε z (0) = ε θ Φ 1 z (0) = ε θ 2 e 0 , Ω ε e -λ ε (0) = 1 ε θ Ω ε (0) = - 2 ε .
This implies n ε (0) = -e 3 , and since

Ω ε 2 = n ε , Φ ε zz , Φ ε z (0) = ε θ 2 e 0 and Φ ε zz (0), e 3 = ε θ-1 . (116)
Comparing ( 115) and (116) yields

Q ε (0) = O(ε) and Q ε z (0), e 3 = O(ε).
However as we pointed out in remark A.4.1, Q ε is loosely defined. We can then evacuate the coefficients of order ε into ϕ ε 0 (which we will do without changing the notations) to obtain:

Q ε (0) = 0 and Q ε z (0), e 3 = 0. ( 117
)
To conclude we write

Q ε ∈ C 3 as Q ε := A ε e 0 + B ε e 0 + C ε e 3 , then (117) yields A ε (0) = B ε (0) = C ε (0) = 0, C ε z (0) = 0. (118)
Evaluating equality (170) (see in the appendix) at z = 0, one has:

R ε z (0) = ε θ (H ε (0) + iV ε (0)) e 0 + iS ε z (0) e 3 . (119)
Estimate (37) then ensures that κ ε := ε θ (H ε (0) + iV ε (0)) is uniformly bounded:

|κ ε | ≤ C(C 0 ).
Step 3: Φ ε is conformal We will linearize the conformality condition:

Φ ε z , Φ ε z = 0.
Injecting (102) in the former yields

Q 2 B ε -P 2 A ε + P QC ε + A ε B ε + (C ε ) 2 2 = 0. ( 120
)
Applying hypothesis 6 and (118) then yields:

z 2 divides B ε . (121)
Conformality also implies ∆ Φ ε , Φ ε z = 0.

Injecting ( 108) and (114) into the former then yields

h ε , ε θ Φ 1 z + Q ε z ε = h ε , ϕ ε 0 + ∆ Φ ε , Ψ ε 0 =: Ψ ε 1 , (122)
with Ψ ε 1 satisfying thanks to (109) and ( 113)

∀κ > 0 Ψ ε 1 ≤ C κ χ 2θ+1-κ . ( 123
) Considering that h ε , ε θ Φ 1 z + Q ε z
ε is a polynomial of degree at most 2θ in z and z, we can state:

h ε , ε θ Φ 1 z + Q ε z ε = 2θ p+q=0 h ε pq ε 2θ-p-q z p z q .
Together (122) and (123) yield:

∀κ > 0 2θ p+q=0 h ε pq ε 2θ-p-q z p z q = O χ 2θ+1-κ .
Applying lemma A.7 then yields:

∀p, q ∀κ > 0 h ε pq = O ε 1-κ . ( 124
)
Step 4: Computing h ε We compute

ε θ Φ 1 z z ε + ε θ Q ε z ε = ε θ - P 2 z ε 2 + B ε z ε e 0 + ε θ Q 2 z ε 2 + A ε z ε e 0 + ε θ P z ε Q z ε + C ε z ε e 3 .
Minimal bubbling for Willmore surfaces 25

Hence

R ε z (0) × ε θ Φ 1 z + Q ε z ε = -2iκ ε ε θ Q 2 2 + A ε z ε e 3 + iκ ε ε θ [P Q + C ε ] z ε e 0 + S ε z (0)ε θ - P 2 2 + B ε z ε e 0 -S ε z (0)ε θ Q 2 2 + A ε z ε e 0 , and
S ε z (0)ε θ Φ 1 z + Q ε z ε = S ε z (0)ε θ - P 2 2 + B ε z ε e 0 + S ε z (0)ε θ Q 2 2 + A ε z ε e 0 + S ε z (0)ε θ [P Q + C ε ] z ε e 3 .
Then

h ε = 2 S ε z (0)ε θ [-P 2 + 2B ε ] z ε + iκ ε ε θ [P Q + C ε ] z ε e 0 + S ε z (0)ε θ [P Q + C ε ] z ε -iκ ε ε θ [Q 2 + 2A ε ] z ε e 3 .
From this we deduce

h ε , ε θ Φ 1 z + Q ε z ε = S ε z (0)ε θ [-P 2 + 2B ε ] z ε + iκ ε ε θ θ[P Q + C ε ] z ε ε θ -P 2 + 2B ε z ε + S ε z (0)ε θ [P Q + C ε ] z ε -iκ ε ε θ [Q 2 + 2A ε ] z ε ε θ [P Q + C ε ] z ε = S ε z (0)ε 2θ |P | 2 |P | 2 + |Q| 2 + 2 P QC ε -2P 2 B ε + 4 |B ε | 2 + |C ε | 2 z ε + iκ ε ε 2θ -P Q |P | 2 + |Q| 2 -P 2 C ε + 2B ε P Q -2P QA ε -C ε Q 2 + 2B ε C ε -2C ε A ε z ε . (125)
Studying ( 125) with (118), (121), (122), (123) and hypothesis 6 in mind, we can write

h ε , ε θ Φ 1 z + Q ε z ε = -iκ ε ε 2θ-1 z + O(r 2 ), (126)
which implies h ε 1,0 = -iκ ε , and in turn thanks to (124):

∀s > 0 κ ε = O ε 1-s . (127)
Then, (122), ( 123) and (125) give us:

h ε , ε θ Φ 1 z + Q ε z ε = S ε z (0)ε 2θ |P | 2 |P | 2 + |Q| 2 + 2 P QC ε -2P 2 B ε + 4 |B ε | 2 + |C ε | 2 z ε + O χ 2θ+1-κ . (128)
A similar process on the remaining polynomial allows us to state

∀κ > 0 S ε z (0) = O ε 1-κ . ( 129
)
Step 5: Conclusion From (106), ( 127) and (128) we deduce:

∀κ > 0 R ε z + |S ε z | ≤ C κ χ 1-κ , ∀p < ∞ ∇ R ε z L p (D) + ∇S ε z L p (D) ≤ C p . (130)
Inequality (119) from [START_REF] Marque | An ε-regularity result with mean curvature control for Willmore immersions and application to minimal bubbling[END_REF] then yields:

∀κ > 0 H ε e λ ε ≤ C κ χ 1-κ . (131)
Letting ε tend to 0 in (131) gives, for any r > 0, thanks to hypothesis 2, the following:

∀κ > 0 H 0 e λ 0 ≤ C κ r 1-κ .
However since Φ 0 is assumed to have a branch point of order θ + 1 at 0, by definition, e λ 0 ∼ Cr θ , which means

∀κ > 0 H 0 ≤ C κ r 1-θ-κ . (132)
By definition of α (see (33)), H 0 r -α . Since α ∈ Z, (132) ensures:

α ≤ θ -1. (133)
This concludes the proof of the desired result.

In the continuity of the previous proof we can improve on a convergence result obtained in [START_REF] Marque | An ε-regularity result with mean curvature control for Willmore immersions and application to minimal bubbling[END_REF]:

Theorem 5.2. Let Φ k : Σ → R 3 be a sequence of Willmore immersions satisfying the hypotheses of theorem 1.1. Assume further that at each concentration point a simple minimal bubble is blown. Then Φ k → Φ 0 C 3,η for all η < 1.

Proof . As before we can reason locally, under hypotheses 1-7, and will continue from (133). Injecting (130) into (108) ensures:

∀κ > 0 ∆ Φ ε ≤ C κ χ θ+1-κ ∀p < ∞ ∇ ∆ Φ ε χ θ L p (D) ≤ C p . (134)
We can then compute

H ε Φ ε z = Φ ε z z × Φ ε z i Φ ε z 2 , ∇ H ε Φ ε z = ∇ Φ ε z z × Φ ε z i Φ ε z 2 + Φ ε z z × ∇ Φ ε z i Φ ε z 2 -∇ Φ ε z , Φ z + Φ z , ∇ Φ z Φ ε z z × Φ ε z i Φ ε z 4 . (135)
Combining (100), ( 114) and (134) yields:

∀κ > 0 H ε Φ ε z χ 1-κ L ∞ (D) ≤ C κ , ∀p < ∞ ∇ H ε Φ ε z L p (D) ≤ C p . (136)
Consequently, injecting (136) into (104) and applying Calderon-Zygmund yields

∀p < ∞ ∇S ε W 2,p (D) + ∇ R ε W 2,p (D) + ∇ Φ ε W 3,p (D) ≤ C(C 0 ). ( 137
)
Which proves theorem 5.2 thanks to classical embeddings.

Remark 5.2.1. We can extend our expansions to the next order. Indeed injecting (136) and ( 130) into (104) yields

∀κ > 0 ∆S ε χ 2-κ L ∞ (D) + ∆ R ε χ 2-κ L ∞ (D) ≤ C κ , ∀p < ∞ ∆∇S ε χ L p (D) + ∆∇ R ε χ L p (D) ≤ C p .
Minimal bubbling for Willmore surfaces 27

Applying corollary A.8 (see in the appendix) then yields

S ε z = S ε z (0) + s ε 1 z + s ε 2 z 2 + σ ε 1 , R ε z = R ε z (0) + r ε 1 z + r ε 2 z 2 + ρ ε 1 , (138)
where the s ε j and the r ε j are uniformly bounded constants and σ ε 1 , ρ ε 1 satisfy:

∀κ > 0 σ ε 1 χ 3-κ + ∇σ ε 1 χ 2-κ + ρ ε 1 χ 3-κ + ∇ ρ ε 1 χ 2-κ ≤ C κ , ∀p < ∞ ∇ 2 σ ε 1 χ L p (D) + ∇ 2 ρ ε 1 χ L p (D) ≤ C p . (139)
Setting

σ ε 0 = s ε 1 z + s ε 2 z 2 + σ ε 1 and ρ ε 0 = r ε 1 z + r ε 2 z 2 + ρ ε 1 yields σ ε 0 χ + |∇σ ε 0 | + ρ ε 0 χ + |∇ ρ ε 0 | ≤ C, ∀p < ∞ ∇ 2 σ ε 0 L p (D) + ∇ 2 ρ ε 0 L p (D) ≤ C p . (140)
We can then do all the reasonings from (108) to (136) for better controls:

H ε Φ ε z χ + S ε z χ + R ε z χ ≤ C. (141)
Injecting this added regularity into the third equation of (104) ensures:

∆ Φ ε χ 3 + ∆∇ Φ χ 2 ≤ C ∀p < ∞ ∆∇ 2 Φ χ L p (D) ≤ C p . (142)
With another application of corollary A.8 (see in the appendix) we can expand Φ ε z in the following manner:

Φ ε z = ε θ Φ 1 z z ε + Q ε z ε + ϕ ε 0 , (143)
with

|ϕ ε 0 | χ θ+1 + |∇ϕ ε 0 | χ θ + ∇ 2 ϕ ε 0 χ θ-1 ≤ C(C 0 ) ∀p < ∞ ∇ 3 ϕ ε 0 χ θ-2 L p (D) ≤ C p (C 0 ). (144)
5.2 Proof of corollary 1.6:

Proof . We only need to show that the inverse of a Chen-Gackstatter torus has a branch point of second residue α = 2. Let Ψ : (C\Z 2 )/Z 2 → R 3 be a parametrization of the Chen-Gackstatter torus, p ∈ R 3 such that d(p, Ψ) > 1, and Φ = ι • ( Ψ -p), the studied inverse. Let b denote its branch point.

It is interesting to notice that the second residue of Φ can be read on its conformal Gauss map Y defined as:

Y Φ = H Φ    Φ | Φ| 2 -1 2 | Φ| 2 +1 2    +     n Φ n Φ , Φ n Φ , Φ     .
Indeed Φ and n Φ are bounded around the branch point, and thus necessarily

Y Φ ∼ b C Φ z -α .
Further, since Φ and Ψ differ by a conformal transform, it has been shown in [START_REF] Marque | Conformal gauss map geometry and application to Willmore surfaces in model spaces[END_REF], that there exists a fixed matrix M ∈ SO(4, 1) such that Y Ψ = M Y Φ . This yields that necessarily Y Ψ ∼ b C Ψ z -α . Hence, considering that Ψ is minimal, we deduce that

Y Ψ =     n Ψ n Ψ , Ψ n Ψ , Ψ     . Since n Ψ is bounded, n Ψ , Ψ ∼ b Cz -α . ( 145
)
We will now use the Enneper-Weierstrass parametrization of Ψ and (145) to compute the second residue of Ψ at its branch point. Chen-Gackstatter is a minimal surface of genus 1 and of Enneper-Weierstrass data centered on the branch point: (f, g) = 2ϕ(z), A ϕz ϕ (z) (see [START_REF] Chen | Elliptische und hyperelliptische Funktionen und vollständige Minimalflächen vom Enneperschen Typ[END_REF]) where ϕ is the Weierstrass elliptic function, of elliptic invariants

g 2 = 60 ∞ m,n=-∞ 1 (m + ni)
4 > 0,

g 3 = 0,
and

A = 3π 2g 2 ∈ R + .
Then, ϕ has the following expansion around 0 (see [START_REF] Apostol | Modular functions and Dirichlet series in number theory[END_REF]):

ϕ(z) = 1 z 2 + O(z 2 ) ϕ z (z) = -2 z 3 + O(z).
Hence we can state that

Ψ = 2 1 z 2 e 0 - 4A
Similarly:

Ψ z × Ψ z = 1 z 2 e 0 - 4A 2 z 4 e 0 - 4A z 3 e 3 + O(1) × 1 z2 e 0 - 4A 2 z4 e 0 - 4A z3 e 3 + O(1) = - 4Ai z 2 z3 e 0 + 32iA 4 r 8 e 3 - 16A 3 i z 4 z3 e 0 - 4Ai z 3 z2 e 0 - 16A 3 i z 3 z4 e 0 + O 1 r 4 = 32iA 4 r 8 - z 2A + z 2 z 8A 3 e 0 - z 2A + z2 z 8A 3 e 0 + e 3 + O(r 4 ) , and
Ψ z 2 = 32A 4 r 8 + 16A 2 r 6 + O(r -4 ) = 32A 4 r 8 1 + r 2 2A 2 + O(r 4 ) , which yields n Ψ = 1 - r 2 2A 2 e 3 - z 2A - z 2 z 8A 3 e 0 - z 2A - z z2 8A 3 e 0 + O(r 4 ). (147)
Combining ( 146) and (147) ensures:

n Ψ , Ψ = 2A 1 z 2 + 1 z2 - 4A 3z 2 - 4A 3z 2 + O 1 r = 2A 3z 2 + 2A 3z 2 + O 1 r . (148)
Considering (148) in light of (145) yields α = 2. Applying theorem 1.5 concludes the proof.

Minimal bubbling for Willmore surfaces 29 5.3 An exploration of the consequences of theorem 1.5

We consider Φ : D → R 3 a true Willmore conformal branched immersion with a single branch point at 0, of multiplicity θ + 1, and of second residue α ≤ θ -1. Then applying theorem 1.4 we can expand Φ around 0 in the following way:

Φ(z) = 2 1 2(θ + 1) e 0 z θ+1 + θ+1-α j=1 A j θ + 1 + j z θ+1+j + Cz θ+1-α zθ+1 + Cz θ+1 zθ+1-α (θ + 1 -α)(θ + 1) e 3 + ξ, (149)
where ξ satisfies ∇ j ξ = O |z| 2(θ+1 -α -j + 1 -ν for all κ > 0 and j ≤ θ + 2 -α.

Moreover, if we do the conformal change of variables Z θ+1 = z θ+1 + Az θ+2 , (149) becomes

Φ(Z) = 2 1 2(θ + 1) e 0 Z θ+1 + A 1 θ + 2 + A e 0 Z θ+2 + CZ θ+1-α Z θ+1 + CZ θ+1 Z θ+1-α (θ + 1 -α)(θ + 1) e 3 + O |Z| θ+3 .
Thus up to doing a conformal change of charts we can assume that A 1 has no component along e 0 , meaning:

A 1 = V 2 e 0 + U e 3 . (150)
Then

Φ z = 1 2 e 0 z θ + θ+1-α j=1 A j z θ+j + C θ + 1 z θ-α zθ+1 + C θ + 1 -α z θ zθ+1-α e 3 + ξ z , (151)
Using Φ conformal, we can expand Φ z , Φ z to the order z 2θ+1 and conclude that

V = A 1 , e 0 = 0.
Then, wishing to expand the Gauss map n we compute:

Φ z 2 = r 2θ 2 + O r 2θ+2 , Φ z × Φ z = 1 2 e 0 + U z θ+1 e 3 + O(r θ+2 ) × 1 2 e 0 + U zθ+1 e 3 + O(r θ+2 ) = ir 2θ 2 e 3 - iU z θ+1 zθ 2 e 0 - iU z θ zθ+1 2 e 0 + O(r 2θ+2 ).
Hence we can write n = e 3 -U z e 0 -U z e 0 + O(r 2 ).

One can differentiate (151), and obtain

Φ zz = θ 2 e 0 z θ-1 + θ+1-α j=1 (θ + j) A j z θ-1+j + (θ -α)C θ + 1 z θ-1-α zθ+1 + θC θ + 1 -α z θ-1 zθ+1-α e 3 + ξ zz .
Taking the scalar product with (152) we find:

Ω 2 = (θ + 1)U z θ -θU z θ + O r θ+1 = U z θ + O(r θ+1 ).
In theorem 4.11 of [START_REF] Michelat | The classification of branched Willmore spheres in the 3-sphere and the 4-sphere[END_REF], A. Michelat and T. Rivière have shown that

A 1 , C e 3 = U C = 0. ( 153
)
This precise equality is found at (4.53), (4.54) and (4.69) (depending on whether θ ≥ 4, θ = 3 or θ = 2) of the aforementioned article and stems from the conformal relation Φ z , Φ z = 0, expanded to the order 3θ + 2 -α (which requires furthering expansion (151) in the way detailed in [START_REF] Michelat | The classification of branched Willmore spheres in the 3-sphere and the 4-sphere[END_REF]) in order to consider the first terms in z. Equality (153) implies notably that either C = 0, and thus that α ≤ θ -2, of U = 0, and thus that |Ω| e -λ = O(r), meaning that 0 is an umbilic point. Besides, theorem 4.11 of [START_REF] Michelat | The classification of branched Willmore spheres in the 3-sphere and the 4-sphere[END_REF] states that Bryant's quartic (see in the introduction, or in [START_REF] Marque | Conformal gauss map geometry and application to Willmore surfaces in model spaces[END_REF]) is then holomorphic across the branched point. Theorem 1.5 has the following corollary:

Corollary 5.3. Let Φ k be a sequence of Willmore immersions of a closed surface Σ satisfying the hypotheses of theorem 1.1. Then at each concentration point p ∈ Σ of multiplicity θ p + 1 on which a simple minimal bubble is blown, either the second residue of Φ ∞ at p satisfies

α p ≤ θ p -2,
or p is an umbilic point for Φ ∞ (Σ).

In both cases Bryant's quartic Q is holomorphic across those branch points.

A Appendix

T |z| 1-2 p -κ L ∞ (D) ≤ C κ µ |z| a L ∞ (D) f L p (D) + u C 1 (∂D) . Additionally if µ ∈ C 1 (D\{0}), a = 0 and ∇µ(z) = O |z| a-1 Then: u zz (z) = P z + |z| a Q with Q ∈ L p (
D) for all p < p and

Q L p (D) ≤ C p µ |z| a L ∞ (D) + ∇µ |z| a-1 L ∞ (D) f L p (D) + u C 1 (∂D) .
In fact Q =

|z| 1-2 p -κ L ∞ (D) ≤ C κ µ |z| a L ∞ (D) f L p (D) + u C 1 (∂D) . Additionally if µ ∈ C 1 (D\{0}), a = 0 and ∇µ(z) = O |z| a-1 Then: u zz (z) = P z + |z| a Q with Q ∈ L p (D) for all p < p and Q L p (D) ≤ C p µ |z| a L ∞ (D) + ∇µ |z| a-1 L ∞ (D) f L p (D) + u C 1 (∂D) . In fact Q = (|z| a T (z)) z |z| a .
To deal with the concentration phenomena, we will need a variant of the last theorem.

Theorem A.3. Let (u ε ) ε>0 ∈ C 2 (D\{0}) solve ∆u ε (z) = χ a f ε (z) in D, with f ε ∈ L p (D) for 2 ≤ p ≤ ∞, a ∈ R and χ := √ ε 2 + r 2 . Then u ε z (z) = P ε (z) + χ a T ε (z) with P ε ∈ C a [X] and T ε = O(χ 1-2 p -κ ) for all κ > 0.
Here a is the upper integral part of a. More precisely one has

T ε χ 1-2 p -κ L ∞ (D) ≤ C κ f ε L p (D) + u ε C 1 (∂D) .
Additionally:

u ε zz (z) = P ε z + χ a Q ε with Q ε ∈ L p (D) for all p < p and Q ε L p (D) ≤ C p f ε L p (D) + u ε C 1 (∂D) .
In fact

Q ε = (χ a T ε (z)) z χ a .
Proof . We first give a small lemma which can help one understand the concrete impact of χ:

Lemma A.4. For all a,b ∈ R + , there exists a constant C a,b such that

ε a r b χ a+b ≤ C a,b .
Similarly for all a ∈ R + , there exists

C a ∈ R * + such that 1 C a ≤ ε a + r a χ a ≤ C a . (154)
Here C a (respectively C a,b ) depends solely on a (and b), and not on ε or r.

We then write ∆u

ε = (ε a + r a ) χ a ε a + r a f ε = (ε a + r a ) f ε , where f ε = χ a ε a +r a f ε satisfies, thanks to (154), f ε L p (D) ≤ C a f ε L p (D). (155)
We can then use Green's formula to write

u ε z (z) = 1 2π ∂D z - x |z -x| 2 ∂ ν u ε (x) -u ε (x)∂ z - x |z -x| 2 dσ(x) - 1 2π D z - x |z -x| 2 (ε a + r a ) f ε (x)dx = 1 2π ∂D z - x |z -x| 2 ∂ ν u ε (x) -u ε (x)∂ z - x |z -x| 2 dσ(x) - 1 2π D z - x |z -x| 2 ε a f ε (x)dx - 1 2π D z - x |z -x| 2 r a f ε (x)dx = I ε 0 (z) + I ε 1 (z) + I ε 2 (z). (156)
We can then successively estimate the three terms as in the proof of theorem A.2 and write

I ε 0 (z) = P ε 0 (z) + z a +1 T ε 0 , (157)
where P ε 0 is a polynomial of degree at most a and whose coefficients are bounded by u ε C 1 (∂D) , and

T ε 0 L ∞ (D) ≤ C u ε C 1 (D) .
Working as in [START_REF] Bernard | Singularity removability at branch points for Willmore surfaces[END_REF] we write

I ε 1 (z) = C ε + ε a T ε 1 (z) I ε 2 (z) = P ε 2 + r a T ε 2 (z) (158)
where C ε is a constant and P ε 2 a polynomial of degree at most a , both bounded by u ε C 1 (∂D) , while

T ε 1 r 1-2 p -κ L ∞ (D) ≤ C κ f ε L p (D) + u ε C 1 (D) T ε 2 r 1-2 p -κ L ∞ (D) ≤ C κ f ε L p (D) + u ε C 1 (D) .
In the end, combining (157) and (158) yields

u ε z = P ε + ε a T ε 1 + r a T ε 3 , (159)
where P ε is a polynomial of degree at most a , T ε 1 is as previously stated and T ε 3 still satisfies

T ε 3 r 1-2 p -κ L ∞ (D) ≤ C κ f ε L p (D) + u ε C 1 (D) .
Minimal bubbling for Willmore surfaces 33

Proceeding similarly then ensures that

u ε zz = P ε z + ε a Q ε 1 + r a Q ε 2 , (160)
where

Q ε 1 = T ε 1, z and Q ε 2 = (r a T ε 3 ) z r a satisfy for all p < p Q ε 1 L p (D) ≤ C p f ε L p (D) + u ε C 1 (∂D) , Q ε 3 L p (D) ≤ C p f ε L p (D) + u ε C 1 (∂D) .
Let us notice that the estimate on Q ε 1 is not stricto sensu derived from the proof of theorem A.2 but from similar classical Calderon-Zygmund estimates.

From (159) we write u ε z = P ε + χ a T ε with

T ε = ε a χ a T ε 1 + r a χ a T ε 3 which then satisfies T ε χ 1-2 p -κ L ∞ (D) ≤ ε a χ a L ∞ (D) T ε 1 χ 1-2 p -κ L ∞ (D) + r a χ a L ∞ (D) T ε 3 χ 1-2 p -κ L ∞ (D) ≤ ε a χ a L ∞ (D) T ε 1 r 1-2 p -κ L ∞ (D) r 1-2 p -κ χ 1-2 p -κ L ∞ (D) + r a χ a L ∞ (D) T ε 3 r 1-2 p -κ L ∞ (D) r 1-2 p -κ χ 1-2 p -κ L ∞ (D) ≤ C κ f ε L p (D) + u ε C 1 (D) , (161)
= P ε z + χ a Q ε with Q ε = ε a χ a Q ε 1 + r a χ a Q ε 3 = (χ a T ε ) z χ a
which then satisfies

Q ε L p (D) ≤ ε a χ a L ∞ (D) Q ε 1 L p (D) + r a χ a L ∞ (D) Q ε 3 L p (D) ≤ C p f ε L p (D) + u ε C 1 (D) , (162)
using lemma A.4.

Remark A.4.1. We must point out that the expansion given by theorem A.3 is by no means unique. Indeed if for instance u z = P ε + χ m T ε , then one could readily write

u z = P ε + ε m+1 + χ m T ε + ε m+1 χ m
with T ε + ε m+1 χ m still satisfying (161). Theorem A.3 can be applied several times to prove an increased regularity on the higher order terms: Lemma A.5. Let u ε ∈ C 2 (D\{0}) such that ∆u ε = χ a f ε , with f ε ∈ L ∞ and ∆ (∇u ε ) = χ a-1 g ε , with g ε ∈ L p . Then u ε z = P ε + µ ε , where P ε is a complex polynomial of degree at most a , and µ ε such that

|µ ε | χ a+1-κ + |∇µ ε | χ a-2 p -κ ≤ C κ f ε L ∞ (D) + g ε L p (D) + u ε C 2 (∂D) ,
and

∇ 2 µ ε χ a-1 L p (D) ≤ C p f ε L ∞ (D) + g ε L p (D) + u ε C 2 (∂D) .
Proof . We apply theorem A.3 twice and decompose u z and (u z ) z :

u ε z = P ε 1 + µ ε 1 (u ε z ) z = P ε 2 + µ ε 2 , (163)
where

|µ ε 1 | χ a+1-κ + |µ ε 2 | χ a-2 p -κ ≤ C κ f ε L ∞ (D) + g ε L p (D) + u ε C 2 (∂D) ∇µ ε 1 χ a L p 1 (D) + ∇µ ε 2 χ a-1 L p 2 (D) ≤ C p 1 ,p 2 f ε L ∞ (D) + g ε L p (D) + u ε C 2 (∂D) , (164)
for all p 1 < ∞ and p 2 < p. We then enjoy two expressions for u zz :

u zz = P ε 1, z + µ ε 1, z = P ε 2 + µ ε 2 . Consequently P ε 1, z -P ε 2 = µ ε 2 -µ ε 1, z ,
which in turn, combined with (164), implies that

D P ε 1, z -P ε 2 χ a s dz ≤ C s f ε L ∞ (D) + g ε L p (D) + u ε C 2 (∂D) ,
for all s < ∞. We decompose

P ε 1, z -P ε 2 = a q=0
p ε q z q , and can state for a given R 0 > 0

D R 0 a q=0 p ε q z q χ a s dz ≤ D a q=0 p ε q z q χ a p ≤ C s f ε L ∞ (D) + g ε L p (D) + u ε C 2 (∂D) .
Changing variables yields

D R 0 a q=0 p ε q ε a-q-2 p z q √ 1 + r 2 a s dz ≤ C s f ε L ∞ (D) + g ε L p (D) + u ε C 2 (∂D) .
And since on D R0 , 1 1+r 2 ≥ 1 1+R 2 0 , we deduce

D R 0 a q=0 p ε q ε a-q-2 p z q s dz ≤ C s,R0 f ε L ∞ (D) + g ε L p (D) + u ε C 2 (∂D) . ( 165
)
Minimal bubbling for Willmore surfaces 35

It is now important to notice that the left-hand term in (165) is in fact a polynomial in R 0 , which is uniformly bounded in ε on compacts of C. All its coefficients are thus uniformly bounded in ε, and straightforward computations then yield:

∀s < ∞ ∀j ≤ a ∀ε > 0 p ε q ε a-q-2 s ≤ C s f ε L ∞ (D) + g ε L p (D) + u ε C 2 (∂D)
which thanks to lemma A.4 translates on P ε 1, z -P ε 2 as

∀s < ∞ P ε 1, z -P ε 2 χ a-2 s ≤ C s f ε L ∞ (D) + g ε L p (D) + u ε C 2 (∂D) , (166)
and

∀s < ∞ P ε 1, z -P ε 2 z χ a-2 s -1 ≤ C s f ε L ∞ (D) + g ε L p (D) + u ε C 2 (∂D) . (167)
Now since µ ε 1, z = µ ε 2 -P ε 1, z -P ε 2 we can combine (164) and (166) to find for all κ > 0

µ ε 1, z χ a-2 p -κ ≤ C κ f ε L ∞ (D) + g ε L p (D) + u ε C 2 (∂D) . (168)
In addition, since µ ε 1, zz = µ ε 2, z -P ε 1, z -P ε 2 z , (164) and ( 167) yield for all p < p:

µ ε 1, zz χ a-1 L p (D) ≤ C p f ε L ∞ (D) + g ε L p (D) + u ε C 2 (∂D) . (169)
Applying similarly theorem A.3 to u ε z yields controls akin to (168) and (169) on the missing terms in the gradient and the Hessian, which concludes the proof.

A cautious reader might have noticed that we have in fact proved the following lemma: We will also use a corresponding result for polynomials in z and z:

Lemma A.7. Let u ∈ N, v ≥ u and P ε = u i+j=0 p ε i,j z i zj such that

∀p < ∞ P ε χ v ≤ C.
Then

∀ν > 0 ∀i + j ≤ u p ε i,j ε v-i-j ≤ C ν .
Applying lemma A.5 several times yields:

Corollary A.8. Let u ε ∈ C 2 (D\{0}) such that, for a ≥ t ∆u ε = χ a f ε 0 , ∆∇u ε = χ a-1 f ε 1 . . . ∆∇ t u ε = χ a-t f ε t 36 N. Marque with f ε j ∈ L ∞ (D) for j ≤ t -1 and f ε t ∈ L p (D). Then

u ε z = P ε + µ ε ,
where P ε is a complex polynomial of degree at most a , and µ ε such that

|µ ε | χ a+1-κ + |∇µ ε | χ a-κ + • • • + |∇ t µ ε | χ a+1-t-2 p -κ ≤ C κ t q=0 f ε q L ∞ (D) + u ε C t+1 (∂D) ,
and

∇ t+1 µ ε χ a-t L p (D) ≤ C p t q=0 f ε q L ∞ (D) + u ε C t+1 (∂D) .
Proof . The proof is a recurrence whose initialization is theorem A.3 and whose heredity is obtained by applying lemma is A.5 to ∇ s u ε .

A.2 Auxiliary formulas

In the following, given Φ a Willmore conformal immersion, and L, S, R defined in (34), we wish to prove:

R z = 2 (H + iV ) Φ z -iS z n, (170)
where V = 1 2 L, n . Indeed, since We first give a version of Gauss-Bonnet formula taking branch points and branched ends into account. We refer the reader to theorem 2.6 in [START_REF] Lamm | Branched Willmore spheres[END_REF]. Proposition A.9. Let Σ be a compact Riemann surface and Φ : Σ → R 3 ∪ {∞} be a branched immersion with a finite number of ends. Let p 1 , . . . , p n be its branch points of respective orders a 1 + 1, . . . , a n + 1 and q 1 , . . . , q m its ends of respective orders b 1 -1, . . . b m -1. We denote χ(Σ) the Euler characteristic of Σ, g the metric induced on Σ by Φ and K its Gauss curvature. Hence Following are a few useful formulas linking the different notions of curvature. The computations are done in appendix A.1 of [START_REF] Marque | An ε-regularity result with mean curvature control for Willmore immersions and application to minimal bubbling[END_REF] .

Proposition A.10. Let Σ be a compact Riemann surface and Φ : Σ → R 3 ∪ {∞} be a branched immersion with a finite number of ends. Let g be the induced metric, n be its Gauss map, H its mean curvature, Å its tracefree second fundamental form and K its Gauss curvature. Then

Figure 1 .Figure 2 .

 12 Figure 1. The transformation of the Bryant's surfaces into a López surface

2 z 4 e 0 - 4A z 3 e 3 +e 3 +

 33 O(r).

A. 1

 1 Weighted Calderon-Zygmund Theorems A.1 and A.2 are taken from Y. Bernard and T. Rivière's [3] (Proposition C.2 and C.3). They were not originally stated with the case p = ∞, but the proof is the same. Theorem A.1. Let u ∈ C 2 (D\{0}) solve ∆u(z) = µ(z)f (z) in D, with f ∈ L p (D) for 2 < p ≤ ∞ and the weight µ satisfying for some a ∈ N |µ(z)| = O (|z| a ) . Then u z (z) = P (z) + |z| a T (z) with P ∈ C a [X] and T = O(|z| 1-2 p -κ ) for all κ > 0. More precisely one has

(

  |z| a T (z)) z |z| a . Remark A.1.1. Theorem A.1 works with a = 0, it is the classic Calderon-Zygmund theorem. Minimal bubbling for Willmore surfaces 31 Theorem A.2. Let u ∈ C 2 (D\{0}) solve ∆u(z) = µ(z)f (z) in D, with f ∈ L p (D) for 2 < p ≤ ∞ and the weight µ satisfying for some a ∈ R + |µ(z)| = O (|z| a ) . Then u z (z) = P (z) + |z| a T (z) with P ∈ C a [X] and T = O(|z| 1-2 p -κ ) for all κ > 0. Here a is the upper integral part of a. More precisely one has T

Lemma A. 6 .

 6 Let u ∈ N, v ≥ u and P ε = u j=0 p ε j z j ∈ C u [X] such that ∀p < ∞ P ε χ v ∈ L p . Then ∀ν > 0 ∀j ≤ u p ε j ε v-j-ν ≤ C ν .

S 3

 3 z = L, Φ z and R z = L × Φ z + 2H Φ z ,we successively compute:R z , Φ z = 0, R z , Φ z = e 2λ 2 2H + i L, n and R z , n = -iS z ,which proves the desired equality. A.Curvature formulas for branched immersions

Σ

  Kdvol g = 2π χ (Σ) +

Σ |∇ n| 2 dvol g = 4 ΣH 2 dvol g - 2 Σg = 2 ΣÅ 2 dvol g + 2 Σ

 42222 KdvolKdvol g .