Olivier Marchal
email: 1.olivier.marchal@univ-st-etienne.fr2.nicolas.orantin@gmail.com

Nicolas Orantin

Isomonodromic deformations of a rational differential system and reconstruction with the topological recursion: the sl 2 case

In this paper, we show that it is always possible to deform a differential equation

by introducing a small formal parameter in such a way that it satisfies the Topological Type properties of Bergère, Borot and Eynard. This is obtained by including the former differential equation in an isomonodromic system and using some homogeneity conditions to introduce . The topological recursion is then proved to provide a formal series expansion of the corresponding tau-function whose coefficients can thus be expressed in terms of intersections of tautological classes in the Deligne-Mumford compactification of the moduli space of surfaces. We present a few examples including any Fuchsian system of sl 2 (C)(x) as well as some elements of Painlevé hierarchies.

Introduction and summary of results

The theory of isomonodromic systems has attracted a renewed attention in the recent years through its interaction with different topics in mathematics and physics such as Conformal Field Theories (CFTs) [START_REF] Harnad | Quantum isomonodromic deformations and the Knizhnik-Zamolodchikov equations[END_REF][START_REF] Iorgov | Isomonodromic Tau-Functions from Liouville Conformal Blocks[END_REF][START_REF] Reshetikhin | The Knizhnik-Zamolodchikov System as a Deformation of the Isomonodromy Problem[END_REF], random matrix models [START_REF] Bertola | The partition function of the two-matrix model as an isomonodromic tau-function[END_REF], N = 2 SUSY theories [START_REF] Cecotti | Ising Model and N=2 Supersymmetric Theories[END_REF], Hitchin integrable systems and moduli spaces of flat connections on a Riemann surface [START_REF] Boalch | Symplectic Manifolds and Isomonodromic Deformations[END_REF][START_REF] Teschner | Quantization of the Hitchin moduli spaces, Liouville theory, and the geometric Langlands correspondence I[END_REF] or Nakajima quiver varieties [START_REF] Hiroe | Moduli spaces of meromorphic connections, quiver varieties, and integrable deformations[END_REF]. In particular, according to the isomonodromic-CFT correspondence, the tau function of an isomonodromic system should be equal to a conformal block of the corresponding theory, i.e. a section of the associated Verlinde bundle seen as a wave function obtained by geometric quantization of a moduli space of connections on the Riemann surface considered. From a physics point of view, this correspondence states that the tau function is a building block for correlation functions of the associated CFT.

On the other hand, the Topological Recursion (TR) [START_REF] Eynard | Invariants of algebraic curves and topological expansion[END_REF], a formalism originally developed for solving random matrix models [START_REF] Chekhov | Free energy topological expansion for the 2-matrix model[END_REF], proved to be another useful tool for studying different problems in mathematics in a unified way, allowing to bring together different fields of mathematics and physics [START_REF] Borot | Lecture notes on topological recursion and geometry[END_REF]. After it was proved that the topological recursion provides an efficient way of computing Gromov-Witten invariants of manifolds with semi-simple cohomology [START_REF] Dunin-Barkowski | Identification of the Givental formula with the spectral curve topological recursion procedure[END_REF][START_REF] Norbury | Gromov-Witten invariants of P 1 and Eynard-Orantin invariants[END_REF], one of the most active and mysterious field of applications lies in the theory of "quantum curves" where the recursion is expected to give a new way to compute the WKB expansion of a solution Ψ(x, ) of a given differential equation with respect to a formal parameter [START_REF] Dijkgraaf | Quantum Curves and D-Modules[END_REF][START_REF] Dumitrescu | Quantum Curves for Hitchin Fibrations and the Eynard-Orantin Theory[END_REF][START_REF] Dumitrescu | The spectral curve of the Eynard-Orantin recursion via the Laplace transform[END_REF][START_REF] Iwaki | Quantum Curve and the First Painlevé Equation[END_REF][START_REF] Norbury | Quantum curves and topological recursion[END_REF][START_REF] Schwarz | Quantum Curves[END_REF]. In particular, if this differential equation comes from the quantization of the A-polynomial of a given knot, the formal series expansion of this solution is conjectured to provide a generating series of some associated knot invariants [START_REF] Borot | All-order asymptotics of hyperbolic knot invariants from non-perturbative topological recursion of A-polynomials[END_REF]. The latter wave function can also be thought of as obtained by geometric quantization of some moduli space of connections on the boundary of a tubular neighborhood of the knot. On the other hand, when the differential equation defined by the quantum curve belongs to some isomonodromic systems, the partition function Z = exp ∞ g=0 2g-2 F g is equal to an associated isomonodromic tau function.

However, not every differential equation gives rise to WKB solutions computed by TR. In [START_REF] Bergère | Determinantal formulae and loop equations[END_REF][START_REF] Bergère | Rational differential systems, loop equations, and application to the q th reductions of KP[END_REF], it was proved that if the -dependent differential system satisfies some conditions, referred to as the Topological Type (TT) property, then the associated partition function has an expansion in whose coefficients can be computed by TR. A good understanding of the TT property is still missing but it could be proved that some interesting examples satisfy it. Among these examples, the simplest isomonodromic systems, i.e. the six Painlevé equations, were proved to be deformed by the introduction of a formal parameter so that they satisfy the TT property [START_REF] Iwaki | Painlevé 2 equation with arbitrary monodromy parameter, topological recursion and determinantal formulas[END_REF][START_REF] Iwaki | Painlevé equations, topological type property and reconstruction by the topological recursion[END_REF]. In the regime → 0, the associated tau functions thus have a formal expansion given by TR. By the results of [START_REF] Bertola | Correlation functions of the KdV hierarchy and applications to intersection numbers over M g,n[END_REF][START_REF] Eynard | Invariants of spectral curves and intersection theory of moduli spaces of complex curves[END_REF], the coefficients of this formal expansion can thus be expressed in terms of intersection of tautological classes in the moduli space of Riemann surfaces generalizing Kontsevich-Witten theorem [START_REF] Kontsevich | Intersection theory on the moduli space of curves and the matrix Airy function[END_REF] obtained in the study of the Airy equation where the partition function is at the same time a KdV tau function and a generating function for Gromov-Witten invariants of the point.

For the understanding of the relation between integrable and isomonodromic systems, TR, knot invariants and CFTs, it is important to understand if this result can be generalized to other isomonodromic systems, namely, if one can deform a rational linear differential system in such a way that it satisfies the TT property. In this article, we prove that it is always possible to deform a rational linear system

d dx Ψ(x) = L(x)Ψ(x) → d dx Ψ(x, ) = L(x, )Ψ(x, )
with L(x) ∈ sl 2 (C)(x), in such a way that -L(x, ) is equal to L(x) for = 1.

-The deformed equation d dx Ψ(x, ) = L(x, )Ψ(x, ) satisfies the TT property. For doing so, as explained in section 3, we consider any generic sl 2 (C) valued rational function L(x) as an element of a coadjoint orbit in a finite dimensional subspace of the sl 2 (C) loop algebra. In other words, after reduction, we consider it as a point in a symplectic space equipped with the corresponding Hitchin integrable system. We remind that this isospectral system has an autonomous Hamiltonian representation before de-autonomizing it to make it an isomonodromic system equipped with a non-autonomous Hamiltonian representation. This is done in section 4 where we also give a prescription to introduce a formal parameter giving rise to a family of isomonodromic systems deforming the original one. Finally, in section 5, we state in Theorem 5.1 our main result which is that the family of systems built in this way satisfies the TT property. Hence, we prove that one can always deform a linear differential system characterized by L(x) ∈ sl 2 (C)(x) in such a way that the associated tau function has an expansion computed by TR. This generalizes the results known up to now on various examples in the context of isomonodromic systems.

Moreover, our result allows to generalize the quantization procedure of [START_REF] Bouchard | Reconstructing WKB from topological recursion[END_REF]. Indeed, to any hyperelliptic curve of genus 0, possibly with interior points in its Newton polytope, we can associate a linear system depending on with the TT property. It means that one can build a quantum curve for any such classical hyperelliptic curve.

In section 6, we present a few examples including the Painlevé equations as well as rank 2 Fuchsian systems and arbitrary sl 2 (C) valued polynomials, the latter including some of the Painlevé hierarchies.

Since our approach uses mainly the Poisson structure of the loop algebra, we believe that it represents a first step towards a generalization to any isomonodromic system. It also allows to make the link with more geometrical approaches of the isomonodromic deformations as explained in section 7. In particular, we believe that the approach developed in this article should extend naturally to L(x) ∈ gl N (C)(x) with N ≥ 2 and to L(x) ∈ g(x) where g is any complex simple Lie subalgebra of gl N (C).

Topological recursion and WKB analysis 2.1 Compatible systems, correlators and loop equations

In this article, we are interested in studying compatible systems of the form

∂ ∂x Ψ(x, t) = L(x, t)Ψ(x, t) ∂ ∂t Ψ(x, t) = A(x, t)Ψ(x, t)
2 This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.

PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0002260

where L(x, t) and A(x, t) belong to a simple Lie algebra g ⊂ M N (C). The general setup maybe extended to an arbitrary semi-simple Lie algebra equipped with a faithful representation ρ used to pushforward quantities in M N (C) (See [START_REF] Belliard | Loop Equations from Differential Systems on Curves[END_REF] for details). Moreover we assume that for all t ∈ C, the functions x → L(x, t) and x → A(x, t) are rational functions of x ∈ Σ = P 1 . Note that loop equations and correlators that will be defined below have been generalized in [START_REF] Belliard | Loop Equations from Differential Systems on Curves[END_REF] for an arbitrary Riemann surface Σ. However for our purpose, we will restrict only to Σ = P 1 in the present article. In particular, since P 1 is a genus 0 Riemann surface, its fundamental group and its corresponding action are trivial and this simplifies some definitions and properties presented in [START_REF] Belliard | Loop Equations from Differential Systems on Curves[END_REF].

To any solution Ψ(x) of a linear differential equation

∂ ∂x Ψ(x) = L(x)Ψ(x)
where L(x) is a rational function of x with values in a Lie algebra g ⊂ M N (C), Eynard et al. [START_REF] Belliard | Integrable differential systems of topological type and reconstruction by the topological recursion[END_REF][START_REF] Belliard | Loop Equations from Differential Systems on Curves[END_REF][START_REF] Bergère | Determinantal formulae and loop equations[END_REF][START_REF] Bergère | Rational differential systems, loop equations, and application to the q th reductions of KP[END_REF] associate a set of elements M (x; E) := Ψ(x)EΨ(x) -1 where E ∈ h is an arbitrary element of the Cartan subalgebra h ⊂ g. Note that since the base field C is algebraically closed and the Lie algebra is finite dimensional then all Cartan subalgebras are conjugate under automorphisms of the Lie algebra, and thus are all isomorphic. These elements satisfy the equation

∂ ∂x M (x, E) = [L(x), M (x, E)] .
Motivated by the correlation functions in random matrix models and CFT's, one can define a set of correlators associated to such a linear differential equation (Cf. [START_REF] Belliard | Integrable differential systems of topological type and reconstruction by the topological recursion[END_REF][START_REF] Belliard | Loop Equations from Differential Systems on Curves[END_REF][START_REF] Bergère | Determinantal formulae and loop equations[END_REF][START_REF] Bergère | Rational differential systems, loop equations, and application to the q th reductions of KP[END_REF]).

Definition 2.1 (Correlation functions).

For n ≥ 1, let the n-point correlators W n : P 1 ⊗ h n → P 1 be defined by

Wn(x1⊗E1, . . . , xn⊗En) :=        Tr [L(x1)M (x1; E1)] dx1 if n = 1 (-1) n n σ∈Sn Tr[M (x σ(1) ;E σ(1) )M (x σ(2) ;E σ(2) )...M (x σ(n) ;E σ(n) )] n i=1 (x σ(i) -x σ(i+1) ) dx1 . . . dxn if n ≥ 2
where S n is the symmetric group of n elements.

In [START_REF] Belliard | Loop Equations from Differential Systems on Curves[END_REF] (Theorem 4.3) it was proved that they satisfy a set of loop equations, or W-constraints. These loop equations greatly simplify in the case of g = sl N (C) presented in [START_REF] Bergère | Rational differential systems, loop equations, and application to the q th reductions of KP[END_REF] and even more in the case where g = sl 2 (C) developed earlier in [START_REF] Bergère | Determinantal formulae and loop equations[END_REF]. Since we will not make direct use of the loop equations in this article, we just refer the reader to Definition 2.6 and Theorem 2.7 of [START_REF] Bergère | Determinantal formulae and loop equations[END_REF] for loop equations arising in g = sl 2 (C).

Let us just mention that generically, this set of loop equations does not have a unique solution. However, restricting the space of solutions by requiring additional properties for the correlators, known as the "Topological Type property", allows to prove uniqueness as well as solve explicitly this set of equations. Note that the correlators, the loop equations and the topological type property are defined for any rational linear differential system ∂ x Ψ(x) = L(x)Ψ(x). In particular there is no need for a second compatible time differential equations to define them. On the other hand, as we show in this paper, the existence of such a compatible equation allows to naturally prove the topological type property which is usually hard to obtain otherwise.

From the point of view of isomonodromic systems, i.e. when the differential equation in x belongs to a compatible system

∂ ∂x Ψ(x, t) = L(x, t)Ψ(x, t) ∂ ∂t Ψ(x, t) = A(x, t)Ψ(x, t)
, one of the main motivations for the introduction of these correlators is the fact that they give access to the associated isomonodromic tau function.

3 This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.

PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0002260

Generalizing the definition of Jimbo, Miwa and Ueno [START_REF] Jimbo | Monodromy preserving deformation of linear ordinary differential equations with rational coefficients I. General theory and τ -function[END_REF][START_REF] Jimbo | Monodromy preserving deformation of linear ordinary differential equations with rational coefficients II[END_REF], Bertola and the first author [START_REF] Bertola | The partition function of the two-matrix model as an isomonodromic tau-function[END_REF] defined a tau function τ BM associated to such a system by identifying its variations with respect to the isomonodromic time t. A similar definition was given by [START_REF] Cafasso | Tau functions as Widom constants[END_REF] for an arbitrary Riemann-Hilbert problem. We believe that these definitions are equivalent whenever the RH problem comes from an isomonodromic system. In general, as we shall see in the following, one could have more than one isomonodromic time and the following formula should be valid for any of them. The existence of a solution to the corresponding set of equations was proven in [START_REF] Bertola | The partition function of the two-matrix model as an isomonodromic tau-function[END_REF].

Let us now recall how one can define such a tau function associated to an isomonodromic system. We shall consider the case g = sl N for N ≥ 2 for this purpose.

Let Ψ(x) be solution to an isomonodormic system and

Ψ(x) = Ψ -(x)x S Ψ + (x) be a Birkhoff factorization where Ψ -(x) ∈ x -1 SL N [[x -1 ]] (resp. Ψ + (x) ∈ SL N [[x]]
) and S ∈ GL N is independent of t. One can then define the form [START_REF] Bertola | The partition function of the two-matrix model as an isomonodromic tau-function[END_REF] ω

BM := Res x→∞ Tr x -S Ψ -(x) -1 dΨ -(x)x S dx dΨ + (x)Ψ + (x) -1
where dΨ + (x) denotes the differential of Ψ + (x) in the space of isomonodromic parameters. It was proved in [START_REF] Bertola | The partition function of the two-matrix model as an isomonodromic tau-function[END_REF] that Theorem 2.1. The differential ω BM is closed and locally defines a tau function by

d ln τ BM = ω BM .
In [START_REF] Bergère | Rational differential systems, loop equations, and application to the q th reductions of KP[END_REF], Bergère, Borot and Eynard remarked that this tau function can be recovered from the knowledge of W 1 by integrating this definition with respect to the isomonodromic time t:

∂ t ln τ BM = x∈C W 1 (x)f (x)dx
where C (resp. f (x)) is a contour (resp. a function) to be specified depending on the system considered.

Topological type property and topological recursion

Let us now assume that the whole isomonodromic system depends on a formal parameter , taking the form

∂ ∂x Ψ(x, t, ) = L(x, t, )Ψ(x, t, ) ∂ ∂t Ψ(x, t, ) = A(x, t, )Ψ(x, t, ) . (2.1)
Under some assumptions, the loop equations then have a unique solution in the space of formal series in and this solution can be computed order by order by the topological recursion procedure of [START_REF] Eynard | Invariants of algebraic curves and topological expansion[END_REF]. Furthermore, it gives a procedure to compute τ BM order by order in up to multiplicative factors independent of t (the BM tau-function is defined up a global multiplicative constant). Using general properties of the topological recursion, the coefficients of this formal series expansion can be expressed in terms of intersection of cohomology classes in the Deligne-Mumford compactification of the moduli space of Riemann surfaces. We first recall the definition of a global classical spectral curve suitable for our purposes and results about the topological recursion developed in [START_REF] Eynard | Invariants of algebraic curves and topological expansion[END_REF]. Note that the topological recursion has known many generalizations after [START_REF] Eynard | Invariants of algebraic curves and topological expansion[END_REF]. In particular, the spectral curve may be replaced by some local data consisting of germs of differential forms [START_REF] Dunin-Barkowski | Identification of the Givental formula with the spectral curve topological recursion procedure[END_REF] but we only need the original simple version in this article.

Definition 2.2 (Classical spectral curve).

A classical spectral curve is composed of the following data. 1. A compact Riemann surface Σ cl equipped with a meromorphic function x : Σ cl → P 1 such that the form dx has simple zeroes at the branchpoints {b i } 1≤i≤m defined by dx(b i ) = 0. The degree of the map

x : Σ cl → P 1 z → x(z)
4 This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.

PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0002260

is denoted by N . Thus, for any generic x 0 ∈ P 1 , we may define the set of preimages x -1 (x 0 ) = {z 1 (x 0 ), . . . , z N (x 0 )} by choosing a labeling of the sheets of the cover. 2. A meromorphic form ω 0,1 on Σ cl holomorphic at the branch points {b i } 1≤i≤m .

3. A symmetric bidifferential ω 0,2 on Σ cl × Σ cl with only double poles on the diagonal without residues.

In the rest of the paper, we denote the set of branchpoints by B = {b i } 1≤i≤m . The number of branchpoints is denoted m = |B| ≥ 1.

We now recall that the topological recursion developed in [START_REF] Eynard | Invariants of algebraic curves and topological expansion[END_REF] may be applied to any classical spectral curve under some genericity assumptions which we shall assume in the following.

Definition 2.3 (Topological Recursion).

For any classical spectral curve (Σ cl , x, ω 0,1 , ω 0,2 ), one may construct recursively symmetric n-forms (ω g,n (z 1 , . . . , z n )) n≥1,g≥0 on Σ cl n called Eynard-Orantin differentials and some symplectic invariants F (g) g≥0 using the topological recursion presented in [START_REF] Eynard | Invariants of algebraic curves and topological expansion[END_REF].

Formulas defining the topological recursion can be found in [START_REF] Eynard | Invariants of algebraic curves and topological expansion[END_REF]. Let us just remind that the recursion is done over the index p = 2g +n starting from ω 0,1 and ω 0,2 that are initial data coming from the definition of the classical spectral curve. The formulas involve residues at the branchpoints and it is known that the Eynard-Orantin differentials produced by the topological recursion are symmetric and may only have poles at the branchpoints (except ω 0,2 ).

In order to apply the topological recursion procedure to the study of the asymptotics of a linear differential equation, the latter must satisfy the so-called Topological type properties. Since we only need to work with reductions of gl N (C), we consider the version of topological type properties for this algebra. Hence one can consider a basis of the corresponding Cartan sub-algebra given by diagonal matrices (e i ) N i=1 with only non-vanishing element [e i ] j,k := δ ij δ jk .

Definition 2.4 (Topological Type Property). A system ∂ x Ψ(x, ) = L(x, )Ψ(x, ), with a N × N matrix L(x, ), is said to satisfy the Topological Type property (TT property for short) if the following conditions are met.

1. Genus 0 property: The Riemann surface Σ cl given by the compactification of the → 0 limit of the characteristic polynomial lim →0 det(yI d -L(x, t)) = 0} ⊂ C 2 has genus 0.

Formal expansion in :

The correlators admit a formal series expansion in .

W 1 (x 1 ⊗ e a1 ) = ∞ k=-1 k W (k) 1 (x 1 ⊗ e a1 ) := ∞ k=-1 k W (k) 1 (z a1 (x 1 )) W n (x 1 ⊗ e a1 , . . . , x n ⊗ e an ) = ∞ k=0 k W (k) n (x 1 ⊗ e a1 , . . . , x n ⊗ e an ) := ∞ k=0 k W (k) n (z a1 (x 1 ), . . . , z an (x n )) , ∀ n ≥ 2,
whose coefficients

W (k) n (z 1 , . . . , z n )dx(z 1 ) . . . dx(z n ) define symmetric n-forms on Σ cl .

Parity property: As functions of , the correlation functions have the parity property

W n (x 1 ⊗ E a1 , . . . , x n ⊗ E an | -) = (-1) n W n (x 1 ⊗ E a1 , . . . , x n ⊗ E an | ) , ∀ n ≥ 1.

Pole property:

The n-forms W (k) n (z 1 , . . . , z n ) with n ≥ 1 and k ≥ 0 and (n, k) / ∈ {(1, 0), (2, 0)} may only have singularities at the branchpoints of the classical spectral curve. Moreover these singularities may only be pole singularities. In other words, they are meromorphic functions on Σ cl with only pole singularities at the branchpoints.

5

This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset. Wg,n (z 1 , . . . , z n ) n-2-2g .

PLEASE CITE THIS ARTICLE AS

Identification of the two points function:

The 2-form W0,2 (z 1 , z 2 )dx(z 1 )dx(z 2 ) defines an admissible symmetric form ω 0,2 for the topological recursion.

Then we have the following theorem [START_REF] Belliard | Integrable differential systems of topological type and reconstruction by the topological recursion[END_REF][START_REF] Bergère | Rational differential systems, loop equations, and application to the q th reductions of KP[END_REF].

Theorem 2.2 (Reconstruction with the Topological Recursion). If the differential system ∂ x Ψ(x) = L(x)Ψ(x) satisfies the TT property and dx has only simple zeros on the classical spectral curve Σ cl then the coefficients of the correlation functions Wg,n (z 1 , . . . , z n ) n≥1,g≥0 identify with the corresponding Eynard-Orantin differentials (ω g,n ) n≥1,g≥0 computed from the application of the topological recursion of [START_REF] Eynard | Invariants of algebraic curves and topological expansion[END_REF] to the classical spectral curve of the system with initial data given by ω 1,0 (z) = W0,1 (z)dx(z) and ω 2,0 (z

1 , z 2 ) = W0,2 (z 1 , z 2 )dx(z 1 )dx(z 2 ), ∀ (z 1 , . . . , z n ) ∈ Σ cl n : Wg,n (z 1 , . . . , z n )dx(z 1 ) . . . dx(z n ) = ω g,n (z 1 , . . . , z n ).
Moreover, if the differential system comes from an isomonodromic system of the form (2.1) then the tau-function ln τ BM admits a formal series expansion in of the form

ln τ BM (t) = +∞ g=0 τ (g) (t) 2g-2
where d dt τ (g) = d dt F (g) and F (g) g≥0 are the so-called symplectic invariant generated by the topological recursion applied to the classical spectral curve.

In the following sections, we shall explain how to build generic sl 2 (C)-isomonodromic system and prove that they satisfy the TT property. This allows to compute a formal -expansion of the tau function and a formal expansion of the correlators by application of the topological recursion. We believe that the construction presented below should extend to sl N (C) and to any simple Lie algebra over C (i.e (A n ) n≥1 , (B n ) n≥1 , (C n ) n≥1 , (D n ) n≥2 and possibly the exceptional cases E 6 , E 7 , E 8 , F 4 and G 2 ) but this is left for future works.

Remark 2.1. sl2 is a rank one reduction of gl 2 obtained by vanishing of the trace. From the identification of the correlation functions with functions on the classical spectral curve, this means that one can work with the usual result for gl 2 by imposing the vanishing of the sum of the values of the functions considered over the two sheets of the classical spectral curve.

In particular, one only needs to prove the TT properties for the element e1 of the gl 2 Cartan subalgebra. Indeed, values for e2 are opposite to tgose for e1 while values for the basis vector e1 -e2 of the sl2 Cartan subalgebra are equal to twice the corresponding values for e1.

SL N isospectral deformation, integrable system and Poisson structure

The construction of isomonodromic systems admitting a WKB solution built by topological recursion relies heavily on the Poisson and symplectic structures on the space of flat connections on a Riemann surface as well as the associated Hitchin integrable system. Let us start by reviewing this structure from the R-matrix point of view [START_REF] Semenov-Tian-Shansky | What is a classical R-matrix?[END_REF]. 6 This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset. PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0002260

Isospectral systems and Poisson structure

Let g be a Lie algebra, and C be an oriented simple closed curve on P 1 . Let us denote by U + (resp. U -) the region outside C (resp. inside C).

One can associate to such a pair a loop space g of smooth maps L : C → g together with a polarization g = g+ ⊕ gwhere g-(resp. g+ ) is the subalgebra of maps admitting an holomorphic extension to U -(resp. to U + ).

One can define an Ad-invariant inner product < •, • >: g × g → C by

∀ (L 1 , L 2 ) ∈ g × g , L 1 , L 2 := 1 2πi x∈C Tr [L 1 (x) • L 2 (x)] dx.
This allows the identification of g with its dual g * in such a way that the g * ± can be identified with g∓ .

The exponentiated group G * acts by coadjoint action through

∀ (f, g) ∈ g * × g , ∀ X ∈ G * , Ad * X (f )(g) = 1 2πi x∈C Tr ([X, f ]g) .
The polarization, which depends on the choice of a curve C, also allows to define a second Lie algebra structure on the infinite dimensional space g in addition to the one inherited from g. Using the classical R-matrix construction, one defines the bracket

∀ (L 1 , L 2 ) ∈ g × g, [L 1 , L 2 ] R := [R(L 1 ), L 2 ] + [L 1 , R(L 2 )]
where

R := 1 2 (P + -P -)
with P ± being the projection operator to U ± . This defines a Lie-Poisson structure on g * through the bracket

∀ (f, g) ∈ g * × g * , {f, g} R (µ) := µ, [df (µ), dg(µ)] R .
Let us denote by I the set of spectral invariants, i.e. the set of Ad * -invariant polynomials on g * . The generalized Adler-Kostant-Symes (AKS) theorem in its R-matrix form [START_REF] Semenov-Tian-Shansky | What is a classical R-matrix?[END_REF] states that the elements of I Poisson commute. In addition, it provides the form of the Hamilton's equations. For H ∈ I, the corresponding equation reads (see for example [START_REF] Harnad | Moment Maps to Loop Algebras, Classical R-Matrix and Integrable Systems[END_REF])

dL dt = [P σ (dH), L].
where

P σ := 1 2 [(1 + σ)P + + (σ -1)P -] (3.1)
for any σ ∈ C and dH denotes the differential of H evaluated at L where the differential of a function f is defined by the first order variation

f (Y + Ad * X Y ) = f (Y ) + Ad * X Y (df ) for any X ∈ G * .
Remark 3.1. These Hamilton's equations are independent of σ. Indeed, since P+ + P-= Id:

Pσ = P0 + σ 2 Id = R + σ 2 Id.
Since dH(L) commutes with L for any element of H ∈ I, one has

dL dt = [Pσ(dH), L] = [R(dH), L] + σ 2 [dH, L] = [R(dH), L]
for any σ.

These equations give a set of isospectral deformations of a given operator L ∈ g, i.e. they imply that the spectrum of L(x) is independent of t.

7

This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.

PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0002260

Remark 3.2. It is important to notice that, even though Hamilton's equations are independent of the value of σ, the auxiliary matrix AH := Pσ(dH) is not! Hence, σ can be used to fine tune the shape of this auxiliary matrix. All the procedure developed here depends explicitly on σ. However, for simplicity of notations, we shall not write the dependance of AH in σ explicitly.

From these equations, one can identify a set of Casimir functions and Hamiltonians of the system considered.

sl N case, finite dimensional subspaces and generators of I

From now on, for simplicity, we shall consider C as a small contour encircling x = ∞. This allows to identify g+ (resp. g-) with elements of g

[[x]] (resp. x -1 g[[x -1 ]]
). We also restrict to the case where g = sl N for N ≥ 2. With this specialization, the ring of invariants I is generated by the coefficients of the expansion of Tr L(x

) k for 2 ≤ k ≤ N , namely ∀ l ∈ Z , ∀ 2 ≤ k ≤ N , h k,l := Res x→∞ x -l-1 Tr L(x) k dx = L(x) k , x -l-1 .
We shall now consider finite dimensional subspaces of g * . Namely, let us fix n ≥ 0, a set of points {a 1 , . . . , a n } together with a 0 := ∞ as well as integers (r ν ) 1≤ν≤n ∈ N n and r 0 ≥ -1. We are finally given a fixed element L 0,r0 ∈ g.

Out of these data, denoting r := n ν=0 r ν , let us define the subspace ĝ * of g * composed of elements with poles of order at most r ν at a ν by -If

r 0 ≥ 0, ĝ * := L(x) := r0 i=0 L 0,i x i + n ν=1 rν i=1 L ν,i (x -a ν ) i , (L ν,k ) ∈ g r -If r 0 = -1, ĝ * := L(x) := n ν=1 rν i=1 L ν,i (x -a ν ) i , (L ν,k ) ∈ g r | - n ν=1 L ν,1 = L 0,-1
where the residue at infinity, -

n ν=1 L ν,1 = L 0,-1 , is fixed.
Let us now define a set of generators of Î, the Ad * -invariant functions on ĝ * . They can be obtained by considering the coefficients of the partial fraction decomposition of each of the terms Tr L(x) k to get

Tr L(x) k = k r0 i=0 H k,0,i x i + n ν=1 k rν i=1 H k,ν,i (x -a ν ) i .
Note that as we shall see, H k,0,kr0 might be vanishing if L 0,r0 is not full rank.

Thus, a set of generators of Î is given by the coefficients of these Ad * -invariant functions defined by

∀ 2 ≤ k ≤ N , 1 ≤ ν ≤ n , 1 ≤ i ≤ kr ν , H k,ν,i := - 1 2πi ∞ µ =ν (x -a µ ) rµ µ =ν (a ν -a µ ) rµ (x -a ν ) i-1 Tr L(x) k dx and ∀ 2 ≤ k ≤ N , 0 ≤ i ≤ kr 0 , H k,0,i := 1 2πi ∞ x -i-1 Tr L(x) k dx.
Indeed,

µ =ν (x-aµ) rµ µ =ν (aν -aµ) rµ (x -a ν ) i-1
Tr L(x) k dx has poles only at ∞ and a ν . Thus one can move the integration contour and get

H k,ν,i = 1 2πi aν µ =ν (x -a µ ) rµ µ =ν (a ν -a µ ) rµ (x -a ν ) i-1 Tr L(x) k dx
8 This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.

PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0002260

= - 1 2πi ∞ µ =ν (x -a µ ) rµ µ =ν (a ν -a µ ) rµ (x -a ν ) i-1 Tr L(x) k dx. (3.2)
The associated Hamilton's equations read

dL(x) dt k,ν,i = [A k,ν,i , L(x)]
where

A k,ν,i =      k x -i-1 (L(x)) k-1 + = k Res x ′ →∞ x ′-i-1 (L(x ′ )) k-1 dx ′ (x ′ -x) if ν = 0 -k µ =ν (x-aµ) rµ µ =ν (aν -aµ) rµ (x -a ν ) i-1 (L(x)) k-1 + otherwise.
For ν = 0, one can move the integration contour to get an expression localized around a ν ,

A k,ν,i (x) = -k Res x ′ →∞ µ =ν (x ′ -a µ ) rµ µ =ν (a ν -a µ ) rµ (x ′ -a ν ) i-1 (L(x ′ )) k-1 dx ′ x ′ -x = k µ =ν (x -a µ ) rµ µ =ν (a ν -a µ ) rµ (x -a ν ) i-1 (L(x)) k-1 +k Res x ′ →aν µ =ν (x ′ -a µ ) rµ µ =ν (a ν -a µ ) rµ (x ′ -a ν ) i-1 (L(x ′ )) k-1 dx ′ x ′ -x = k µ =ν (x -a µ ) rµ µ =ν (a ν -a µ ) rµ (x -a ν ) i-1 (L(x)) k-1 + k (x -a ν ) i-1 (L(x)) k-1 -,aν
where the notation [f (x)] -,aν stands for the polar part of f (x) at x = a ν . In addition, one can add to A k,ν,i any matrix commuting with L(x) without changing Hamilton's equations. Hence, we may remove the first contribution in the previous formula and get Hamilton's equations of the form

dL(x) dt k,ν,i = Ãk,ν,i , L(x) where Ãk,ν,i = k (x -a ν ) i-1 (L(x)) k-1 -,aν
is the projection of (xa ν ) i-1 (L(x)) k-1 to its polar part at a ν .

Remark 3.3. For convenience, in some cases, it is useful to shift this auxiliary matrix by a matrix proportional to L(x), i.e. taking σ = 0 in eq. (3.1). This is for example the case for recovering Painlevé IV equation.

Casimirs and symplectic leaves

In order to describe the underlying integrable system, let us now restrict to symplectic leaves of ĝ * which are given by coadjoint orbits by first identifying a set of Casimir functions.

Lemma 3.1. For any

ν = 0, H k,ν,i is a Casimir for r ν (k -1) + 1 ≤ i ≤ r ν k while H k,0,i is a Casimir for r 0 (k -1) ≤ i ≤ r 0 k
Proof. This follows from a direct computation. Indeed, for ν = 0, (L(x))

k-1 being a polynomial of degree r ν (k -1) in (xa ν ) -1 , A k,ν,i is non-vanishing only for 1 ≤ i ≤ r ν (k -1). In the same way, A k,0,i gives non-vanishing contributions for 0 ≤ i ≤ r 0 (k -1) -1. 9 This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.

PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0002260

By AKS theorem, the functions H k,ν,i are in involution. Thanks to the preceding lemma, one can check that the number of non-Casimir Hamiltonians is equal to

N k=2 r 0 (k -1) + n ν=1 r ν (k -1) = r N (N -1) 2
which matches with half the dimension of a generic symplectic leaf as expected for this integrable system. Note that this is also true in the case r 0 = -1 since the residues L ν,1 are no longer independent.

Remark 3.4. The point at infinity seems to play a specific role in this presentation. It indeed does since the contour C separates it from all other points. A different choice of contour would give rise to a different loop space and a different presentation of the same problem. It is expected that all the results should be independent of this choice of contour. This might be a manifestation of some invariance under the mapping class group of the surface considered.

Reduction

The symplectic leafs described above are too large for our purpose. Because of the invariance of the correlators of the topological recursion under conjugation of L(x) by a constant loop in SL N , we would like to work with a symplectic manifold whose points represent coadjoint orbits under the SL N action. By fixing the value of L 0,r0 , we have already fixed part of the SL N action but not all of it. Note also that fixing L 0,r0 could have been achieved by symplectic reduction by considering L 0,r0 as a moment map.

Indeed, according to section 3.2, the N -1 flows corresponding to Hamiltonians H k,0,r0(k-1)-1 k∈ 2,N lead to auxiliary matrices A k,0,r0(k-1)-1 k∈ 2,N independent of x and commuting with L 0,r0 . The corresponding Hamiltonian flows generate conjugations by elements of the stabilizer Stab L0,r 0 of L 0,r0 . Once again, extra care has to be taken for the cases r 0 ∈ {-1, 0} but the same conclusion holds in these cases.

Working modulo the coadjoint action of SL N , one may fix the values of the corresponding Hamiltonians by symplectic reduction.

This reduction leads to a reduced symplectic space ĝred := ĝ * Stab L0,r 0 of dimension

D red = D -2(N -1) = (rN -2)(N -1) def = 2g
where g is equal to the genus of the spectral curve det(y -L(x)) = 0 ⊂ C 2 , i.e. to the dimension of its Jacobian, giving rise to the usual Hitchin Hamiltonian system. This reduced space is obtained by fixing a level set for the Hamiltonians H k,0,r0(k-1)-1 k∈ 2,N and modding out by the action of the stabilizer Stab L0,r 0 of L 0,r0 . The symplectic reduction implies that the remaining auxiliary matrices should have a contribution coming from the projection to these orbits leading to a shift

A k,l → A k,l + N i=2 β i A i,0,r0(i-1)-1 for some constants (β 2 , . . . , β N ) ∈ C N -1 .

Spectral Darboux coordinates in the sl 2 case

When restricting to such an orbit, one is left with a Hitchin integrable system of dimension 2g which is a fibration above a base parametrized by the spectral invariants H k,ν,l with fiber equal to the Jacobian of the spectral curve considered (i.e. a point in the fiber is a line bundle over the spectral curve, giving the Abelian property of the fibers).

In the present article, restricting to the sl 2 (C) case, we shall not need this geometric picture since one can make all the computations explicit. However, the knowledge of this integrable system still 10 This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.

PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0002260

provides a guiding principle for our work and is very useful for the generalization of this procedure to higher rank algebras as well as base curves of different topologies.

It is thus useful to remind how one can equip a generic symplectic leaf with spectral Darboux coordinates [START_REF] Adams | Isospectral Hamiltonian Flows in Finite and Infinite Dimensions[END_REF][START_REF] Adams | Spectral Darboux Coordinates and Liouville-Arnold Integration in Loop Algebras[END_REF] allowing to identify our examples with well-known Lax pairs found in the literature. In the following, we shall often use a representation of sl 2 by 2 × 2 matrices with a basis given by σ 3 := 1 0 0 -1 = e 1e 2 , σ + := 0 1 0 0 and σ -:= 0 0 1 0 .

These Darboux coordinates are easy to describe in the sl 2 (C) case and we shall restrict to this case in the present section. Considering L 0,r0 such that the component (L 0,r0 ) 2,1 is vanishing, [L(x)] 2,1 has g zeroes located at x = q i , i ∈ 1, g . At these values, L(x) takes the form

L(q i ) = p i [L(q i )] 1,2 0 -p i
where p i = [L(q i )] 1,1 is an eigenvalue of L(q i ). The set (p i , q i ) g i=1 form a system of Darboux coordinates on the symplectic leaf going through L(x).

Since p i is an eigenvalue of L(q i ), (q i , p i ) satisfy det(p i -L(q i )) = 0 and such a pair defines a point on the spectral curve. Fixing a point in the Hitchin base, i.e. fixing the spectral curve, the set (p i , q i ) g i=1 defines g points on the spectral curve, i.e. a point in its Jacobian which is the fiber of the integrable system above the corresponding point in the base.

Note that the element L(x) can be recovered by Lagrange interpolation. Since (p i , q i ) g i=1 provide Darboux coordinates, they satisfy the evolution equations

∀ i ∈ 1, g , ∂qi ∂t k,ν,l = ∂H k,ν,l (p,q) ∂pi ∂pi ∂t k,ν,l = - ∂H k,ν,l (p,q) ∂qi (3.3)
for any triple (k, ν, l) such that H k,ν,l is not a Casimir.

Representative of a reduced orbit in the sl 2 case

In the following, it will sometimes be convenient to represent a reduced orbit by one of its representatives. Since all the results presented in this article are invariant under the action of Stab L0,r 0 , we may choose any representative we need. For the derivation of the topological type property, one shall make the following choice.

1. In the case where L 0,r0 = σ 3 , we fix the value of an additional element L ν,k to be of the form

L ν,k = U 1 V -U .
2. In the case where (L 0,r0 ) = σ + , we fix the value of an additional element L ν,k to be of the form

L ν,k = 0 V 1 0 .
Remark 3.5. The previous choices of representatives are possible only if the second term L ν,k has a non-zero off-diagonal value in the σ 3 case (in which case after a potential transposition leaving σ 3 invariant, we may always set it to the upper-right entry) or a non-zero lower-left entry in the σ + case. However, in these degenerate cases, the degree of the Poisson manifold decreases and thus the generic theory presented in this paper fails to apply directly. Thus, in the rest of the paper, we will only consider cases where the orbits admit a representative of the form presented above.

Remark 3.6. In practice, when r 0 ≥ 1, one often chooses to fix the shape of L 0,r0-1 in this way. In this article, this choice will be particularly useful in the study of the potential pole singularity at

x = ∞ presented in appendix A.
11 This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.

PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0002260 4 Building an -deformed isomonodromic system from an isospectral one

In the preceding section, we reminded how to build a symplectic leaf going through a given point L(x) whose dimension equals the genus of the associated spectral curve such that the corresponding Hamiltonian flows are isospectral. In this section, we show that one can associate to any such isospectral system a (non-unique) isomonodromic system by replacing the autonomous Hamiltonians by non-autonomous ones. We also prove that one can introduce a formal parameter rescaling Hamilton's equations around the point = 0 in the base of the integrable system where the fiber is completely degenerate.

Isomonodromic systems from non-autonomous Hamiltonians

Let us remind how to build an isomonodromic system by de-autonomization of an isospectral one. If a time t preserves the spectrum of a matrix L, the compatibility of the system

L(t)ψ(t) = µψ ∂ ∂t ψ(t) = A(t)ψ(t) takes the Lax form ∂L ∂t = [A(t), L] (4.1)
which is the form of the Hamiltonian flows studied up to now. We are actually interested in the compatibility of a system of the form

∂ ∂x Ψ(x, t) = L(x, t)Ψ(x, t) ∂ ∂t Ψ(x, t) = A(x, t)Ψ(x, t) which reads ∂L ∂t - ∂A ∂x = [A, L] . (4.2)
One can trade the former compatibility condition eq. (4.1) to the latter eq. (4.2) by replacing an autonomous system by a non-autonomous one. Namely, let us assume that one makes the matrix L(x, t, a) depend on a parameter a ∈ I such that

∂L ∂a a=t = ∂A ∂x .
Then the former compatibility condition includes both the Hamiltonian flow with respect to t and the variation in t coming from the a = t condition and gives the isomonodromic compatibility condition thanks to the relation between A and L. Let us now prove that one can build such an isomonodromic system in the N = 2 case.

Isomonodromic systems with linear auxiliary matrices

To derive an isomonodromic system admitting the topological type property for N = 2, it is enough to have one isomonodromic time leading to an auxiliary matrix of the form

A(x) = M x + B p(x) where p(x) ∈ C[x] is a C-valued polynomial in x while (M, B) ∈ sl 2 (C) × sl 2 (C) do not depend on x.
In this section, we prove that it is always possible to derive such an isomonodromic deformation with σ = 0 in eq. (3.1). We will see in some of the examples of section 6 that it is sometimes more convenient to consider other values of σ but for the purpose of the article we only need to prove that there exists at least one suitable isomonodromic deformation.

For this purpose, one can consider different possible deformations of an element L ∈ ĝred depending on its pole structure.

12 This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.

PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0002260

1. If there exists ν = 0 such that r ν = 1, one can see that

A 2,ν,1 (x) = 2 L ν,1 (x -a ν ) .
Hence, 1 2

∂A 2,ν,1 (x) ∂x = - L ν,1 (x -a ν ) 2 = ∂L(x) ∂a ν .
One thus gets an isomonodromic system with identifying t

= 1 2 t 2,ν,1 = a ν leading to    ∂ ∂x Ψ(x, t) = L(x, a ν = t)Ψ(x, t) ∂ ∂t Ψ(x, t) = Lν,1 (x-t) + N k=2 β k A k,0,r0(k-1)-1 Ψ(x, t)
where one has taken into account the reduction for defining the auxiliary matrix.

2. If there exists ν = 0 such that r ν ≥ 2, one can define

A 2,ν,rν 2 = L rν (x -a ν )
in such a way that 1 2

∂A 2,ν,rν ∂x = - L rν (x -a ν ) 2 .
This allows to define an isomonodromic system through the identification t 3. Finally, if r 0 ≥ 2 for ν = 0, one has

= 1 2 t 2,ν,rν      ∂ ∂x Ψ(x, t) = L(x) -t Lr ν (x-aν ) 2 Ψ(x, t) ∂ ∂t Ψ(x, t) = Lr ν (x-aν ) + N k=2 β k A k,0,r0(k-1)-1 Ψ(x, t) .
A 2,0,r0-2 (x) 2 = L 0,r0-1 + L 0,r0 x leading to 1 2 ∂A 2,0,r0-2 (x) ∂x = L 0,r0 .
One can thus define an isomonodromic system by identifying t

= 1 2 t 2,0,r0-2    ∂ ∂x Ψ(x, t) = [L(x) + tL 0,r0 ] Ψ(x, t) ∂ ∂t Ψ(x, t) = L 0,r0-1 + L 0,r0 x + N k=1 β k A k,0,r0(k-1)-1 Ψ(x, t) (4.3)
With these three possibilities, one exhausts all possible cases for N = 2 apart from the cases n = 0 and r 0 ∈ {0, 1}. These cases are irrelevant from the integrable systems perspective developed here since the corresponding dimension of the reduced phase space D red = (r 0 N -2)(N -1) is non positive. Moreover, the case (n, r 0 ) = (0, 0) is equivalent to L(x) = L 0,0 independent of x for which there exist explicit solutions Ψ(x) and no need for this theory. Cases for (n, r 0 ) = (0, 1), i.e. L(x) = L 0,0 + L 0,1 x correspond to either the Airy (if rank(L 0,1 ) = 1) or the Hermite-Weber (if rank(L 0,1 ) = 2) systems that are well-known to satisfy the TT property (as a consequence of [START_REF] Bouchard | Reconstructing WKB from topological recursion[END_REF] since the corresponding Newton polytopes of the classical spectral curves have no interior points, or as specific limits of the Painlevé I or II Lax systems as explained in [START_REF] Iwaki | Painlevé 2 equation with arbitrary monodromy parameter, topological recursion and determinantal formulas[END_REF]). Note that for N ≥ 3, cases with n = 0 and r 0 = 1 would provide a positive dimension of the reduced phase space and thus would require a similar analysis as the one developed above. [START_REF] Bouchard | Reconstructing WKB from topological recursion[END_REF] This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.

PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0002260

Theorem 4.1. The integrable system defined on the coadjoint orbit through any sl 2 valued rational function L(x) can be deformed into an isomonodromic system

∂ ∂x Ψ(x, t) = L(x, t)Ψ(x, t) ∂ ∂t Ψ(x, t) = A(x, t)Ψ(x, t)
where A(x) = Mx+B p(x) with p ∈ C[X] and (M, B) ∈ (sl 2 (C)) 2 and L(x, t = 0) = L(x). Remark 4.2. It is worth noticing that the different cases considered are non-exclusive. If an element L(x) falls into more than one of the cases presented above, all the corresponding matrices A are automatically compatible and one may choose to pick any of them for its purposes. In our case, we will only need to take a single matrix to prove the topological type property, and all other compatible matrices will not be considered in our procedure.

Introduction of .

One shall now introduce a formal parameter in order to replace Hamilton's equation for the autonomous system by ∂L(x,

) ∂t = [A(x, ), L(x, )]
and thus the compatibility condition for the non-autonomous one by

∂L(x, t, ) ∂t - ∂A(x, t, ) ∂x = [A(x, t, ), L(x, t, )] .
One can easily obtain such a deformation by finding a rescaling of the variable x, the Hamiltonians H k,ν,i and the times t k,ν,i by some dq factor, with

d q ∈ Q, such that t • H = O( 2 ), L 0,r0 is kept independent of

and each component of L(x) as well as Tr [L(x)]

2 are homogenous in . For this purpose, one rescales x → dx x, H k,ν,i → d k,ν,i H k,ν,i and t k,ν,i → dt k,ν,i t k,ν,i . Theorem 4.2. For any L ∈ sl 2 giving rise to an isomonodromic system as in section 4.2, there exists a set of exponents (d x ,d k,ν,i ,d t k,ν,i ) such that the compatibility equations are transformed into

∂L(x, t, ) ∂t - ∂A(x, t, ) ∂x = [A(x, t, ), L(x, t, )]
by a rescaling of t → dt t and

(x, aν , H 2,0,i , H 2,ν,i ) →      1 r 0 +1 x, 1 r 0 +1 aν , 2r 0 -i r 0 +1 H 2,0,i , 2r 0 +i r 0 +1 H 2,ν,i if rank(L 0,r 0 ) = 2 2 2r 0 +1 x, 2 2r 0 +1 aν , 2(2r 0 -i-1) 2r 0 +1 H 2,0,i , 2(2r 0 +i-1) 2r 0 +1 H 2,ν,i if rank(L 0,r 0 ) = 1 .
Proof. The proof follows from the explicit construction of such deformations in the different cases introduced for building the isomonodromic system from the isospectral one. These are not the only possible ways of introducing but it shows that it is always possible to find at least one. One shall only give the explicit computation in the case r 0 ≥ 0. The case r 0 = -1 can be treated in a similar way.

Case 1.

Let us first consider the case ν = 0 and r ν = 1. In this case, according to section 4.2, one considers t = a ν as isomonodromic time. It has the same degree as x. The associated Hamiltonian is H 2,ν,1 . We first compute d 2,ν,1 . Its value depends on the leading order L 0,r0 . Indeed, the degree in x of Tr L(x) 2 depends on the rank of L 0,r0 . For this purpose, one shall consider the cases where L 0,r0 as rank 1 or 2 separately. The degree of Tr L(x) 2 is equal to d x (2r 0 -1) in the first case and 2r 0 d x in the second case.

The homogeneity constraint for Tr L(x) 2 imposes that

d 2,ν,1 = d x (2r 0 + 1) if rank(L 0,r0 ) = 2 d x (2r 0 ) if rank(L 0,r0 ) = 1 .
14 This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.

PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0002260

The fact that the degrees of t and H sum up to 2 is thus equivalent to 2 = d x + d 2,ν,1 which leads to

d x = 1 r0+1 if rank(L 0,r0 ) = 2 2 2r0+1 if rank(L 0,r0 ) = 1 .
One finally gets the associated deformed system through the rescaling

(t, x, aν , H 2,0,i , H 2,ν,i ) →      1 r 0 +1 t, 1 r 0 +1 x, 1 r 0 +1 aν , 2r 0 -i r 0 +1 H 2,0,i , 2r 0 +i r 0 +1 H 2,ν,i if rank(L 0,r 0 ) = 2 2 2r 0 +1 t, 2 2r 0 +1 x, 2 2r 0 +1 aν , 2(2r 0 -i-1) 2r 0 +1 H 2,0,i , 2(2r 0 +i-1) 2r 0 +1 H 2,ν,i if rank(L 0,r 0 ) = 1
Note that all exponents are rational numbers.

Case 2:

Let us now consider the case when r ν ≥ 2 for some ν = 0. The general procedure defines t = t 2,ν,rν as isomonodromic time. From the homogeneity of L(x) and with similar computations as above, one finds that

d t2,ν,r ν = (2 -r ν )d x .
Then the homogeneity of Tr L(x) 2 gives

d 2,ν,i = d x (2r 0 + i) if rank(L 0,r0 ) = 2 d x (2r 0 -1 + i) if rank(L 0,r0 ) = 1 . d t2,ν,r ν + d 2,ν,rν = 2 then implies d x = 1 r0+1 if rank(L 0,r0 ) = 2 2 2r0+1 if rank(L 0,r0 ) = 1
.

One gets the associated systems by rescaling

(t, x, aν , H 2,0,i , H 2,ν,i ) →      2-rν r 0 +1 t, 1 r 0 +1 x, 1 r 0 +1 aν , 2r 0 -i r 0 +1 H 2,0,i , 2r 0 +i r 0 +1 H 2,ν,i if rank(L 0,r 0 ) = 2 2(2-rν ) 2r 0 +1 t, 2 2r 0 +1 x, 2 2r 0 +1 aν , 2(2r 0 -i-1) 2r 0 +1 H 2,0,i , 2(2r 0 +i-1) 2r 0 +1 H 2,ν,i if rank(L 0,r 0 ) = 1
Note again that all exponents are rational numbers.

Case 3. Let us finally consider the case

L(x) = r0 i=0 L 0,i x i + n ν=1 rν i=1 L ν,i (x -a ν ) i with r 0 ≥ 2.
The isomonodromic time t = 1 2 t 2,0,r0-2 . The homogeneity condition for Tr L(x) 2 implies that

d 2,0,i = d x (2r 0 -i) if rank(L 0,r0 ) = 2 d x (2r 0 -1 -i) if rank(L 0,r0 ) = 1 .
On the other hand, from the homogeneity of L(x), one sees that t is of order x r0 leading to

d t2,0,r 0 -2 = r 0 d x . The requirement t.H = O 2 thus reads 2 = d t2,0,r 0 -2 + d 2,0,r0-2 = d x (2r 0 + 2) if rank(L 0,r0 ) = 2 d x (2r 0 + 1) if rank(L 0,r0 ) = 1 ,
i.e.

d x = 1 r0+1 if rank(L 0,r0 ) = 2 2 2r0+1 if rank(L 0,r0 ) = 1 .
15 This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.

PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0002260

One gets the associated systems by rescaling

(t, x, aν , H 2,0,i , H 2,ν,i ) →    r 0 r 0 +1 t, 1 r 0 +1 x, 1 r 0 +1 aν , 2r 0 -i r 0 +1 H 2,0,i , 2r 0 +i r 0 +1 H 2,ν,i if rank(L 0,r 0 ) = 2 2r 0 2r 0 +1 t, 2 2r 0 +1 x, 2 2r 0 +1 aν , 4r 0 -2i-2 2r 0 +1 H 2,0,i , 4r 0 +2i-2 2r 0 +1 H 2,ν,i if rank(L 0,r 0 ) = 1 .
Remark 4.3. If we have proved that one can always -deform such an isomonodromic system, there exist many other ways to do it. Examples will be given by some of the Painlevé equations. In any case, the procedure for determining the introduction of always follows from some homogeneity conditions as well as the constraint for the order of the product tH.

In the introduction of , one can consider x as a dynamical variable, i.e. changing the value of changes the base curve, or keep it independent of . We always considered the first possibility here but we shall see in examples that another choice can be made.

After performing this non-autonomous deformation and introducing , the dynamics of the system is described by

∀ i ∈ 1, g , ∂qi( ,t) ∂t k,ν,l = ∂H k,ν,l (p,q,t) ∂pi ∂pi( ,t) ∂t k,ν,l = - ∂H k,ν,l (p,q,t) ∂qi (4.4)
It differs from equations 3.3 by the factor and the fact that the Hamiltonians H k,ν,l (p, q, t) are non-autonomous now, i.e. they may explicitly depend on the times t.

It is also interesting to remark that this rescaling implies a rescaling of the quadratic differential

Tr(L(x) 2 ) (dx) 2 → 2 Tr(L(x) 2 ) (dx)
2 used to define the base of the corresponding Hitchin integrable system.

Let us now study how the differential equation in x gets deformed by introducing . This depends on the rank of L 0,r0 . We use the representative of the orbit considered in remark 3.7 of section 3.6.

-Let us first consider the case L 0,r0 = σ 3 . After rescaling, the differential equation in x reads

dx ∂ ∂x Ψ(x, t, ) = -r0dx [L(x, t)] 1,1 -[L(x, t)] 2,2 + -dx(r0+1) [L(x, t)] 1,2 + -dx(r0-1) [L(x, t)] 2,1 Ψ(x, t, ) (4.5)
where we made explicit any dependence in . Remember that in all cases, d x = 1 r0+1 . It is natural to gauge transform (this gauge transformation acts by conjugation on the Lax matrices and does not change the compatibility conditions) the system by

Ψ(x, t, ) → dx 2 0 0 -dx 2 Ψ(x, t, )
and multiply equation (4.5) by r0dx to get the -deformed equation:

∂ ∂x Ψ(x, t, ) = L(x, p, q, t)Ψ(x, t, )
where the Lax matrix depends on only through the solutions (p, q) of Hamilton's equations.

Let us also notice that the rescaling can be expressed at the level of the spectral Darboux coordinates as

∀ i ∈ 1, g : p i → r0dx p i q i → dx q i .
-Let us now consider the case L 0,r0 = σ + . The rescaled differential equation reads

dx ∂ ∂x Ψ(x, t, ) = - 2r 0 -1 2 dx [L(x, t)] 1,1 -[L(x, t)] 2,2 + -r0dx [L(x, t)] 1,2
16 This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.

PLEASE CITE THIS ARTICLE AS

∂ ∂x Ψ(x, t, ) = L(x, p, q, t)Ψ(x, t, )
which is the same as above.

The spectral Darboux coordinates are rescaled according to

∀ i ∈ 1, g : p i → 2r 0 -1 2r 0 +1 p i q i → 2 2r 0 +1 q i .

Topological recursion for isomonodromic systems and quantum curve

For any g = sl 2 (C) valued rational function L(x), we have explained how to build an isomonodromic system

∂ ∂x Ψ(x, t) = L(x, t)Ψ(x, t) ∂ ∂t Ψ(x, t) = A(x, t)Ψ(x, t)
such that L(x, t = 0) = L(x) and the compatibility condition

∂L ∂t - ∂A ∂x = [A, L]
has a non-autonomous Hamiltonian formulation. Moreover, we could introduce a parameter in order to replace the above system by

∂ ∂x Ψ(x, t, ) = L(x, t, )Ψ(x, t, ) ∂ ∂t Ψ(x, t, ) = A(x, t, )Ψ(x, t, ) (5.1)
with the compatibility condition

∂L(x, t, ) ∂t - ∂A(x, t, ) ∂x = [A(x, t, ), L(x, t, )] (5.2)
where lim →0 L(x, t, ) ∈ ĝ is non-vanishing. We can now state the main result of this article.

Theorem 5.1. The -dependent isomonodromic systems built in section 4.3 satisfy the topological type property with ω 0,2 given by the unique Bergman kernel on the classical spectral curve obtained as the compactification of

Σ cl := (x, y) ∈ C 2 | lim →0 det (y -L(x, t, )) = 0 .
Proof. The proof is technical and consists in checking one by one the conditions of the topological type property. This is done in the various sections of appendix A.

As a corollary, this implies that the corresponding isomonodromic tau function can be computed by topological recursion. 17 This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.

PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0002260 Corollary 5.1. The tau function τ BM (t, ) of the -dependent isomonodromic systems of section 4.3 reads

τ BM (t, ) = α exp ∞ g=0 2g-2 F g (t)
where α is a constant independent of t and the coefficients F g (t) are computed by topological recursion. This also implies that τ BM (t, ) is a generating function of intersections of tautological classes in the Deligne-Mumford compactification of the moduli space of Riemann surfaces.

Proof. Corollary 5.1 is a direct consequence of the fact that the system satisfies the Topogical Type property (Theorem 5.1) using Theorem 2.2.

The above result can be interpreted as a generalization of Kontsevich-Witten theorem [START_REF] Kontsevich | Intersection theory on the moduli space of curves and the matrix Airy function[END_REF] by bringing a correspondence between isomonodromic tau functions and generating functions of correlators of CohFT's. It also provides a set of Virasoro constraints annihilating this tau function.

In [START_REF] Bouchard | Reconstructing WKB from topological recursion[END_REF], it was proved that one can construct a differential equation

d dx Ψ(x, ) = L(x, )Ψ(x, )
with the topological type property whenever the Newton polytope of the classical spectral curve P (x, y) = lim →0 det(y -L(x, )) = 0 does not have any interior point. In particular, this procedure consists in associating a linear differential operator d dx -L(x, ) to a classical spectral curve P (x, y) = 0. In the literature, this linear differential operator is often referred to as a "quantum curve" associated to a classical spectral curve P (x, y) = 0. Up to now, the only examples of quantum curves obtained from a classical spectral curve of genus 0 with interior points in its Newton polytope were given by the Painlevé equations [START_REF] Iwaki | Painlevé equations, topological type property and reconstruction by the topological recursion[END_REF].

However, theorem 5.1 can be used to provide a map that associates a quantum curve to any genus 0 hyperelliptic classical spectral curve. Indeed, let us consider a genus 0 hyperelliptic classical spectral curve presented in the form

y 2 = Q(x)
where Q(x) is a rational function. From the pole structure of Q(x) we may present this equation as the characteristic polynomial of a sl 2 valued function L(x) with fixed poles. One can then apply the result of the present work to define a -deformation of L(x) so that one gets a differential operator d dx -L(x, ) satisfying the topological type property, hence producing a quantum curve whose symbol in the → 0 limit is precisely y 2

-Q(x) = 0.
It is worth noticing that, in this setup, the spectral curve of L(x, ) has non-vanishing genus for = 0 and the classical spectral curve corresponds to a fully degenerate curve in this family.

Examples

Painlevé examples

Let us apply the formalism to the case of Painlevé equations, finding the Hamiltonians and the auxiliary systems. All these cases correspond to sl 2 (C) cases with r = 2 giving rise to D = 4 and a reduced space of dimension 2.

Painlevé I

Let us discuss in details the example leading to the Painlevé I equation, showing how our procedure produces the deformation considered in [START_REF] Iwaki | Painlevé equations, topological type property and reconstruction by the topological recursion[END_REF]. For this purpose, let us consider

ĝ := L(x) := 2 i=0 L 0,i x i , (L 0,i ) 0≤i≤2 ∈ sl 2 (C) 3 /SL 2 18
This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset. fixed. This is an example where , n = 0, r 0 = 2 and L 0,r0 does not have rank 2. Thus, the characteristic polynomial of L(x) has degree 3 in x.

PLEASE CITE THIS ARTICLE AS

The non-Casimir Hamiltonians are given by

H 2,0,0 = Res x→∞ x -1 Tr [L(x)] 2 dx and H 2,0,1 = Res x→∞ x -2 Tr [L(x)] 2 dx
with the associated auxiliary matrices

A 2,0,0 2 = x -1 L(x) + = 0 1 0 0 x + L 0,1 and A 2,0,1 2 = x -2 L(x) + = 0 1 0 0 .
Let us now explicitly express a point in a given symplectic leaf. Denoting

∀ i ∈ 1, 2 , L 0,i := u i v i w i -u i ,
the coefficients of the characteristic polynomial read

       H 2,0,0 = 2(u 2 0 + v 0 w 0 ) H 2,0,1 = 2(2u 1 u 0 + v 1 w 0 + w 1 v 0 ) H 2,0,2 = 2(w 0 + u 2 1 + v 1 w 1 ) H 2,0,3 = 2w 1 .
One shall now restrict to a symplectic leaf by setting H 2,0,3 = 2α 3 = 0 and H 2,0,2 = 2α 2 and map it to a point in the reduced space by fixing the value of H 2,0,1 = 2α 1 and u 1 = 0. A point in this reduced space hence reads

L(x) = 0 1 0 0 x 2 + 0 v 1 α 3 0 x + u 0 v 2 1 + α1 α3 -v 1 α2 α3 α 2 -v 1 α 3 -u 0 .
We can now identify the isomonodromic time t with t 2,0,0 leading to

L(x) = 0 1 0 0 x 2 + 0 v 1 α 3 0 x + u 0 v 2 1 + α1 α3 -v 1 α2 α3 + 2t α 2 -v 1 α 3 -u 0 . One indeed has ∂ L(x, t) ∂t 2,0,0 = ∂ Ã2,0,0 ∂x where Ã2,0,0 2 = x -1 L(x, t) + = 0 1 0 0 x + β
where β is a matrix independent of x.

The spectral Darboux coordinates (p, q) are given by

q = v 1 -α2 α3 p = u 0 leading to the representation L(x) = p x 2 + q + α2 α3 x + q q + α2 α3 + α1 α3 + 2t α 3 (x -q) -p .
One recovers the Lax representation of Painlevé I used in [START_REF] Iwaki | Painlevé equations, topological type property and reconstruction by the topological recursion[END_REF] by setting α 1 = α 2 = 0 and α 3 = 1. The Hamiltonian of the system is given by H 2,0,0 .

19 This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.

PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0002260

Hence, our isomonodromic system takes the form

       ∂ ∂x Ψ = 0 1 0 0 x 2 + L 0,1 x + L 0,0 + 2t 2,0,0 0 1 0 0 Ψ ∂ ∂t2,0,0 Ψ = 0 1 0 0 x + β Ψ .
We shall now detail the introduction of the parameter . From the computation of

Tr( L(x, t) 2 ) = x 3 H 2,0,3 + x 2 H 2,0,2 + xH 2,0,1 + H 2,0,0
one can read the homogeneity degree of the Hamiltonians with respect to x. If x has degree d x and imposing H 2,0,3 of degree 0 (this assumption can be released easily but leads to more complicated expressions), one gets that H 2,0,i has degree

d 2,0,i = (3 -i)d x .
From equation (3.2), one can compute the value of the Casimirs in terms of the entries of the matrix leading to d t2,0,0 = 2d x . The condition d t2,0,0 + d 2,0,0 = 2 imposes d x = 2 5 as expected for the Painlevé I case from [START_REF] Iwaki | Painlevé equations, topological type property and reconstruction by the topological recursion[END_REF].

Painlevé II

Let us consider again n = 0 and r 0 = 2 but now fix L 0,2 = σ 3 . This implies that the characteristic polynomial of L(x) has degree 4 in opposition to the Painlevé I case above.

As in the preceding section, the non-Casimir Hamiltonian are given again by

H 2,0,0 = Res x→∞ x -1 Tr [L(x)] 2 dx and H 2,0,1 = Res x→∞ x -2 Tr [L(x)] 2 dx
with the associated auxiliary matrices

A 2,0,0 2 = x -1 L(x) + = σ 3 x + L 0,1 and A 2,0,1 2 = x -2 L(x) + = σ 3 .
We can compute the spectral Darboux coordinates for this system. For simplicity let us consider a symplectic leaf of the form (H 2,0,3 , H 2,0,2 , H 2,0,1 ) = (0, α 2 , α 1 ). Considering a representative of the reduced orbit as before, one has a Lax matrix of the form

L(x) = σ 3 x 2 + 0 v 1 1 0 x + u 0 v 0 w 0 -u 0 with 2u 0 + v 1 = α 2 and v 0 + v 1 w 0 = α 1 . One obtains the spectral Darboux coordinates q = -w 0 p = q 2 + u 0 giving L(x) = σ 3 x 2 + 0 α 2 -2(p -q 2 ) 1 0 x + (p -q 2 ) α 1 + q α 2 -2(p -q 2 ) -q -(p -q 2
) .

Remark that p = p-q 2 gives an alternative Darboux coordinate dual to q. In terms of the coordinates (p, q), one gets

L(x) = σ 3 x 2 + 0 α 2 -2p 1 0 x + p α 1 + q [α 2 -2p] -q -p .
Following the general procedure, one can identify the isomonodromic time with t 2,0,0 and consider L(x, t) := σ 3 x 2 + L 0,1 x + L 0,0 + 2tσ 3 20 This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.

PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0002260

leading to the isomonodromic system

∂ ∂x Ψ = σ 3 x 2 + L 0,1 x + L 0,0 + 2t 2,0,0 σ 3 Ψ ∂ ∂t2,0,0 Ψ = [σ 3 x + β] Ψ .
The explicit expression of a point in the reduced space can be performed as before and allows to introduce the parameter in the same way. One finds the homogeneity degrees d 2,0,0 = 4d x and d t2,0,0 = 2d x .

This implies that the rescaling should follow from

d x = 2 4 + 2 = 1 3 .
A Lax representation of Painlevé II is recovered by setting H 2,0,1 = 0, H 2,0,2 = t, H 2,0,3 = θ. The Hamiltonian of the system is given by H 2,0,0 .

Painlevé III

In order to recover Painlevé III equation, let us consider the case where n = 1, r 0 = 0, a 1 = 0, r 1 = 2 and fix L 0 = σ 3 to get an element of the form

L(x) = σ 3 + L 1 x -1 + L 2 x -2 .
The non-Casimir Hamiltonians are

H 2,1,1 = -Res x→0 Tr [L(x)] 2 dx and H 2,1,2 = -Res x→0 x Tr [L(x)] 2 dx
with the associated auxiliary matrices

A 2,1,1 2 = [L(x)] + = -L 0 and A 2,1,2 2 = [xL(x)] + = -L 0 x -L 1 .
The Hamiltonian of the system is H 2,1,2 while H 2,1,1 should be fixed by symplectic reduction.

In order to recover a known representation of Painlevé III, one should consider some other auxiliary matrices by shifting A 2,1,2 by a matrix proportional to L(x). We shall then recover the Lax representation considered in [29]. For this purpose, one defines the auxiliary matrix

A(x) := A 2,1,2 + 2xL(x) = xL 0 + L 1 -L 2 x -1 obtained from (3.1) with σ = 2. Hence ∂A(x) ∂x = 1 2 L 0 + L 2 x -2 .
This motivates the introduction of a time t through a rescaling L 0 → tL 0 and L 2 → tL 2 leading to the isomonodromic system

∂ ∂x Ψ = tσ 3 + L 1 x -1 + tL 2 x -2 Ψ ∂ ∂t Ψ = xtσ 3 + L 1 -tL 2 x -1 Ψ .
Let us now introduce . We cannot apply the general result of the preceding section since we have considered σ = 0. However, we can take d x = 0 which immediately imposes that t should be of degree 1 and thus H as well and we recover the deformation of [START_REF] Iwaki | Painlevé equations, topological type property and reconstruction by the topological recursion[END_REF].

21

This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset. PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0002260

Painlevé IV

Let us consider the case where n = 1, r 0 = 1, r 0 = 1 and L 0,1 = σ 3 ,

L(x) = L 0,1 x + L 0,0 + L 1,1 x -a 1
The non-Casimir invariant functions are

H 2,0,0 = Res x→∞ x -1 Tr [L(x)] 2 dx and H 2,1,1 = -Res x→a1 Tr [L(x)] 2 dx
with the associated auxiliary matrices

A 2,0,0 2 = x -1 L(x) + = L 0,1 and A 2,1,1 2 = -[L(x)] -= -L 1,1 (x -a 1 ) -1 .
The Hamiltonian is H 2,1,1 while H 2,0,0 should be fixed by symplectic reduction. Following the general procedure for n ≥ 1, the isomonodromic time t 2,1,1 should be identified with the position a 1 of the second pole. This leads to the isomonodromic system

∂ ∂x Ψ = L 0,1 x + L 0,0 + L1,1 x-t2,1,1 Ψ ∂ ∂t2,1,1 Ψ = - L1,1 x-t2,1,1 Ψ .
The introduction of follows from the general procedure with d 2,1,1 = 3 leading to

d x = 1 2 .

Painlevé V

Let us consider the case where n = 2, r 0 = 0, r 1 = 1, r 2 = 1, a 1 = 0 and a 2 = 1,

L(x) = L 0,0 + L 1,1 x -1 + L 2,1 x -1
with L 0,0 = σ 3 . We consider the non-Casimir functions

H = -Res x→∞ x Tr [L(x)] 2 dx and a = Res x→∞ Tr [L(x)] 2 dx
with the associated auxiliary matrices

A 2 = -[xL(x)] + = -xL 0,0 -L 1,1 and A ′ 2 = [L(x)] + = L 0,0 .
In order to recover Painlevé V equation, following [START_REF] Harnad | R-matrix construction of electromagnetic models for the Painlevé transcendents[END_REF][START_REF] Its | On Some Hamiltonian Properties of the Isomonodromic Tau Functions[END_REF], one shall consider the Hamiltonian

H * := H - a 2 4 .
One recovers a Lax representation of Painlevé V by considering σ = 1 in (3.1) leading to the auxiliary matrix

A 1 = [xL(x) + aL(x)] + .
As in Painlevé III case, one obtains an isomonodromic system by rescaling L 0,0 → tL 0,0 and H * → H * t . One shall now introduce . For this purpose, one can see that H * has the same degree as t 2 if x is independent of leading to the rescaling described in [START_REF] Iwaki | Painlevé equations, topological type property and reconstruction by the topological recursion[END_REF].

22

This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset. PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0002260

Painlevé VI

Let us consider the case where n = 3, r 0 = -1, a 1 = 0, a 2 = 1 and r 1 = r 2 = r 3 = 1, i.e. the simplest Fuchsian system.

As in any Fuchsian system, one finds an isomonodromic system by identifying t = a 3 (or any of the simple poles) and considering the Hamiltonian H 2,3,1 .

From the general procedure to introduce , one has that x and t have the same homogeneity degree. We present below the case of general Fuchsian systems that includes Painlevé VI.

Fuchsian systems

For N = 2, one can consider a Fuchsian system with an arbitrary number of poles n ≥ 3 such that

L(x) = n i=1 L i,1 x -a i .
We shall consider here the case where one fixes L 0,-1 = -n i=1 L i,1 = σ 3 , other cases being easily obtained in a similar way. One can assume that a 1 = 0 and a 2 = 1. The remaining (a j ) 3≤j≤n can be identified with isomonodromic times: t i = a i+2 , associated with the Hamiltonians

H i := H 2,i+2,1 with i ∈ 1, n -2 . Let us denote L i,1 = u i σ 3 + v i σ + + w i σ -.
One now introduces by rescaling

u i → u i and v i → 2 v i
while leaving all other quantities invariant. This implies

H i → 2 H i and t i → t i .
Since we are in a case where r 0 = -1 which was not explicitly detailed in the introduction of , let us use this opportunity to study it more carefully. The full isomonodromic system reads

     d dx Ψ = n i=1 -1 uiσ3+ -2 viσ++wiσ- x-ai Ψ d dai+2 Ψ = -1 uiσ3+ -2 viσ++wiσ- x-ai + -1 c i σ 3 Ψ
where c i is a constant independent of . After multiplying by and performing the gauge transformation

Ψ → 1 2 0 0 -1 2 Ψ, one gets the system      d dx Ψ = n i=1 uiσ3+viσ++wiσ- x-ai Ψ d dai+2 Ψ = uiσ3+viσ++wiσ- x-ai + c i σ 3 Ψ .

Second elements of the Painlevé hierarchies

Let us consider some simple higher dimensional examples by remaining in the g = sl 2 (C) case. We consider examples leading to the second elements of some of the Painlevé hierarchies.

23

This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset. PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0002260

Painlevé II 2 hierarchy

Let us first consider n = 0, r 0 = 3 and L 0,3 = σ 3 ,

L(x) := 3 i=0 L 0,i x i .
One considers the non-Casimir invariant functions

H 2,0,0 = Res x→∞ x -1 Tr L(x) 2 dx , H 2,0,1 = Res x→∞ x -2 Tr L(x) 2 dx and H 2,0,2 = Res x→∞ x -3 Tr L(x) 2 dx
together with the auxiliary matrices

A 2,0,0 2 = x -1 L(x) + = σ 3 x 2 + L 0,2 x + L 0,1 , A 2,0,1 2 = x -2 L(x) + = σ 3 x + L 0,2 , A 2,0,2 2 = x -3 L(x) + = σ 3 .
One shall proceed by symplectic reduction by fixing the value of H 2,0,2 . One can make it an isomonodromic system by considering.

L(x) := L(x) + 2t 2,0,1 σ 3 + 2t 2,0,0 [2σ 3 x + L 0,2 ] .
For completeness, we introduced here the two isomonodromic times corresponding to this four dimensional symplectic space even though we only need one to apply our procedure. The isomonodromic system is a Lax representation of (P II 2 ) 2 [START_REF] Chiba | Multi-poisson approach to the Painlevé equations: From the isospectral deformation to the isomonodromic deformation[END_REF][START_REF] Chiba | The first, second and fourth Painlevé equations on weighted projective spaces[END_REF]. To introduce , one sees that H 2,0,0 (resp. H 2,0,1 ) has homogenous degree d 2,0,0 = 6d x (resp. d 2,0,1 = 5d x ) with respect to x while d t2,0,0 = 2d x and d t2,0,1 = 3d x . Hence, one shall consider

d x = 2 2 + 6 = 2 3 + 5 = 1 4
for the rescaling.

Painlevé IV hierarchy

Let us consider n = 1, r 0 = 2 and r 1 = 1 together with L 0,2 = σ 3 ,

L(x) := 2 i=0 L 0,i x i + L 1,1 x -1 .
One considers the non-Casimir invariant functions

H 2,0,0 = Res x→∞ x -1 Tr L(x) 2 dx , H 2,0,1 = Res x→∞ x -2 Tr L(x) 2 dx and H 2,1,1 = Res x→∞ Tr L(x) 2 dx
together with the auxiliary matrices

A 2,0,0 2 = x -1 L(x) + = σ 3 x + L 0,1 , A 2,0,1 2 = x -2 L(x) + = σ 3 , A 2,1,1 2 = [L(x)] + = σ 3 x 2 + L 0,1 x + L 0,0 .
One shall fix the value of H 2,0,1 for obtaining the reduced space. One can make it an isomonodromic system by considering

L(x) := L(x) + 2t 2,0,0 σ 3 + 2t 2,1,1 [2σ 3 x + L 0,1 ] + 2t 2 2,1,1 σ 3
The isomonodromic systems is a Lax representation of (P IV) 2 . The main isomonodromic time is t 2,0,0 in this case. [START_REF] Eynard | Invariants of algebraic curves and topological expansion[END_REF] This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.

PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0002260

Remark 6.1. The appearance of the t 2 2,1,1 term comes from the fact that dA2,1,1/dx does depend on t2,1,1 linearly. This order 2 term appears for correcting this. In general, one could have higher order terms used for correcting such issues in higher degree examples (such as higher order elements in this hierarchy).

Once again, one only needs the introduction of the isomonodromic deformation t2,0,0 for our purpose but we introduce the other one for making the connection with the integrable hierarchy.

To introduce , one sees that H 2,0,0 (resp. H 2,1,1 ) has homogenous degree d 2,0,0 = 4d x (resp. d 2,1,1 = 5d x ) with respect to x while d t2,0,0 = 2d x and d t2,1,1 = d x . Hence, one shall consider

d x = 2 2 + 4 = 2 1 + 5 = 1 3
for the rescaling.

Painlevé I hierarchy

Let us consider n = 0, r 0 = 3 and L 0,3 = σ + ,

L(x) := 3 i=0 L 0,i x i .
One considers the non-Casimir invariant functions

H 2,0,0 = Res x→∞ x -1 Tr L(x) 2 dx , H 2,0,1 = Res x→∞ x -2 Tr L(x) 2 dx and H 2,0,2 = Res x→∞ x -3 Tr L(x) 2 dx
together with the auxiliary matrices

A 2,0,0 2 = x -1 L(x) + = σ + x 2 + L 0,2 x + L 0,1 , A 2,0,1 2 = x -2 L(x) + = σ + x + L 0,2 , A 2,0,2 2 = x -3 L(x) + = σ + .
One shall fix the value of H 2,0,2 for obtaining the reduced space. One can make it an isomonodromic system by considering

L(x) := L(x) + 2t 2,0,1 σ + + 2t 2,0,0 [2σ + x + L 0,2 ] .
The isomonodromic systems is a Lax representation of (P I) 2 .

To introduce , one sees that H 2,0,0 (resp. H 2,0,1 ) has homogenous degree d 2,0,0 = 5 (resp. d 2,0,1 = 4d x ) with respect to x while d t2,0,0 = 2d x and d t2,0,1 = 3d x . Hence, one shall consider

d x = 2 2 + 5 = 2 3 + 4 = 2 7
for the rescaling.

Painlevé II 1 hierarchy

Let us now take n = 2, r 0 = 2, r 1 = 1 and fix L 0,2 = σ + ,

L(x) := 2 i=0 L 0,i x i + L 1,1 x -1 .
One considers the non-Casimir invariant functions

H 2,0,0 = Res x→∞ x -1 Tr L(x) 2 dx , H 2,0,1 = Res x→∞ x -2 Tr L(x) 2 dx and H 2,1,1 = Res x→∞ Tr L(x) 2 dx
25 This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.

PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0002260 together with the auxiliary matrices

A 2,0,0 2 = x -1 L(x) + = σ + x + L 0,1 , A 2,0,1 2 = x -2 L(x) + = σ + , A 2,1,1 2 = [L(x)] + = σ + x 2 + L 0,1 x + L 0,0 .
One shall fix the value of H 2,0,1 for defining the reduced space.

One can make it an isomonodromic system by considering

L(x) := L(x) + 2t 2,0,0 σ + + 2t 2,1,1 [2σ + x + L 0,1 ] + 2t 2 2
,1,1 σ + which gives a Lax representation of (P II 1 ) 2 . The main isomonodromic time is t 2,0,0 in this case.

To introduce , one sees that H 2,0,0 (resp. H 2,1,1 ) has homogenous degree d 2,0,0 = 3d x (resp. d 2,1,1 = 4d x ) with respect to x while d t2,0,0 = 2d x and d t2,1,1 = d x . Hence, one shall consider

d x = 2 2 + 3 = 2 1 + 4 = 2 5
for the rescaling.

Polynomial of arbitrary order

Let us show how our procedure applies to arbitrary degree polynomials as the ones coming from integrables hierarchies. We shall consider the case n = 0 and arbitrary values of r 0 with L 0,r0 = σ 3 allowing to treat elements of the Painlevé II 2 hierarchy

L(x) = σ 3 x r0 + r0-1 i=0 L 0,i x i .
Other values of L 0,r0 or values of n can be treated in a similar way. The functions

H 2,0,i := Res x→∞ x -i-1 Tr L(x) 2 dx
provide Casimir functions for i ≥ r 0 and non-Casimir for 0 ≤ i ≤ r 0 -1 with associated auxiliary functions

∀ 0 ≤ i ≤ r 0 -1 , A 2,0,i 2 = x -i-1 L(x) + .
We shall fix H 2,0,r0-1 for defining the reduced space which has dimension 2r 0 -1 and is equipped with the Hamiltonians H 2,0,i for i ∈ 0, r 0 -2 .

We can now define the associated isomonodromic system by introducing the times t = (t 2,0,i ) r0-2 i=0 by defining

L(x, t) = r0 k=0 L0,k (t)x k
where L0,k (t), k = 0, . . . , r 0 , are obtained as the unique solution of the set of equations

∀ 0 ≤ k ≤ r 0 , ∀ 0 ≤ k + i + 2 ≤ r 0 , ∂ L0,k (t) ∂t 2,0,i = 2(k + 1) L0,k+i+2
with initial condition L0,k (0) = L 0,k .

One can now introduce following the general procedure leading to

d x = 1 r 0 + 1
26 This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.

PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0002260

This article explains how the topological recursion can be used to express a tau function of isomonodromic deformations of SL 2 connections over the Riemann sphere in some particular regimes. It is natural to ask whether this procedure can be generalized to any meromorphic G-connection over an arbitrary Riemann surface for a reductive group G. It seems very likely to be the case and we plan to address this problem in a future work. From the topological recursion side, the TT properties have been defined for such a general setup in [START_REF] Belliard | Integrable differential systems of topological type and reconstruction by the topological recursion[END_REF]. On the other hand, the Hitchin integrable system is well known to be described by a similar R-matrix construction. It seems also possible to deform it to an isomonodromic system by de-autonomization and introduce a parameter through the same steps as the ones followed in the present article. The question is to know whether we can do it in a way which gives rise to the TT properties. If the properties 2, 3, 5 and 6 of Definition 2.4 should be easy to obtain, the genus 0 property as well as the pole property will certainly prove to be harder to get, even though the → 0 limit of Lax equations are likely to allow to obtain property 1, we cannot rely on a simple explicit matrix representation of the corresponding Lie algebra.

On the other hand, it should be noted that recent progress have been made for classifying such isomonodromic systems by establishing isomorphisms between moduli spaces of connections and Nakajima quiver varieties [START_REF] Hiroe | Moduli spaces of meromorphic connections, quiver varieties, and integrable deformations[END_REF][START_REF] Woodhouse | The symplectic and twistor geometry of the general isomonodromic deformation problem[END_REF]. This allows to associate a quiver to an isomonodromic system of the type studied here. These quivers have symmetries that preserve the differential of the associated isomonodromic tau function, hence the free energies defined in our paper. In the simplest case, such symmetries give rise to the well-known Harnad duality [START_REF] Harnad | Dual isomonodromic deformations and moment maps to loop algebras[END_REF]. In particular, it means that even if our result applies only to rank 2 systems, they also imply that the tau function of Harnad's dual system can also be computed by topological recursion, giving access to higher rank systems. It would be very nice to understand better the action of the symmetries acting on quivers in terms of the topological recursion. It is highly plausible that it would shed some light on the mysterious symplectic invariance property of the partition functions computed by topological recursion.

Another natural generalization of our work would be to work with classical spectral curves embedded into C * × C * instead of C 2 . This is usually obtained by considering exponentiated variables as in the study of Gromov-Witten invariants of local CY three-folds. There exists a corresponding integrable system defined by Goncharov and Kenyon [START_REF] Goncharov | Dimers and cluster integrable systems[END_REF] which considers the space of connections on graphs on a torus. This replaces the usual Hitchin system by trading differential operators by finite difference operators. We believe that one should be able to apply our procedure to such a system to obtain a partition function computed by the topological recursion on this new type of classical spectral curve.

Even further is the possibility to quantize the whole picture (and all the generalizations described above) by replacing the classical R-matrix by a quantum R-matrix construction. Such a quantization of isomonodromic systems gives rise to the Knizhnik-Zamolodchikov-Bernard equations [START_REF] Harnad | Quantum isomonodromic deformations and the Knizhnik-Zamolodchikov equations[END_REF][START_REF] Reshetikhin | The Knizhnik-Zamolodchikov System as a Deformation of the Isomonodromy Problem[END_REF]. This allows us to interpret the result in terms of sections of the associated Verlinde bundle, i.e. in terms of some conformal blocks of an associated CFT. From the topological recursion perspective, we believe that it corresponds to the so-called β-deformations taking as input a non-commutative version of the classical spectral curve.

Finally, all our computations give rise to formal series in , where can be thought of as a local coordinate around a genus 0 singular point in the base of Hitchin's integrable system. It would be nice to be able to consider the expansion around any fiber of the Hitchin system. However, we do not have any hope to be able to do so without introducing some non-perturbative effects breaking the -expansion properties. There exists a conjecture in the literature [START_REF] Borot | Geometry of spectral curves and all order dispersive integrable system[END_REF] for such a non-perturbative generalization. This conjecture was mainly motivated by some random matrix model examples. Since isomonodromic systems are generalizations of random matrix models and admit some Riemann-Hilbert representations, we could hope to use this approach to give a more firm ground to the conjecture if not proving it. A first step in this direction seems to have been reached by Iwaki who could prove this conjecture in the case of Painlevé I equation by building a two parameter family of solutions as discrete Fourier transforms of tau functions built by topological recursion [START_REF] Iwaki | 2-parameter τ -function for the first Painlevé equation -topological recursion and direct monodromy problem via exact WKB analysis[END_REF].

A Topological type property and isomonodromic systems

In this appendix, we prove theorem 5.1 by checking that the conditions 1-6 of the topological type property given in definition 2.3 are fulfilled by systems considered in the present paper.

A.1 Proof of condition 2 : Formal expansion in

There are several ways to justify the existence of a formal expansion in . The first one used in [START_REF] Bergère | Determinantal formulae and loop equations[END_REF][START_REF] Bergère | Rational differential systems, loop equations, and application to the q th reductions of KP[END_REF][START_REF] Iwaki | Painlevé 2 equation with arbitrary monodromy parameter, topological recursion and determinantal formulas[END_REF][START_REF] Iwaki | Painlevé equations, topological type property and reconstruction by the topological recursion[END_REF] is to consider formal WKB solutions of the system of equations (5.1). In this context, it is equivalent to assume that

Ψ(x, t, ) = exp ∞ k=-1 Ψ (k) (x, t) k = ∞ k=0 Ψ(k) (x, t) k exp 1 Ψ (-1) (x, t) . (A.1) Since L(x, t, ) = d dx Ψ(x, t, ) Ψ -1 (x, t, ) and A(x, t, ) = d dt Ψ(x, t, ) Ψ -1 (x, t,
), one automatically gets that L(x, t, ), A(x, t, ) and M (x, t, ) admits a formal series expansion in of the form:

L(x, t, ) = ∞ k=0 L (k) (x, t) k , A(x, t, ) = ∞ k=0 A (k) (x, t) k , M (x, t, ) = ∞ k=0 M (k) (x, t) k . (A.2)
However, according to the general geometric construction presented in this article, we propose an alternative way to obtain the same result. Indeed, in this paper, we introduce the parameter using rescaling of parameters and therefore it is necessary to justify the existence of formal expansions (A.2) within this formalism. First of all, let us mention that after rescaling, the values of the Casimirs defining the symplectic leaf as well as the value of the moment map used for the reduction are independent of . In our construction, the dependance in only comes from the values of the Darboux coordinates (p i ( ), q i ( )) g i=1 . Given the non-autonomous Hamiltonian H and the associated isomonodromic time t built in section 4.2, the -deformed Hamilton's equations are given by:

∀ i ∈ 1, g , ∂q i ( ) ∂t = ∂ ∂p i H(p, q, t) and - ∂p i ( ) ∂t = ∂ ∂q i H(p, q, t). (A.3)
30 This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.

PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0002260

We get that (p i ( ), q i ( )) g i=1 admit a formal series expansion in . Note that standard analytic results prove that they admit a Taylor expansion at any order thus determining the coefficients of the formal expansion. However, the convergence of the full series are not guaranteed. Moreover, for = 0, (q i ( = 0), p i ( = 0)) g i=1 are solutions of the algebraic equations ∂ ∂pi H i (p, q, t) = 0 and ∂ ∂qi H i (p, q, t) = 0. Sub-leading coefficients of the formal series expansions are then determined recursively using (A.3). Finally, since (p i ( ), q i ( )) g i=1 admit a formal series expansion in and since the -dependence only comes from these Darboux coordinates, we get that L(x, t, ) has a formal expansion given by (A.2). Indeed, by definition of the spectral Darboux coordinates of section 3.5 and the reconstruction of the element L(x, t) through Lagrange interpolation, L(x, t) has components that are rational functions of (p i ( ), q i ( )) g i=1 . Eventually, by definition, it implies that the auxiliary matrix A(x, t, ) also has a formal expansion given by (A.2).

Since ∂ x Ψ(x, t, ) = L(x, t, )Ψ(x, t, ) with L(x, t, ) having a formal series expansion in , it implies that Ψ(x, t, ) admit a formal WKB expansion given by equation (A.1). It also implies that log Ψ admits a formal series expansion in which is equivalent to say that Ψ admit a formal WKB expansion.We then conclude using its definition that M (x, t, ) admits a formal series expansion in . Note that this result is also a direct consequence of the equations ∂ x M (x, t, ) = [L(x, t, ), M (x, t, )] and ∂ t M (x, t, ) = [A(x, t, ), M (x, t, )] from which we may find recursively the coefficients of the formal series (see section A.4.1 for explicit formulas in the sl 2 (C) case). Finally, we conclude from definition 2.1 that the correlation functions admit an expansion in .

Remark A.1. We stress again that these series expansions, as well as the WKB expansion for Ψ(x, t, ) are only considered as formal series in . In particular, the issue of convergence of these series remains and is well-known to be a difficult question involving Stokes sectors, Borel summability and tedious analytical considerations [START_REF] Iwaki | Exact WKB analysis and cluster algebras[END_REF][START_REF] Iwaki | Exact WKB Analysis and Cluster Algebras II: Simple Poles, Orbifold Points, and Generalized Cluster Algebras[END_REF]. However, since our main goal is to prove the formal reconstruction of these expansions via the topological recursion (order by order in powers of ), we may consider only formal series in leaving out the issues of convergence of these series.

A.2 Proof of condition 3 : Parity of the formal series expansion of the correlation functions

Transposing (5.2) and (5.1), one finds out that, for a given matrix Γ ∈ G (independent of x and t), we may define L Γ (x, t, ) := ΓL(x, t, ) t Γ -1 , A Γ := ΓA(x, t, ) t Γ -1 and Ψ Γ (x, t, ) = Ψ(x, t, ) t Γ -1 .

These matrices satisfy

d dx Ψ Γ (x, t, ) = Ψ Γ (x, t, )L Γ (x, t, ) d dt Ψ Γ (x, t, ) = Ψ Γ (x, t, )A Γ (x, t, )
and are related via the corresponding compatibility equation

- dL Γ dt + dA Γ dx = [A Γ , L Γ ] ,
where the relation between A Γ and L Γ is the same as the one between A and L.

As an element of a reduced space, L Γ (x, t, ) is given by Darboux coordinates (p Γ i ( , t), q Γ i ( , t)) g i=1 satisfying the evolution equations

∀ i ∈ 1, g ,    ∂q Γ i ( ,t) ∂t k,ν,l = - ∂H k,ν,l (p Γ ,q Γ ,t) ∂p Γ i ∂p Γ i ( ,t) ∂t k,ν,l = ∂H k,ν,l (p Γ ,q Γ ,t) ∂q Γ i
since the transposition changes the sign of the Poisson bracket. Note that H k,ν,l (p Γ , q Γ , t) are the same Hamiltonians as the one computed from L(x, t, ) as function of the Darboux coordinates. Hence, these evolution equations are the same as the one for the Darboux coordinates describing the evolution of L(x, t, -). In addition, because they are given by the same polynomial equations, the values of (p i (t, -), q i (t, -))

g i=1 and p Γ i (t, ), q Γ i (t, ) g i=1 coincide for = 0.
31 This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.

PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0002260

Hence, for any Γ, p Γ i (t, ), q Γ i (t, ) g i=1 = (p i (t, -), q i (t, -)) g i=1 and L Γ (x, t, ) coincides with a point in the orbit of L(x, t, -) meaning that ∃ Γ ∈ SL L0,r 0 , L Γ (x, t, ) = L(x, t, -). This is a sufficient condition for the parity property required for the TT property. Let us be more specific and construct explicitly this matrix Γ for g = sl 2 .

A.2.1 Specific form of the Γ matrix in the sl 2 case

Depending on its rank, we can generically consider two possible leading orders L 0,r0 , namely L 0,r0 = σ 3 or L 0,r0 = σ + . In both cases, one can obtain information about the matrix Γ required to obtain the parity property.

Indeed, by definition, the matrix Γ satisfies the identity ΓL t = L † Γ where L † (x, t, ) = L(x, t, -). Since L 0,r0 is chosen independently of , t and x, we may project the last identity at x r0 and get ΓL t 0,r0 = L 0,r0 Γ (because L † 0,r0 = L 0,r0 since it is chosen fixed by the SL 2 action). -Case L 0,r0 = σ 3 : In this case, ΓL t

0,r0 = L 0,r0 Γ is equivalent to Γ 1,2 = Γ 2,1 = 0. Since Γ ∈ G = SL 2 , we end up with Γ(t, ) = α(t, ) 0 0 1 α(t, )
where α(t, ) is an appropriate function of t and to be determined. Note that this general form is in agreement with the Painlevé cases (except Painlevé I for which L 0,r0 = σ + that is studied below) studied in [START_REF] Iwaki | Painlevé equations, topological type property and reconstruction by the topological recursion[END_REF].

-Case L 0,r0 = σ + : In this case, the SL 2 action also implies that L 0,r0-1 = C + l(t, )σ + where C is a constant (i.e. independent of t, and x) matrix. The identity ΓL t 0,r0 = L 0,r0 Γ is equivalent to Γ 2,1 = Γ 1,2 and Γ 2,2 = 0. Since Γ ∈ G = SL 2 and Γ is defined up to a global sign (that does not change the value of the determinant), we get that we may always choose Γ 1,2 = Γ 2,1 = i. In other words,

Γ(t, ) = α(t, ) i i 0 . (A.4)
Using the additional information on L 0,r0-1 and projecting the identity ΓL t = L † Γ at order x r0-1 gives ΓC t -CΓ = l(t, -)σ + Γl(t, )Γσ -.

The l.h.s. and r.h.s. of the previous equation may be computed using (A.4)

ΓC t -CΓ = 0 α(t, )C 2,1 -2iC 1,1 2iC 1,1 -α(t, )C 2,1 0
l(t, -)σ + Γ -l(t, )Γσ -= i(l(t, -) -l(t, )) 0 0 0 .
Hence, we end up with the fact that l(t, ) = l(t, -) which is equivalent to say that l(t, ) is an even function of and α(t, )C 2,1 = 2iC 1,1 . Thus, in the case when C 2,1 = 0 we get

Γ = 2i C1,1 C2,1 i i 0 .
Note that this general form as well as the fact that l(t, ) is an even function of is verified in the case of Painlevé I developed in details in [START_REF] Iwaki | Painlevé equations, topological type property and reconstruction by the topological recursion[END_REF] where l(t, ) = q(t, ) is an even function of while (C 1,1 , C 2,1 ) = (0, 4) and Γ = 0 i i 0 .

In the case C 2,1 = 0, we only get C 1,1 = 0 so that C is proportional to σ + and thus may be absorbed in the l(t, ) factor in the definition of L 0,r0-1 . Thus, L 0,r0-1 is proportional to σ + and in particular det L(x, t, ) x→∞ = O(x 2r0-2 ), i.e. that L(x, t, ) is twice degenerate at infinity. 32 This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.

PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0002260

In this highly degenerate cases, the missing factor α(t, ) of Γ should be studied specifically but the previous study still proves that l(t, ) is always an even function of . Note that this degenerate case cannot happen for (ν, r 0 ) = (0, 1) since it would give L(x) = (x + c)σ + and thus det L(x) = 0.

A.3 Proof of property 1 : Genus 0 property

For a given t, the classical spectral curve is defined as det(yI 2 -L(x, t, )) =0 = 0 = y 2 + det L (0) (x, t).

Since det L (0) (x, t) is a rational function of x with poles at x ∈ {a 0 , a 1 , . . . , a n }, this defines a Riemann surface with possible punctures at those points. In order to prove the topological type property, we must prove that this always defines a Riemann surface of genus 0. In the case of sl 2 (C), the situation is particularly simple. Indeed, taking order 0 into the compatibility equations leads to L (0) , A (0) = 0.

In other words, all leading orders in of the differential systems must commute. In particular, this is equivalent to the fact that their eigenspaces coincide. Since they are elements of sl 2 (C), there are only two possible cases.

-0 is an eigenvalue of L (0) , and therefore it must be an eigenvalue of order 2. This means that L (0) and A (0) are proportional to σ + (remember that we took the leading order L 0,r0 equal to either σ 3 or σ + ). Projecting (5.1) at order 0 using (A.1) we get Ψ(0) ∂ x Ψ (-1) = L (0) Ψ(0) . Since Ψ(0) is invertible (Ψ ∈ G = SL 2 (C)) and Ψ (-1) is diagonal while L (0) is proportional to σ + , we end up with a contradiction. -L (0) admits two opposite non-zero eigenvalues ±s 0 . This means that each eigenspace is of dimension 1. Therefore, all matrices A (0) are also diagonal in this basis with two opposite eigenvalues (since they belong to sl 2 (C)). Thus, the matrices L (0) and A (0) are proportional to each other. In conclusion, the previous discussion leads to the fact that ∃ α(x, t) ∈ C / A (0) (x, t) = α(x, t)L (0) (x, t) where α(x, t) ∈ C * .

At the level of determinants, we end up with det A (0) (x, t) = α(x, t) 2 det L (0) (x, t), so that the spectral curve is given by

y 2 + det A (0) (x, t) α(x, t) 2 = 0.
From section 4.2, we know that there exists an auxiliary matrix A(x, t, ) of the form A(x, t, ) = Â(x, t, )

β(x) (A.5)
where β(x) = 1 or β(x) = xa i and Â(x, t, ) is a sl 2 valued degree 1 polynomial in x. This leads to a spectral curve given by y 2 + det Â(0) (x, t) β(x) 2 α(x, t) 2 = 0. Thus, we end up with the following theorem: Theorem A.1 (Genus 0). The family of classical spectral curves Σ t is of genus 0.

√ -det A (0) 1 2 (A (0) )1,2 √ -det A (0)        M (k) (x, t) 1,1 M (k) (x, t) 1,2 M (k) (x, t) 2,1   =              ∂ t M (k-1) (x, t) 1,1 - k-1 j=0
A (k-j) (x, t), M (j) (x, t) 1,1

∂ t M (k-1) (x, t) 1,2 - k-1 j=0
A (k-j) (x, t), M (j) (x, t) 1,2 ∂ t M (k-1) (x, t) 2,1 -k-1 j=0 A (k-j) (x, t), M (j) (x, t) 2,1 k-1 j=1 M (j) (x, t) 1,1 M (k-j) (x, t) 1,1 + M (j) (x, t) 1,2 M (k-j) (x, t) 2,1

             def =      RHS (k) 1 (x, t) RHS

RHS

(k) 3

RHS (k) 4        def =              β(x)∂ t M (k-1) (x, t) 1,1 - k-1 j=0
Â(k-j) (x, t), M (j) (x, t)

1,1 β(x)∂ t M (k-1) (x, t) 1,2 - k-1 j=0
Â(k-j) (x, t), M (j) (x, t)

1,2 β(x)∂ t M (k-1) (x, t) 2,1 - k-1 j=0
Â(k-j) (x, t), M (j) (x, t)

2,1 k-1 j=1 M (j) (x, t) 1,1 M (k-j) (x, t) 1,1 + M (j) (x, t) 1,2 M (k-j) (x, t) 2,1              . (A.8)
Since the quantity β(x) is no longer involved in denominators of equations (A.7) and (A.8) we get that the induction does not create any singularity at the zeros of β(x) and thus that the only possible singularities of M (k) (x, t) are at the branchpoints (i.e. at the zeros of det A (0) (x, t)) or at x = ∞). Note that the conclusion remains valid if β(x) is any arbitrary polynomial in x (and not necessarily of degree at most 1) since it simplifies in the computations. However, in this case, the control of the singularity at x = ∞ presented below may require some improvements.

35 This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset. In this section we prove that the correlation functions are regular at x = ∞. As we will see below, the proof differs depending on the order of det Â(0) (x, t) at x → ∞.

PLEASE CITE THIS ARTICLE AS

• Let us first discuss the case when det Â(0) (x, t) x→∞ = αx 2 + O(x) with α = 0. In other words det Â(0) (x, t) is a polynomial of degree 2 in x. Since entries Â(0) (x) 1≤i,j≤2 are at most linear in x, formula (A.6) immediately shows that M (0) (x) x→∞ = O(1) and therefore is regular at x = ∞. Then, an easy induction using (A.7) and the fact that β(x) is at most linear, proves that for all k ≥ 0: M (k) (x, t) x→∞ = O(1) and thus that the matrices M (k) k≥0 are regular at x = ∞. Indeed, by induction all numerators involved in (A.7) are at most quadratic in x and thus dividing by det Â(0) (x, t) leads to at most constant terms at infinity. Eventually, using the definition of the correlation functions (2.1), we conclude that the correlation functions are regular at x = ∞.

• Let us now discuss the case when det Â(0) (x, t) x→∞ = αx + O(1) with α = 0. Note that section 4.2 implies that this case may only happen for ν = 0, r 0 ≥ 2 and if L 0,r0 does not have full rank. In the sl 2 case, this is equivalent to say that L 0,r0 is of rank 1. Using the SL 2 action presented in section 3.6, we may always get to the situation L 0,r0 = σ + and L 0,r0-1 (t, ) = C + l(t, )L 0,r0 where C is a constant matrix (in the sense independent of x, t and ) . A typical example of this situation is the Painlevé I Lax pair (Cf. section 6.1.1). Formula (4.3) shows that we must have A (0) (x, t, ) = L 0,r0 x + L 0,r0-1 (t, ) + β 2 (t, )A 2,0,r0-1 . Moreover, since A 2,0,r0-1 = 2 [x -r0 L(x)] + = L 0,r0 , we end up with A(x, t, ) = L 0,r0 x + C + (l(t, ) + β 2 (t, )) L 0,r0 with L 0,r0 = σ + . (A.9) Formula (A.6) then immediately shows that we have

M (0) (x, t) = L 0,r0 x 1 2 + 1 2 I 2 + O x -1 2 .
We may now prove by induction that M (k) (x, t)

x→∞ = O x -1 2
for k ≥ 1. Indeed for k ≥ 1 we have:

β(x)∂ t M (k-1) (x, t) = O x 1 2 since β(x) x→∞ = O(x) A (k) , M (0) = O x 1 2 , ∀ k ≥ 1 A (k-j) , M (j) = O x 1 2 , ∀ k ≥ 1 , 1 ≤ j ≤ k -1 RHS (k) i = O x 1 2
, ∀ i ∈ {1, 2, 3} RHS (k) 4

= O x -1 (A.10)

Note that the first identity is valid for k = 1 since

∂ t M (0) (x, t) = O x -1 2
. The second and third identity hold since (A.9) gives A (k) (x, t) = (l(t, ) + β 2 (t, ))

(k) L 0,r0 ∝ L 0,r0 , so that the first two terms L 0,r0 x 1 2 and 1 2 I 2 do not contribute to the commutator A (k) , M (0) . The fourth identity is a direct consequence of the first two using definition (A.8) while the fifth one is obvious from the induction assumption and the fact that M (0) is not involved in the sum defining RHS

Remark 4 . 1 .

 41 The condition rν ≥ 2 is necessary in order to introduce the deformation t Lr ν (x-aν )2 without changing the value of rν. The preceding deformation allows to consider the simple pole case rν = 1.

  DOI:10.1063/5.0002260+ -dx(r0-1) [L(x, t)] 2,1 Ψ(x,

34

 34 This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset. PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.00022601 2 Â(0) (x, t) 1,1 RHS (k) 3 (x, t) -Â(0) (x, t) 2,1 -det Â(0) (x, t) RHS (k) 4 (x, t) , M (k) (x, t) 2,2 = -M (k) (x, t)

A. 4 . 2

 42 DOI:10.1063/5.0002260 Study at x = ∞

  , inserting (A.10) in (A.7) and the fact that Â(0) (x, t) = O(x) and det Â(0) (x, t) = αx+O[START_REF] Adams | Isospectral Hamiltonian Flows in Finite and Infinite Dimensions[END_REF] with α = 0 shows that M (k) = O x -1 2 thus concluding the induction. Finally, the definition of the correlation functions (2.1) (involving denominators1 xi-x σ(i) )shows that the correlation functions are regular at infinity (even ifM (0) (x, t) = O x 1 2 is not). 36This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.PLEASE CITE THIS ARTICLE ASDOI:10.1063/5.0002260

DOI:10.1063/5.0002260 5. Leading order in : For

  n ≥ 1, the correlators satisfy W n = O( n-2 ) as functions of . In other words, for any n ≥ 2, the coefficients W

	(k) n combining properties 1, 2 and 5 is equivalent to say that we have are identically zero. Note that 0≤k≤n-1
	Wn (z 1 , . . . , z n ) =	∞
		g=0

This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.PLEASE CITE THIS ARTICLE ASDOI:10.1063/5.0002260

Acknowledgements: Both authors are grateful to Kohei Iwaki for sharing with them a preliminary version of [34]. This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.

PLEASE CITE THIS ARTICLE AS

DOI:10.1063/5.0002260 N.O. would like to thank Mattia Cafasso, John Harnad and Volodya Roubstov for helpful discussions on the subject as well as Université de Lyon for its hospitality and the Center for Quantum Geometry of Moduli Spaces for its support. O.M. would like to thank Université Lyon 1, Université Jean Monnet and Institut Camille Jordan for material support. This work was supported by the LABEX MILYON (ANR-10-LABX-0070) of Université de Lyon, within the program « Investissements d'Avenir » (ANR-11-IDEX-0007) operated by the French National Research Agency (ANR).

A.4 Proof of property 4 : Pole property

A.4.1 Computation of M (k) k≥0

In the sl 2 (C) case, the correlation functions are defined as polynomials of the matrix M (x, t) = Ψ(x, t)e 1 Ψ(x, t) -1 . We remind here that even if e 1 does not belong to sl 2 , we can still use it for proving the TT properties as explained in remark 2.1. It is a projection of rank 1 satisfying M 2 = M with Tr M = 1. It admits a formal -expansion of the form

It satisfies the differential system

In order to prove that the correlation functions admit singularities only at the branchpoints, we will prove this property for M . The strategy follows the one already used in [START_REF] Iwaki | Painlevé equations, topological type property and reconstruction by the topological recursion[END_REF] for Painlevé equations. Using det M (0) = 0 and the time-differential equations, we get (See [START_REF] Iwaki | Painlevé equations, topological type property and reconstruction by the topological recursion[END_REF] for details)

(A.6) Thus, it is obvious that x → M (0) (x, t) may only be singular at the branchpoints with only at most pole singularity there (note in particular the factor β(x) in (A.5) always simplifies in the previous formula so that M (0) (x, t) is regular at the poles of A (0) (x, t)). Next, following [START_REF] Iwaki | Painlevé equations, topological type property and reconstruction by the topological recursion[END_REF], we may determine the matrices M (k) k≥1 recursively. However, we need to extend the strategy used in [START_REF] Iwaki | Painlevé equations, topological type property and reconstruction by the topological recursion[END_REF] in order to make it sufficient for our general case. The recursion provided in [START_REF] Iwaki | Painlevé equations, topological type property and reconstruction by the topological recursion[END_REF] is given by projecting the relation:

] and det M = 0 at order k . In [START_REF] Iwaki | Painlevé equations, topological type property and reconstruction by the topological recursion[END_REF], only the entries (1, 1) and (1, 2) of the first identity were used since the entries (2, 1) and (2, 2) are not independent from them. However, this strategy generates formulas involving

that may be problematic in the generic case (it was not necessary for the Lax pairs considered in [START_REF] Iwaki | Painlevé equations, topological type property and reconstruction by the topological recursion[END_REF] since the top right entries of the matrices A (0) considered there were simple) because they create fake poles at the zeros of A (0) 1,2 . In order to avoid such problems, one simply has to remember that the entry (2, 1) may also be used and provide alternative formulas. Following this strategy we obtain the over-determined system of equations for (

Thus we find the recursion

Generically, the last recursion may create singularities at the zeros of det A (0) (x, t), at x = ∞ and at singularities of the matrices A (k) (x, t) k≥0 . In our present case, we know that we have

where β(x) = 1 or β(x) = xa i and Â(x, t, ) is a sl 2 valued polynomial in x of degree at most 1. Thus, since det

, we may reduce the previous identities to

A.5 Proof of property 6 : Identification of the two point function

In order to prove the reconstruction by the topological recursion, we need to verify that ω 0,1 and ω 0,2 defined by the determinantal formulas identify with the data of the classical spectral curve. Using (2.1), this is equivalent to checking that

The first equation is automatically satisfied since

√

det A (0) A (0) (Cf. (A.6)) and the fact that y = -

as well as the fact that for B ∈ sl 2 : Tr B 2 = -2 det B.

• Let us verify that the second equation holds when Â(0) (x, t) = A 1 (t)x + A 0 (t) with det A 1 = 0. In this case we first observe that

Thus, we get det Â(0

where the two branchpoints satisfy

Thus, we have

a straightforward computation from (A.11) gives

• Let us verify that the second equation holds when Â(0) (x, t) = A 1 (t)x + A 0 (t) with det A 1 = 0. We get det Â(0) (x, t) = -Tr(A 0 A 1 )x + det A 0 = -Tr(A 0 A 1 )(xa) where the only branchpoint satisfies a = det A0 Tr(A0A1) . Note that Tr(A 0 A 1 ) = 0, otherwise the classical spectral curve would be trivial (no branchpoints). We have

37 This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.

PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0002260

Using the parametrization x(z) = a + z 2 giving x ′ (z) = 2z and x(z)a = z, we end up with

so that Tr(M (0) (x(z 1 ), t)M (0) (x(z 2 ), t))

We conclude that in both situations (one or two branchpoints), the existence of a linear auxiliary matrix presented in section 4.2 is sufficient to prove the identification of the initial data of the topological recursion (classical spectral curve and Bergman kernel) with their corresponding quantities given by the determinantal formulas.

A.6 Proof of property 5 : Leading order in

This follows from the other properties using the double recursion introduced in [START_REF] Iwaki | Painlevé equations, topological type property and reconstruction by the topological recursion[END_REF].