Laurent Ménard
email: laurent.menard@normalesup.org

Percolation probability and critical exponents for site percolation on the UIPT

Keywords: MSC 2010 Classification: 05A15, 05A16, 05C12, 05C30, 60C05, 60D05, 60K35

We derive three critical exponents for Bernoulli site percolation on the Uniform Infinite Planar Triangulation (UIPT). First we compute explicitly the probability that the root cluster is infinite. As a consequence, we show that the off-critical exponent for site percolation on the UIPT is β = 1/2. Then we establish an integral formula for the generating function of the number of vertices in the root cluster. We use this formula to prove that, at criticality, the probability that the root cluster has at least n vertices decays like n -1/7 . Finally, we also derive an expression for the law of the perimeter of the root cluster and use it to establish that, at criticality, the probability that the perimeter of the root cluster is equal to n decays like n -4/3 . Among these three exponents, only the last one was previously known. Our main tools are the so-called gasket decomposition of percolation clusters, generic properties of random Boltzmann maps, as well as analytic combinatorics.

Introduction

Percolation on random planar maps has been studied intensively since the pioneering work of Angel [START_REF] Angel | Growth and percolation on the uniform infinite planar triangulation[END_REF]. The main feature of random planar maps making this study so fruitful is the spatial Markov property. It can be used with two different approaches. The first approach is to perform an exploration process of percolation interfaces with the so-called peeling process. This is the approach developed by Angel [START_REF] Angel | Growth and percolation on the uniform infinite planar triangulation[END_REF] to prove that the threshold for Bernoulli site percolation on the Uniform Infinite Planar Triangulation (UIPT, the limit in law of large uniform random triangulations for the local topology [START_REF] Angel | Uniform infinite planar triangulations[END_REF]) is 1/2. This approach has been later used by several authors to study other models of percolation on maps, see for example [START_REF] Angel | Percolations on random maps I: half-plane models[END_REF][START_REF] Budd | Simple peeling of planar maps with application to site percolation[END_REF][START_REF] Gorny | The geometry of a critical percolation cluster on the UIPT[END_REF][START_REF] Ménard | Percolation on uniform infinite planar maps[END_REF][START_REF] Richier | Universal aspects of critical percolation on random half-planar maps[END_REF] and the references therein. The second approach is more global and consists in decomposing the map into the cluster of the root vertex and pieces filling the faces of this cluster. Such a decomposition is often called the Gasket decomposition, see for instance the works of Borot, Bouttier, Duplantier and Guitter [START_REF] Borot | Nesting statistics in the O(n) loop model on random planar maps[END_REF][START_REF] Borot | A recursive approach to the O(n) model on random maps via nested loops[END_REF][START_REF] Borot | Loop models on random maps via nested loops: the case of domain symmetry breaking and application to the Potts model[END_REF]. This second approach has been used very recently to study percolation on random finite triangulations by Bernardi, Curien and Miermont [START_REF] Bernardi | A Boltzmann approach to percolation on random triangulations[END_REF], following the previous work by Curien and Kortchemski [START_REF] Curien | Percolation on random triangulations and stable looptrees[END_REF]; and on other natural models of random finite planar maps by Curien and Richier [START_REF] Curien | Duality of random planar maps via percolation[END_REF]. This work builds on the article by Bernardi-Curien-Miermont [START_REF] Bernardi | A Boltzmann approach to percolation on random triangulations[END_REF] to study site percolation on the UIPT (of type 1, with loops and multiple edges allowed). Our first main result is an explicit formula for the probability that percolation from the root occurs.

Theorem 1. Let P p

∞ denote the law of the type 1 UIPT, with vertices colored black with probability p and white with probability 1 -p, and conditioned on the event where the root edge has both end vertices colored black. Let C denote the site percolation cluster of the root vertex under P p ∞ . Then, for every p ∈ [0, 1], we have

P p ∞ (|C| = ∞) = 2 √ 2p -1 √ 3 -cos 3 2 3 arccos( √ p) cos 2 3 arccos( √ p) 3/2 + p(2p -1) p (2 √ 3 -3(1 -p)) 1 p⩾1/2 .
In particular, the critical exponent is β = 1/2: as p → 1/2 + one has To the best of our knowledge, this is the first formula of this type and the first calculation of the critical exponent β for percolation on the UIPT. A similar formula was obtained for percolation on the Uniform Infinite Half-Plane Triangulation (UIHPT) by Angel and Curien [START_REF] Angel | Percolations on random maps I: half-plane models[END_REF] using the peeling process. They also calculate explicitly the probability that percolation from the root occurs and then obtain β = 1 in the UIHPT setting. The relation between these two exponents remains quite mysterious, and we do not know of any strategy to obtain the exponent β without first computing explicitly the probability of percolation.

P p ∞ (|C| = ∞) = 3 1/4 15 26 1 + 4 √ 3 3 p - 1 2 + O p - 1 2 .
Our second main result gathers estimates for the volume and perimeter of critical percolation clusters.

Theorem 2. With the notation of Theorem 1 we have for an explicit constant κ > 0:

P 1/2 ∞ (|V (C)| ⩾ n) ∼ n→∞ κ n -1/7 .
Furthermore, let ∂C denote the root face of C, that is the face of C containing the root face of the UIPT. Then there is an explicit constant κ ′ > 0 such that

P 1/2 ∞ (|V (∂C)| = n) ∼ n→∞ κ ′ n -4/3 .
The perimeter exponent 4/3 was established by Curien and Kortchemski [START_REF] Curien | Percolation on random triangulations and stable looptrees[END_REF] using the gasket decomposition but with a different approach than the present work. The exponent 1/7 was conjectured in [START_REF] Gorny | The geometry of a critical percolation cluster on the UIPT[END_REF] using heuristics for the peeling process, and the present work is the first time it is established rigorously. Previous works that computed volume exponents for critical percolation models on infinite random planar maps did so for cluster hulls (part of the maps separated from infinity by a percolation interface). The reason being that all these works use variations around the peeling process, which is particularly well suited to study percolation interfaces -and therefore hulls -but not so useful to study the geometry of the clusters themselves. For instance, let H denote the hull of the root cluster C in the UIPT; i.e. H is the complement of the only infinite connected component of T ∞ \ C, where T ∞ is the UIPT with a critical site percolation configuration. Gorny, Maurel-Segala and Singh [START_REF] Gorny | The geometry of a critical percolation cluster on the UIPT[END_REF] proved that

P 1/2 ∞ (|V (H)| ⩾ n) ≍ n -1/8 and P 1/2 ∞ (|V (∂H)| ⩾ n) ≍ n -1/6
, where u n ≍ v n means that u n /v n is bounded.

Main ingredients and organization of the paper. For every site percolated finite rooted planar triangulation t, we denote by v • (t) and v • (t) its number of white (resp. black) vertices, and by e(t) its number of edges. For p ∈ (0, 1) and for t > 0 small enough, we can consider a random finite percolated triangulation t whose law gives a probability proportional to t e(t) p v•(t) (1 -p) v•(t) to every finite triangulation t. Let us denote by Z the partition function of this model and by P p the corresponding probability: Z(p, t) := t t e(t) p v•(t) (1 -p) v•(t) < +∞, and P p (t = t) = t e(t) p v•(t) (1 -p) v•(t) Z(p, t) .

The partition function Z and its generalizations to triangulations with a boundary (with additional parameters counting respectively the number of boundary vertices and boundary edges) are studied in Section 2. In particular, we establish a rational parametrization for the generalized partition function. This parametrization significantly simplifies the study of its analytic properties started in [START_REF] Bernardi | A Boltzmann approach to percolation on random triangulations[END_REF].

The gasket decomposition resides in the following statement. There exists an explicit sequence of positive numbers (q k (p, t)) k⩾1 such that, for every finite non atomic map m, the probability that the root cluster C(t) of the random triangulation t is equal to m is given by

P p (C(t) = m) = f ∈Faces(m) q deg(f ) (p, t) Z(p, t) .
In this sense, the root cluster is a Boltzmann random planar map associated to the weight sequence (q k (p, t)) k⩾1 . The properties of such random maps depend on the asymptotic behavior in k of the weight sequence, which can be inferred from the generating function of the weights. In our case, the weight generating function is closely related to the generalized partition function of percolated triangulations counted by edges, boundary edges, and boundary vertices. A crucial consequence of this relation is that it allows to calculate explicitly the so-called pointed disk generating function of the root cluster in terms of the singularities of the generalized partition function of percolated triangulations in Equations ( 12) and [START_REF] Borot | Nesting statistics in the O(n) loop model on random maps of arbitrary topologies[END_REF]. See Section 3 for details.

To study the origin cluster of the UIPT, we can condition the random triangulation t to have 3n edges. By continuity for the local topology, this gives

P p (C(t) = m|e(t) = 3n) = [t 3n ] f ∈Faces(m) q deg(f ) (p, t) [t 3n ]Z(p, t) → n→∞ P p ∞ (C = m). (1)
With a careful study of the dependency in t of the weight sequence (q k (p, t)) k⩾1 performed in Section 6.1, we are able to compute the above limit, giving the law of the root cluster in the UIPT on the event where it is finite.

We are then able to establish an integral formula for the sum of the limiting probability over every finite map, which is the probability that the root cluster of the UIPT is finite. See Proposition 7. In particular, this calculation uses the explicit universal form of the pointed disk generating functions and cylinder generating functions of Boltzmann maps. With the help of our rational parametrizations, we can then calculate explicitely our integral formula for the probability that the root cluster is finite and establish Theorem 1. This is performed in Section 4.

The formula we obtain for the limit ( 1) is quite easy to sum over maps with the same perimeter. As a consequence we are also able to calculate explicitely the law of the perimeter of the root cluster and obtain the second statement of Theorem 2 on the perimeter of the root cluster. This is done at the beginning of Section 5.2.

Finally, to compute the tail probability of the number of vertices in the root cluster at criticality (p = 1/2), we establish an integral formula for the generating function E p ∞ g |V (C)| in Section 5.2. See in particular identity (29), which also originates from our explicit formula for the limit [START_REF]Maple companion file[END_REF]. The expression we get involves two quantities. The first is the generating series derived from the asymptotics of the coefficients in t of the weights (q k (p, t)) k⩾1 studied in Section 6.1 for which we have an explicit parametric expression. The second quantity involves the pointed disk generating functions of Boltzmann maps with modified weight sequence (g (k-2)/2 q k (p, t)) k⩾1 . We can analyze the behavior as g → 1 -of these modified pointed disk generating functions using the bivariate generating functions associated to the Bouttier-Di Francesco-Guitter bijection [START_REF] Bouttier | Planar maps as labeled mobiles[END_REF] presented in Section 5.1 and their singularities obtained in Section 6.2. We put all that together in Section 5.2 to study the singular behavior as g →

1 -of E 1/2 ∞ g |V (C)| and prove Theorem 2.
On the robustness of our approach. We believe that it should be possible to adapt the strategy of this paper to other models of random planar maps and other statistical models. Indeed, what is needed first is the gasket decomposition; which exists for example for percolation on other models of maps [START_REF] Curien | Duality of random planar maps via percolation[END_REF], or for the O(n) model on maps [START_REF] Borot | Nesting statistics in the O(n) loop model on random planar maps[END_REF][START_REF] Borot | A recursive approach to the O(n) model on random maps via nested loops[END_REF]. We then need information on the generalized generating series of maps with a boundary and their singularities. In the present work we have explicit parametric expressions to simplify calculations, but we really only need to identify the nature of the singularities for the critical exponents. Any model for which such information is available should fall into the scope of our method.

In another article [START_REF] Albenque | Geometric properties of spin clusters in random triangulations coupled with an ising model[END_REF], we derive several critical exponents for the sign clusters in finite and infinite planar triangulations coupled with an Ising model. In particular, we establish in [START_REF] Albenque | Geometric properties of spin clusters in random triangulations coupled with an ising model[END_REF] counterparts for the main results of the present work. Theorem 1 and its counterpart for Ising model are unrelated, but Theorem 2 is a particular case of its Ising version at infinite temperature.

Links with other critical exponents for percolation on the UIPT. The volume exponent 1/7 of Theorem 2 suggests that the scaling limit of a large critical percolation cluster in the UIPT should be of quantum dimension 7/8. There is also very strong evidence (see for instance the works of Bernardi-Holden-Sun [START_REF] Bernardi | Percolation on triangulations: a bijective path to Liouville quantum gravity[END_REF] and Holden-Sun [START_REF] Bernardi | Percolation on triangulations: a bijective path to Liouville quantum gravity[END_REF]) that in some sense, this scaling limit is a CLE 6 on an independent pure Liouville Quantum Gravity surface. The quantum dimension 7/8 of the scaling limit of the root cluster and the dimension 91/48 of the gasket of a CLE 6 agree with the KPZ [START_REF] Knizhnik | Fractal structure of 2D-quantum gravity[END_REF] relation:

1 - 1 2 91 48 = 2 3 1 - 7 8 2 + 1 3 1 - 7 8 .
In a similar fashion, the perimeter exponent 4/3 of Theorem 2 agrees with the KPZ relation and the dimension 7 4 of a SLE 6 curve [6] The value of the critical exponent β agrees heuristically with known quantities of the UIPT and Kesten's scaling relations. Indeed, Bernardi-Holden-Sun [START_REF] Bernardi | Percolation on triangulations: a bijective path to Liouville quantum gravity[END_REF] and Holden-Sun [START_REF] Holden | Convergence of uniform triangulations under the cardy embedding[END_REF] established that the the number of percolation pivotal points in a random triangulation of size n is of order n 1/4 . In this sense, the quantum dimension of the set of pivotal points of critical percolation on the UIPT should be 1 (the map itself having dimension 4). This dimension can also be predicted with the KPZ relation and the dimension 3/4 of critical percolation pivotal points in Euclidean geometry [START_REF] Garban | Pivotal, cluster, and interface measures for critical planar percolation[END_REF]. On the other hand, the one arm exponent α 1 of critical percolation on the UIPT should be 1/2 from the quantum dimension 7/8 of the clusters. Kesten's scaling relation then states that the dimension of pivotal points should be α 1 /β, giving 1 with the exponent of Theorem 1.

Generating series

Let T be the generating series of rooted triangulations with a (not necessarily simple) boundary counted by edges (variable t), boundary length (parameter y), and boundary vertices (parameter p). That is, we define: Using the quadratic method, the authors of [START_REF] Bernardi | A Boltzmann approach to percolation on random triangulations[END_REF] establish the following algebraic equation for T 1 ≡ T 1 (p, t) that will be our starting point:

T (p,
64 T 3 1 t 5 -27 p 3 t 5 -96 T 2 1 p t 4 + 30 T 1 p 2 t 3 + p 3 t 2 + T 2 1 p t -T 1 p 2 = 0.
Up to a multiplicative constant p, the series T 1 is simply the generating series of triangulations of the 1-gon counted by edges and admits a proper rational parametrization: Lemma 3. Let U be the unique power series in t 3 having constant term 0 and satisfying

t 3 = ŵ(U ) := 1 2 U (1 -U )(1 -2U ).
The series t T 1 (p, t) seen as a series in t 3 admits the following proper rational parametrization in U :

t T 1 (p, t) = T1 (p, U ) := p U 1 -3U 1 -2U .
Furthermore, the series U (t 3 ) has a unique dominant singularity at t 3 = (t c ) 3 := √ 3 36 with the following singular behavior:

U (t 3 ) = 3 - √ 3 6 - √ 2 6 1 - t t c 3 1/2 + √ 3 54 1 - t t c 3 - 5 √ 2 648 1 - t t c 3 3/2 + O 1 - t t c 3 2 . ( 3
)
Proof. This result could be qualified as folklore since the series T 1 is just the generating series of triangulations of the 1-gon multiplied by p. The rational parametrization given is also the classical one. The fact that U (t 3 ) is unique comes from the Lagrangean form of the equation that defines it. We can also see that the algebraic equation satisfied by t T 1 (t) has a unique solution that is a power series with constant term 0. By composition and unicity we can see that indeed tT 1 (t) = T1 (p, U (t 3 )) as power series in t since both verify the algebraic equation and have constant term 0. The singular behavior of U is also very classical and without difficulties. See the Maple companion file [START_REF]Maple companion file[END_REF] for details.

Injecting the parametrization of t 3 and t T 1 by U in the equation for T then allows to establish a proper rational parametrization for T : Lemma 4. Recall the definition of the power series U ≡ U (t 3 ) from Lemma 3. Let V ≡ V (p, U, y) be the unique power series in Q[p, U ] y ⊂ Q[p] t 3 , y with constant term in y equal to 0 satisfying:

y = ŷ(p, U, V ) := 2 V (2 -4 U -V ) 4 p U (1 -U ) (1 -2U ) + 2 U (1 -3U ) V + 2 (1 -3U )V 2 -V 3 .
The series T (p, t, ty) seen as a series in t 3 and y is algebraic and admits the following proper rational parametrization in U and V

T (p, t, ty) = T (p, U, V ) := 4 p U (1 -U ) (1 -2U ) + 2 U (1 -3U ) V + 2 (1 -3U )V 2 -V 3 4 U (1 -U ) (1 -2U ) .
In addition, for any p ∈ (0, 1] and any fixed t ∈ (0, t c ], the series V (p, U (t Proof. All computations are available in the Maple companion file [START_REF]Maple companion file[END_REF].

The fact that V is uniquely defined as a power series comes from the Lagrangean form of the equation V = y × R(p, U, V ) with R a rational fraction such that R(p, U, 0) ̸ = 0. This form also implies by inductive calculation of the coefficients in y of V that they are all rational in p and U . Similarly, the equation ( 2) verified by T (p, t, ty) takes the form

p(T (p, t, ty) -p)(T (p, t, ty) -(1 -p)) = y T (p, t, ty) tT 1 (p, t) -p(T (p, t, ty) -p) -y 2 p t 3 T (p, t, ty) 2 .
Here again, by inductive calculation of the coefficients, we can see that this last equation has a unique solution that is a power series in t 3 and y with constant term in y equal to p. Note that since the series tT 1 (p, t) is algebraic, this equation also ensures that T (p, t, ty) is algebraic. By composition, T p, U (t 3 ), V (p, U (t 3 ), y) is a power series in t 3 and y with constant term in y equal to p. We can verify that it satisfies the same algebraic equation as T (p, t, ty), and therefore the two power series are identical. Now, fix t 3 ∈ (0, t c ]. The function V → ŷ(p, U (t 3 ), V ) has poles and stationary points that we can locate. Let us start with the poles. The denominator of ŷ(p, U (t 3 ), V ) is a polynomial of degree three in V . It is positive with positive derivative at V = 0, and changes signs between 1 -2U and 2(1 -2U ). Since the coefficient of V 3 is -1, this leaves two possibilites for the poles of ŷ(p, U (t 3 ), V ): there is always a pole between 1 -2U and 2(1 -2U ), and either there is no additional real pole or there are two negative poles.

The stationary points of ŷ(p, U (t 3 ), V ) are the roots of the polynomial 2 , where U stands for U (t 3 ). By computing the values of the polyomial at 0, 1 -2U and 2(1 -2U ), we can see that it has 4 real roots

-2V 4 + (-16U + 8)V 3 - 4(3U -1)(3U -2)V 2 -16U p(2U -1)(U -1)V -16U p(U -1)(2U -1)
V -(p, U ) < 0 < V + (p, U ) ⩽ 1 -2U ⩽ V l (p, U ) < 2(1 -2U ) < V r (p, U ) (the
case of a double root at 1-2U only happens when U = U c and is treated separately in the Maple file). It is also easy to see that ∂ V ŷ(p, U (t 3 ), 0) > 0, and therefore y + (p, t)

:= ŷ(p, U (t 3 ), V + (p, U (t 3 ))) > 0. We can define y -(p, t) := ŷ(p, U (t 3 ), V -(p, U (t 3 ))) < 0 when ŷ(p, U (t 3 ), V ) has no pole between V -(p, U (t 3 ))
and 0, and y -(p, t) := -∞ when it has such a pole. By singular inversion, the inverse function

V (p, U (t 3 ), y) of ŷ(p, U (t 3 ), V ) is analytic in the domain C \ ((-∞, y -(p, t)] ∪ [y + (p, t), +∞)
) and singular at both points y ± (p, t) when they are finite.

By composition y → T (p, t, ty) is analytic in the same domain. Since it has nonnegative coefficients, we know that it is singular at its radius of convergence. Therefore y + (p, t) is its radius of convergence and y -(p, t) ⩽ -y + (p, t). Checking that y -(p, t) < -y + (p, t) just with our current material is cumbersome but could be done. However, we do not need to go through this since this inequality is a direct consequence of equations ( 13) established with no computations at the end of Section 3. Finally, to see that y + (p, t) > 1, we first note that for p fixed, y + (p, t) is non increasing with respect to |t| and then check y + (p, t c ) > 1 by direct computation.

The root cluster as a Boltzmann map

In this whole section, p ∈ (0, 1) and t ∈ (0, t c ] are fixed.

The weight sequence from [7] with the edge parameter

In [7, Section 2.2] it was established that for a random site-percolated triangulation (conditioned on the event where both ends of the root edge are colored in black), the cluster of the root is a random Boltzmann map with weight sequence q(p, t) = (q k (p, t)) k⩾1 given by

q k (p, t) = 1 p   (pt) 3/2 δ {k=3} + (pt 3 ) k/2 l⩾0 k + l -1 k -1 [y l ]T (1 -p, t, ty)   , (4)
for k ⩾ 1 (see [START_REF] Bernardi | A Boltzmann approach to percolation on random triangulations[END_REF]Equation (9)]). We briefly recall here what this statment means and how to obtain it. For every percolated rooted planar triangulation t, we denote respectively by v • (t) and v • (t) the number of white (resp. black) vertices. Let T perc be the set of all percolated rooted triangulations with both end vertices of the root edge black. For every t ∈ T perc , we denote by C(t) the cluster of its root edge viewed as a planar rooted map (whose root is the root edge of the triangulation). The fact that the root cluster is a random Boltzmann map comes from the following identity (see Equation ( 7) of [START_REF] Bernardi | A Boltzmann approach to percolation on random triangulations[END_REF]): for every non atomic map m one has

t∈T perc : C(t)=m t e(t) p v•(t) (1 -p) v•(t) = p 2 • f ∈F (m) q deg(f ) (p, t).
(

) 5
This identity stems from the classical gasket decomposition [START_REF] Borot | A recursive approach to the O(n) model on random maps via nested loops[END_REF] sometimes called island decomposition [START_REF] Bernardi | A Boltzmann approach to percolation on random triangulations[END_REF] and we briefly explain how it is obtained in the following lines (see Figure 2 for an illustration of these arguments). Fix m a non atomic map and t a site-percolated triangulation with root cluster m colored black. For such a triangulation t, each face of its root cluster is filled with a triangulation with an arbitrary boundary of white vertices, and a necklace of triangles with no additional vertices between this triangulation with a boundary and the cluster. For each cluster face of degree k ⩾ 1 filled with a triangulation with perimeter l ⩾ 1, there are k+l-1 k-1 different possible necklaces, and each of these necklaces require an additional k + l edges. Taking into account the case where the cluster face has degree 3 and can be part of the triangulation, this writes:

t∈T perc : C(t)=m t e(t) p v•(t) (1 -p) v•(t) = p |V (m)| t |e(m)| f ∈F (m)   1 deg(f )=3 + l⩾0 deg(f ) + l -1 deg(f ) -1 t deg(f )+l [y l ]T (1 -p, t, y)   , = p 2 • f ∈F (m) 1 p (pt) deg(f )/2   1 deg(f )=3 + l⩾0 deg(f ) + l -1 deg(f ) -1 t deg(f )+l [y l ]T (1 -p, t, y)   ,
where we used Euler's relation

|V (m)| -2 = f ∈F (m) (deg(f )/2 -1)
in the last display. The expression [START_REF] Angel | Uniform infinite planar triangulations[END_REF] follows easily and we refer the reader to [7, Section 2.2] for additional details. Note for future reference that for any k ⩾ 1, the weight q k (p, t) is 1/p times the generating series of all triangulations of the k-gon with a weight t per inner edge, a weight √ pt per boundary edge, and a weight 1 -p per inner vertex adjacent to the boundary.

We can define the partition function of our percolated triangulations by:

Z(p, t) = 1 p 2 • t∈T perc t e(t) p v•(t) (1 -p) v•(t) . ( 6
)
From identity [START_REF] Angel | Uniform infinite planar triangulations[END_REF], denoting by M the set of all non atomic rooted planar maps, we see that

Z(p, t) = m∈M f ∈F (m) q deg(f ) (p, t). ( 7
)
Notice that this sum is finite when p ∈ (0, 1) and t ∈ (0, t c ], meaning that the weight sequence q(p, t) is admissible in the sense of of [START_REF] Miermont | An invariance principle for random planar maps[END_REF]. We will need the asymptotic behavior of the coefficients in t 3 of the series Z(p, t):

Proposition 5. Fix p ∈ (0, 1), we have

[t 3n ]Z(p, t) ∼ n→∞ √ 2 3p -3 + 2 √ 3 2p √ π t -3n c n -5/2 .
Proof. By opening the root edge of the triangulations appearing in the sum (6), we can see that, for p ∈ (0, 1) and |t| ⩽ t c , the partition function Z(p, t) is given by

Z(p, t) = 1 p 2 t T 2 (p, t) = 1 p 2 t 3 t 2 T 2 (p, t) = 1 t 3 tT 1 (p, t) 1 + 1 -p p 3 tT 1 (p, t) . ( 8
)
From Lemma 3, The function tT 1 (t) seen as a series in t 3 has a unique dominant singularity at t 3 c and we can obtain its asymptotic expansion at t 3 c from the expansion of U . As a consequence, the function Z(p, t) also has a unique dominant singularity at t 3 c ans we can easily obtain the following asymptotic expansion:

Z(p, t) = 3 √ 3 7 -4 √ 3 + p 4p - √ 3 3p -27 + 16 √ 3 4p 1 -t 3 /t 3 c + 2 √ 2 3p -3 + 2 √ 3 3p 1 -(t/t c ) 3 ) 3/2 + O 1 -t 3 /t 3 c ) 2 .
The asymptotic behavior of [t 3n ]Z(p, t) then follows from the classical transfer Theorem [20, Theorem VI.4].

Generating series of the weights and of the associated Boltzmann maps

The generating series of the weight sequence q(p, t) is straightforward to compute. As shown by the following lines, the expression we obtain is valid for every p ∈ (0, 1), t in the greater domain of analycity of U (t 3 ), and z close enough to 0. More precisely we need pt 3 z ∈ C \ [1, +∞) and

1 -pt 3 z -1
must belong to the domain of analycity in y of T (t, 1 -p, ty) which was studied in the previous Section. Under these conditions for p, t and z, the computation of the weight sequence generating function is as follows:

F q(p,t) (z) : = k⩾1 q k (p, t)z k , = 1 p (pt) 3/2 z 3 + 1 p l⩾0 [y l ]T (1 -p, t, ty) k⩾1 k + l -1 k -1 (pt 3 ) k/2 z k , = 1 p (pt) 3/2 z 3 + 1 p pt 3 z l⩾0 [y l ]T (1 -p, t, ty) 1 -pt 3 z -l-1 , = 1 p (pt) 3/2 z 3 + 1 p pt 3 z 1 -pt 3 z T 1 -p, t, t 1 -pt 3 z . ( 9
)
We will need expressions for the pointed and unpointed disk generating functions associated to the weight sequence q(p, t). For every l ⩾ 0, let M l denote the set of all rooted planar maps with root face of degree l (for l = 0, this set contains only the atomic map). The unpointed disk generating function is defined as follows for |z| large enough:

W q(p,t) (z) := l⩾0   m∈M l f ∈F (m)\{fr} q deg(f ) (p, t)   z -l-1 ,
where we denote the root face of a planar map by f r . From our discussion establishing (5), we can compute the coefficients of these series. Indeed, from Equation ( 5) and the fact that q l is 1/p times the generating series of triangulations with a boundary of perimeter l counted with a weight t per inner edge, a weight √ pt per boundary edge, and a weight 1 -p per inner vertex adjacent to the boundary, we have:

m∈M l t∈T perc : C(t)=m t e(t) p v•(t) (1 -p) v•(t) = p √ p t -l q l (p, t) t∈T l t e(t) p vout(t) .
Comparing with the right hand side of ( 5), this gives for every l ⩾ 1:

m∈M l f ∈F (m)\{fr} q deg(f ) (p, t) = 1 p √ pt l [y l ]T (p, t, y), ( 10
)
and thus

W q(p,t) (z) = 1 p z T p, t, t pt 3 z . ( 11
)
The pointed disk generating function is defined similarly:

W • q(p,t) (z) = l⩾0   m∈M l |V (m)| f ∈F (m)\{fr} q deg(f ) (p, t)   z -l-1 .
It has the following universal form: For p ∈ (0, 1) and and t ∈ (0, t c ] fixed, there exists real numbers

c + (p, t) > 2 and c -(p, t) ∈ (-c + (p, t), c + (p, t)) (the lower bound is excluded since our maps are not bipartite) such that for z ∈ C \ [c -(p, t), c + (p, t)] one has W • q(p,t) (z) = 1 (z -c + (p, t))(z -c -(p, t)) . ( 12
)
This expression, sometimes called the one-cut Lemma, appears in several references. See for example [START_REF] Budd | The peeling process of infinite Boltzmann planar maps[END_REF]Proposition 12] or [START_REF] Borot | More on the O(n) model on random maps via nested loops: loops with bending energy[END_REF]Section 6.1]. It is also established in these articles that the two disk generating functions W and W • have the same domain of analycity. Comparing our two expressions [START_REF] Borot | Loop models on random maps via nested loops: the case of domain symmetry breaking and application to the Potts model[END_REF] and ( 12), we see that as a consequence

c ± (p, t) = 1 p t 3 y ± (p, t) , ( 13
)
where y ± (p, t) are the respective positive and negative singularities in y of the series T (p, t, ty) defined in Lemma 4. Note that, as mentioned in the proof of Lemma 4, this directly implies that y -(p, t) < -y + (p, t).

Percolation probability

Recall that for any triangulation with a site percolation configuration, C denotes the percolation cluster of its root vertex. Also recall that P p ∞ denotes the law of the UIPT with vertices colored independently black with probability p and white with probability 1 -p, conditioned on the event where the root edge has both end vertices colored black. We want to identeify the law of C under P p ∞ on the event where it is finite. To express that law, we will need some additional notation.

Fix p ∈ (0, 1) and recall the definition of the power series V (p, U, y) of Lemma 4. Define V c ≡ V c (p, z) as the following power series in z:

V c (p, z) = V 1 -p, U (t 3 c ), 1 1 -z . ( 14
)
From Lemma 4, we know that this series is analytic on C \ 1 - tc) , +∞ . Now let ∆(p, z) be the power series in z defined by:

1 y + (1-p,
∆(p, z) = ∆(p, V c ) := 3 V c (2 √ 3 -3V c ) 9V 3 c -9( √ 3 + 1)V 2 c + 3(3 + 2 √ 3)V c -2(1 -p) √ 3 (3(p -1) + 2 √ 3) √ 3 -3V c 3 9V 3 c -9V 2 c √ 3 + 4(1 -p) √ 3 . ( 15
)
We will see in Lemma 8 that z → ∆(p, z) is also analytic on

C \ 1 - 1 y + (1 -p, t c ) , +∞ = C \ 1 y + (p, t c ) , +∞ , (16)
and in Lemma 9 that it is the generating series of positive numbers (δ k (p)) k⩾1 verifying

[t 3n ]q k (p, t) pt 3 -k [t 3n ]Z(p, t) → n→∞ δ k (p) (17)
for every k ⩾ 1.

We can now identify the law of the root cluster on the event where it is finite: Proposition 6. For every p ∈ (0, 1) and every non atomic rooted finite map m we have

P p ∞ (C = m) =   f ∈F (m) q deg(f ) (p, t c )   • f ∈F (m) (pt 3 c ) deg(f )/2 δ deg(f ) (p) q deg(f ) (p, t c ) .
Proof. For p ∈ (0, 1) and n ⩾ 1, we denote by P p n the law of a uniform triangulation of the sphere with 3n edges and vertices colored independently black with probability p and white with probability 1 -p, conditioned on the event where both end vertices of the root edge are colored black. From equation [START_REF] Bernardi | A Boltzmann approach to percolation on random triangulations[END_REF], we can write that for every finite non atomic rooted map m we have

P p n (C = m) = [t 3n ] f ∈F (m) q deg(f ) (p, t) [t 3n ]Z(p, t) . ( 18
)
Since the event {C = m} is continuous for the local topology, the limit as n → ∞ of the previous display gives access to the annealed law of the percolation cluster of the root in the UIPT.

As can be seen from Proposition 5 and Lemma 9, the asymptotic behaviors of [t 3n ]Z(p, t) and of [t 3n ](pt 3 ) -k/2 q k (p, t) for every fixed k ⩾ 1 are all of the form Cst • t -3n c n -5/2 . It is then very classical to establish the limit

lim n→∞ [t 3n ] f ∈F (m) q deg(f ) (p, t) [t 3n ]Z(p, t) = lim n→∞ [t 3n ](pt 3 ) e(m) f ∈F (m) (pt 3 ) -deg(f )/2 q deg(f ) (p, t) [t 3n ]Z(p, t) , = f ∈F (m) (pt 3 c ) deg(f )/2 δ deg(f ) (p) f ′ ∈F (m),f ′ ̸ =f q deg(f ) (p, t c ),
proving the proposition.

Our next result is an integral formula for the probability that the root cluster C of the UIPT is infinite: Proposition 7. Recall the definition of the quantities c + (p, t c ) and c + (p, t c ) given in equations [START_REF] Borot | More on the O(n) model on random maps via nested loops: loops with bending energy[END_REF] and [START_REF] Borot | Nesting statistics in the O(n) loop model on random maps of arbitrary topologies[END_REF]. For every p ∈ (0, 1) one has

P p ∞ (|C| < ∞) = 1 2π c + (p,tc) c -(p,tc) dz z ∆ p, pt 3 c z z + c + (p, t c ) + c -(p, t c ) 2 (c + (p, t c ) -z)(z -c -(p, t c )), (19)
where the function ∆(p, z) is defined in Equation [START_REF] Budd | The peeling process of infinite Boltzmann planar maps[END_REF].

Proof. We start with the formula established in Proposition 6. Summing over every finite non atomic map gives for every p ∈ (0, 1):

P p ∞ (|C| < ∞) = m∈M f ∈F (m) (pt 3 c ) deg(f )/2 δ deg(f ) (p) f ′ ∈F (m)\{f } q deg(f ′ ) (p, t c ).
By opening the root edges of the maps m involved in the previous display, the sum transforms into a sum over maps in M 2 , the set of all rooted planar maps with root face of degree 2:

P p ∞ (|C| < ∞) = m∈M 2 f ∈F (m)\{fr} (pt 3 c ) deg(f )/2 δ deg(f ) (p) f ′ ∈F (m)\{f,fr} q deg(f ′ ) (p, t c ),
where f r denotes the root face of a map. By rooting the non root face f involved in the last sum (so that f lies on the right hand side of this additional root) we can transform the sum into a sum over maps of the cylinder. For l 1 , l 2 ⩾ 1, let M (l 1 ,l 2 ) denote the set of all planar maps with two distinct marked rooted faces f 1 , f 2 of respective degree l 1 , l 2 (rooted so that the corresponding marked faces lie on the right hand side of their respective root). Taking into account the deg(f ) possible roots for the face f in the sum we have:

P p ∞ (|C| < ∞) = k⩾1 1 k m∈M (2,k) (pt 3 c ) k/2 δ k (p) f ′ ∈F (m)\{f 1 ,f 2 } q deg(f ′ ) (p, t c ). (20)
Setting for every k ⩾ 0

φ k (p) = 1 k m∈M (2,k) f ′ ∈F (m)\{f 1 ,f 2 } q deg(f ′ ) (p, t c ),
and denote the generating series of these numbers by

Φ(p, z) = k⩾1 φ k (p)z k ,
the sum in [START_REF] Ph | Analytic combinatorics[END_REF] takes the form of the Hadamard product of ∆(p, pt 3 c z) and Φ(p, z) evaluated at z = 1:

P p ∞ (|C| < ∞) = k⩾1 (pt 3 c ) k/2 δ k (p) • φ k (p) = ∆ p, pt 3 c z ⊙ Φ(p, z)| z=1 .
We will compute this Hadamard product with the help of its contour integral representation:

P p ∞ (|C| < ∞) = 1 2iπ γ dz z ∆ p, pt 3 c z Φ(p, 1/z),
where the contour γ must lie in a region enclosing 0 where both functions z → ∆ p, pt 3 c z and z → Φ(p, 1/z) are analytic (see for example [20, Section VI.10.2]).

To compute this integral, we first have to compute Φ. To this end, we introduce the cylinder generating functions of Boltzmann planar maps. It is defined for |z 1 | and |z 2 | large enough by

W cyl q(p,t) (z 1 , z 2 ) = l 1 ,l 2 ⩾0 m∈M (l 1 ,l 2 ) f ′ ∈F (m)\{f 1 ,f 2 } q deg(f ′ ) (p, t) z l 1 +1 1 z l 2 +1 2 .
In a similar fashion than the pointed disk generating function W • of Section 2, this series has a universal form involving the two constants c ± (p, t):

W cyl q(p,t) (z 1 , z 2 ) = 1 2(z 1 -z 2 ) 2 W • q(p,t) (z 1 )W • q(p,t) (z 2 ) z 1 z 2 - c + (p, t) + c -(p, t) 2 (z 1 + z 2 ) + c + (p, t)c -(p, t) -1 .
This formula appears in [START_REF] Borot | Nesting statistics in the O(n) loop model on random maps of arbitrary topologies[END_REF], [START_REF] Eynard | Counting surfaces[END_REF] and a multitude of other references. We refer to [START_REF] Albenque | Geometric properties of spin clusters in random triangulations coupled with an ising model[END_REF] for a proof and a review of the literature on this matter.

Of particular interest to us will be the generating function where the first root face has degree 2:

[z -3 1 ]W cyl q(p,t) (z 1 , z) = k⩾0 z -k-1 m∈M (2,k) f ′ ∈F (m)\{f 1 ,f 2 } q deg(f ′ ) (p, t), = -z + W • q(p,t) (z) z 2 - c + (p, t) + c -(p, t) 2 z + 1 4 c + (p, t)c -(p, t) - 1 8 c 2 + (p, t) + c 2 -(p, t) .
The antiderivative of this function has a simple expression giving the identity:

Φ(p, 1/z) = k⩾0 1 k z -k m∈M (2,k) f ′ ∈F (m)\{f 1 ,f 2 } q deg(f ′ ) (p, t c ), = z 2 2 - z + c + (p,tc)+c -(p,tc) 2 2 (z -c + (p, t c ))(z -c -(p, t c )).
From this expression, we see that the function

z → Φ(p, z) is analytic on C \ [c -(p, t c ), c + (p, t c )]. From Lemma 8, we know that z → ∆(p, pt 3 c z) is analytic on C \ [c + (p, t c ), +∞).
We cannot directly pick an appropriate contour γ to compute our Hadamard product of series evaluated at 1, however, if we take w ∈ (0, 1) we have

∆ p, pt 3 c z ⊙ Φ(p, z)| z=w = 1 2iπ γ(w) dz z ∆ p, w pt 3 c z Φ(p, 1/z), (21)
where the contour γ(w) encloses the interval [c -(p, t c ), c + (p, t c )] and crosses the positive real line at some point inside the interval (c + (p, t c ), c + (p, t c ) + ε) for some ε > 0. Note that the in the last display, the function z → ∆(p, w pt 3 c z)/z is well defined and continuous at z = 0. Using the fact that ∆(p, w pt 3 c z) is analytic inside γ(w) and deforming the contour gives, for every w ∈ (0, 1),

∆ p, pt 3 c z ⊙ Φ(p, z)| z=w = 1 2iπ γ(w) dz z ∆ p, w pt 3 c z z 2 2 - z + c + (p,tc)+c -(p,tc) 2 2 (z -c + (p, t c ))(z -c -(p, t c )) = - 1 4iπ γ(w) dz z ∆ p, w pt 3 c z z + c + (p, t c ) + c -(p, t c ) 2 (z -c + (p, t c ))(z -c -(p, t c )) = 1 2π c + (p,tc) c -(p,tc) dz z ∆ p, w pt 3 c z z + c + (p, t c ) + c -(p, t c ) 2 (c + (p, t c ) -z)(z -c -(p, t c )).
Taking the limit as w → 1 finaly gives the proposition

We are now ready to prove our first main result:

Proof of Theorem 1. The proof basically consists on computing the integral [START_REF] Eynard | Counting surfaces[END_REF], which still requires some work. The change of variables z = (1 -1/ȳ)(pt 3 c ) -1/2 gives:

P ν ∞ (|C| < ∞) = 1 2π p t 3 c y + (p,tc) y + (p,tc)-1 y -(p,tc) y -(p,tc)-1 dȳ ȳ(ȳ -1) ∆ p, V (1 -p, U (t 3 c ), ȳ) ȳ -1 ȳ + 1 2y -(p, t c ) + 1 2y + (p, t c ) × 1 y + (p, t c ) - ȳ -1 ȳ ȳ -1 ȳ - 1 y -(p, t c ) ,
where ∆ is defined in Equation ( 15) and V (1 -p, U, ȳ) is defined in Lemma 4. To calculate this new integral, we want to do the change of variables ȳ = ŷ(1 -p, U (t 3 c ), 2

√ 3/3 -V )
, where ŷ is defined in Lemma 4 (as we will see shortly, this change of variables instead of simply ȳ = ŷ(1 -p, U (t 3 c ), V ) simplifies calculations a bit more).

Recall that, from Lemma 4 and its proof, the function

V → ŷ(1 -p, U (t 3 c ), V ) is an increasing bijection from [V -(1 -p, U (t 3 c )), V + (1 -p, U (t 3 c ))] onto [y -(1 -p, t c ), y + (1 -p, t c )],
and is analytic on

(V -(1 -p, U (t 3 c )), V + (1 -p, U (t 3 c ))
). In view of this, for our change of variables, we want 2

√ 3/3 -V to be in a sub-interval of [V -(1 -p, U (t 3 c )), V + (1 -p, U (t 3 c ))
]. This will also enable us to use the rational parametrization of the series ∆ of Lemma 9. The details of the calculations that follow are available in the Maple file [START_REF]Maple companion file[END_REF].

We have to solve for V the two equations

ŷ 1 -p, U (t 3 c ), 2 √ 3 3 -V = y ± (p, t c ) y ± (p, t c ) -1 .
There is an interesting symmetry to exploit: for every p ∈ (0, 1) and every V ∈ C one has

ŷ p, U (t 3 c ), V = ŷ(1 -p, U (t 3 c ), 2 √ 3 3 -V ) ŷ(1 -p, U (t 3 c ), 2 √ 3 3 -V ) -1
, therefore, we want to solve for V the two equations

ŷ p, U (t 3 c ), V = y ± (p, t c ). (22)
Since ŷ p, U (t 3 c ), V is a rational fraction in V of degree 2 for its numerator and 3 for its denominator, the solutions of equation ( 22) are the roots of a polynomial of degree 3. By definition, one of the root of this polynomial is V ± (p, U (t 3 c )), and it will even be a double root except for V -(p, U (t 3 c )) when it is a negative pole of ŷ. Fortunately, we can compute explicitly these values. Indeed, the stationary points of V → ŷ(p, U (t 3 c ), V ) are the roots of the polynomial

9V 3 -9V 2 √ 3 + 4p √ 3 V - √ 3/3 .
The roots of the polynomial of degree 3 are given by

V m (p) = - √ 3 3 √ p cos arccos( √ p) 3 < 0, V l (p) = √ 3 3 √ p cos arccos( √ p) 3 -π 3 ∈ [0, 2 √ 3/3], V r (p) = 2 √ 3 3 -V m (p) ⩾ 2 √ 3/3.
From there, we see that

V -(p, U (t 3 c )) = V m (p)
when ŷ has no negative pole, which is the case when p > 1 2 -5 √ 3 18 ∼ 0.018. Since we are interested in p ⩾ 1/2, the case when ŷ has negative poles will not bother us. The positive singularity is given by

V + (p, U (t 3 c )) = √ 3 3 when p ⩾ 1/2, V l (p) when p ⩽ 1/2.
We can compute the third root of equation ( 22) from the constant term of the polynomial, which is always

-2 √ 3p 9 . When p > 1 2 -5 √ 3
18 , this root is given by

V i ± (p) = 2 √ 3p 9V ± (p, U (t 3 c )) 2 . Since V ± (p, U (t 3 c )) ⩽ √ 3/3, we can see that 2 √ 3/3 -V ± (p, U (t 3 c )) is always larger than √ 3/3, and therefore V + (1 -p, U (t 3 c )). From the explicit expression of V ± (p, U (t 3 c )), we can also check that 2 √ 3/3 -V i ± (p) is always between 0 and V + (1 -p, U (t 3 c )).
The correct bounds for our change of variables are then V i ± (p) since they are the only solutions inside the interval

[V -(1 -p, U (t 3 c )), V + (1 - p, U (t 3 c ))]. Note that we have 0 < V i + (p) < V i -(p).
This discussion also provides the following nice factorizations:

ŷ(1 -p, U (t 3 c ), 2 √ 3 2 -V )) -1 ŷ(1 -p, U (t 3 c ), 2 √ 3 2 -V )) - 1 y ± (p, t c ) = 1 ŷ(p, U (t 3 c ), V )) - 1 y ± (p, t c ) , = V -V i ± (p) V -V ± (p, U (t 3 c )) 2 2 V V -2 √ 3 3
.

With the explicit expression for ∆ in terms of V given in Lemma 9, it is easy to check:

∆ (p, V ) ∂ V ŷ(1 -p, U (t 3 c ), V ) ŷ(1 -p, U (t 3 c ), V ) (ŷ(1 -p, U (t 3 c ), V ) -1) = 1 3 (2 √ 3 -3(1 -p)) V - √ 3 3 2 .
Using the change of variables ȳ = ŷ(1 -p, U (t 3 c ), 2 √ 3/3 -V ) we finally get:

P ν ∞ (|C| < ∞) = 1 12 p t 3 c (2 √ 3 -3(1 -p)) π V i -(p) V i + (p) dV (V i -(p) -V )(V -V i + (p)) × (V -V + (p, U (t 3 c )))(V -V -(p, U (t 3 c ))) V 2 √ 3 3 -V V - √ 3 3 2 × 1 ŷ(p, U (t 3 c ), V ) + 1 2 1 y + (p, t c ) + 1 y -(p, t c )
.

Embarrassingly, we were not able to simplify directly this integral (it is equal to 1 when p < 1/2 since we know that the critical percolation threshold is 1/2 !). Fortunately, when p > 1/2, Maple is able to give an expression in terms of V ± (p, U (t 3 c )). Using basic trigonometric identities, we can further simplify this expression into the one given in the theorem after injecting the expressions for V ± (p, U (t 3 c )).

Cluster volume and perimeter

Admissibility equations and volume-modified weight sequence

We review in this Section additional background on Boltzmann planar maps that we will use in the proof of Theorem 2. We refer the reader to the references [START_REF] Bernardi | A Boltzmann approach to percolation on random triangulations[END_REF][START_REF] Budd | The peeling process of infinite Boltzmann planar maps[END_REF][START_REF] Miermont | An invariance principle for random planar maps[END_REF] for details.

For p ∈ (0, 1) and t ∈ (0, t c ], consider the two following bivariate power series in (z 1 , z 2 ):

f • (p, t; z 1 , z 2 ) = k,k ′ ⩾0 z k 1 z k ′ 2 2k + k ′ + 1 k + 1, k, k ′ q 2+2k+k ′ (p, t), f ⋄ (p, t; z 1 , z 2 ) = k,k ′ ⩾0 z k 1 z k ′ 2 2k + k ′ k, k, k ′ q 1+2k+k ′ (p, t).
These two functions are linked with Boltzmann maps with weight sequence (q k (p, t c )) k⩾1 by the Bouttier-Di Francesco-Guitter bijection [START_REF] Bouttier | Planar maps as labeled mobiles[END_REF] and will be instrumental in the remaining of this work. In our case, we can compute alternative expressions for these functions that will be more amenable to analysis, see Lemma 11.

Recall the parameters c + (p, t) and c -(p, t) defined in Equation ( 13). Since the weight sequence q(p, t) is admissible, the two functions defined above are well defined at least in the domain

|z 1 | ⩽ z + (p, t) and |z 2 | ⩽ z ⋄ (p, t)
where z + (p, t) and z ⋄ (p, t) are positive real numbers defined by

c ± (p, t) = z ⋄ (p, t) ± 2 z + (p, t). ( 23
)
Furthermore, from Proposition 4.2 and Lemma 4.4 of [START_REF] Bernardi | A Boltzmann approach to percolation on random triangulations[END_REF], (z + (p, t), z ⋄ (p, t)) is the minimal solution of the system of equations

   f • (p, t; z + (p, t), z ⋄ (p, t)) = 1 - 1 z + (p,t)) , f ⋄ (p, t; z + (p, t), z ⋄ (p, t)) = z ⋄ (p, t). ( 24
)
In addition, when p = 1/2 and t = t c , the weight sequence is critical (see Theorem 1.1 of [START_REF] Bernardi | A Boltzmann approach to percolation on random triangulations[END_REF]), and

(z + (1/2, t c ), z ⋄ (1/2, t c ))
is the unique solution of the system of equations ( 24) such that

(∂ z 2 + √ z 1 ∂ z 1 ) f ⋄ (p, t; z + (p, t), z ⋄ (p, t)) = 1. ( 25
)
For g ∈ (0, 1], let (z + (p, t; g), z ⋄ (p, t; g)) be the unique solution in (0, z + (p, t)] × (0, z ⋄ (p, t)] of the system of equations

   f • (p, t; z + (p, t; g), z ⋄ (p, t; g)) = 1 - g z + (p,t;g)) , f ⋄ (p, t; z + (p, t; g), z ⋄ (p, t; g)) = z ⋄ (p, t; g). ( 26
) Define c ± (p, t; g) = 1 √ g z ⋄ (p, t; g) ± 2 z + (p, t; g) .
From e.g. [START_REF] Bernardi | A Boltzmann approach to percolation on random triangulations[END_REF]Equations (4.3) and (4.4)], the pointed disk generating function of Boltzmann maps with modified weight sequence q(p, t; g) := g (k-2)/2 q k (p, t) k⩾1 defined as in ( 10) is given by

W • q(p,t;g) (z) = 1 (z -c + (p, t; g))(z -c -(p, t; g)) . ( 27
)

Proof of Theorem 2

We start with the perimeter exponent as it is the easiest of the two. From Equations ( 5) and ( 10), we can write

P p n (|V (∂C)| = k) = [t 3n ] q k (p, t) • [z -(k+1) ]W q(p,t) (z) [t 3n ]Z(p, t) = 1 p [t 3n ] pt 3 -k q k (p, t) • [y k ]T (p, t, ty) [t 3n ]Z(p, t) .
Applying Lemmas 9 and 10 directly give the limit of the previous display as n → ∞. Indeed, the series pt 3 -k q k (p, t) and [y k ]T (p, t, ty) have the same unique dominant singularity in t 3 and it is easy to get an asymptotic exansion at t 3 c of the product at this singularity from the two separate expansions. By continuity for the local topology, this limit is the probability that |V (∂C)| = k in the UIPT :

P p ∞ (|V (∂C)| = k) = 1 p δ k (p) t k c T k (p, t c ) + pt 3 c -k q k (p, t c ) θ k (p) .
Using the asymptotics (38), (40), ( 42) and (44) derived in Section 6.3 gives that for p = 1/2,

P 1/2 ∞ (|V (∂C)| = k) ∼ k→∞ κ ′ k -4/3 with κ ′ = -8 1 Γ(4/3) 8 3 5 6 351 + 2 3
1 3 117 3 5/6 2Γ(-2/3) ≃ 0.454, (28)
proving the second statement of Theorem 2.

We now turn on the first statement of the Theorem on the volume of the root cluster. Although we are only interested in the case p = 1/2, we start the proof with a generic p ∈ (0, 1) as it will be easier to follow. From proposition 6 we have for every g ⩽ 1

E p ∞ g |V (C)| = m∈M g |V (m)| f ∈F (m) (pt 3 c ) deg(f )/2 δ deg(f ) (p) f ′ ∈F (m)\{f } q deg(f ′ ) (p, t c ).
Applying Euler's formula then gives

E p ∞ g |V (C)| = g m∈M f ∈F (m) (g pt 3 c ) deg(f )/2 δ deg(f ) (p) f ′ ∈F (m)\{f } g (deg(f ′ )-2)/2 q deg(f ′ ) (p, t c ).
Using the exact same line of reasoning as in the proof of Theorem 1, we can compute this sum as the Hadamard product evaluated at z = 1 of the functions z → ∆(p, gpt 3 c ) and the function z → Φ(p, z) where c + (p, t c ) and c -(p, t c ) are replaced by the two constants c + (p, t c ; g) and c -(p, t c ; g) associated to the pointed disk generating function

((z -c + (p, t c ; g))(z -c -(p, t c ; g) -z)) -1/2 of Boltzmann maps with weight sequence q(p, t c ; g) = g (k-2)/2 q k (p, t c ) k⩾1
introduced in Section 5.1. Indeed, the first function should be obvious, and the second comes from the antiderivative of cylinder generating functions associated to the weight sequence q(p, t c ; g) instead of the weight sequence q(p, t c ; 1), which has the same universal form. Of course, it remains to calculate these two constants, which are less explicit than their counterparts c + (p, t c ) and c -(p, t c ). Nevertheless, after mimicking the part of the proof of Theorem 1 leading to [START_REF] Eynard | Counting surfaces[END_REF], we arrive at the formula

E p ∞ g |V (C)| = g 2π c + (p,tc;g) c -(p,tc;g) dz z ∆ p, g pt 3 c z z + c + (p, t c ; g) + c -(p, t c ; g) 2 × (c + (p, t c ; g) -z)(z -c -(p, t c ; g)). (29)
We want to study the asymptotic behavior of the integral (29) when g → 1 -. To do so, we must first study the dependency in g of the two constants c ± (p, t c ; g). We do not have a simple formula for the unpointed or pointed disk generating function of the Boltzmann maps as was the previously the case. However, the two constants can be studied via the solution the system of equations [START_REF] Miermont | An invariance principle for random planar maps[END_REF]. Let us denote z + (g) = z + (1/2, t c ; g) and z ⋄ (g) = z ⋄ (1/2, t c ; g) the solution of this system. We can calculate an expansion as g → 1 -of these two quantities with Lemma 12. Indeed, the development of f ⋄ gives

z ⋄ -z ⋄ (g) = ∂ z 1 f ⋄ (z + , z ⋄ ) z + -z + (g) + 1 - √ z + ∂ z 1 f ⋄ (z + , z ⋄ ) (z ⋄ -z ⋄ (g)) -κ ⋄ (z + -z + (g)) + √ z + (z ⋄ -z ⋄ (g)) 7/6 + o (z + -z + (g)) + √ z + (z ⋄ -z ⋄ (g)) 7/6
. This yields in turn

z + -z + (g) = √ z + (z ⋄ -z ⋄ (g)) + κ ⋄ ∂ z 1 f ⋄ (z + , z ⋄ ) (z + -z + (g)) + √ z + (z ⋄ -z ⋄ (g)) 7/6 + o (z + -z + (g)) + √ z + (z ⋄ -z ⋄ (g)) 7/6
.

Therefore we have

z ⋄ -z ⋄ (g) = 1 √ z + z + -z + (g) - κ ⋄ √ z + ∂ z 1 f ⋄ (z + , z ⋄ ) 2(z + -z + (g)) 7/6 + o z + -z + (g) 7/6 . ( 30
)
Now the development of f • gives

1 z + - g z + (g) = - 1 (z + ) 2 z + -z + (g) + κ ⋄ √ z + + κ • 2(z + -z + (g)) 7/6 + o 2(z + -z + (g)) 7/6
.

Expanding the left hand side of this identity then gives

1 -g z + + 1 -g (z + ) 2 z + -z + (g) = κ ⋄ √ z + + κ • 2(z + -z + (g)) 7/6 + o 2(z + -z + (g)) 7/6
, which in turns gives

z + (g) = z + - 1 2 κ ⋄ √ z + + κ • z + -6/7 (1 -g) 6/7 + o (1 -g) 6/7 .
Plugging this expression in equation (30) gives

z ⋄ (g) = z ⋄ - 1 2 √ z + κ ⋄ √ z + + κ • z + -6/7 (1 -g) 6/7 + o (1 -g) 6/7 .
We finally get an expansion for c ± (g):

c + (g) = z ⋄ (g) + 2 z + (g) = c + (1) - 1 √ z + κ ⋄ √ z + + κ • z + -6/7 (1 -g) 6/7 + o (1 -g) 6/7 , c -(g) = z ⋄ (g) -2 z + (g) = c -(1) + o (1 -g) 6/7 .
The change of variable z = ϕ(g, ξ)

:= c -(g) + (c + (g) -c -(g)) ξ in Equation (29) gives E p ∞ g |V (C)| = g(c + (g) -c -(g)) 2 2π 1 0 dξ ξ(1 -ξ) ∆ 1/2, g t 3 c /2 ϕ(g, ξ) ϕ(g, ξ) ϕ(g, ξ) + c + (g) + c -(g) 2 , = (c + (1) -c -(1)) 2 2π 1 0 dξ ξ(1 -ξ) ∆ 1/2, g t 3 c /2 ϕ(g, ξ) √ g ϕ(g, ξ) ϕ(1, ξ) + c + (1) + c -(1) 2 + O (1 -g) 6/7
where we used the developments of c ± (g) to obtain the second equality. In a similar fashion than in the proof of Lemma 12, the dominant singular term of the integral comes for the singularity at g = 1

and ξ = 1 of the term

∆ 1/2, √ g t 3 c /2 ϕ(g,ξ) √ g ϕ(g,ξ)
. The asymptotic expansion (41) established in Section 6.3 gives

∆ (1/2, z) z ∼ z→1/2 16 3 5 6 351 + 4 3
1 3 117 (1 -2z) -4/3 .
Applying the same techniques as in the proof of Lemma 12, we see that the main singular term of the expansion of E p ∞ g |V (C)| comes from the integral

κ π 1 0 dξ ξ(1 -ξ) 1 -2 g t 3 c /2 ϕ(g, ξ) -4/3 = κ 8 1 -2 g t 3 c /2c -(g) 4/3 2 F 1 4 3 , 3 2 ; 3; 2 g t 3 c /2 (c + (g) -c -(g)) 1 -2 g t 3 c /2c -(g) , with κ = (c + (1) -c -(1)) 2 ϕ(1, 1) + c + (1)+c -(1) 2 t 3 c /2 2 4 3 5 6

351

+ 3

1 3
117 .

Using the singular expansion of the hypergeometric function at 1 and the developments of c ± (g), we finally get

E p ∞ g |V (C)| = 1- κ 8 1 -2 t 3 c /2c -(1) 4/3 36 √ 3 Γ 5 6 Γ 2 3 π 3 2 2 t 3 c /2 1 -2 t 3 c /2c -(1) 1 √ z + κ ⋄ √ z + + κ • z + -6/7 (1 -g) 6/7 1/6 + o (1 -g) 1/7 .
This expansion is unfortunately not enough to extract the asymptotic behavior of the probabilities

P 1/2 ∞ (|V (C)| = n) as n → ∞.
Indeed, the generating series of these probabilities is E p ∞ g |V (C)| , but this function could have singularities of modulus 1 other than 1 contributing to the asymptotic. However, we do not have this problem for tail probabilities. For every n, denote

p n = P 1/2 ∞ (|V (C)| ⩾ n). A simple computation gives n⩾0 p n g n = 1 -g E p ∞ g |V (C)| 1 -g ∼ g→1 - κ 1 (1 -g) -6/7 ,
with

κ 1 = κ 8 1 -2 t 3 c /2c -(1) 4/3 36 √ 3 Γ 5 6 Γ 2 3 π 3 2 2 t 3 c /2 1 -2 t 3 c /2c -(1) 1 √ z + 1/6 κ ⋄ √ z + + κ • z + -8/7 .
From there, a classical Tauberian Theorem (see e.g. Theorem VI.13 of [START_REF] Ph | Analytic combinatorics[END_REF] and the following discussion) establishes the asymptotic behavior of p n given in the Theorem with 6 Technical Lemmas

κ = κ 1 Γ(8/7) = 63 3 1 3 + 4 3

Dependency in t of the weights

Recall that V c be the power series in z defined in [START_REF] Bouttier | Planar maps as labeled mobiles[END_REF] by

V c (p, z) = V (1 -p, U (t 3 c ), 1/(1 -z))
, where V is defined in Lemma 4. Also recall that ∆(p, z) is the power series in z defined in [START_REF] Budd | The peeling process of infinite Boltzmann planar maps[END_REF] by:

∆(p, z) = ∆(p, V c (p, z)) := 3 V c (2 √ 3 -3V c ) 9V 3 c -9( √ 3 + 1)V 2 c + 3(3 + 2 √ 3)V c -2(1 -p) √ 3 (3(p -1) + 2 √ 3) √ 3 -3V c 3 9V 3 c -9V 2 c √ 3 + 4(1 -p) √ 3 = k⩾0 δ k (p) z k .
Now that we know that qk (p, t) has a unique dominant singularity at t 3 c , it will follow from the general form of Puiseux expansions of algebraic series near their singularities (see [START_REF] Ph | Analytic combinatorics[END_REF]Theorem VII.7,p. 498]) that [t 3n ]q k (p, t) has the same asymptotic behavior as [t 3n ]Z(p, t) if we can prove that there exists two positive constants c, c ′ , depending on k and p such that for n large enough:

c • [t 3n ]Z(p, t) ⩽ [t 3n ]q k (p, t) ⩽ c ′ • [t 3n ]Z(p, t).
The upper bound is easily obtained by putting a triangulation of the 1-gon with a white boundary vertex inside a cycle of k black vertices and summing over every possible necklace between the two. The lower bound is obtained similarly starting from cycle of k black vertices by putting a single white vertex inside the cycle and an edge joining this additional vertex to every boundary vertex, and putting an arbitrary triangulation with white boundary together with a matching necklace on the outside of the cycle. To sum up, we have proved that for every k, qk (p, t) seen as a series in t 3 has a unique dominant singularity at t 3 c and that its asymptotic expansion at t 3 c is of the form

qk (p, t) = qk (p, t c ) -ãk (p) 1 - t 3 t 3 c + bk (p) 1 - t 3 t 3 c 3/2 + o 1 - t 3 t 3 c 3/2 . ( 33
)
This finishes the proof proof of the first statement [START_REF] Curien | Percolation on random triangulations and stable looptrees[END_REF] of the Lemma and we now have to identify the generating series of the numbers δ k (p) = bk (p)/κ(p), where κ(p) is the coefficient of the term

1 -t 3 t 3 c 3/2
in the asymptotic expansion of Z(p, t) calculated in Proposition 5.

Using the rational parametrization of Lemma 4, we can find a rational parametrization for F .

Indeed, if V = V (1 -p, U, y) is the power series in Q[p, U ] y ⊂ Q[p] t 3 , y defined in Lemma 4, we have F (p, t, z) = 1 p ŷ 1 -p, U (t 3 ), V 1 -p, U (t 3 ), 1/(1 -z) × T 1 -p, U (t 3 ), V 1 -p, U (t 3 ), 1/(1 -z) ,
where ŷ and T are rational fractions defined in Lemma 4. For z fixed such that |1/(1-z)| < y + (1-p, t c ) (which includes a neighborhood of 0 for z since y + (1 -p, t c ) > 1, the series F (p, t, z) seen as a series in t 3 has non negative coefficients and has radius of convergence t 3 c . In addition, this implies that (t 3 , z) → F (p, t, z) is analytic in the larger domain D(0, t 3 c ) × D(0, y + (1 -p, t c )/(y + (1 -p, t c ) -1)). We will produce an asymptotic expansion of F (p, t, z) near t 3 c using our rational parametrization. To do so, we start by computing the asymptotic expansion of V p, U (t 3 ), 1/(1 -z) near t 3 c , with z fixed.

First, writing

1 1 -z = ŷ 1 -p, U (t 3 ), V 1 -p, U (t 3 ), 1/(1 -z) = ŷ 1 -p, U (t 3 c ), V 1 -p, U (t 3 c ), 1/(1 -z)
we get an algebraic equation between

V z c = V 1 -p, U (t 3 c ), 1/(1 -z) , V 1 -p, U (t 3
), 1/(1 -z) and U (t 3 ). Plugging the asymptotic expansion (3) of U (t 3 ) in this equation, we obtain an asymptotic expansion for V p, U (t 3 ), 1/(1 -z) of the form:

V 1 -p, U (t 3 ), 1/(1 -z) = V z c -a 1 (V z c ) 1 - t 3 t 3 c 1/2 + a 2 (V z c ) 1 - t 3 t 3 c + a 3 (V z c ) 1 - t 3 t 3 c 3/2 + o 1 - t 3 t 3 c 3/2
, where the a i 's are explicit rational fractions whose expressions are given in the Maple companion file [START_REF]Maple companion file[END_REF]. Injecting in turn the asymptotic expansions in t 3 of U and V in T and ŷ we find an asymptotic expansion for F of the form:

F (p, t, z) = F (p, t c , z) + Ã(p, V z c ) 1 - t 3 t 3 c + B(p, V z c ) 1 - t 3 t 3 c + o 1 - t 3 t 3 c 3/2
, where à and B are explicit rational functions, and are analytic on the disk D(0, y + (1 -p, t c )/(y + (1p, t c ) -1) (this is obvious from their expressions: z → V z c is analytic in this region, and the poles of à and B fall outside it, see the Maple companion file [START_REF]Maple companion file[END_REF]). Note that the error term in the previous expansion is a priori not uniform in z. To ensure that à and B are the respective generating series of the numbers ãk and bk , we see that, as power series in (z, t 3 ) we have

Ã(p, V z c ) = lim t→tc F (p, t, z) -F (p, t c , z) • (1 -t 3 /t 3 c ) -1 , and B(p, V z c ) = lim t→tc F (p, t, z) -F (p, t c , z) -Ã(p, V z c ) 1 - t 3 t 3 c • (1 -t 3 /t 3 c ) -3/2 .
Combined with the analycity properties of these series, this ensures that Ã(p, V z c ) and B(p, V z c ) are indeed the generating series of the numbers ãk and bk .

Finally, the generating series of the numbers δ k is then given by

∆(p, z) = z κ(p) (1 -z) B(p, V z c ) = ŷ 1 -p, U (t 3 c ), V z c -1 κ(p) B(p, V z c ),
giving the expression of the Lemma. See the Maple file [START_REF]Maple companion file[END_REF] for detailed computations.

Applying the same proof to the function T (p, t, ty) instead of F (p, t, z) defined in Equation (32) allows to establish the asymptotic behavior in n of [t 3n ]t k T k (p, t) for every k ⩾ 0. We do not reproduce the proof as it will be almost exactly the same as the proof of Lemma 9, with no additional difficulties but with the function

1 -p y F 1 -p, t, 1 - 1 y instead of F (p, t, z).
The statement is as follows.

Lemma 10. Fix p ∈ (0, 1) and k ⩾ 1. One has

[t 3n ]t k T k (p, t) [t 3n ]Z(p, t) → n→∞ θ k (p) (34)
where the generating series of the numbers θ k (p) is given by

Θ(p, y) := k⩾0 θ k (p) y k = 1 -p y ∆ 1 -p, 1 - 1 y , ( 35
)
which is analytic on C \ [y + (p, t c ), +∞).

BDFG functions

We start with an integral formula for the functions f • and f ⋄ .

Lemma 11. Fix p ∈ (0, 1) and t ∈ (0,

t c ]. Then for 0 < z 1 ⩽ z + (p, t) and 0 < z 2 ⩽ z ⋄ (p, t) we have f • (p, t; z 1 , z 2 ) = pt 3 2z 2 + t 3 π 1- √ pt 3 (z 2 +2 √ z 1 ) -1 1- √ pt 3 (z 2 -2 √ z 1 ) -1 dz 1 -pt 3 z 2 -1 z 2pt 3 z 1 • T (1 -p, t, tz) 1 -z 1 -pt 3 (z 2 + 2 √ z 1 ) z 1 -pt 3 (z 2 -2 √ z 1 ) -1 , f ⋄ (p, t; z 1 , z 2 ) = pt 3 (2z 2 + z 2 1 ) + t 3 /p π 1- √ pt 3 (z 2 +2 √ z 1 ) -1 1- √ pt 3 (z 2 -2 √ z 1 ) -1 dz • T (1 -p, t, tz) 1 -z 1 -pt 3 (z 2 + 2 √ z 1 ) z 1 -pt 3 (z 2 -2 √ z 1 ) -1 .
Proof. Fix p ∈ (0, 1) and t ∈ (0, t c ]. Replacing the weights by their expression (4) in the definitions of the two functions gives

f • (p, t; z 1 , z 2 ) = 1 p 2(pt) 3/2 z 2 + pt 3 p k,k ′ ,l⩾0 2k + k ′ + l + 1 k + 1, k, k ′ , l (pt 3 z 1 ) k ( pt 3 z 2 ) k ′ [y l ]T (1 -p, t, ty), f ⋄ (p, t; z 1 , z 2 ) = 1 p (pt) 3/2 (2z 2 + z 2 1 ) + pt 3 p k,k ′ ,l⩾0 2k + k ′ + l k, k, k ′ , l (pt 3 z 1 ) k ( pt 3 z 2 ) k ′ [y l ]T (1 -p, t,

ty).

We can express these two functions as Hadamard products. Indeed, define the trivariate power series in (z 1 , z 2 ; z):

h • (z 1 , z 2 ; z) := k,k ′ ,l⩾0 2k + k ′ + l + 1 k + 1, k, k ′ , l (pt 3 z 1 ) k ( pt 3 z 2 ) k ′ z l , = 1 2pt 3 z 1     1 -z -pt 3 z 2 1 -z -pt 3 (z 2 + 2 √ z 1 ) 1 -z -pt 3 (z 2 -2 √ z 1 ) -1     ; and h ⋄ (z 1 , z 2 ; z) := k,k ′ ,l⩾0 2k + k ′ + l k, k, k ′ , l (pt 3 z 1 ) k ( pt 3 z 2 ) k ′ z l , = 1 1 -z -pt 3 (z 2 + 2 √ z 1 ) 1 -z -pt 3 (z 2 -2 √ z 1
) .

Then we have

f • (p, t; z 1 , z 2 ) = pt 3 2z 2 + t 3 T (1 -p, t, tz) ⊙ h • (z 1 , z 2 ; z)| z=1 , f ⋄ (p, t; z 1 , z 2 ) = pt 3 (2z 2 + z 2 1 ) + pt 3 p T (1 -p, t, tz) ⊙ h ⋄ (z 1 , z 2 ; z)| z=1 .
We can calculate these two Hadamard products as contour integrals in a similar fashion than in the proof of Theorem 1 where we established [START_REF] Garban | Pivotal, cluster, and interface measures for critical planar percolation[END_REF]. For (z 1 , z 2 ) ∈ (0, z + (p, t)] × (0, z ⋄ (p, t)] the functions h • and h ⋄ are analytic in the domain |z| < 1 -pt 3 (z 2 + 2 √ z 1 ) that contains the ,t) . This last domain is not empty since y + (p, t) > 1 from Lemma 4. Therefore, we can represent the Hadamard product as a contour integral similar to [START_REF] Garban | Pivotal, cluster, and interface measures for critical planar percolation[END_REF] if

domain |z| < 1 -pt 3 c + (p, t) = 1 -1 y + (p
1 -1 y + (p,t) -1
⩽ y + (1 -p, t). We can check that this is the case since we have explicit formulas for y + (p, t c ) and y + (p, t) is increasing in t (see the Maple file [START_REF]Maple companion file[END_REF] for details). Therefore, for a contour γ enclosing 0 and a point in the interval 1 - We obtain the expressions of the Lemma after simplifications and taking contours converging to the

cut 1 -pt 3 (z 2 -2 √ z 1 ) -1 , 1 -pt 3 (z 2 + 2 √ z 1 ) -1
.

When p = 1/2 and t = t c , we can compute explicitly an asymptotic expansion of f • and f ⋄ at the point (z + (p, t), z ⋄ (p, t)) that will be used in the proof of Theorem 2. and we have the following asymptotic expansions at (z + , z ⋄ ) -:

f • (z 1 , z 2 ) = 1 - 1 z + + 1 z + 1 z + - √ z + ∂ z 1 f ⋄ (z + , z ⋄ ) z 1 -z + + ∂ z 1 f ⋄ (z + , z ⋄ ) (z 2 -z ⋄ ) + κ • (z + -z 1 ) + √ z + (z ⋄ -z 2 ) 7/6 + o (z + -z 1 ) + √ z + (z ⋄ -z 2 ) 7/6 , f ⋄ (z 1 , z 2 ) = z ⋄ + ∂ z 1 f ⋄ (z + , z ⋄ ) z 1 -z + + 1 - √ z + ∂ z 1 f ⋄ (z + , z ⋄ ) (z 2 -z ⋄ ) + κ ⋄ (z + -z 1 ) + √ z + (z ⋄ -z 2 ) 7/6 + o (z + -z 1 ) + √ z + (z ⋄ -z 2 ) 7/6
, where The singular parts of the expansions of f • and f ⋄ come from the singularities of the form (1 -y/2) α in the development of T (1/2, t c , t c y) at y = 2 for α ∈ {2/3, 4/3} (it is straightforward to check that the linear term (1 -y/2) contributes only to non singular parts in the expansion). Indeed, for α ∈ {2/3, 4/3}, set

κ ⋄ = 4 2
I α (z 1 , z 2 ) = 1 π 1 0 dξ ξ -1/2 (1 -ξ) -1/2 1 - ϕ(z 1 , z 2 ; ξ) 2 α , = 1 - 1 2(1 -t 3 c /2(z 2 -2 √ z 1 )) α 1 π 1 0 dξ ξ -1/2 (1 -ξ) -1/2 •     1 -ξ 1- √ t 3 c /2(z 2 -2 √ z 1 ) 1- √ t 3 c /2(z 2 +2 √ z 1 ) -1 1 -2 t 3 c /2(z 2 -2 √ z 1 )     α , = 1 - 1 2(1 -t 3 c /2(z 2 -2 √ z 1 )) α 2 F 1     -α, 1 2
; 1;

1- √ t 3 c /2(z 2 -2 √ z 1 ) 1- √ t 3 c /2(z 2 +2 √ z 1 ) -1 1 -2 t 3 c /2(z 2 -2 √ z 1 )     ,
where we used Euler's integral representation of the hypergeometric function 2 F 1 in the last line. This last equality is valid when 0 < z 1 ⩽ z + and 0 < z 2 ⩽ z ⋄ since in this case t 3 c /2(z 2 + 2 √ z 1 ) ⩽ 1/2 and the variable in the hypergeometric function is in (0, 1]. Furthermore, using the values of z + and z ⋄ , a simple computation done in the Maple file [START_REF]Maple companion file[END_REF] gives: The first singular term of I 4/3 (z 1 , z 2 ) is similar, but with exponent 11/6 instead of 7/6. This means that the first singular term in the development of f ⋄ (1/2, z 1 , z 2 ) is from -3 5/6 2 t 3 c /p I 2/3 (z 1 , z 2 ) and we have

1- √ t 3 c /2(z 2 -2 √ z 1 ) 1- √ t 3 c /2(z 2 +2 √ z 1 )
f ⋄ (z 1 , z 2 ) = f ⋄ (z + , z ⋄ ) + ∇f ⋄ (z + , z ⋄ ) • z 1 -z + , z 2 -z ⋄ + 3 5/6 2 2t 3 c κ (z + -z 1 ) + √ z + (z ⋄ -z 2 ) 7/6 + o (z + -z 1 ) + √ z + (z ⋄ -z 2 )
7/6

.

The statement for f ⋄ follows using the fact that f ⋄ (z + , z ⋄ ) = z ⋄ and the criticality equation [START_REF] Ménard | Percolation on uniform infinite planar maps[END_REF].

The expansion for f Note that these two singular expansions and the corresponding asymptotics were established in [START_REF] Bernardi | A Boltzmann approach to percolation on random triangulations[END_REF] using different, and more involved techniques. The rational parametrization that we have simplifies this analysis a lot.

Using ∆, we get the following expansion:

∆ (1/

Figure 1 :

 1 Figure 1: Plot of the probability P p ∞ (|C| = ∞) for p ∈ [0, 1].

Figure 2 :

 2 Figure 2: Gasket decomposition: A face of the root cluster (black vertices in the figure) is filled with a triangulation with arbitrary white boundary (upper right figure, grey vertices can be either black or white) and a necklace (lower right figure) triangulating the region between the two.

Lemma 12 . 3 32

 123 Write f • (z 1 , z 2 ) = f • (1/2, t c ; z 1 , z 2 ), f ⋄ (z 1 , z 2 ) = f ⋄ (1/2, t c ; z 1 , z 2 ), z + = z + (1/2, t c ) and z ⋄ = z ⋄ (1/2, t c ). Then z + = 27 √ and z ⋄ = 3 1/4 √ 2 4

  The respective values of z + and z ⋄ are computed from[START_REF] Borot | Nesting statistics in the O(n) loop model on random maps of arbitrary topologies[END_REF] and the explicit values of the singularities y + (1/2, t c ) = 2 and y -(1/2, t c ) = -4 of the function y → T (1/2, t c , t c y) (see the Maple

- 1 1 -2 t 3 c / 2 (z 2 - 2 √ z 1 ) 2 ,

 132212 -z 1 ) + √ z + (z ⋄ -z 2 ) + O (z + -z 1 ) 2 + (z ⋄ -z 2 )Using the standard asymptotic development of hypergeometric functions at 1 we see that the first singular term in the development ofI 2/3 is -κ (z + -z 1 ) + √ z + (z ⋄ -z 2 )

  3 ), y) and T (p, t, ty) seen as series in y both have radius of convergence y + (p, t) > 1 where it is singular. Furthermore, both series can be analytically continued in the domain C \ ((-∞, y -(p, t)] ∪ [y + (p, t), +∞)), where y -(p, t) is on the negative real line (-∞ included) and is such that y -(p, t) < -y + (p, t).

  , y + (1 -p, t) , we havef • (p, t; z 1 , z 2 ) = pt 3 2z 2 + -p, t, tz) h • (z 1 , z 2 , 1/z), f ⋄ (p; z 1 , z 2 ) = pt 3 (2z 2 + z 2 1 ) +

		1	-1
	y + (p,t)
	t 3 2iπ γ	dz z T (1 t 3 c /p
		2iπ

γ dz z T (1 -p, t, tz) h ⋄ (z 1 , z 2 , 1/z).

  • (1/2, z 1 , z 2 ) is obtained similarly by replacing φ by

	and as a consequence
	t 3 c /2	-k	q k (1/2, t c ) ∼ k→∞	-2 3 5/6 Γ(-2/3)	2 -k k -5/3 .			(40)
	1 -pt 3 c z 2 -1 y 2pt 3 c z 1 -1 -t 3 c /2z 2 -1 T (1/2, t c , t c y) 2 t 3 c z 1 √ 3 2 -3 5/6 2	1 -	y 2	2/3	+	√ 2 3	1 -	1 2t 3 c z 1	1 -	y 2	-3 1/6 1 -	y 2	4/3	.

  2, z) ∼

	z→1/2	32 3 351 5 6	+	1 3 117 8 3	(1 -2z) -4/3 ,	(41)
	giving
	δ k (1/2) ∼ k→∞	1 Γ(4/3)	32 3 351 5 6	+	1 3 117 8 3	2 -k k 1/3 .	(42)
	Finally, using the expression (35), we get
	Θ (1/2, y) ∼ y→2		8 3 351 5 6	+	1 3 117 2 3	(1 -y/2) -4/3 ,	(43)
	giving
	θ k (1/2) ∼ k→∞	1 Γ(4/3)		8 3 351 5 6	+	1 3 117 2 3	2 k k 1/3 .	(44)

Acknowledgments. The author thanks Pierre Nolin for stimulating discussions on critical exponents and two anonymous referees for their careful reading of this article. This work is partially supported by the ANR grant ProGraM (Projet-ANR-19-CE40-0025) and the Labex MME-DII (ANR11-LBX-0023-01).

Lemma 8. Fix p ∈ (0, 1). The series ∆(p, z) is analytic in the domain

Proof. By definition, the series z → V c (p, z) is analytic on C \ 1 - (1-p,tc) , +∞ . In addition, if z is in this domain, the denominator of ∆ does not vanish (see the discussion in the proof of Theorem 1 after Equation [START_REF] Gorny | The geometry of a critical percolation cluster on the UIPT[END_REF]). The only statement left to prove is the equality

) .

Since we have an explicit expression for y + (p, t c ) (again, see the discussion in the proof of Theorem 1 after Equation ( 22)), it is a straightforward verification.

Lemma 9. Fix p ∈ (0, 1). For every k ⩾ 1 one has

Proof. In the whole proof, p ∈ (0, 1) is fixed. All calculations are available in the Maple companion file [START_REF]Maple companion file[END_REF]. We start by proving that, for every k ⩾ 1, the series

seen as a series in t 3 is algebraic and has a unique dominant singularity at t 3 c . In view of (9), the generating series of these modified weights is given by

Injecting this into the algebraic equation verified by T , we get an algebraic equation verified by F :

where Pol 1 and Pol 2 are explicit polynomials. The form of this equation has the following consequences. First, using Lemma 4, the series F (p, t, z) is algebraic. Second, it is the unique solution of this equation with constant term in z equal to T (1 -p, t, t)/p, and we can compute its coefficients in z inductively. These coefficients are the modified weights qk (p, t) and their expressions are then rational fractions in p, t 3 , tT 1 (p, t) and T (1 -p, t, t), whose denominator (up to a factor p) are the k-th power of

A quick glance at equation ( 2) with y = 1 shows that the last display is the derivative of the algebraic equation verified by T (1 -p, t, t) with respect to T . Therefore this quantity can only be 0 at singularities of T (1 -p, t, t), which leaves only t 3 c according to Lemma 4. As a consequence, we just proved that the series q(p, t) are all algebraic series in t 3 and all have a unique dominant singularity at t 3 c .

file [START_REF]Maple companion file[END_REF] for details):

.

Now, the change of variable

in the expressions of Lemma 11 gives

We will see that, when p = 1/2 and t = t c , the main term in the expansion of these functions stems from the singular behavior of T at y + = 2.

From Lemma 13 and the discussion that follows, we know that the function y → T (1/2, t c , t c y) is analytic for y ∈ [0, 2) and has the following expansion as y → 2 -:

.

The function

is twice differentiable on [0, 2). In addition, for ξ ∈ (0, 1) and if (z 1 , z 2 ) ∈ [0, z + ] × [0, z ⋄ ], the quantity ϕ(z 1 , z 2 ; ξ) varies in a subset of (0, (1 -1/y + (1/2, t c )) -1 ) = (0, 2). We then have

As a consequence, the function

is twice differentiable on [0, z + ] × [0, z ⋄ ] and has the following asymptotic expansion when (z 1 , z 2 ) → (z + , z ⋄ ):

The first singular term in the development of f • is then the one from -

.

The statement for f • the follows from f

Perimeter asymptotics at criticality

We gather in this Section asymptotics of several quantities appearing in this paper when p = 1/2 and t = t c . They are all consequences of the following Lemma characterizing the singularities of the function y → V (1/2, U (t 3 c ), y) defined in Lemma 4.

Lemma 13. The function

In addition it has the following asymptotic expansion in a slit neighborhood of 2:

. (36)

Proof. From Lemma 4, we already know that y → V (1/2, U (t 3 c ), y) is analytic on the domain C \ ((-∞, y -(1/2, t c )] ∪ [y + (1/2, t c ), +∞)). We also know that y ± (1/2, t c ) are the values of ŷ(1/2, U (t 3 c ), V ) at V = V ± (1/2, t c ) the two stationary points of ŷ enclosing 0. We can easily compute the corresponding values in the Maple file [START_REF]Maple companion file[END_REF]. The expansion is then easily obtained by singular inversion. This Lemma combined with the rational expressions in terms of V that we have allow to immediately compute asymptotic expansions. Indeed, the expressions T (1/2, U (t 3 c ), V ) and ∆(1/2, V ) defined respectively in Lemma 4 and Lemma 9 are singular only when V is singular, and we can get an asymptotic expansion at their unique dominant singularity by plugging the expansion of V in their expression. As usual, calculations are available in the Maple companion file [START_REF]Maple companion file[END_REF].

We get the following expansion for T :


Théo Lacour

Ruzanna Hayrapetyan

Isabelle Séverin

Isabelle Mouche

Francis Finot

Marie-Christine Chagnon

Could CTA be a reliable alternative method to identify non-genotoxic substances in Food Contact Material extract ?

à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

Introduction Cell Transformation Assay (CTA) is able to detect non-genotoxic carcinogen agents but is it able to identify Non-Intentionally Added Substances (NIAS) ?

In this work the objective is to minaturize the classical CTA. In a second step, the aim is to etablish a more specific and sensitive labeling of transformed foci.

Finally we would like to identify transcriptomics markers impacted by a non-genotoxic cancerigen agent.

The CTA using Bhas 42 cells, is a 3T3 mouse fibroblast cell line transfected with several copies of the Ras gene allowing a cell transformation process. Indeed the process of cellular transformation is possible because it is a cell already initiated. This process is described by the formation of foci characterized by :

-deep basophilic staining -spindle-shaped cells -multilayer growth (piling up of cells)

-random cell orientation -invasive cell growth into the background monolayer.

Transformed foci :

Non-transformed foci :

Could CTA be a reliable alternative method to identify non-genotoxic substances in Food Contact Material extract ? Identification of the foci is done by an experienced experimenter. This staining allows an easier identification of the foci because it is realized under fluorescence.

In parallel, a quantitative examination can be performed. The fluorescent properties of BMVC allow a spectrophotometric determination with excitation at 435 nm and emission at 560 nm.

Unpaired t test with Welch's correction, **** : P value<0,0001

Results

A transformics assay allows the examination of the transcriptome to target which steps of the cell cycle are impacted by the test product.

Previous publications illustrate its interest:

•HAPs (2, 4-diaminotoluene, benzo(a)pyrene, 2-acetylaminoflourene, or 3methycholanthrene) : ➢14 genes regulated in common for the 24h and 32 days groups which shows the predictivity of this test

•38 chemicals (tumor promoters, non-tumor promoters, genotoxic carcinogens and food additives) : ➢The cells are treated for 48 hours, resulting in only a list of genes commonly known to be involved in human cancers.

Conclusion

•The miniaturized 96-well CTA identifies genotoxic and non-genotoxic carcinogens using fewer cells, less medium and, most importantly, less test material.

•Labeling of BMVC on Bhas 42 cells combined with CTA would allow for more accurate and sensitive quantification, in addition to allowing for a non-examiner-dependent quantification.

•The transcriptome analysis allowed to identify the impacted cell cycle steps, such as immune response, cell adhesion, apoptosis and cell proliferation.

Outlook

❖Automation of rinsing and identification of foci using a microplate washer and a microplate reader.

❖Use BMVC labeling to identify the potential non-genotoxic carcinogenicity of NIAS. Non-intentionally added substances (NIAS) are chemical compounds that are present in food contact materials (FCMs) and could therefore migrate into food, but they are not added for a technical reason during the production process. Often their presence is not known by the consumer and not even by the producer.

T.

  Lacour 1,2 and R. Hayrapetyan 3 , I. Severin 3 , I. Mouche 1 , F. Finot 1 , M.-C. Chagnon 3 1 GenEvolutioN, Porcheville, France 2 Ecole Pratique des Hautes Etudes, Paris, France 3 Derttech "Packtox", NUTOX, INSERM U1231, Dijon, France

  ❖Reduction of SVF variability by using a serum substitute, non-fetal serum or by reducing the proportion of serum in the medium. ❖Identification of non-genotoxic (epigenetic) damage with hypermethylation or hypomethylation as well as with different DNA compactions • Maeshima H, Ohno K, Tanaka-Azuma Y, Nakano S, Yamada T. Identification of tumor promotion marker genes for predicting tumor promoting potential of chemicals in BALB/c 3T3 cells. Toxicol In Vitro. 2009 Feb;23(1):148-57. doi: 10.1016/j.tiv.2008.10.005. Epub 2008 Nov 1. PMID: 19000923. • Mascolo MG, Perdichizzi S, Vaccari M, Rotondo F, Zanzi C, Grilli S, Paparella M, Jacobs MN, Colacci A. The transformics assay: first steps for the development of an integrated approach to investigate the malignant cell transformation in vitro. Carcinogenesis. 2018 Jul 3;39(7):955-967. doi: 10.1093/carcin/bgy037. Erratum in: Carcinogenesis. 2018 Jul 3;39(7):968. PMID: 29554273; PMCID: PMC6031005. • Rohrbeck A, Salinas G, Maaser K, Linge J, Salovaara S, Corvi R, Borlak J. Toxicogenomics applied to in vitro carcinogenicity testing with Balb/c 3T3 cells revealed a gene signature predictive of chemical carcinogens. Toxicol Sci. 2010 Nov;118(1):31-41. doi: 10.1093/toxsci/kfq246. Epub 2010 Aug 16. PMID: 20713471. LACOUR, M.Sc Student, GenEvolutioN Innovation Award 2018, theo.lacour@genevolution.fr