Martin Traizet

OPENING NODES ON HOROSPHERE PACKINGS

Introduction

In this paper, a horosphere packing will be a finite, connected set of distinct horospheres in hyperbolic space H 3 , such that any two horospheres are either disjoint or tangent. We will denote n the number of horospheres and m the number of tangency points.

Horospheres have constant mean curvature equal to 1 (CMC-1 for short). In [START_REF] Pacard | Attaching handles to Bryant surfaces[END_REF], Pacard and Pimentel have constructed complete, embedded CMC-1 surfaces in hyperbolic space by desingularization of a horosphere packing. Heuristically, the horospheres are slightly shrinked and a suitable small catenoid is glued at each tangency point (see Figure 1).

On the other hand, CMC-1 surfaces in hyperbolic space admit a Weierstrass-type representation in term of meromorphic data, as discovered by Bryant [START_REF] Bryant | Surfaces of mean curvature one in hyperbolic space[END_REF]. In particular, they have a meromorphic Gauss map G and a holomorphic quadratic differential Q. Our goal in this paper is to revisit the result of Pacard and Pimentel using Bryant's representation. We prove: Theorem 1. Given a packing of n horospheres with m tangency points, there exists a smooth family (M s ) 0<s<ε of complete, embedded CMC-1 surfaces in hyperbolic space H 3 , such that M s converges when s → 0 to the given horosphere packing. The surfaces M s have genus m -n + 1 and n catenoid-cousin-type ends. They have finite total curvature.

We prove this theorem using Bryant's representation. Our point of view on Riemann surfaces is that of opening nodes. The construction follows the general strategy developped by the author in [START_REF] Traizet | Adding handles to Riemann minimal examples[END_REF][START_REF] Traizet | An embedded minimal surface with no symmetries[END_REF][START_REF] Traizet | On the genus of triply periodic minimal surfaces[END_REF] to construct minimal surfaces in euclidean space R 3 by opening nodes, using the Weierstrass representation. The new point is that instead of having to solve a period problem, which is homology invariant, we have to solve a monodromy problem, which is homotopy invariant hence more difficult.

The surfaces we construct will in fact depend on n complex parameters c 1 , • • • , c n , which are the limit points of the horospheres on the ideal boundary of H 3 (identified with the Riemann sphere) and n positive real parameters ξ 1 , • • • , ξ n which represent the speed at which each horosphere is "deflated" to accomodate the catenoidal necks, where we think of s as the time parameter. We can normalize ξ 1 = 1, so our family depends on 3n real parameters, which is the expected dimension for the space of CMC-1 surfaces with n ends.

One may ask what are the possible topologies that one can achieve with this construction. In other words, given the number of ends n, what are the possible genera g? Since one can always ignore some tangency points (see Remark 3), this boils down to the following question: Question 1. What is the maximum number of tangency points that a packing of n horospheres may have?

For n ≥ 3, Pacard and Pimentel explain how to construct horosphere packings with up to 3n -6 tangency points, which gives genus up to 2n -5. It turns out that one can do much better! We investigate Question 1 in Section 2. We prove that the number of tangency points is always less than 5n, and we give examples where it grows like 4n. Still, the question remains open.

The rest of the paper is organized as follows. In Section 3, we explain the principle of Bryant representation of CMC-1 surfaces in hyperbolic space, and recall standard material about monodromy and opening nodes. In Section 4, we construct a family of candidates for the holomorphic data of the CMC-1 immersions we want to construct. The monodromy problem is solved in Section 5 using the implicit function theorem. In Section 6, we study the geometry of the surfaces we have constructed and prove that they are embedded. Finally, an appendix contains several results of independent interest that we use in this construction.

1.1. Related works. Several authors have used using Bryant representation to study CMC-1 surfaces in hyperbolic space: Benoît Daniel [START_REF] Daniel | Surfaces de Bryant dans H 3 de type fini[END_REF][START_REF] Daniel | Flux for Bryant surfaces and applications to embedded ends of finite total curvature[END_REF][START_REF] Daniel | Minimal disks bounded by three straight lines in Euclidean space and trinoids in hyperbolic space[END_REF], Ricardo Sa Earp and Eric Toubiana [START_REF] Sa | On the geometry of constant mean curvature one surfaces in hyperbolic space[END_REF], Wayne Rossman, Masaaki Umehara and Kotaro Yamada [START_REF] Wayne Rossman | Irreducible constant mean curvature 1 surfaces in hyperbolic space with positive genus[END_REF], Masaaki Umehara and Kotaro Yamada [START_REF] Umehara | A parametrization of the Weierstrass formulae and perturbation of complete minimal surfaces in R 3 into the hyperbolic 3-space[END_REF][START_REF] Umehara | Complete surfaces of constant mean curvature-1 in the hyperbolic 3-space[END_REF][START_REF] Umehara | A duality on CMC-1 surfaces in hyperbolic space and a hyperbolic analogue of the Osserman inequality[END_REF]. Many examples have been constructed which are inspired from known minimal surfaces in R 3 .

Pascal Collin, Laurent Hauswirth and Harold Rosenberg [START_REF] Pascal Collin | The geometry of finite topology Bryant surfaces[END_REF] have proved the fundamental result that a properly embedded CMC-1 surface in H 3 of finite topology must have finite total curvature, and the Gauss map extends meromorphically to the conformal compactification.

Maximum number of tangency points of a horosphere packing

We start with the 2-dimensional case, namely horocycle packings in the hyperbolic plane H 2 , because in this case we know the exact answer.

Theorem 2. For n ≥ 2, the maximum number of tangency points that a packing of n horocycles may have is 2n -3.

Proof: We work in the disk model of H 2 . Let us first see that the bound 2n -3 can be achieved. Start with two tangent horocycles. Given two tangent horocycles, there exists two other horocycles which are tangent to both of them, by Apolonius three circle theorem (with the ideal boundary of H 2 as the third circle). Choose one of them and iterate the process. The number of tangency points is increased by 2 at each step, so this gives, for each n ≥ 2, a packing of n horocycles with 2n -3 tangency points.

Next let us prove that one cannot do better. Consider an arbitrary packing of n horocycles with m tangency points. We construct a planar graph with n + 1 vertices and m+n edges as follows: The vertices are the n points p 1 , • • • , p n where the horocycles touch the unit circle (ideal boundary of H 2 ), plus the point at ∞. We connect two vertices p i and p j by the geodesic joining them if the corresponding horocycles are tangent. We also connect each point p i to ∞ by a radial arc. This gives a planar graph. Moreover, each of its faces has at least 3 edges in its boundary. By a standard inequality for planar graph, the number of edges is at most 3N -6 where N is the number of vertices. This gives m + n ≤ 3(n + 1) -6 which gives the result.

2

Let us now consider the 3-dimensional case of horosphere packings in H 3 . Start with three mutually tangent horospheres. Given three mutually tangent horospheres, there exists two other horospheres which are tangent to all three (by the 3-dimensional version of Apolonius theorem). Choose one of them and iterate the process. This gives, for each n ≥ 3, a packing of n horospheres with m = 3n -6 tangency points. This is the construction described by Pacard and Pimentel in [START_REF] Pacard | Attaching handles to Bryant surfaces[END_REF]. A packing obtained by this process is called an apolonian packing.

Here is how to do better. We work in the half-space model of H 3 . Let M be the equilateral lattice in the horizontal plane generated by the vectors (1, 0) and (cos π 3 , sin π 3 ). Consider some radius R. Let M R = M ∩ D(0, R). For each point p ∈ M R , consider a euclidean sphere S p of radius 1 2 centered at (p, 1 2 ). Stack on top of that the horizontal plane x 3 = 1, which is a horosphere in the half-space model and is tangent to all horospheres S p . Let n(R) be the total number of horospheres and m(R) the number of tangency points. For example, for R = 1, we get n(1) = 8 and m(1) = 19, so this is already better than an apolonian packing. For R = 2, we get n(2) = 20 and m(2) = 61.

Let us estimate the ratio m(R)/n(R) for large values of R. Each sphere S p is tangent to at most 6 other spheres, so there are at most 3(n(R) -1) tangency points between them. Therefore, adding the n(R) -1 tangency points with the horizontal horosphere, m(R) ≤ 4(n(R) -1).

On the other hand, if a sphere S p has less than 6 tangency points with the other spheres,

then p ∈ M R \ D(0, R -1). The cardinal of this set is O(R). Since n(R) = O(R 2 ), this gives m(R) ≥ 4n(R) -O( n(R)). Hence lim R→∞ m(R) n(R) = 4.
In general, we have the following bound.

Theorem 3. For a packing of n ≥ 5 horospheres, the number m of tangency points satisfies m ≤ 5n -16.

Proof. First consider a packing of n = 5 horospheres. Then the five horospheres cannot be pairwise tangent to each other. Indeed, if this is the case, one can assume (by an isometry) that one of the horosphere is the horizontal plane x 3 = 1. The statement then boils down to the fact that four circles of the same radius in the plane cannot be pairwise tangent to each other. This implies that m ≤ 9 for n = 5.

Consider then a packing of n > 5 horospheres. We work in the upper-half space model of H 3 . We may choose the model so that all horospheres have distinct euclidean radius. (Indeed, generically, an isometry will change the ratio of the euclidean radii of two horospheres.) By the following lemma, the horosphere which has the smallest euclidean radius has at most 5 tangency points. Removing that horosphere, Theorem 3 follows by induction.

2 Lemma 1. Consider three horospheres S 1 , S 2 , S 3 such that S 1 , S 2 are tangent, S 1 , S 3 are tangent, and S 2 , S 3 are either disjoint or tangent. Let R i be the euclidean radius of S i and (p i , R i ) its center. Assume that R 1 < R 2 ≤ R 3 . Then the angle θ 1 of the triangle p 1 p 2 p 3 at p 1 satisfies θ 1 > π 3 . Proof: If S i and S j are disjoint then the distance between their centers is greater than

R i + R j , hence ||p i -p j || 2 + (R i -R j ) 2 > (R i + R j ) 2 which gives ||p i -p j || 2 > 4R i R j .
If S i and S j are tangent, these inequalities are equalities. Trigonometry in the triangle p 1 p 2 p 3 gives

cos θ 1 = ||p 1 -p 2 || 2 + ||p 1 -p 3 || 2 -||p 2 -p 3 || 2 2||p 1 -p 2 || ||p 1 -p 3 || ≤ R 1 R 2 + R 1 R 3 -R 2 R 3 2R 1 √ R 2 R 3 < 1 2 .
Let m max (n) be the maximum number of tangency points that a packing of n horospheres may have. Collecting the above results:

Corollary 1. lim sup n→∞ m max (n) n ∈ [4, 5].
3. Background material x, x = -1, x 0 > 0} with the induced metric, is the Minkowski model of H 3 . Bryant identifies L 4 with the space of 2 × 2 hermitian matrices by identifying (x 0 , x 1 , x 2 , x 3 ) with the matrix

X = x 0 + x 3 x 1 + ix 2 x 1 -ix 2 x 0 -x 3 .
Then H 3 becomes the set of positive definite hermitian matrices with determinant 1.

3.1.2. Bryant representation. Recall that a holomorphic map F : Σ → SL(2, C) is null if det(F -1 dF ) = 0.
Theorem 4 (Bryant [START_REF] Bryant | Surfaces of mean curvature one in hyperbolic space[END_REF]). Let Σ be a simply connected Riemann surface. Let F : Σ → SL(2, C) be a holomorphic null immersion. Then F F * : Σ → H 3 is a smooth conformal CMC-1 immersion. Conversely, if f : Σ → H 3 is a conformal CMC-1 immersion, there exists a holomorphic null immersion F : Σ → SL(2, C) such that f = F F * . Moreover, F is unique up to right multiplication by a constant matrix H ∈ SU (2).

If Σ is not simply connected, then F is only well defined on the universal cover of Σ and will have SU (2)-valued monodromy.

3.1.3. Global meromorphic data. Assume we are given a CMC-1 immersion f : Σ → H 3 . By Theorem 4, we can write locally f = F F * where F is a null holomorphic, SL(2, C)valued map. Consider the matrix of holomorphic 1-forms

A(z) = (dF (z))F (z) -1 ∈ sl(2, C)
where sl(2, C) is the Lie algebra of 2 × 2 matrices whose trace is zero. Then A is well defined on Σ: replacing F by F H does not change A. Define G = A 11 A 21 and Ω = A 21 , where A ij denote the coefficients of A. Then G is a meromorphic function and Ω is a holomorphic 1-form, both globally defined on Σ (except in the exceptional case where A 21 ≡ 0). Since the trace and determinant of A are zero, we have

(1) A = G -G 2 1 -G Ω
The function G is the Gauss map introduced by Bryant, and Q = Ω dG is the Hopf quadratic differential. We will see the geometric meaning of the Gauss map G in the next section. We call (G, Ω) the meromorphic data for the immersion f .

Conversely, here is a recipe to construct CM C-1 immersions in H 3 . Start with a Riemann surface Σ, a meromorphic function G and a holomorphic 1-form Ω on Σ, such that Ω has a zero at each pole of G with twice the multiplicity, and has no other zeros. Define the matrix A by [START_REF] Bryant | Surfaces of mean curvature one in hyperbolic space[END_REF]. Then A is a holomorphic, sl(2, C)-valued 1-form on Σ which does not vanish. Solve the linear differential system [START_REF] Pascal Collin | The geometry of finite topology Bryant surfaces[END_REF] dF (z) = A(z)F (z)

with initial data F (z 0 ) = F 0 ∈ SL(2, C).
The solution is a multi-valued holomorphic null immersion F : Σ → SL(2, C). (It is an immersion because A(z) = 0). If its monodromy happens to be in SU (2), then the immersion f = F F * : Σ → H 3 is well defined and has CMC-1.

Remark 1. Many authors consider instead the matrix F -1 dF and write

(3) F -1 dF = g -g 2 1 -g ω
where g is a meromorphic function and ω is a holomorphic 1-form. It turns out that (g, ω) is the Weierstrass data for the corresponding minimal surface in R 3 via the Lawson correspondence [START_REF] Lawson | Complete minimal surfaces in S 3[END_REF], which is why (3) has become so popular. The problem is that unless Σ is simply connected, F -1 dF is not globally defined on Σ. This is not surprising, since the Lawson correspondence is local. So this point of view is not appropriate if we are to construct high genus examples.

3.1.4. Bryant representation in the half-space model. A more familiar model of H 3 is the half-space x 3 > 0 in R 3 , with conformal metric x -2 3 (dx 2 1 +dx 2 2 +dx 2 3 ). The following results are proved in [START_REF] Rosenberg | The global theory of minimal surfaces in flat spaces[END_REF]. An orientation preserving isometry from the Minkowski model to the half-space model is given by

Φ(x 0 , x 1 , x 2 , x 3 ) = x 1 x 0 -x 3 , x 2 x 0 -x 3 , 1 x 0 -x 3 .
The immersion Φ • f is given in the half-space model by (4)

x 1 + ix 2 = F 11 F 21 + F 12 F 22 |F 21 | 2 + |F 22 | 2 x 3 = 1 |F 21 | 2 + |F 22 | 2
where F ij denote the coefficients of the matrix F . The ideal boundary of H 3 in this model is C ∪ {∞}. The Gauss map G has the following geometric interpretation: the normal geodesic ray originated from Φ • f (z) (in the direction of the mean curvature vector) hits the ideal boundary at the point G(z). where H ∈ SL(2, C) and X ∈ L 4 . The action preserves H 3 and its kernel is {±I 2 }, so we recognize P SL(2, C) as the group of direct isometries of H 3 . The action of SL(2, C) on H 3 extends to the ideal boundary as homographic transformation of the Riemann sphere, namely

H • z = H 11 z + H 12 H 21 z + H 22 .
If f : Σ → H 3 is a conformal CMC-1 immersion with Gauss map G and null holomorphic map F , then H • f has Gauss map H • G and null holomorphic map HF .

3.1.6. Horospheres. The Gauss map is constant on a horosphere, and that constant is the limit point of the horosphere on the ideal boundary of H 3 in the half-space model. (This follows from the geometric interpretation of the Gauss map in the half-space model.) If the horosphere is not a horizontal plane, then its meromorphic data is Σ = C, G = c, Ω = λ dz where c and λ are complex constants. The constant λ has no geometrical meaning and depends on the chosen conformal parametrization of the horosphere. The matrix A is given by

A = λ c -c 2 1 -c .
If the horosphere is a horizontal plane then G ≡ ∞ and we are in the exceptional case where the meromorphic data (G, Ω) is not defined. In this case, one has

A = 0 λ 0 0
for some complex number λ. In any case, since the matrix A is constant, the solution to (2) is F (z) = exp(zA)F 0 .

3.2. Linear differential systems. In this section, we recall standard facts about linear differential systems on a Riemann surface and setup some notations. A basic reference is [START_REF] Teschl | Ordinary differential equations and dynamical systems[END_REF]. Let Σ be a Riemann surface and A an n × n matrix of holomorphic 1-forms on Σ.

We consider the first order linear differential system on Σ

(5)

dY (z) = A(z)Y (z) 3.2.1.
Local theory: the principal solution. Assume that Σ is simply connected. Given

z 0 ∈ Σ, (5) has a unique solution Y : Σ → GL(n, C) such that Y (z 0 ) = I n . Following [14], we write Y (z) = Π(z, z 0 ). The map Π : Σ × Σ → GL(n, C
) is holomorphic in both variables and is called the principal solution. It satisfies

Π(z 3 , z 2 )Π(z 2 , z 1 ) = Π(z 3 , z 1 ). Given Y 0 ∈ GL(n, C), the solution Y such that Y (z 0 ) = Y 0 is given by Y (z) = Π(z, z 0 )Y 0 .
3.2.2. Global theory. Now assume that Σ is not simply connected. Then the principal solution Π(z, z 0 ) is not well defined: it depends on the homotopy class of the path from z 0 to z. If γ : [0, 1] → Σ is a path from z 0 to z, the solution Y of ( 5) such that Y (z 0 ) = I n , which exists in a simply connected neighborhood of γ(0), can be analytically continued along γ. Its value at γ(1) will be denoted Π(γ). When the path γ is clear from the context, we will still use the notation Π(z, z 0 ).

Let γ 1 , γ 2 : [0, 1] → Σ be two paths such that γ 1 (1) = γ 2 (0). We denote by γ 2 • γ 1 the path obtained by composing γ 1 and γ 2 . (The usual notation is γ 1 • γ 2 but in this context, it is more convenient, and customary, to reverse the order). Then ( 6)

Π(γ 2 • γ 1 ) = Π(γ 2 )Π(γ 1 ).
In particular, Π is a morphism from the fundamental group π 1 (Σ, z 0 ) to GL(n, C). (With the usual notation for the product in the fundamental group, it would be an antimorphism.)

3.2.3. Monodromy. The monodromy of a solution is usually defined as follows. Let Y be a solution of ( 5) in a simply connected neighorhood

U of z 0 . Let γ ∈ π 1 (Σ, z 0 ). Analytic continuation of Y along γ gives another solution of (5) in U , which is denoted γ • Y . There exists a matrix M ∈ GL(n, C) such that γ • Y = Y M . The matrix M is called the monodromy of Y along γ and is denoted M γ (Y ).
In term of the principal solution, one has

(7) M γ (Y ) = Y (z 0 ) -1 Π(γ)Y (z 0 ).

3.3.

Opening nodes. We recall the standard construction of opening nodes. Consider n copies of the Riemann sphere

C = C ∪ {∞}, labelled C 1 , • • • , C n . Consider 2m distinct points p 1 , • • • , p m , q 1 , • • • , q m in the disjoint union C 1 ∪ • • • ∪ C n . Identify p i with q i for 1 ≤ i ≤ m.
This defines a Rieman surface with nodes which we denote Σ 0 . We assume Σ 0 is connected.

To open nodes, consider local complex coordinates v i : V i → D(0, 1) in a neighborhood of p i and w i :

W i → D(0, 1) in a neighborhood of q i , with v i (p i ) = 0 and w i (q i ) = 0. We assume that the neighborhoods V 1 , • • • , V m , W 1 , • • • , W m are disjoint in C 1 ∪ • • • ∪ C n . Consider, for each 1 ≤ i ≤ m, a complexe parameter t i with |t i | < 1. If t i = 0, identify p i with q i as above. If t i = 0, remove the disks |v i | ≤ |t i | and |w i | ≤ |t i |. Identify the point z ∈ V i with the point z ∈ W i such that v i (z)w i (z ) = t i .
This creates a Riemann surface, possibly with nodes, which we denote Σ t , where t = (t 1 , • • • , t n ). When all t i are non zero, Σ t is a genuine Riemann surface of genus g = m -n + 1.

Observe that if t

i = 0, the circle |v i | = 1 is homologous in Σ to the circle |w i | = 1 with the opposite orientation. Consequently, if ω is a holomorphic 1-form on Σ t , ( 8
)
|v i |=1 ω = - |w i |=1 ω.
This makes the following definition natural.

Definition 1 (Bers).

A regular differential on Σ t is a holomorphic 1-form, which is allowed to have simples poles at p i and q i if t i = 0, with opposite residues.

By a theorem of Fay [START_REF] Fay | Theta Functions on Riemann Surfaces[END_REF], the space of regular differentials on Σ t has dimension g and admits a basis which depends holomorphically on t in a neighborhood of 0. For 1 ≤ i ≤ n, let J + i and J - i be the set of indices j such that p j ∈ C i and q j ∈ C i , respectively. As a consequence of Fay's theorem, one has: Theorem 5. For t in a neighborhood of 0, and for a = (a j ) 1≤j≤m ∈ C m satisfying (9)

j∈J + i a j - j∈J - i a j = 0 for 1 ≤ i ≤ n
there exists a unique regular differential ω = ω t,a on Σ t such that (10)

|v j |=1 ω = a j for 1 ≤ j ≤ m.
Moreover, ω t,a depends holomorphically on t (away from the nodes).

Proof: from Cauchy theorem in C i and (8), we see that ( 9) is necessary for ω to exist. If t = 0, the map ω → ( |v j |=1 ω) 1≤j≤m is an isomorphism from the space of regular differentials on Σ 0 to the space of vectors a ∈ C m satisfying [START_REF] Pacard | Attaching handles to Bryant surfaces[END_REF]. (This follows from the fact that a holomorphic 1-form on the Riemann sphere with simple poles is entirely determined by its residues.) Using Fay's theorem, this remains true for t in a neighborhood of 0. 2 Fay's proof is rather abstract and non-constructive. For an elementary proof of Theorem 5 based on the contraction mapping principle, see [START_REF] Traizet | Opening infinitely many nodes[END_REF]. One has a similar result for meromorphic 1-forms with poles at some points r 1 , • • • , r k distinct from the nodes: Theorem 6. For t in a neighborhood of 0, one can define a regular meromorphic differential ω on Σ t with poles at the points r 1 , • • • , r k by prescribing its principal part at each pole and its periods as in [START_REF] Rosenberg | The global theory of minimal surfaces in flat spaces[END_REF], replacing (9) by the restriction coming from the residue theorem, namely

(11) j∈J + i a j - j∈J - i a j + 2πi r j ∈C i Res r j ω = 0 for 1 ≤ i ≤ n.
Recall that the principal part of a meromorphic 1-form ω at a pole r is its equivalence class under the relation: ω ∼ ω if ω -ω is holomorphic in a neighborhood of r. The analogue of Fay's theorem for meromorphic differentials with simple poles is proved by Masur in [START_REF] Masur | The extension of the Weil-Petersson metric to the boundary of Teichmuller space[END_REF]. For a proof of Theorem 6 in the case of poles of arbitrary order, see [START_REF] Traizet | Opening infinitely many nodes[END_REF]. We will also need the following result to compute the partial derivatives of ω with respect to t.

Theorem 7. The partial derivative ∂ ∂t i ω t,a at t = 0 has two double poles at p i and q i , with principal parts

-dv i v 2 i Res q i ω 0,a w i at p i -dw i w 2 i Res p i ω 0,a v i at q i
and has vanishing periods on all circles

|v j | = 1.
This is proved in [START_REF] Traizet | On the genus of triply periodic minimal surfaces[END_REF], Lemma 3. See also [START_REF] Traizet | Opening infinitely many nodes[END_REF], Remark 5.6.

4.

The meromorphic data (Σ, G, Ω)

4.1.
Notations. The horospheres of our given horosphere packing are denoted S 1 , • • • , S n . We define the following sets, which we use to index various quantities:

I = {(i, j) ∈ [[1, n]
] 2 : i < j and S i , S j are tangent}

J + i = {j : (i, j) ∈ I} J - i = {j : (j, i) ∈ I} J i = J + i ∪ J - i .
Without loss of generality, we may assume (by applying an isometry) that each horosphere S i is not a horizontal plane in the half-space model. We fix a conformal parametrisation f i : C → S i and let G = c i Ω = λ i dz be its meromorphic data (see Section 3.1.6). For j ∈ J i , we denote p ij the point in C such that f i (p ij ) is the point S i ∩ S j . Without loss of generality, we may assume (by changing the parametrization f i ) that for each i, the disks D(p ij , 1) for j ∈ J i and D(0, 1) are pairwise disjoint. 4.2. The Riemann surface Σ. We consider n copies of the complex plane, denoted C 1 , • • • , C n , and m copies of the Riemann sphere C ∪ {∞}, denoted C ij for (i, j) ∈ I. We think of p ij as a point in C i . The points 0, 1 and ∞ in C ij will be denoted respectively 0 ij , 1 ij and ∞ ij . Heuristically, the reader should think of C 1 , • • • , C n as the parametrization domain for the horospheres S 1 , • • • , S n , and C ij \ {0, ∞} for (i, j) ∈ I as the parametrization domain for the catenoidal necks connecting them.

We define a Riemann surface Σ 0 with 2m nodes by identifying p ij with 0 ij and p ji with ∞ ij , for (i, j) ∈ I and call it Σ 0 . To open nodes, we consider the natural coordinates

v ij = z -p ij in a neighborhood of p ij in C i , w ij = z in a neighborhood of 0 ij in C ij , v ji = z -p ji in a neighborhood of p ji in C j and w ji = 1 z in a neighborhood of ∞ ij in C ij .
We open nodes as explained in Section 3.3, introducing a complex parameter t ij to open the node p ij ∼ 0 ij and another parameter t ji to open the node p ji ∼ ∞ ij . Let t = (t ij , t ji ) (i,j)∈I be the collection of these parameters. This defines a Riemann surface (possibly with nodes) which we denote Σ, or Σ t when we need to emphasize the dependance on the parameter t. We denote Σ the compactification of Σ obtained by adding the points ∞ 1 , • • • , ∞ n , where ∞ i denotes the point at infinity in C i .

For (i, j) ∈ I, we define α ij as the homology class of the circle |z -p ij | = 1 in C i , with the positive orientation. This is homologous, in Σ, to the unit circle in C ij , with the negative orientation, and to the circle |z -p ji | = 1 in C j , also with the negative orientation.

4.3.

The Gauss map G. Here are our requirements on the Gauss map G. At ∞ i , it should take the value c i . It should have a simple pole in each Riemann sphere C ij (because on each catenoidal neck, we expect a point where the mean curvature vector is vertical pointing up, in the half-space model). We choose the identification of the Riemann sphere with C ∪ {∞} so that this pole is z = 1. The following proposition tells us that these requirements completely determine G.

Proposition 1. For t small enough, there exists a unique meromorphic function G = G t on Σ t with the following properties :

• G t has m simple poles at the points

1 ij for (i, j) ∈ I, • G t (∞ i ) = c i for i = 1, • • • , n.
Moreover, G t depends holomorphically on t (away from the nodes and its poles) and at t = 0, we have

(12) G 0 (z) = c i in C i c j + c j -c i z -1 in C ij
Proof: We first define the differential µ = dG and we recover G by integration. By Theorem 6, there exists a unique meromorphic differential µ t on Σ t which has m double poles at 1 ij for (i, j) ∈ I with principal part

r ij dz (z -1) 2
and has vanishing period on all cycles α ij . Here the m complex numbers r ij are free parameters. At t = 0, we have

µ 0 = 0 in C i r ij dz (z-1) 2 in C ij .
Lemma 2. For t in a neighborhood of 0, there exist unique values of the parameters r ij such that

(13) ∞ j ∞ i µ t = c j -c i for (i, j) ∈ I.
Moreover, each r ij depends holomorphically on t, and when t = 0, r ij = c i -c j .

Proof: First observe that ( 13) is a linear system of m linear equations with m unknowns r ij , (i, j) ∈ I. Hence it suffices to prove that the system (13) is invertible when t = 0. In [START_REF] Taylor | Introduction to Differential Equations[END_REF], it is understood that the path from ∞ i to ∞ j goes through C ij . There is no canonical way to choose this path, but all choices are homologous modulo α ij . Since µ has no period on α ij , ∞ j ∞ i µ t is a well defined holomorphic function of t. Moreover, by Lemma 4 in [START_REF] Traizet | On the genus of triply periodic minimal surfaces[END_REF], this function extends holomorphically at t = 0 with value

∞ j ∞ i µ 0 = p ij ∞ i 0 + ∞ ij 0 ij r ij dz (z -1) 2 + ∞ j p ji 0 = -r ij .
The Lemma follows.

2

Returning to the proof of Proposition 1, we define the function G t on Σ t by

G t (z) = c 1 + z ∞ 1 µ t .
By Lemma 2, and the fact that µ t has no residues and no period on the cycles α ij , G t is well defined on Σ t (meaning that the integral does not depend on the path from ∞ 1 to z) and has all desired properties. 2

4.4. The holomorphic differential Ω. Here are our requirements on the holomorphic differential Ω. It should have a double pole at ∞, with leading term λ i dz, just like the horosphere S i . It also needs a double zero at each pole of G. We define Ω by prescribing poles, principal parts and periods, using Theorem 6. Then we adjust the parameter t so that Ω has the required zeros.

Definition 2. Consider m complex parameters a ij , (i, j) ∈ I and let a = (a ij ) (i,j)∈I .

We define Ω = Ω t,a as the unique meromorphic 1-form on Σ t with n double poles at

∞ 1 , • • • , ∞ n , with principal part λ i dz + j∈J + i a ij dz z - j∈J - i a ji dz z at ∞ i
and periods

α ij Ω = 2πi a ij for (i, j) ∈ I.
It depends holomorphically (away from its poles and the nodes) on t. Moreover, at t = 0 we have

(14) Ω 0,a =        λ i dz + j∈J + i a ij dz z -p ij - j∈J - i a ji dz z -p ji in C i -a ij dz z in C ij
Note that the residue of the prescribed principal part at ∞ i is

Res ∞ i Ω = - j∈J + i a ij + j∈J - i a ji
so Equation ( 11) holds.

Proposition 2. For a in a neighborhood of 0, there exists a unique value t(a), depending holomorphically on a, such that Ω t(a),a has a double zero at each pole of G, and has no other zeros in Σ (provided all parameters a ij are non-zero). Moreover, for each (i, j) ∈ I, we have

(15) a ij = 0 ⇒ t ij (a) = t ji (a) = 0 (16) ∂t ij (a) ∂a ij | a=0 = -1 2λ i (17) ∂t ji (a) ∂a ij | a=0 = 1 2λ j .
Proof: let (i, j) ∈ I. Using Theorem 7 and ( 14), we compute the partial derivatives of

Ω t,a in C ij at (t, a) = (0, 0). ∂Ω t,a ∂t ij = -dz z 2 Res p ij Ω 0,0 z -p ij = -λ i dz z 2 (18) ∂Ω t,a ∂t ji = -dw ji w 2 ji Res p ji Ω 0,0 z -p ji = λ j dz (recall w ji = 1 z ) (19) ∂Ω t,a ∂a ij = -dz z (20)
The partial derivatives of Ω t,a in C ij with respect to all other parameters t k and a k are zero. Write Ω t,a = f ij (t, a, z)dz in C ij . We want to solve f ij (t, a, 1) = f ij (t, a, 1) = 0. The Jacobian matrix of (f ij (t, a, 1), f ij (t, a, 1)) with respect to (t ij , t ji ) is

-λ i λ j 2λ i 0 .
The existence of the solution t(a) then follows from the implicit function theorem applied to the map (f ij (t, a, 1), f ij (t, a, 1)) (i,j)∈I whose Jacobian has block diagonal form. Next consider some (i, j) ∈ I and assume that a ij = 0. If t ij = t ji = 0, then Ω t,a has two simple poles at 0 ij and ∞ ij , with residue ±a ij = 0, hence is holomorphic in C ij , so Ω t,a ≡ 0 in C ij . So [START_REF] Traizet | Adding handles to Riemann minimal examples[END_REF] follows from uniqueness in the implicit function theorem. Equations ( 16) and ( 17) are obtained by differentiating f ij (t(a), a, 1) = f ij (t(a), a) = 0 with respect to a ij , using ( 18), ( 19) and [START_REF] Umehara | Complete surfaces of constant mean curvature-1 in the hyperbolic 3-space[END_REF]. Finally, if all parameters a ij are non-zero, then ( 16) and [START_REF] Traizet | On the genus of triply periodic minimal surfaces[END_REF] imply that all parameters t ij and t ji are non-zero, so Σ is a genuine compact Riemann surface of genus g = m -n + 1. Ω is a meromorphic 1-form with n doubles poles, so its number of zeros, counting multiplicity, is 2n + 2g -2 = 2m. This ensures that Ω has no other zeros than the double zeros it has at the m poles of G. 2 4.5. Partial derivatives with respect to the parameter a ij . Proposition 3. For (i, j) ∈ I we have at a = 0

(21) ∂G t(a) ∂a ij =          c j -c i 2λ i 1 z -p ij in C i c j -c i 2λ j 1 z -p ji in C j 0 in C k for k = i, j. (22) ∂Ω t(a),a ∂a ij =                  dz z -p ij in C i -dz z -p ji in C j (1 -z) 2 2z 2 dz in C ij 0 in C k for k = i, j and in C k for (k, ) = (i, j).
Proof: Recall that µ = dG. By Theorem 7,

∂µ t ∂t ij = - dz (z -p ij ) 2 Res 0 ij µ 0 = -(c i -c j ) dz (z -p ij ) 2 in C i .
Hence, by the chain rule and ( 16),

(23) ∂µ t(a) ∂a ij = c i -c j 2λ i dz (z -p ij ) 2 in C i .
Integrating, we obtain the first line of [START_REF] Umehara | A duality on CMC-1 surfaces in hyperbolic space and a hyperbolic analogue of the Osserman inequality[END_REF]. The proof of the second line is entirely similar. Regarding (22), we have (using again Theorem 7)

∂Ω t,a ∂t ij = 0 ∂Ω t,a ∂a ij = dz z -p ij in C i .
The first line of (22) follows from the chain rule. The proof of the second line is similar. Using Equations ( 18), ( 19), [START_REF] Umehara | Complete surfaces of constant mean curvature-1 in the hyperbolic 3-space[END_REF] and the chain rule, we have

∂Ω t(a),a ∂a ij = -λ i dz z 2 × -1 2λ i + λ j dz × 1 2λ j - dz z = (1 -z) 2 2z 2 dz in C ij .

The monodromy problem

5.1. Formulation of the problem. We consider the matrix A = A a defined by [START_REF] Bryant | Surfaces of mean curvature one in hyperbolic space[END_REF] with G = G t(a) and Ω = Ω t(a),a . Each coefficient of A is a holomorphic differential on Σ = Σ t(a) . Let 0 i be the point z = 0 in C i . Let F : Σ → SL(2, C) be the solution of dF = AF with initial condition F (0 1 ) = M 1 , where M 1 ∈ SL(2, C) is a matrix we can prescribe. (Observe that F (z) ∈ SL(2, C) because A(z) ∈ sl(2, C).) The solution F is of course only well defined on the universal cover of Σ. We need to adjust the parameters so that F has SU (2)-valued monodromy, so f = F F * is well defined on Σ. Taking 0 1 as a base point for the fundamental group and using [START_REF] Lawson | Complete minimal surfaces in S 3[END_REF], this is equivalent to

(24) ∀γ ∈ π 1 (Σ, 0 1 ), M -1 1 Π(γ)M 1 ∈ SU (2)
where Π denotes the principal solution of dF = AF on Σ (see Section 3.2.1). Instead of using a set of generators of π 1 (Σ, 0 1 ), which would involve in a complicated way the "combinatorics" of our given horosphere packing, we reformulate the monodromy problem in a more "local" way as follows. For (i, j) ∈ I, let γ ij ∈ π 1 (Σ, 0 i ) be a loop in C i with base point 0 i which goes around p ij and does not encircle any other node. We also define Γ ji as a path connecting 0 i to 0 j through C ij (to be defined more precisely later on).

Proposition 4. Given n matrices M 1 , • • • , M n in SL(2, C
), assume that for all (i, j) ∈ I:

(25) M -1 i Π(γ ij )M i ∈ SU (2) (26) M -1 j Π(Γ ji )M i ∈ SU (2)
Then (24) is satisfied. Moreover, the solution F of dF = AF with initial condition 25) and (26) hold for all (i, j) ∈ I.

F (0 1 ) = M 1 satisfies F (0 i ) ∈ M i × SU (2), hence f (0 i ) = M i M * i for 1 ≤ i ≤ n. Proof: Assume that (
Let Γ ij = Γ -1 ji . Using (6) M -1 i Π(Γ ij )M j = M -1 j Π(Γ ji )M i -1 ∈ SU (2). Define γ ji = Γ ji γ ij Γ ij ∈ π 1 (Σ, 0 j ). Using (6) again, M -1 j Π(γ ji )M j = M -1 j Π(Γ ji )M i M -1 i Π(γ ij )M i M -1 i Π(Γ ij )M j ∈ SU (2)
. In other words, (25) and (26) also hold for (j, i) ∈ I. Let C i be C i minus the disks D(p ij , 1) for j ∈ J i . The fundamental group π 1 (C i , 0 i ) is the free group with generators γ ij for j ∈ J i . Hence (25) implies that

(27) ∀δ ∈ π 1 (C i , 0 i ) M -1 i Π(δ)M i ∈ SU (2). Any element γ ∈ π 1 (Σ, 0 1 ) is homotopic to a product of the form δ k Γ i k i k-1 δ k-1 Γ i k-1 i k-2 δ k-2 • • • δ 2 Γ i 2 i 1 δ 1
where k ∈ N * , i 1 = i k = 1 and δ j ∈ π 1 (C i j , 0 i j ) for 1 ≤ j ≤ k. So (24) follows from ( 26) and (27). 2

Choice of the matrices

M 1 , • • • M n .
By the last statement of Proposition 4, choosing the matrices M i amounts to prescribe the image of the points 0 1 , • • • , 0 n . Recall from the introduction that we want to "deflate" the horosphere S i at speed ξ i . This suggests the following choice. Consider n fixed, positive numbers ξ 1 , • • • , ξ n and a real parameter s. (These are the same parameters as in the introduction). Let O i = f i (0) ∈ S i , where f i is the chosen conformal parametrization of the horosphere

S i . Choose M i (s) ∈ SL(2, C) so that s ∈ [0, ∞) → M i (s)M i (s) * (in the Minkowski model)
is the parametrization at speed ξ i of the geodesic ray normal to the horosphere S i at the point O i (in the direction of the mean curvature vector). The matrix M i (s) is unique up to right multiplication by an element in SU (2), which is clearly irrelevant for the monodromy problem.

Main result.

To solve the monodromy problem, we need to adjust the complex parameters p ij and p ji for (i, j) ∈ I. We will denote p 0 ij and p 0 ji the value of these parameters corresponding to the given horosphere packing (namely, such that

f i (p 0 ij ) = f j (p 0 ji ) = S i ∩ S j
). The matrix of holomorphic 1-forms A depends holomorphically on the parameters a = (a ij ) (i,j)∈I and p = (p ij , p ji ) (i,j)∈I and will be denoted A a,p . The principal solution of dF = A a,p F will be denoted Π a,p . Our goal is to prove: Proposition 5 (solution of the monodromy problem). For s > 0 small enough, there exists unique values a(s) and p(s) such that Equations (25) and (26), with M i = M i (s) and Π = Π a(s),p(s) , are satisfied for all (i, j) ∈ I. Moreover, a(s) and p(s) are smooth functions of s for s = 0, and extend continuously at s = 0 with value a ij (0) = 0 and p ij (0) = p 0 ij . 5.4. Choice of an isometry. From now on, we fix a couple (i, j) ∈ I. We have in mind to solve Equations ( 25) and (26). The computations will be simplified by applying a well chosen isometry. Let h be an orientation preserving isometry of H 3 such that h(S i ) is the horosphere x 3 = 1 and h(S i ∩ S j ) = (0, 0, 1), in the half-space model. (This isometry is unique up to composition by a rotation around the vertical axis.) Since the horospheres S i and S j are tangent, h(S j ) is the sphere of radius 1 2 centered at (0, 0, 1 2 ). The limit points of h(S i ) and h(S j ) are respectively ∞ and 0.

As explained in Section 3.1.5, the isometry h corresponds to a matrix H ∈ SL(2, C), unique up to sign. We have H • c i = ∞ and H • c j = 0, where the dot means the action by homography on the Riemann sphere. So H may be written in the form

H = 1 √ c j -c i ρ -ρc j ρ -1 -ρ -1 c i
where ρ is some complex number.

We use hats to denote the action of the isometry h on various objects: F = HF solves d F = A F where A = HAH -1 . An elementary computation gives

(28) A = 1 c j -c i (G -c i )(G -c j ) -ρ 2 (G -c j ) 2 ρ -2 (G -c i ) 2 -(G -c i )(G -c j ) Ω.
The principal solution of Y = AY is Π = HΠH -1 . Equations ( 25) and ( 26) are equivalent to

(29) M -1 i Π(γ ij ) M i ∈ SU (2) (30) M -1 j Π(Γ ji ) M i ∈ SU (2)
where M i = HM i and M j = HM j .

Remark 2. All these quantities: H, ρ, A, Π, M i , M j actually depend on both indices i and j because the chosen isometry h does. However, since i and j are fixed until the very end of Section 5.9, this dependence will not be written to make notations lighter.

5.5. Computation of the matrices M i (s) and M j (s). Consider the matrix

Ξ(s) = e s/2 0 0 e -s/2 ∈ SL(2, C).
Then s → Φ(Ξ(s)Ξ(s) * ) is the parametrization at unit speed of the positive vertical axis (oriented upwards) in the half-space model. (Here Φ is the isometry from the Minkowski model to the half-space model given in Section 3.1.4). The horosphere S i = h(S i ) is parametrized by f i = h • f i . We need to compute the corresponding null holomorphic map F i . By substitution of G = c i and Ω = λ i dz in (28), we obtain (31)

A i = 0 λ i 0 0 λ i = ρ 2 λ i (c i -c j )
By our choice of the isometry h, we have

F i (p 0 ij ) = I 2 . Hence (32) F i (z) = exp((z -p 0 ij ) A i ).
Then by our choice of the matrix M i (s) in Section 5.2, (33)

M i (s) = exp(-p 0 ij A i )Ξ(ξ i s)
In the same way, the horosphere S j has null holomorphic map F j given by (34)

F j (z) = exp((z -p 0 ji ) A j ). ( 35
) A j = 0 0 λ j 0 λ j = ρ -2 λ j (c j -c i )
and recalling that the mean curvature vector of S j at (0, 0, 1) points down, (36) M j (s) = exp(-p 0 ji A j )Ξ(-ξ j s) 5.6. Partial derivatives of the matrix A. We return to the matrix A a,p given by ( 28) with G = G t(a) and Ω = Ω t(a),a . Taking the derivative of (28) and using Equations ( 12), ( 14), ( 21) and ( 22), we obtain the following result after simplification: Proposition 6. At a = 0, we have:

(37) A 0,p = A i dz ∂ A a,p ∂a ij = c i -c j 2 1 0 0 -1 dz z -p ij in C i (38) A 0,p = A j dz ∂ A a,p ∂a ij = c j -c i 2 1 0 0 -1 dz z -p ji in C j (39) A 0,p = 0 ∂ A a,p ∂a ij = c j -c i 2 z -ρ 2 ρ -2 z 2 -z dz z 2 in C ij
where the matrices A i and A j are given by (31) and (35).

Expansion of Π(γ ij ).

Proposition 7. We have the following expansion

(40) Π a,p (γ ij ) = I 2 + a ij πi(c i -c j ) exp(-p ij A i ) 1 0 0 -1 exp(p ij A i ) + O(||a|| 2 ).
Proof: when a = 0, the principal solution in C i is given by Π 0,p (z, 0 i ) = exp(z A i ), which is well defined, so Π 0,p (γ ij ) = I 2 . By Proposition 9 in Appendix A, Equation (37) and the residue theorem:

∂ Π a,p (γ ij ) ∂a ij = c i -c j 2 γ ij exp(-z A i ) 1 0 0 -1 exp(z A i ) dz z -p ij = πi(c i -c j ) exp(-p ij A i ) 1 0 0 -1 exp(p ij A i ).
For (k, l) = (i, j), the partial derivative of A with respect to a k is holomorphic at p ij . Hence the partial derivative of Π(γ ij ) with respect to a k is zero. Equation (40) follows. 2

Expansion of Π(Γ ji ).

Recall that exp is a local diffeomorphism from a neighborhood of 0 in M 2 (C) to a neighborhood of I 2 in GL(2, C). We denote the inverse diffeomorphism by log. Also, for a matrix M and a complex number λ, we define M λ = exp(λ log M ).

Proposition 8. We have

(41) Π a,p (Γ ji ) × Π a,p (γ ij ) -2 2πi log a ij = exp(-p ji A j ) exp(p ij A i ) + O(a)
where O(a) is a well defined holomorphic function of (a, p) which vanishes at a = 0.

Proof: Let us fix the value of all parameters except a ij . We must first see that the left-hand side of (41) is a well defined function of a ij for a ij = 0, which of course means that Π(Γ ji ) itself is not. To see this, we have to define precisely the path Γ ji from 0 i to 0 j . For t ij and t ji non zero, we define Γ ji as the composition of the following two paths:

• The following path from 0 i to p ij + t ij in C i : a path from 0 i to p ij + 1 in C i minus all unit disks around the nodes (depending continuously on the parameter p ij in a neighborhood of p 0 ij ), composed with the spiral parametrized by

x → p ij + (t ij ) x for x ∈ [0, 1].
• The following path from p ji + t ji to 0 j in C j : the spiral parametrized by x → p ji + (t ji ) 1-x for x ∈ [0, 1], composed with a path from p ji + 1 to 0 j .

We can compose these two paths because the points p ij +t ij and p ji +t ji are both identified with 1 ij when opening nodes. Observe that we need a determination of the arguments of t ij and t ji to define the spirals. In other words, if we take t ij and t ji to live in the universal cover of the punctured unit disk (so arg t ij and arg t ji are well defined), Γ ji depends continuously on t ij and t ji . Now if the argument of a ij is increased by 2π, then the arguments of t ij and t ji are increased by the same amount. Hence the homotopy class of Γ ji is multiplied on the right by (γ ij ) 2 and Π(Γ ji ) is multiplied on the right by Π(γ ij ) 2 . Consequently, the left hand side of (41) is unchanged, so is a well defined holomorphic function of a ij for a ij = 0.

Next we prove that the left hand side of (41) is uniformly bounded. Because it is a well defined function of a ij , we can assume that arg a ij ∈ [-2π, 2π]. Using (6), we write

(42) Π(Γ ji ) = Π(0 j , p ji + 1) Π(p ji + 1, p ji + t ji ) Π(p ij + t ij , p ij + 1) Π(p ij + 1, 0 i ).
Since the path from 0 i to p ij + 1 stays in a fixed compact set of C i minus the nodes, where A is uniformly bounded, the fourth factor in (42) is uniformly bounded. We estimate the third factor using Proposition 10 from Appendix B. For this, we need an integral estimate of || A|| on the circle of center p ij and radius |t ij |/2.

I claim that A = O(a ij ) in compact subsets of C ij \{0
, ∞} (even if the other parameters a k are non-zero). Indeed, A depends holomorphically on a ij , and if a ij = 0, then by [START_REF] Traizet | Adding handles to Riemann minimal examples[END_REF], t ij = t ji = 0, so Ω = 0 and A = 0 in C ij . Also, by [START_REF] Traizet | An embedded minimal surface with no symmetries[END_REF],

a ij = O(t ij ). Consequently, since A is a matrix-valued 1-form, |z-p ij |= |t ij | 2 || A|| = |v ij |= |t ij | 2 || A|| = |w ij |=2 || A|| ≤ C|t ij |
for some uniform constant C. By Proposition 10 in Appendix B, Π(p ij + t ij , p ij + 1) is uniformly bounded. The first and second factors in (42) are estimated in the exact same way. We conclude that Π(Γ ji ) is uniformly bounded (although not well defined -but we assumed that arg a ij ∈ [-2π, 2π]). The left-hand side of ( 41) is now a bounded, well defined holomorphic function of (a, p) on the set a ij = 0. By Riemann extension theorem (in several variables), it extends holomorphically at a ij = 0.

To compute its value at a = 0, assume that all parameters a k for (k, ) = (i, j) are zero. By (37), A a,p -A i = O(a ij ) in compact subsets of C i minus the nodes. By Point (2) of Proposition 10 (with A = A i ), we obtain

|| Π a,p (p ij + t ij , 0 i ) -Π i (p ij + t ij , 0 i )|| ≤ C |a ij log |a ij || where Π i is the principal solution of Y = A i Y in C i , namely Π i (z, 0 i ) = exp(z A i ). This gives lim a ij →0 Π a,p (p ij + t ij , 0 i ) = exp(p ij A i ).
Arguing in the same way, we obtain

lim a ij →0
Π a,p (0 j , p ji + t ji ) = exp(-p ji A j ).

Proposition 8 follows. 2 5.9. Solution of the monodromy problem. We are now ready to prove Proposition 5. The unitary group SU (2) is not a complex manifold so we have to leave the realm of holomorphic functions. We introduce a small positive real number τ and have in mind to apply the implicit function theorem at τ = 0. We write

a ij = τ b ij c i -c j
where b ij is a complex number in a neighborhood of a non-zero central value b 0 ij . The computation will be simplified by knowing a priori the order of each parameter as a function of τ . The correct orders are s = -τ log τ p ij = p 0 ij + sq ij where q ij is a complex parameter in a neighborhood of 0. One issue here is that the function τ → τ log τ does not extend as a differentiable function at τ = 0. We solve this problem by writing τ = e -1/t 2 where t is a real parameter in a neighborhood of 0. Both τ and τ log τ extend smoothly at t = 0, and all parameters are smooth functions of t. Let b = (b ij ) (i,j)∈I and q = (q ij , q ji ) (i,j)∈I .

Recall that exp maps the Lie algebras sl(2, C) and su(2, C) to the Lie groups SL(2, C) and SU (2, C), respectively. We define

P ij = P ij (t, b, q) = log M i (s) -1 Π a,p (γ ij ) M i (s) ∈ sl(2, C) Q ij = Q ij (t, b, q) = log M j (s) -1 Π a,p (Γ ji ) M i (s) ∈ sl(2, C).
We want to solve P ij ∈ su(2, C) and Q ij ∈ su(2, C). We compute P ij using Proposition 7:

(43) Π(γ ij ) = I 2 + πiτ b ij exp(-p ij A i ) 1 0 0 -1 exp(p ij A i ) + O(τ 2 ) using (40)
.

P ij = πiτ b ij Ξ(-ξ i s) exp(-sq ij A i ) 1 0 0 -1 exp(sq ij A i )Ξ(ξ i s) + O(τ 2 ) using (33)
(44)

P ij = πiτ b ij 1 2 λ i sq ij 0 -1 + O(τ 2 )
We compute Q ij using Proposition 8:

Π(Γ ji ) = exp(-p ji A j ) exp(p ij A i ) Π(γ ij ) 2 2πi log a ij + O(τ ) using (41). Π(γ ij ) 2 2πi log a ij = I 2 -sb ij exp(-p ij A i ) 1 0 0 -1 exp(p ij A i ) + O(τ ) using (43). M j (s) -1 Π(Γ ji ) M i (s) = Ξ(ξ j s) exp(-sq ji A j ) 1 -sb ij 0 0 1 + sb ij exp(sq ij A i )Ξ(ξ i s)+O(τ ). ( 45
) Q ij = s ξ i +ξ j 2 -b ij λ i q ij -λ j q ji b ij - ξ i +ξ j 2 + O(τ ).
For a matrix M = (M k ) 1≤k, ≤2 ∈ sl(2, C), we have

M ∈ su(2, C) ⇔ Re(M 11 ) = 0 and M 12 + M 21 = 0.
Define the function F = (F ij ) (i,j)∈I for t = 0 by

F ij (t, b, q) = 1 τ Re(P 11 ij ) , 1 τ s (P 12 ij + P 21 ij ) , 1 s Re(Q 11 ij ) , 1 s (Q 12 ij + Q 21 ij ) .
We want to solve F(t, b, q) = 0. By (44) and (45), F extends smoothly at t = 0, with

F ij (0, b, q) = -πIm(b ij ) , 2πi b ij λ i q ij , ξ i +ξ j 2 -Re(b ij ) , λ i q ij -λ j q ji . Taking b 0 ij = ξ i +ξ j 2
> 0, we have F(0, b 0 , 0) = 0. It is straightforward that the partial differential of F with respect to the variables (b, q) at (0, b 0 , 0) is an isomorphism. By the implicit function theorem, for t in a neighborhood of 0, there exists (b(t), q(t)) depending smoothly on t such that F(t, b(t), q(t)) = 0. Proposition 5 is proved, and the monodromy problem is solved. 2

Embeddedness

Here is what we have achieved so far. For each small enough value of the parameter t > 0, we have constructed a null holomorphic map F which has SU (2)-valued monodromy. All parameters are now smooth functions of t > 0. To ease notation, the dependence on t will not be written. Let f : Σ → H 3 be the CMC-1 immersion associated to F . It remains to prove that f (Σ) is embedded. We work in the half-space model, so f (z) = Φ(F (z)F (z) * ) is given by formula [START_REF] Daniel | Flux for Bryant surfaces and applications to embedded ends of finite total curvature[END_REF]. Fix a small number ε > 0. We consider the following disjoint domains in Σ:

C ε i = {z ∈ C i : ∀j ∈ J i , |z -p 0 ij | > ε} C ε ij = {z C ij : ε < |z| < 1 ε }.
The complement of these domains in Σ are annuli which we call transition regions. We also fix some large number R and define C ε,R i = C i ∩ D(0, R).

6.1. Geometry of the image of C ε,R i . Fix some i, 1 ≤ i ≤ n, and consider an isometry h such that h(S i ) is the horosphere x 3 = 1 and h maps f i (0) to the point (0, 0, 1). The isometry h is represented by a matrix H ∈ SL(2, C) which has the form

H = 1 √ c -c i ρ -ρc ρ -1 -ρ -1 c i
where ρ, c are some complex numbers (ρ not the same as in Section 5.4). As in Section 5.4, we use hats to denote the action of h, so f = h • f , F = HF and so on. Equations ( 28) and (31) hold true, with c in place of c j . By construction, f (0 i ) parametrizes the vertical axis at speed ξ i as s varies, so F (0 i ) = Ξ(ξ i s), up to right multiplication by SU (2). We have

A = A i + O(τ ). Since C ε,R i is a fixed compact domain, F (z) = exp(z A i )Ξ(ξ i s) + O(τ ) for z ∈ C ε,R i .
From this, we conclude that f (C ε,R i ) converges smoothly to (a subdomain of) the horosphere x 3 = 1 as t → 0. Moreover, from (4), we get

x 3 (z) = e ξ i s + O(τ ) for z ∈ C ε,R
i so for t small enough, the image of C ε,R i lies above the horosphere x 3 = 1.

6.2. Geometry of the end at ∞ i . Next we prove that the image of |z| > R in C i is embedded. I claim that for t > 0 small enough, the Gauss map G has multiplicity 1 at ∞ i . This is delicate because G is constant when t = 0. We work in the local coordinate w = 1 z in a neighborhood of ∞ i and write G(w) = G(1/w). From (21), we obtain ∂ G (w)

∂a ij = c j -c i 2λ i 1 (1 -p ij w) 2 for j ∈ J + i ∂ G (w) ∂a ji = c i -c j 2λ i 1 (1 -p ji w) 2 for j ∈ J - i . G (0) = j∈J + i ∂ G (0) ∂a ij a ij + j∈J - i ∂ G (0) ∂a ji a ji + O(||a|| 2 ) = j∈J + i c j -c i 2λ i τ b ij c i -c j + j∈J - i c i -c j 2λ i τ b ji c j -c i + O(τ 2 ) = - τ ζ i 2λ i + o(τ ) where ζ i = 1 2 j∈J i (ξ i + ξ j ) > 0.
Hence for t > 0 small enough, G (0) = 0, so the Gauss map has multiplicity one at the end. To study the geometry of the end, we consider again the isometry h introduced in Section 6.1. Then G = H • G has a simple pole at w = 0 with residue

Res w=0 G = Res w=0 ρ 2 G -c G -c i = ρ 2 c i -c G (0) -2 λ i τ ζ i where λ i = ρ 2 λ i (c i -c). From (28) we obtain G 2 Ω = -A 21 -λ i dz = λ i dw w 2 .
By Theorem 8 in Appendix C (with α = -2 λ i τ ζ i and α 2 β = λ i ), there exists a uniform positive (independent of t) such that the image of 0 < |w| < is the vertical graph x 3 = u(x 1 , x 2 ) of a function u. Moreover, at infinity we have

log u(x 1 , x 2 ) (τ ζ i + o(τ )) log x 2 1 + x 2 2
so x 3 > 1 on the end. Replacing R by -1 if necessary, we obtain that f (C ε i ) is embedded. and moreover lies above the horosphere x 3 = 1 (using the maximum principle). In other words, f (C ε i ) lies on the mean-convex side of the horosphere S i . This guarantees that the images f (C ε i ) for 1 ≤ i ≤ n are disjoint. Remark 3. From this, we conclude that we can always ignore a tangency point by simply removing the corresponding couple (i, j) from I, and still obtain an embedded surface for t > 0.

6.3. Geometry of the catenoidal necks. Fix a couple (i, j) ∈ I. Consider again the isometry h introduced in Section 5.4, which maps the horosphere S i to the horosphere x 3 = 1 and the horosphere S j to the sphere of radius 1 2 centered at (0, 0, 1 2 ). In this section, we prove that after a blowup of ratio 1/τ , the image f (C ε ij ) converges to a vertical catenoid.

By computations similar to the computation of Q ij in Section 5.9, we have ( 46)

F (1 ij ) = Π(p ij + t ij , O i ) F (0 i ) = I 2 + s 2 ξ i -b ij 2 λ i q ij 0 b ij -ξ i + O(τ ).
By (39), we have in

C ε ij (47) A(z) = τ A(z) + O(τ 2 ) with A(z) = - b ij 2 z -ρ 2 ρ -2 z 2 -z dz z 2 . Let F (z) = 1 τ ( F (z) -F (1 ij )).
Then using ( 46) and (47),

d F (z) = 1 τ A(z) F (z) = ( A(z) + O(τ ))(I 2 + O(s)) = A(z) + O(s).
Since C ε ij is a fixed compact set, we obtain by integration

F (z) = - b ij 2 log z ρ 2 (z -1 -1) ρ -2 (z -1) -log z + O(s). Write x k (z) = x k (1 ij ) + τ x k (z) for 1 ≤ k ≤ 3. Using (4) and F (z) = I 2 + O(s), we obtain ( x 1 (z) + i x 2 (z), x 3 (z)) = F 12 (z) + F 21 (z), -2 Re F 22 (z) + O(s). lim t→0 ( x 1 (z) + i x 2 (z), x 3 (z)) = - ξ i + ξ j 4 ρ 2 (z -1 -1) + ρ -2 (z -1), 2 log |z| .
This is the parametrization of a vertical catenoid of necksize

ξ i +ξ j 2 .
This means that after a blowup of ratio 1 τ at f (1 ij ), the image of C ε ij converges smoothly to a catenoid. Also observe that the image of the circle |z| = ε lies above the image of |z| = 1 ε . Finally, (46) gives

x 3 (1 ij ) = 1 + s ξ i -ξ j 2 + O(τ ).
Hence the catenoidal neck lies below the image of C ε i .

6.4. Geometry of the transition regions. Fix (i, j) ∈ I and let U ij be the annulus in Σ bounded by the circles |z -

p ij | = ε in C i and |z| = ε in C ij .
We consider again the isometry h introduced in Section 5.4. Let us prove that the mean curvature vector of f is almost vertical in U ij . Given the geometric interpretation of the Gauss map given in Section 3.1.4, an elementary computation shows that the angle θ(z) between the mean curvature vector at f (z) and the vertical axis is related to the Gauss map G(z) by

(48) sin θ(z) 1 + cos θ(z) = x 3 (z) |G(z) -x 1 (z) -ix 2 (z)| .
The function G -1 is holomorphic in the annulus U ij and is bounded by Cε on the boundary circles, for some uniform constant C. By the maximum principle, G -1 is bounded by Cε in U ij . The norm of the holomorphic map F (z) -I 2 is bounded by Cε on the boundary of U ij , so is bounded by Cε in U ij by the maximum principle. Hence the function x 1 + i x 2 is uniformly bounded in U ij , and the height x 3 satisfies | x 3 -1| ≤ Cε in U ij . Using (48), we obtain sin θ(z)

1 + cos θ(z) ≤ Cε in U ij .
Hence by choosing ε small enough, we can ensure that θ(z) < π 2 . This implies that f (U ij ) is locally a vertical graph. Since we have already seen that it is a graph on the boundary circles, it is globally a graph. Moreover, by the maximum principle, it lies above the lowest point of the top boundary component of the catenoidal neck f (C ε ij ). The image of the annulus bounded by the circles |z -p ji | = ε in C j and |z| = 1 ε in C ij is studied in the same way, using an isometry which maps the horosphere S j to the horosphere x 3 = 1. This proves that f (Σ) is embedded and concludes the proof of Theorem 1. Proof. We use the letter C for uniform constants, depending only on c and ρ but not on t. By Cauchy theorem,

A(z) = 1 2πi |w|=ρ A(w) w -z dw - 1 2πi |w|=ρ -1 t
A(w) w -z dw.

Hence for t ≤ |z| ≤ 1,

||A(z)|| ≤ 1 2π c ρ -1 + ct t(1 -ρ -1 ) ≤ C.
We can connect 1 and z (in the universal cover of Ω) by a path of length less than 1 + c . The first Point of Proposition 10 follows from Gromwall inequality (Lemma 2.7 in [START_REF] Teschl | Ordinary differential equations and dynamical systems[END_REF]). Using Cauchy formula in the same way, we obtain for t ≤ |z| ≤ 1

||A(z) -A(z)|| ≤ 1 2π ct ρ -1 + ct |z| -ρ -1 t ≤ Ct |z| .
By the variation of constants formula (Theorem 3.12 in [START_REF] Teschl | Ordinary differential equations and dynamical systems[END_REF]) where U (z) ∈ GL(2, C) is well defined, holomorphic in D(0, 1) and satisfies U (0) = I 2 , and Y 0 ∈ GL(2, C) is a constant matrix. The monodromy of F on the unit circle γ is M γ ( F ) = M γ (F ) = Y -1 0 e 2πiA 0 Y 0 . Since f is well defined, its monodromy M γ (F ) belongs to SU (2) so its eigenvalues are complex number of modulus 1. This implies that the eigenvalues λ 1 , λ 2 of A 0 are real numbers, so αβ is real. To compute z A 0 , we write A 0 = P DP -1 with

Y (z) = Y (z) + Y (z)
P = 1 ∆ 1/4 -λ 2 αλ 1 α -1 λ 1 -λ 2 ∈ SL(2, C), D = λ 1 0 0 λ 2 .
By standard matrix theory, we can write P -1 Y 0 = T H where T is upper triangular and H ∈ SU (2). Then (51) gives First of all, our hypothesis imply the following bounds:

(53) 1 2 ≤ ∆ ≤ 3 2 , - 3
2 ≤ λ 2 ≤ -1 2 , |α -1 λ 1 | ≤ c and 1 2c ≤ |αλ 1 | ≤ 3c 2 .
Hence the matrix P is uniformly bounded. From (52), we obtain 1 = det F (z) = z det(U (z))z -1 det T.

Hence det U (z) is constant, and since U (0) = I 2 , we obtain det U (z) = det T = 1. The monodromy of F is given by M γ (F ) = H -1 T -1 e 2πiD T H = H -1 e 2πiλ 1 T 12 T 22 (e 2πiλ 1 -e 2πiλ 2 ) 0 e 2πiλ 2 H ∈ SU (2).

Since λ 1 -λ 2 is not an integer, this implies that T 12 = 0, so the matrix T is diagonal. Then T and z D commute. Equation (52) implies that U (z)P T is uniformly bounded on the unit circle. Since U (z)P T is holomorphic, it is uniformly bounded in the unit disk by the maximum principle. Taking z = 0, we obtain that T is uniformly bounded, hence U (z) is uniformly bounded in the unit disk. Expanding the product in (52), we obtain

1

 1

Figure 1 .

 1 Figure 1. Left: a packing of three horospheres with three tangency points in the ball model of H 3 . Right: a CMC-1 surface of genus one with three ends.

3. 1 . 5 .

 15 Isometries. The Lie group SL(2, C) acts isometrically on L 4 by the representation H • X = HXH *

2 course,

 2 Y (z) is only well defined in the universal cover of Ω: the value of Y (z) depends on the determination of arg z.) Proposition 10.(1) Assume that for some constant c, |z|=ρ ||A|| ≤ c and |z|=ρ -1 t ||A|| ≤ ct. Then for t ≤ |z| ≤ 1 and | arg z| ≤ c , ||Y (z)|| is bounded by a constant depending only on c, c and ρ. (2) Let A be another matrix-valued map satisfying the same hypotheses as A and let Y (z) be the solution of Y = A Y with initial condition Y (1) = I n . Assume moreover that |z|=ρ ||A -A|| ≤ ct. Then for t ≤ |z| ≤ 1 and | arg(z)| ≤ c , ||Y (z) -Y (z)|| ≤ Ct| log t| for some constant C depending only on c, c and ρ.

z 1 Y

 1 (w) -1 ( A(w) -A(w)) Y (w)dw.This gives|| Y (z) -Y (z)|| ≤ Ctz 1 |dw| |w| ≤ Ct (| arg z| + | log t|) .

  (52)F (z) = N (z)U (z)P z D P -1 Y 0 = N (z)U (z)P z D T HAssume that 1 c ≤ |α 2 β| ≤ c and ||F (z)|| ≤ c on the unit circle, for some real number c. We need a uniform bound (depending only on c) of U (z) in the unit disk. The theory of fuchsian systems gives us a bound of U (z) by construction, but this bound is not uniform as αβ → 0 (because we are approaching the resonant case). To obtain a uniform bound, we must use the fact that the monodromy of F is in SU (2).

4 z λ 1 (z αλ 1 λ 2 ( 1 +

 4121 -T 11 λ 2 + O(z)) z λ 1 (T 22 αλ 1 + O(z)) z 1+λ 1 (T 11 α -1 λ 2 + O(z)) z 1+λ 2 (-T 22 λ 2 + O(z)) H.where O(z) is holomorphic and uniformly bounded. Using (4) and the bounds (53), we obtain(x 1 + ix 2 )(z) = -1 O(z) + |z| ∆ O(1)) x 3 (z) = 1 |z| 2+2λ 2 √ ∆ |T 22 λ 2 | 2 (1 + O(z) + |z| ∆ O(1))where O(z) and O(1) are real analytic functions that have uniformly bounded derivatives and O(z) vanishes at the origin. The conclusions of Theorem 8 follow.

Consider a domain Ω ⊂ C and a point z 0 ∈ Ω. Let A λ (z) ∈ GL(n, C) be a family of matrices depending holomorphically on (λ, z) for z ∈ Ω and λ in a neighborhood of 0. Let Π λ denote the principal solution of Y = A λ Y in Ω. Proposition 9. For any γ ∈ π 1 (Ω, z 0 ),

By the variation of constants formula (Theorem 3.12 in [START_REF] Teschl | Ordinary differential equations and dynamical systems[END_REF])

Taking z = γ(1), the result follows. 2

Appendix B. Uniform estimates of the solution of Y = AY in an annulus

In this section, we consider the annulus Ω ⊂ C defined by ρ -1 t < |z| < ρ, where ρ > 1 is some fixed number and t is a small positive parameter. We are aiming for estimates which are uniform with respect to t.

Theorem 8. Let f : D * (0, 1) → H 3 be a conformal, CMC-1 immersion of the punctured closed unit disk. Assume that the Gauss map G has a simple pole at 0, with residue α, and the holomorphic differential Ω is holomorphic at 0, with Ω(0) = βdz. Assume that 0 < |αβ| ≤ 1 8 . Then there exists ε > 0 such that in the half-space model, f (D * (0, ε)) is a vertical graph x 3 = u(x 1 , x 2 ) of a (positive) function u over an exterior domain in the plane. The number ε only depends on a bound on |α 2 β| ±1 and ||F (z)|| on the unit circle. Moreover, αβ is real and at infinity, the function u has the following asymptotic behavior:

Remark 4. In this paper, we are interested in the case where αβ → 0 and we need a uniform positive ε. The conclusions of Theorem 8 remain true without the hypothesis |αβ| ≤ 1 8 but the proof is more involved, as the fuchsian system can be resonant. In particular, one can prove that αβ is always a real number in (-1 4 , ∞). Proof: We use the theory of fuchsian systems to compute F (z) such that f = F F * in the punctured disk. The system F = AF is fuchsian provided the matrix A(z) has a simple pole at 0, which is not the case here (it has a double pole). To circumvent this problem, we introduce the matrix where

Now the matrix A has a simple pole at 0, with residue

and the system (50) is fuchsian. The eigenvalues of A 0 are

where ∆ = 1 + 4αβ.

The system (50) is called resonant if λ 1 -λ 2 = √ ∆ is a non-zero integer. It follows from our hypothesis that 1 2 ≤ |∆| ≤ 3 2 and ∆ = 1, so the system is non-resonant. By the standard theory of fuchsian systems (Proposition 11.2 in [START_REF] Taylor | Introduction to Differential Equations[END_REF]), the solution of (50) has the form (51) F (z) = U (z)z A 0 Y 0