Olivier Marchal
email: 1olivier.marchal@univ-st-etienne.fr

Mohamad Alameddine
email: 2mohamad.alameddine@univ-st-etienne.fr

Hamiltonian representation of isomonodromic deformations of twisted rational connections: The Painlevé 1 hierarchy 1

Keywords: 6, 2 Symmetric Darboux coordinates, isomonodromic times, coordinates

In this paper, we build the Hamiltonian system and the corresponding Lax pairs associated to a twisted connection in gl 2 (C) admitting an irregular and ramified pole at infinity of arbitrary degree, hence corresponding to the Painlevé 1 hierarchy. We provide explicit formulas for these Lax pairs and Hamiltonians in terms of the irregular times and standard 2g Darboux coordinates associated to the twisted connection. Furthermore, we obtain a map that reduces the space of irregular times to only g non-trivial isomonodromic deformations. In addition, we perform a symplectic change of Darboux coordinates to obtain a set of symmetric Darboux coordinates in which Hamiltonians and Lax pairs are polynomial. Finally, we apply our general theory to the first cases of the hierarchy: the Airy case (g = 0), the Painlevé 1 case (g = 1) and the next two elements of the Painlevé 1 hierarchy.

3 Classical spectral curve and connection with topological recursion [START_REF] Dumitrescu | Quantization of spectral curves for meromorphic Higgs bundles through topological recursion[END_REF] 1 Introduction and summary of the results

Isomonodromic deformations have been studied since the beginning of the twentieth century [START_REF] Picard | Mémoire sur la théorie des fonctions algébriques de deux variables[END_REF][START_REF] Fuchs | Sur quelques équations différentielles linéaires du second ordre[END_REF][START_REF] Garnier | Solution du problème de riemann pour les systèmes différentiels linéaires du second ordre[END_REF][START_REF] Painlevé | Sur les équations différentielles du second ordre et d'ordre supérieur dont l'intégrale générale est uniforme[END_REF][START_REF] Gambier | Sur les équations différentielles du second ordre et du premier degré dont l'intégrale générale est à points critiques fixes[END_REF][START_REF] Schlesinger | Über eine klasse von differentialsystemen beliebiger ordnung mit festen kritischen punkten[END_REF] and is still currently an active domain in modern mathematics. If the initial restriction to Fuchsian singularities is now very well understood, many questions in the case of irregular singularities remain open. If a geometrical understanding of the Hamiltonian representation of the isomonodromic equations for a generic meromorphic connection in gl d (C) with d ≥ 2 is now well understood [START_REF] Hurtubise | On the geometry of isomonodromic deformations[END_REF][START_REF] Bertola | Hamiltonian structure of rational isomonodromic deformation systems[END_REF], the explicit expression for the Hamiltonians and Lax pairs was derived on a case by case basis until some very recent results. In [START_REF] Gaiur | Isomonodromic deformations: Confluence, reduction & quantisation[END_REF], the authors obtained explicit expressions for the Hamiltonians using confluences of isomonodromic deformations of Fuchsian systems. In particular, this method may only obtain results for deformations obtained by confluences of simple poles limiting its range of application for d ≥ 3. Independently, in [START_REF] Marchal | Hamiltonian representation of isomonodromic deformations of general rational connections on[END_REF], the authors proposed a generic construction and some explicit formulas for the Lax pairs and Hamiltonians associated to meromorphic connections in gl 2 (C) such that all leading orders at each pole are assumed to be diagonalizable. 1 In addition, they proposed an explicit map from the geometric set of irregular times (defined in [START_REF] Boalch | Symplectic manifolds and isomonodromic deformations[END_REF][START_REF] Boalch | Simply-laced isomonodromy systems[END_REF][START_REF] Boalch | Twisted local wild mapping class groups: configuration spaces, fission trees and complex braids[END_REF]) to a smaller set of isomonodromic times complemented by a set of trivial times and showed that the Darboux coordinates are independent of the trivial times. Thus, it is natural to wonder if the method of [START_REF] Marchal | Hamiltonian representation of isomonodromic deformations of general rational connections on[END_REF] can be extended to the case where a meromorphic connection in gl 2 (C) exhibits a pole whose leading order cannot be diagonalized.

The main purpose of this article is to provide a positive answer to this question and to obtain some explicit expressions for both the Hamiltonian system and the Lax pairs in the so-called "twisted case", i.e. for meromorphic connections in gl 2 (C) such that the leading order of the connection at a pole is non-diagonalizable. Such poles are also referred to as "ramified poles" in the literature. In [START_REF] Marchal | Hamiltonian representation of isomonodromic deformations of general rational connections on[END_REF] results indicated that the formulas are independent at each pole so that we focus, without loss of generality, to the case of only one ramified pole at infinity in this paper (the position of the pole playing no role). In the end, combining the results of the present paper and those of [START_REF] Marchal | Hamiltonian representation of isomonodromic deformations of general rational connections on[END_REF] completes the study of all meromorphic connections in gl 2 (C).

Let us emphasize that the "twisted case" requires a specific and non-trivial analysis. Indeed, the underlying geometry, in particular the definition of the irregular times at a ramified pole, is more difficult and less understood than the non-ramified case. Main results in this area are [START_REF] Boalch | Symplectic manifolds and isomonodromic deformations[END_REF][START_REF] Biquard | Wild non-abelian hodge theory on curves[END_REF][START_REF] Boalch | Simply-laced isomonodromy systems[END_REF][START_REF] Boalch | Twisted local wild mapping class groups: configuration spaces, fission trees and complex braids[END_REF] and shall be used throughout the article. One of the main difference in the twisted case is the necessity to introduce a ramified cover around each pole with some associated local coordinates in order to be able to "diagonalize" the singular part of the connection around the pole. Moreover, the definition of irregular times differs since for example the eigenvalues of the leading order of the Lax matrix are necessarily the same (because the matrix is assumed to be non-diagonalizable) so that the dimension of the space of irregular times and associated deformations drastically change. All these important changes require a detailed analysis of the twisted case that we propose in the present paper.

In particular, our main results that can be seen as a summary and plan of the article are:

• For any isomonodromic deformation, characterized by a vector α in the tangent space, we provide an explicit gauge transformation between the geometric Lax pair L(λ), Ãα (λ) and the companion-like Lax pair (L(λ), A α (λ)) in terms of apparent singularities (q i ) 1≤i≤g , their dual coordinates (p i ) 1≤i≤g and the irregular times t in Proposition 2.2.

• A general expression of the companion-like Lax pair (L(λ), A α (λ)) in terms of the Darboux coordinates (q i , p i ) 1≤i≤g is given in Propositions 2.4, 4.2 and 4.3, complemented with equation [START_REF] Bergère | Determinantal formulae and loop equations[END_REF][START_REF] Bertola | Hamiltonian structure of rational isomonodromic deformation systems[END_REF]. These results follow from the local asymptotics at infinity of the wave matrix obtained in Proposition 2.3 following the geometric construction of the twisted meromorphic connection.

• Explicit expressions of the evolutions of the Darboux coordinates relatively to irregular times and a proof that these evolutions are indeed Hamiltonian with an explicit expression of the latter are given in Theorem 5.1.

• A symplectic change from Darboux coordinates (q i , p i ) 1≤i≤g to the set of symmetric Darboux coordinates (Q i , P i ) 1≤i≤g for which the geometric Lax pair L(λ), Ãα (λ) and the Hamiltonians are polynomial is provided in Definition 6.2. The explicit polynomial expressions are given in Theorem 6.1 and Propositions 6.3 and 6.4.

• A natural reduction of the space of deformations relatively to irregular times (of dimension 2g + 4) to a subspace of non-trivial isomonodromic deformations (of dimension g) complemented by a space of trivial deformations is detailed in Section 7. Note in particular that this map is explicit both at the level of the tangent space (Definition 7.1) and at the level of times (Definition 7.2). The terminology "trivial deformations" stands for the fact that the evolutions of the (shifted) Darboux coordinates (q i , pi ) 1≤i≤g relatively to these times are proven trivial in Theorem 7.2.

• Simpler formulas in terms of the Darboux coordinates for the companion-like Lax pair (L(λ), A α (λ)) and Hamiltonians, after a canonical choice of the trivial times (given in Definition 8.1), are provided in Proposition 8.2 and Theorem 8.1.

• Some simpler polynomial expressions in terms of the symmetric Darboux coordinates of the geometric Lax pair L(λ), Ãα (λ) and Hamiltonians, after a canonical choice of the trivial times (defined in Definition 8.1), are provided in Proposition 8.1 and Theorem 8.1.

• The connection with the quantization of classical spectral curves via the topological recursion of [START_REF] Eynard | Invariants of algebraic curves and topological expansion[END_REF] is presented as a by-product in Section 3.

2 Twisted meromorphic connections at infinity

Twisted meromorphic connections and irregular times

The space of gl 2 (C) meromorphic connections has been studied from many different perspectives.

In the present article, we shall mainly follow the point of view of the Montréal group [START_REF] Adams | Darboux coordinates and Liouville-Arnold integration in loop algebras[END_REF][START_REF] Adams | Darboux coordinates on coadjoint orbits of lie algebras[END_REF] together with some insight from the work of P. Boalch [START_REF] Boalch | Symplectic manifolds and isomonodromic deformations[END_REF]. Let us first define the space we shall be studying.

Definition 2.1 (Space of meromorphic connections with a pole at infinity). Let r ∞ ≥ 3 be a given integer. We shall consider

F ∞,r∞ := L(λ) = r∞-1 k=1 L[∞,k] λ k-1 / { L[∞,k] } ∈ (gl 2 (C)) r∞-1 /GL 2 (2-1)
where GL 2 acts simultaneously by conjugation on all the coefficients { L[∞,k] } 1≤k≤r∞-1 . The corresponding meromorphic connection is defined by

d Ψ = L(λ)dλ Ψ ⇔ ∂ λ Ψ = L(λ) Ψ (2-2)
where Ψ is referred to as the wave matrix.

The space F ∞,r∞ corresponds to the space of meromorphic connections with a pole at infinity whose order is prescribed by the integer r ∞ . Since r ∞ ≥ 3, the pole at infinity is said to be irregular. The generic case where the leading order L[∞,r∞-1] is diagonalizable has been studied in [START_REF] Marchal | Hamiltonian representation of isomonodromic deformations of general rational connections on[END_REF] where a complete construction of the associated Hamiltonian systems is provided. In this article, we shall deal with the so-called "twisted" case of [START_REF] Boalch | Twisted local wild mapping class groups: configuration spaces, fission trees and complex braids[END_REF]. It corresponds to the case where the leading order L[∞,r∞-1] is assumed to be non-diagonalizable. In the literature, this case is also referred to as "ramified" at infinity. We introduce the following definition. Definition 2.2 (Set of twisted meromorphic connections at infinity). Let r ∞ ≥ 3 be a given integer. We shall consider the subset of F∞,r∞ defined by

F∞,r∞ = L(λ) = r∞-1 k=1 L[∞,k] λ k-1 / { L[∞,k] } ∈ (gl 2 (C)) r∞-1 and L[∞,r∞-1] is not diagonalizable /GL 2 (2-3)
F ∞,r∞ can be given a Poisson structure inherited from the Poisson structure of a corresponding loop algebra and this space has been intensively studied from the point of view of isospectral and isomonodromic deformations. Following P. Boalch's works, we can use the Poisson structure on F R,r∞ in order to describe it as a bundle whose fibers are symplectic leaves obtained by fixing the irregular type and monodromies of L(λ) ∈ F∞,r∞ . The general theory for the non-twisted or twisted case has been described in [START_REF] Boalch | Symplectic manifolds and isomonodromic deformations[END_REF][START_REF] Boalch | Simply-laced isomonodromy systems[END_REF][START_REF] Boalch | Twisted local wild mapping class groups: configuration spaces, fission trees and complex braids[END_REF]. Let us briefly review this perspective in the twisted setting and use it to define local coordinates on F∞,r∞ trivializing the fibration.

The main difference when dealing with the twisted case at infinity is that one needs to introduce a two-sheeted cover above infinity and define the local coordinate at infinity by

z ∞ (λ) = λ -1 2 (2-4)
The general theory of [START_REF] Boalch | Twisted local wild mapping class groups: configuration spaces, fission trees and complex braids[END_REF] implies the following proposition. . For any given L(λ) in an orbit of F∞,r∞ , there exists a local gauge matrix G ∞ (z) around ∞ such that

G ∞ (z) = G ∞,-1 z + G ∞,0 + ∞ k=1 G ∞,k z -k with G ∞,-1 of rank 1 (2-5)
and

• Ψ ∞ (z) def = G ∞ (z)
Ψ is a formal fundamental solution, also known as a Turritin-Levelt fundamental form (or Birkhoff factorization):

Ψ ∞ (λ) = Ψ (reg) ∞ (z) diag exp - 2r∞-2 k=1 t ∞,k k z k + 1 2 ln z , exp - 2r∞-2 k=1 (-1) k t ∞,k k z k + 1 2 ln z = Ψ (reg) ∞ (z) diag exp - 2r∞-2 k=1 t ∞,k k λ k 2 + 1 4 ln λ , exp - 2r∞-2 k=1 (-1) k t ∞,k k λ k 2 + 1 4 ln λ (2-6)
where Ψ

(reg)

∞ (z) ∈ GL 2 [[z -1 ]] is holomorphic at z = ∞. • The associated Lax matrix L ∞ = G ∞ LG -1 ∞ + (∂ λ G ∞ )G -1 ∞ has a diagonal singular part at ∞: L ∞ (λ) = diag - 1 2 2r∞-2 k=1 (-1) k t ∞,k z k-2 + 1 4z 2 , - 1 2 2r∞-2 k=1 (-1) k t ∞,k z k-2 + 1 4z 2 + O(1) = diag - 1 2 2r∞-2 k=1 t ∞,k λ k 2 -1 + 1 4λ , - 1 2 2r∞-2 k=1 (-1) k t ∞,k λ k 2 -1 + 1 4λ + O(1)
(2-7)

The complex numbers (t ∞,k ) 1≤k≤2r∞-2 define the "irregular times" at infinity that we shall

denote t = {(t ∞,k ) 1≤k≤2r∞-2 } the irregular type of L ∈ F∞,r∞ .
Remark 2.1. The coefficients t = (t ∞,k ) 1≤k≤2r∞-2 are also referred to as "spectral times" or "KP times" in part of the literature.

Note that the diagonal part could be expressed as diag(t

∞ (1) ,k , t ∞ (2) ,k ) with t ∞ (2) ,k = (-1) k t ∞ (1)
,k . This immediately follows from the fact that Tr L = Tr L ∞ + O(1) and by definition Tr L may only involve even powers of z. In the same way, we also have the following remark. [START_REF] Adams | Darboux coordinates on coadjoint orbits of lie algebras[END_REF][START_REF] Bergère | Rational differential systems, loop equations, and application to the q th reductions of KP[END_REF][START_REF] Bergère | Determinantal formulae and loop equations[END_REF][START_REF] Bertola | Hamiltonian structure of rational isomonodromic deformation systems[END_REF][START_REF] Bertola | The partition function of the two-matrix model as an isomonodromic τ function[END_REF][START_REF] Biquard | Wild non-abelian hodge theory on curves[END_REF][START_REF] Boalch | Symplectic manifolds and isomonodromic deformations[END_REF][START_REF] Boalch | Simply-laced isomonodromy systems[END_REF][START_REF] Boalch | Twisted local wild mapping class groups: configuration spaces, fission trees and complex braids[END_REF][START_REF] Chekhov | Hermitean matrix model free energy: Feynman graph technique for all genera[END_REF][START_REF] Chekhov | Hermitian matrix model free energy: Feynman graph technique for all genera[END_REF][START_REF] Chekhov | Matrix eigenvalue model: Feynman graph technique for all genera[END_REF]. Similarly, note that the determinant det L ∞ only involves even powers of z and satisfies

det L ∞ = 1 4 t 2 ∞,2r∞-2 λ 2r∞-4 + 1 4 t 2 ∞,2r∞-3 + 1 2 t ∞,2r∞-2 t ∞,2r∞-4 λ 2r∞-5 + O λ 2r∞-6 (2-8) so that from det L ∞ = det( L + G -1 ∞ ∂ λ G ∞ ) the coefficients of L[r∞-1]
as well as the upper-right coefficient of L[r∞-1] (that is necessarily 1) shall be fixed in (2-13).

Choice of representative normalized at infinity

Fixing the irregular type of L(λ) does not fix it uniquely. In fact, the space M∞,r∞,t := L(λ) ∈ F∞,r∞ / L(λ) has irregular type t

(2-9) is a symplectic manifold of dimension dim MR,r,t = 2r ∞ -6 = 2g (2-10) where g := r ∞ -3 (2-11)
is the genus of the spectral curve defined by det(yI 2 -L(λ)) = 0.

For any value of the irregular times, the Montréal group introduced a set of local Darboux coordinates (q i , p i ) 1≤i≤g on M∞,r∞,t . Indeed, in each orbit in F∞,r∞ , the global action of GL 2 (C) implies that we may choose the leading coefficient L[r∞-1] as a lower triangular matrix with identical coefficients on the diagonal (which is the standard form for a non-diagonalizable matrix of size 2). Furthermore, the remaining action allows to fix the coefficients on the diagonal of the subleading order L[r∞-2] at equal values. Combining this choice with Remark 2.2, we obtain the existence of a unique element for which L(λ) is of the form

L(λ) = -1 2 t ∞,2r∞-2 0 1 4 (t ∞,2r∞-3 ) 2 -1 2 t ∞,2r∞-2 λ r∞-2 + -1 2 t ∞,2r∞-4 1 X -1 2 t ∞,2r∞-4 λ r∞-3 +O λ r∞-4 .
(2-12) One may thus identify M∞,r∞,t with the space of such representatives

M∞,r∞,t ∼ L(λ) = r∞-1 k=1 L[∞,k] λ k-1 / L[∞,r∞-1] = -1 2 t ∞,2r∞-2 0 1 4 (t ∞,2r∞-3 ) 2 -1 2 t ∞,2r∞-2 and L[∞,r∞-2] = -1 2 t ∞,2r∞-4 1 δ ∞ -1 2 t ∞,2r∞-4 , δ ∞ ∈ C . (2-13)
In the following, we shall use the notation L(λ) whenever we consider such a representative and call it a representative "normalized at infinity".

Darboux coordinates

The work of the Montréal group implies that the space M∞,r∞,t is a symplectic manifold of dimension 2r ∞ -6 = 2g. Consequently, one may define a set of Darboux coordinates (q i , p i ) 1≤i≤g on M∞,r∞,t that we shall present in this section. Let L(λ) ∈ M∞,r∞,t be a representative of the form described above in eq. [START_REF] Adams | Darboux coordinates on coadjoint orbits of lie algebras[END_REF][START_REF] Bergère | Rational differential systems, loop equations, and application to the q th reductions of KP[END_REF][START_REF] Bergère | Determinantal formulae and loop equations[END_REF][START_REF] Bertola | Hamiltonian structure of rational isomonodromic deformation systems[END_REF][START_REF] Bertola | The partition function of the two-matrix model as an isomonodromic τ function[END_REF][START_REF] Biquard | Wild non-abelian hodge theory on curves[END_REF][START_REF] Boalch | Symplectic manifolds and isomonodromic deformations[END_REF][START_REF] Boalch | Simply-laced isomonodromy systems[END_REF][START_REF] Boalch | Twisted local wild mapping class groups: configuration spaces, fission trees and complex braids[END_REF][START_REF] Chekhov | Hermitean matrix model free energy: Feynman graph technique for all genera[END_REF][START_REF] Chekhov | Hermitian matrix model free energy: Feynman graph technique for all genera[END_REF][START_REF] Chekhov | Matrix eigenvalue model: Feynman graph technique for all genera[END_REF]. By definition, the entry L(λ)

1,2
is a monic polynomial function of λ of degree r ∞ -3 = g. We thus define (q i ) 1≤i≤g as the g zeroes of L(λ)

1,2 ∀ i ∈ 1, g : L(q i ) 1,2 = 0. (2-14)
This defines half of the spectral Darboux coordinates. The second half is obtained by evaluating the entry L(λ)

1,1

at λ = q i , ∀ i ∈ 1, g : p i := L(q i ) 1,1 . (2-15)
Let us remark that, for any i ∈ 1, g , the pair (q i , p i ) is by definition a point on the spectral curve defined by det(yI 2 -L(λ)) = 0. In other words, we have

∀ i ∈ 1, g : det(p i I 2 -L(q i )) = 0. (2-16)
As in the non-twisted case,the previous construction provides a local description of the space FR,r as a trivial bundle F∞,r∞ → B where the base B = t is the set of irregular times. The fiber above a point t ∈ B is M∞,r∞,t that we equip with spectral Darboux coordinates (q i , p i ) 1≤i≤g .

The space B is a space of isomonodromic deformations meaning that any vector field ∂ t ∈ T t B gives rise to a deformation of L(λ) preserving its generalized monodromy data. There exist different equivalent ways to characterize the property of being an isomonodromic vector field. The one that we shall use in this article is the existence of a compatible system of the form

∂ λ Ψ(λ, t) = L(λ) Ψ(λ, t) ∂ t Ψ(λ, t) = Ãt (λ) Ψ(λ, t) (2-17)
where Ãt (λ) is a polynomial function of λ with a pole at infinity lower or equal to r ∞ -2 (the order of the pole at infinity of L). Equations [START_REF] Adams | Darboux coordinates on coadjoint orbits of lie algebras[END_REF][START_REF] Bergère | Rational differential systems, loop equations, and application to the q th reductions of KP[END_REF][START_REF] Bergère | Determinantal formulae and loop equations[END_REF][START_REF] Bertola | Hamiltonian structure of rational isomonodromic deformation systems[END_REF][START_REF] Bertola | The partition function of the two-matrix model as an isomonodromic τ function[END_REF][START_REF] Biquard | Wild non-abelian hodge theory on curves[END_REF][START_REF] Boalch | Symplectic manifolds and isomonodromic deformations[END_REF][START_REF] Boalch | Simply-laced isomonodromy systems[END_REF][START_REF] Boalch | Twisted local wild mapping class groups: configuration spaces, fission trees and complex braids[END_REF][START_REF] Chekhov | Hermitean matrix model free energy: Feynman graph technique for all genera[END_REF][START_REF] Chekhov | Hermitian matrix model free energy: Feynman graph technique for all genera[END_REF][START_REF] Chekhov | Matrix eigenvalue model: Feynman graph technique for all genera[END_REF][START_REF] Dumitrescu | Quantization of spectral curves for meromorphic Higgs bundles through topological recursion[END_REF][START_REF] Eynard | Quantization of classical spectral curves via topological recursion[END_REF][START_REF] Eynard | Invariants of algebraic curves and topological expansion[END_REF][START_REF] Eynard | Topological recursion in random matrices and enumerative geometry[END_REF] are referred to as a Lax pair whose compatibility condition is

∂ λ Ãt (λ) -∂ t L(λ) + L(λ), Ãt (λ) = 0. (2-18)

Scalar differential equation and companion gauge

Let us now consider an orbit in FR,r and a representative L(λ) of this orbit normalized at infinity as above. Let Ψ(λ) be a wave matrix solution to the linear system

∂ λ Ψ(λ) = L(λ) Ψ(λ). (2-19)
The differential system ∂ λ Ψ(λ) = L(λ) Ψ(λ) may be rewritten into a scalar differential equation for Ψ1,1 that is equivalent to a companion like matrix system. More precisely, defining

Ψ(λ) = G(λ) Ψ(λ) with G(λ) = 1 0 L1,1 L1,2 (2-20)
we get that Ψ is a solution of the companion-like system

∂ λ Ψ(λ) = L(λ)Ψ(λ) with L(λ) = 0 1 L 2,1 L 2,2 (2-21)
given by

L 2,1 = -det L + ∂ λ L1,1 -L1,1 ∂ λ L1,2 L1,2 , L 2,2 = Tr L + ∂ λ L1,2 L1,2 . (2-22)
Note in particular that the first line of Ψ and Ψ is obviously the same:

Ψ 1,1 = Ψ1,1 def = ψ 1 and Ψ 1,2 = Ψ1,2 def = ψ 2 so that we immediately get Ψ(λ) = Ψ1,1 (λ) Ψ1,2 (λ) ∂ λ Ψ1,1 (λ) ∂ λ Ψ1,2 (λ) = ψ 1 (λ) ψ 2 (λ) ∂ λ ψ 1 (λ) ∂ λ ψ 2 (λ)
.

(2-23)

The companion-like system (2-21) is equivalent to say that ψ 1 and ψ 2 satisfy the linear ODE:

[∂ λ ] 2 -L 2,2 (λ)∂ λ -L 2,1 (λ) ψ i = 0 (2-24)
which is sometimes referred to as the "quantum curve".

Introduction of a scaling parameter

In order to make the connection with formal -transseries appearing in the quantization of classical spectral curves via topological recursion of [START_REF] Eynard | Invariants of algebraic curves and topological expansion[END_REF], we shall also introduce a formal parameter by a simple rescaling of the irregular times.

t ∞,k → k 2 -1 t ∞,k , ∀ k ∈ 1, 2r ∞ -2 , λ → -1 λ (2-25)
This very simple rescaling implies that the differential system reads

∂ λ Ψ(λ, ) = L(λ, ) Ψ(λ, ). (2-26)
However, for readers uneasy with this additional parameter, we stress here that may be fixed to 1 in the rest of the paper except for Section 3.

Explicit expressions of the gauge transformation

Using the Darboux coordinates (q i , p i ) 1≤i≤g and the irregular times t, one may obtain the explicit expression of the gauge transformation relating Ψ and Ψ. In order to do so, we shall introduce an intermediate wave matrix Ψ for the following proposition.

Proposition 2.2. The matrices Ψ and Ψ are related by the gauge transformations

Ψ(λ, ) = G 1 (λ, ) Ψ(λ, ) with G 1 (λ, ) = 1 0 1 2 t ∞,2r∞-2 λ + g 0 1 Ψ(λ, ) = J(λ, )Ψ(λ, ) with J(λ, ) =    1 0 Q(λ, ) g j=1 (λ-q j ) 1 g j=1 (λ-q j )    (2-27)
where Q is the unique polynomial in λ of degree g -1 such that (with the convention that empty products are set to 1)

Q(q i , ) = -p i , ∀ i ∈ 1, g (2-28)
i.e.

Q(λ, ) = - g i=1 p i j =i λ -q j q i -q j (2-29)
and the coefficient g 0 is given by

g 0 = 1 2 t ∞,2r∞-4 + 1 2 t ∞,r∞-2 g j=1 q j (2-30)
Proof. The proof consists in observing that

G(λ) = (G 1 (λ, )J(λ, )) -1 =   0 1 -Q(λ) -1 2 t ∞,2r∞-2 λ + g 0 g j=1 (λ -q j ) g j=1 (λ -q j )   (2-31)
recovers the matrix [START_REF] Adams | Darboux coordinates on coadjoint orbits of lie algebras[END_REF][START_REF] Bergère | Rational differential systems, loop equations, and application to the q th reductions of KP[END_REF][START_REF] Bergère | Determinantal formulae and loop equations[END_REF][START_REF] Bertola | Hamiltonian structure of rational isomonodromic deformation systems[END_REF][START_REF] Bertola | The partition function of the two-matrix model as an isomonodromic τ function[END_REF][START_REF] Biquard | Wild non-abelian hodge theory on curves[END_REF][START_REF] Boalch | Symplectic manifolds and isomonodromic deformations[END_REF][START_REF] Boalch | Simply-laced isomonodromy systems[END_REF][START_REF] Boalch | Twisted local wild mapping class groups: configuration spaces, fission trees and complex braids[END_REF][START_REF] Chekhov | Hermitean matrix model free energy: Feynman graph technique for all genera[END_REF][START_REF] Chekhov | Hermitian matrix model free energy: Feynman graph technique for all genera[END_REF][START_REF] Chekhov | Matrix eigenvalue model: Feynman graph technique for all genera[END_REF][START_REF] Dumitrescu | Quantization of spectral curves for meromorphic Higgs bundles through topological recursion[END_REF][START_REF] Eynard | Quantization of classical spectral curves via topological recursion[END_REF][START_REF] Eynard | Invariants of algebraic curves and topological expansion[END_REF][START_REF] Eynard | Topological recursion in random matrices and enumerative geometry[END_REF][START_REF] Fuchs | Sur quelques équations différentielles linéaires du second ordre[END_REF][START_REF] Gaiur | Isomonodromic deformations: Confluence, reduction & quantisation[END_REF][START_REF] Gambier | Sur les équations différentielles du second ordre et du premier degré dont l'intégrale générale est à points critiques fixes[END_REF]. Indeed, we first have that G2,2 (λ) = L1,2 (λ) and by definition L1,2 is a monic polynomial of degree g = r ∞ -3 with zeroes given by (q i ) 1≤i≤g . Similarly, the entry G2,1 (λ) is polynomial in λ of degree g + 1 = r ∞ -2. Moreover, it satisfies G2,1 (q i ) = p i for all i ∈ 1, g because of [START_REF] Adams | Darboux coordinates on coadjoint orbits of lie algebras[END_REF][START_REF] Bergère | Rational differential systems, loop equations, and application to the q th reductions of KP[END_REF][START_REF] Bergère | Determinantal formulae and loop equations[END_REF][START_REF] Bertola | Hamiltonian structure of rational isomonodromic deformation systems[END_REF][START_REF] Bertola | The partition function of the two-matrix model as an isomonodromic τ function[END_REF][START_REF] Biquard | Wild non-abelian hodge theory on curves[END_REF][START_REF] Boalch | Symplectic manifolds and isomonodromic deformations[END_REF][START_REF] Boalch | Simply-laced isomonodromy systems[END_REF][START_REF] Boalch | Twisted local wild mapping class groups: configuration spaces, fission trees and complex braids[END_REF][START_REF] Chekhov | Hermitean matrix model free energy: Feynman graph technique for all genera[END_REF][START_REF] Chekhov | Hermitian matrix model free energy: Feynman graph technique for all genera[END_REF][START_REF] Chekhov | Matrix eigenvalue model: Feynman graph technique for all genera[END_REF][START_REF] Dumitrescu | Quantization of spectral curves for meromorphic Higgs bundles through topological recursion[END_REF][START_REF] Eynard | Quantization of classical spectral curves via topological recursion[END_REF]. Finally its leading coefficients at infinity are

G2,1 (λ) = - 1 2 t ∞,2r∞-2 λ r∞-2 +   1 2 t ∞,r∞-2 g j=1 q j -g 0   λ r∞-3 + O λ r∞-4
(2-32) so that taking

g 0 = 1 2 t ∞,2r∞-4 + 1 2 t ∞,r∞-2 g j=1 q j (2-33) provides G2,1 (λ) = -1 2 t ∞,2r∞-2 λ r∞-2 -1 2 t ∞,2r∞-4 λ r∞-3 + O λ r∞-4
. Hence, with this choice of g 0 , it is equal to L1,1 (λ). Consequently, G(λ) recovers the matrix G(λ) of equation [START_REF] Adams | Darboux coordinates on coadjoint orbits of lie algebras[END_REF][START_REF] Bergère | Rational differential systems, loop equations, and application to the q th reductions of KP[END_REF][START_REF] Bergère | Determinantal formulae and loop equations[END_REF][START_REF] Bertola | Hamiltonian structure of rational isomonodromic deformation systems[END_REF][START_REF] Bertola | The partition function of the two-matrix model as an isomonodromic τ function[END_REF][START_REF] Biquard | Wild non-abelian hodge theory on curves[END_REF][START_REF] Boalch | Symplectic manifolds and isomonodromic deformations[END_REF][START_REF] Boalch | Simply-laced isomonodromy systems[END_REF][START_REF] Boalch | Twisted local wild mapping class groups: configuration spaces, fission trees and complex braids[END_REF][START_REF] Chekhov | Hermitean matrix model free energy: Feynman graph technique for all genera[END_REF][START_REF] Chekhov | Hermitian matrix model free energy: Feynman graph technique for all genera[END_REF][START_REF] Chekhov | Matrix eigenvalue model: Feynman graph technique for all genera[END_REF][START_REF] Dumitrescu | Quantization of spectral curves for meromorphic Higgs bundles through topological recursion[END_REF][START_REF] Eynard | Quantization of classical spectral curves via topological recursion[END_REF][START_REF] Eynard | Invariants of algebraic curves and topological expansion[END_REF][START_REF] Eynard | Topological recursion in random matrices and enumerative geometry[END_REF][START_REF] Fuchs | Sur quelques équations différentielles linéaires du second ordre[END_REF][START_REF] Gaiur | Isomonodromic deformations: Confluence, reduction & quantisation[END_REF][START_REF] Gambier | Sur les équations différentielles du second ordre et du premier degré dont l'intégrale générale est à points critiques fixes[END_REF] ending the proof.

Remark 2.3. By definition, the matrix Ψ(λ, ) satisfies the Lax system:

∂ λ Ψ(λ, ) = Ľ(λ, ) Ψ(λ, ) ∂ t Ψ(λ, ) = Ǎt (λ, ) Ψ(λ, ) (2-34)
for any irregular time t ∈ t. In particular, the corresponding Lax matrices Ľ(λ, ) and Ǎt (λ, ) are given by

Ľ(λ, ) = J(λ, )L(λ, )J -1 (λ, ) + (∂ λ J(λ, ))J -1 (λ, ) Ǎt (λ, ) = J(λ, )A t (λ, )J -1 (λ, ) + ∂ t (J(λ, ))J -1 (λ, ) (2-35)
and are polynomial functions of λ with no singularities at λ ∈ {q 1 , . . . , q g }.

Note that by definition, the entries of Ľ are related to those of L by

Ľ1,1 (λ, ) = -Q(λ, ), Ľ1,2 (λ, ) = g j=1 (λ -q j ), Ľ2,2 (λ, ) = L 2,2 (λ, ) + Q(λ, ) - g j=1 λ -q j , Ľ2,1 (λ, ) = g j=1 (λ -q j ) ∂Q(λ, ) ∂λ + L 2,1 (λ, ) g j=1 (λ -q j ) -L 2,2 (λ, ) Q(λ, ) g j=1 (λ -q j ) - Q(λ, ) 2 g j=1 (λ -q j ) (2-36)
Similarly, the entries of L are related to those of Ľ by

L1,1 (λ, ) = Ľ1,1 (λ, ) - 1 2 t ∞,2r∞-2 λ + g 0 Ľ1,2 (λ, ) L1,2 (λ, ) = Ľ1,2 (λ, ) L2,1 (λ, ) = Ľ2,1 (λ, ) - 1 2 t ∞,2r∞-2 λ + g 0 2 Ľ1,2 (λ, ) + 1 2 t ∞,2r∞-2 λ + g 0 Ľ1,1 (λ, ) -Ľ2,2 (λ, ) + 1 2 t ∞,2r∞-2 L2,2 (λ, ) = Ľ2,2 (λ, ) + 1 2 t ∞,2r∞-2 λ + g 0 Ľ1,2 (λ, ) (2-37)

Wronksians and asymptotics of the wave functions

Combining the gauge transformations G ∞ , G 1 and J, we obtain the following proposition.

Proposition 2.3. The scalar wave functions ψ 1 = Ψ 1,1 = Ψ1,1 = Ψ1,1 and ψ 2 = Ψ 1,2 = Ψ1,2 = Ψ1,2 have the following expansions around ∞.

ψ 1 (λ) λ→∞ = exp - 1 2r∞-2 k=1 t ∞,k k λ k 2 - 1 4 ln λ + O(1) , ψ 2 (λ) λ→∞ = exp - 1 2r∞-2 k=1 (-1) k t ∞,k k λ k 2 - 1 4 ln λ + O(1) . (2-38)
Proof. The proof is done in Appendix A.

For convenience, we shall also define the Wronskians associated to the Lax systems and provide their explicit expressions that follow from the previous proposition: Definition 2.3 (Wronskians). Let us define W (λ, ) = det Ψ(λ, ), W (λ, ) = det Ψ(λ, ) and W (λ, ) = det Ψ(λ, ) the Wronskians associated to the corresponding wave matrices. They are given by

W (λ) = W0 exp 1 λ 0 P1 (s)ds , W (λ) = W (λ) = W0 exp 1 λ 0 P1 (s)ds W (λ) = W 0 g i=1 (λ -q i ) exp 1 λ 0 P1 (s)ds . (2-39)
where W 0 and W0 are unknown constants (in the sense independent of λ) and where we have defined

P (1) ∞,k = -t ∞,2k+2 , ∀ k ∈ 0, r ∞ -2 , (2-40)
and regrouped them into the polynomial P1 :

P1 (λ) = r∞-2 j=0 P (1) ∞,j λ j = - r∞-2 j=0 t ∞,2j+2 λ j . (2-41)
Proof. The proof starts with

W (λ) = (ψ 1 (λ)∂ λ ψ 2 (λ) -ψ 2 (λ)∂ λ ψ 1 (λ)).
From the general construction and Proposition 2.3, W (λ) exp -1 λ 0 P1 (s)ds is a polynomial function of λ of degree r ∞ -3 = g. Moreover, since L 2,2 (λ) = ∂ λ W (λ) W (λ) , we get that the zeroes of W (λ) are simple poles of L 2,2 (λ). Since L 2,2 (λ) may only have poles at infinity or at (q i ) 1≤i≤g , we get that W (λ) exp -1 λ 0 P1 (s)ds = W 0 g i=1 (λ -q i ) for some constant W 0 (i.e. independent of λ). Formulas for W (λ) and W (λ) follow from W (λ) and the gauge transformations using the determinants of G 1 and J.

Explicit expression for the Lax matrix L

In this section we shall provide an explicit expression for the matrix L(λ) in terms of irregular times and Darboux coordinates. Only g coefficients of the matrix shall remain undetermined at this stage. These coefficients will be put in one-to-one correspondence with the upcoming Hamiltonians. In order to write down the Lax matrix in a compact form, we shall introduce the following definition.

Definition 2.4. We define the following quantities:

P (2) ∞,k = 1 4 2r∞-2 j=2k-2r∞+6 (-1) j t ∞,j t ∞,2k-j+4 , ∀ k ∈ r ∞ -2, 2r ∞ -4 P (2) ∞,r∞-3 = 1 4 2r∞-3 j=1 (-1) j t ∞,j t ∞,2r∞-j-2 (2-42)
and regroup them into the polynomial function P2 :

P2 (λ) = 2r∞-4 k=r∞-3 P (2) ∞,k λ k (2-43)
We shall also define

P2 (λ) = P2 (λ) - r∞-4 k=0 H ∞,k λ k (2-44)
where the g = r ∞ -3 coefficients (H ∞,k ) 0≤k≤r∞-4 remain undetermined at this stage.

Using the previous definition, we obtain the following proposition.

Proposition 2.4. The Lax matrix L(λ, ) is given by

L(λ, ) = 0 1 L 2,1 (λ, ) L 2,2 (λ, ) (2-45) with L 2,2 (λ, ) = P1 (λ) + g j=1 λ -q j L 2,1 (λ, ) = -P2 (λ) + r∞-4 k=0 H ∞,k λ k - g j=1 p j λ -q j (2-46)
Coefficients (H ∞,k ) 0≤k≤r∞-4 shall be determined later in Proposition 5.1.

Proof. The proof is based on the fact that the entries of L are rational functions of λ with poles only at ∞ or at apparent singularities (q i ) 1≤i≤g . Using the knowledge of the asymptotics expansion at ∞ provides the result. This is detailed in Appendix B.

Classical spectral curve and connection with topological recursion

Before turning to deformations relatively to the irregular times, let us briefly mention the connection of the present setup with the classical spectral curve and the topological recursion. This section being independent of the others, we stress that readers with no interest in topological recursion or in WKB expansions may skip the content of this section.

Let us first recall how one may obtain the classical spectral curve from a Lax system. When dealing with a Lax system of the form

∂ λ Ψ(λ, ) = L(λ, )Ψ(λ, ), (3-1)
it is standard to define the "classical spectral curve" as lim

→0 det(yI 2 -L(λ, )) = 0. It is important
to note that the classical spectral curve is unaffected by the gauge transformations Ψ(λ, ) → G(λ, )Ψ(λ, ) with G(λ, ) regular in . Indeed, the conjugation of the Lax matrix does not change the characteristic polynomial and the additional term (∂ λ G)G -1 disappears in the limit → 0. In particular, in our setup, it means that one may compute the classical spectral curve using either L, Ľ or L:

lim →0 det(yI 2 -L(λ, )) = lim →0 det(yI 2 -Ľ(λ, )) = lim →0 det(yI 2 -L(λ, )). ( 3-2)
In our case, the general expression of the matrix L(λ, ) implies that the classical spectral curve is y 2 -P1 (λ, = 0)y + P2 (λ, = 0) = 0.

(3-3)

It defines a Riemann surface Σ of genus g = r ∞ -3 whose coefficients are determined by (2-40) and Definition 2.4. Note that only g coefficients remained undetermined at this stage (i.e. (H ∞,k ) 0≤k≤r∞-4 ) that can be mapped with the so-called filling fractions ( i ) 1≤i≤g naturally associated to the Riemann surface. Moreover, the present twisted case corresponds to the case where infinity is a ramification point of the Riemann surface. In other words, the twisted case happens when a pole of the connection is also a ramification point of the underlying classical spectral curve. The asymptotic expansions of the differential form ydx at each pole is in direct relation with the asymptotics of the wave functions (2-38) since we have

y 1 (z) z→∞ = - 1 2 2r∞-2 k=1 t ∞,k x(z) k 2 -1 - 1 4x(z) + O x(z) -3 2 y 2 (z) z→∞ = - 1 2 2r∞-2 k=1 (-1) k t ∞,k x(z) k 2 -1 - 1 4x(z) + O x(z) -3 2 (3-4)
where y 1 (z) and y 2 (z) are the expressions of y(z) in both sheets.

Let us now discuss the connection of the present work with the Chekhov-Eynard-Orantin topological recursion [START_REF] Chekhov | Hermitean matrix model free energy: Feynman graph technique for all genera[END_REF][START_REF] Chekhov | Hermitian matrix model free energy: Feynman graph technique for all genera[END_REF][START_REF] Chekhov | Matrix eigenvalue model: Feynman graph technique for all genera[END_REF][START_REF] Eynard | Invariants of algebraic curves and topological expansion[END_REF][START_REF] Eynard | Topological recursion in random matrices and enumerative geometry[END_REF] as given in [START_REF] Marchal | Hamiltonian representation of isomonodromic deformations of general rational connections on[END_REF]. Recent works [START_REF] Marchal | Isomonodromic deformations of a rational differential system and reconstruction with the topological recursion: the sl 2 case[END_REF][START_REF] Eynard | Quantization of classical spectral curves via topological recursion[END_REF] have shown how to quantize the classical spectral curve using topological recursion. Indeed, applying the topological recursion to the classical spectral curve (3-3) generates Eynard-Orantin differentials (ω h,n ) h≥0,n≥0 that can be regrouped into formal -transseries to define formal wave functions ψ TR 1 , ψ TR 2 that satisfy a quantum curve, i.e. a linear ODE of degree 2 with pole singularities at infinity and apparent singularities at λ = q i and whose → 0 limit recovers the classical spectral curve. The construction presented in [START_REF] Marchal | Isomonodromic deformations of a rational differential system and reconstruction with the topological recursion: the sl 2 case[END_REF][START_REF] Eynard | Quantization of classical spectral curves via topological recursion[END_REF] implies that this ODE is the same as the one defined by the Lax matrix L(λ, ) of the present paper so that we get

Ψ(λ, ) = C ψ TR 1 (λ, ) ψ TR 2 (λ, ) ∂ψ TR 1 (λ, ) ∂ψ TR 2 (λ, ) (3-5)
where C is a constant (independent of λ) matrix. In other words, the topological recursion reconstructs our wave functions ψ 1 and ψ 2 making the classical spectral curve the only necessary object to build the full Lax system. However, the price to pay in this perspective is the mandatory introduction of the formal parameter to define the formal -transseries and then ψ TR 1 , ψ TR

2

. As explained in Section 2.5, this formal parameter can be removed by proper rescaling at the level of the Lax system but it is unclear how the topological recursion wave functions may be defined after this rescaling, since there is no more formal parameter to define the series. This issue is in deep relation with the analytical meaning that might be given to the formal -transseries. In particular, it is presently unclear how to resum analytically the -transseries to obtain non-formal identities and current works are in progress to tackle this problem.

General isomonodromic deformations and auxiliary matrices 4.1 Definition of general isomonodromic deformations

The previous sections provide a natural set of parameters for which we may consider deformations, namely the irregular times (t ∞,k ) 1≤k≤2r∞-2 . In order to study deformations relatively to these parameters we introduce the following definition. Definition 4.1. We define the following general deformation operators.

L α = 2r∞-2 k=1 α ∞,k ∂ t ∞,k (4-1)
where we define the vector α

∈ C 2r∞-2 = C 2g+4 by α = 2r∞-2 k=1 α ∞,k e k . (4-2)
Deformations defined by Definition 4.1 shall be seen as general isomonodromic deformations in F∞,r∞ . Associated to a vector α are general auxiliary Lax matrices Ãα (λ), Ǎα (λ) and A α (λ) defined by

Ãα (λ) = L α [ Ψ(λ)] Ψ-1 (λ) ⇔ L α [ Ψ(λ)] = Ãα (λ) Ψ(λ) Ǎα (λ) = L α [ Ψ(λ)] Ψ-1 (λ) ⇔ L α [ Ψ(λ)] = Ǎα (λ) Ψ(λ) A α (λ) = L α [Ψ(λ)]Ψ -1 (λ) ⇔ L α [Ψ(λ)] = A α (λ)Ψ(λ) (4-3)
In particular, Ãα (λ) and Ǎα (λ) are polynomial functions of λ while A(λ) may also have additional poles at {q 1 , . . . , q g }. Note that (L(λ), A α (λ)), Ľ(λ), Ǎα (λ) and L(λ), Ãα (λ) provide equivalent Lax pairs but expressed in three different gauges. The corresponding compatibility equations are

L α [L] = [A α , L] + ∂ λ A α L α [ Ľ] = [ Ǎα , Ľ] + ∂ λ Ǎα L α [ L] = [ Ãα , L] + ∂ λ Ãα . (4-4)
We shall now use the asymptotic expansions of the wave matrices in order to obtain information on the general form of the auxiliary matrices. Then, we shall use the compatibility equations in order to determine the evolutions of the Darboux coordinates under general isomonodromic deformations and prove that these evolutions are Hamiltonian as performed in a similar way for the non-twisted case in [START_REF] Marchal | Hamiltonian representation of isomonodromic deformations of general rational connections on[END_REF].

General form of the auxiliary matrix A α (λ, )

Using compatibility conditions one may easily obtain two of the entries of A α (λ). Indeed, since L is a companion-like matrix, compatibility equations (4-4) imply that

[A α (λ)] 2,1 = ∂ λ [A α (λ)] 1,1 + [A α (λ)] 1,2 L 2,1 (λ), [A α (λ)] 2,2 = ∂ λ [A α (λ)] 1,2 + [A α (λ)] 1,1 + [A α (λ)] 1,2 L 2,2 (λ), (4-5)
so that only the first line of A α (λ) remains unknown at this stage. The other two entries of the compatibility equation (4-4) leads to

L α [L 2,1 (λ)] = 2 ∂ 2 [A α (λ)] 1,1 ∂λ 2 + 2 L 2,1 (λ) ∂ λ [A α (λ)] 1,2 + [A α (λ)] 1,2 ∂ λ L 2,1 (λ) -L 2,2 (λ) ∂ λ [A α (λ)] 1,1 , L α [L 2,2 (λ)] = 2 ∂ 2 [A α (λ)] 1,2 ∂λ 2 + 2 ∂ λ [A α (λ)] 1,1 + L 2,2 (λ) ∂ λ [A α (λ)] 1,2 + [A α (λ)] 1,2 ∂ λ L 2,2 (λ) (4-6)
that shall be used later to determine the evolution equations for (q i , p i ) 1≤i≤n . Before studying the compatibility equations, let us observe that the asymptotic expansions of the wave matrix Ψ at infinity allows to determine the general form of the auxiliary matrix A α (λ, ). This leads us to the following results.

Proposition 4.1. The asymptotic expansion of entry [A α (λ)] 1,2 at infinity is given by

∀ M ≥ 1 : [A α (λ)] 1,2 λ→∞ = M i=-1 ν (α) ∞,i λ i + O λ -M -1 . (4-7) Moreover, coefficients ν (α) ∞,k -1≤k≤r∞-3
are determined by

M ∞       ν (α) ∞,-1 ν (α) ∞,0 . . . ν (α) ∞,r∞-3       =       2α ∞,2r∞-3 (2r∞-3) 2α ∞,2r∞-5 (2r∞-5) . . . 2α ∞,1 1       (4-8)
where M ∞ is a lower triangular Toeplitz matrix of size (r ∞ -1) × (r ∞ -1) independent of the deformation α:

M ∞ =          t ∞,2r∞-3 0 . . . 0 t ∞,2r∞-5 t ∞,2r∞-3 0 . . . . . . . . . . . . . . . . . . t ∞,3 . . . . . . 0 t ∞,1 t ∞,3 . . . t ∞,2r∞-3          . (4-9)
Proof. The proof is presented in Appendix C.

The previous proposition may be used to determine the general form of the entry

[A α (λ)] 1,2 . Proposition 4.2. Entry [A α (λ)] 1,2 is given by [A α (λ)] 1,2 = ν (α) ∞,-1 λ + ν (α) ∞,0 + g j=1 µ (α) j λ -q j . (4-10) Coefficients µ (α) j 1≤j≤g
are determined by the linear system

V ∞       µ (α) 1 . . . . . . µ (α) g       =       ν (α) ∞,1 ν (α) ∞,2 . . . ν (α) ∞,r∞-3       (4-11)
where

V ∞ is a (r ∞ -3) × g matrix V ∞ =         1 1 . . . . . . 1 q 1 q 2 . . . . . . q g . . . . . . . . . . . . q r∞-4 1 q r∞-4 2 . . . . . . q r∞-4 g         , (4-12)
Proof. We know that [A α (λ)] 1,2 is rational in λ with only simple poles in {q 1 , . . . , q g } and a pole at infinity. Proposition 4.1 provides the asymptotics at infinity so that (4-10) holds. Moreover the expansion at infinity of

g j=1 µ (α) j λ -q j = ∞ k=1 g j=1 µ (α) j q k-1 j λ -k (4-13)
identifies with (4-7) only with [START_REF] Bergère | Determinantal formulae and loop equations[END_REF][START_REF] Bertola | Hamiltonian structure of rational isomonodromic deformation systems[END_REF][START_REF] Bertola | The partition function of the two-matrix model as an isomonodromic τ function[END_REF][START_REF] Biquard | Wild non-abelian hodge theory on curves[END_REF][START_REF] Boalch | Symplectic manifolds and isomonodromic deformations[END_REF][START_REF] Boalch | Simply-laced isomonodromy systems[END_REF][START_REF] Boalch | Twisted local wild mapping class groups: configuration spaces, fission trees and complex braids[END_REF][START_REF] Chekhov | Hermitean matrix model free energy: Feynman graph technique for all genera[END_REF].

Note that we may determine coefficients ν

(α) ∞,k k≥r∞-2
by the fact that

g j=1 µ (α) j i =j (λ -q i ) = [A α (λ)] 1,2 -ν (α) ∞,-1 λ -ν (α) ∞,0 g i=1 (λ -q i ) = ∞ k=1 ν (α) ∞,k λ -k g i=1 (λ -q i ) = ∞ k=1 ν (α) ∞,k λ -k g i=0 (-1) g-i e g-i ({q 1 , . . . , q g })λ i (4-14)
where the l.h.s. is a polynomial in λ and (e k ({q 1 , . . . , q g })) 0≤k≤g are the elementary symmetric polynomials. Thus, for all m ≥ 1:

0 = g+m k=m (-1) g+m-k ν (α) ∞,k e g+m-k ({q 1 , . . . , q g }) = ν (α) ∞,g+m + g+m-1 k=m (-1) g+m-k ν (α) ∞,k e g+m-k ({q 1 , . . . , q g }) (4-15)
so that we obtain the recursive relations

∀ m ≥ 1 : ν (α) ∞,g+m = g+m-1 k=m (-1) g+m-1-k ν (α) ∞,k e g+m-k ({q 1 , . . . , q g }) (4-16)
In particular, we get for m = 1:

ν (α) ∞,r∞-2 = g k=1 (-1) g-k ν (α) ∞,k e g+1-k ({q 1 , . . . , q g }) (4-17)
Let us now perform similar computation for [A α (λ)] 1,1 . We obtain the following proposition.

Proposition 4.3. The entry [A α (λ)] 1,1 is given by [A α (λ)] 1,1 = r∞-1 i=0 c (α) ∞,i λ i + g j=1 ρ (α) j λ -q j . (4-18) with ∀ j ∈ 1, n : ρ (α) j = -µ (α) j p j (4-19) Coefficients c (α) ∞,k 1≤k≤r∞-1 are determined by M ∞          c (α) ∞,r∞-1 . . . c (α) ∞,k . . . c (α) ∞,1          =             α ∞,2r∞-3 2r∞-3 t ∞,2r∞-2 - α ∞,2r∞-2 2r∞-2 t ∞,2r∞-3 . . . r∞-1 m=k α ∞,2k+2r∞-2m-3 2k+2r∞-2m-3 t ∞,2m - α 2k+2r∞-2m-2 2k+2r∞-2m-2 t ∞,2m-1 . . . r∞-1 m=1 α ∞,2r∞-2m-1 2r∞-2m-1 t ∞,2m - α 2r∞-2m 2r∞-2m t ∞,2m-1             (4-20)
with the matrix M ∞ given by (4-9).

Proof. The proof is done in Appendix D.

Note that c (α)

∞,0 is not determined but will play no role in the rest of the paper. In the previous propositions, one may easily observe that M ∞ , ν

(α) ∞,k -1≤k≤r∞-3 and c (α) ∞,k 1≤k≤r∞-1
are independent of the Darboux coordinates and depend only on irregular times and the deformation α. On the contrary, µ

(α) j 1≤j≤g
and V ∞ depend on the Darboux coordinates.

General Hamiltonian evolutions

The previous sections provide the general form of the matrices L(λ, ) and A α (λ) through Propositions 2.4, 4.1, 4.2, 4.3 and equation (4-5). As we shall see below, inserting the previous knowledge into the compatibility equations (4-6) provides the evolutions of the Darboux coordinates.

The first step is to look at order (λ -

q j ) -2 in L α [L 2,2 (λ)].
We obtain, for all j ∈ 1, g :

L α [q j ] = 2µ (α) j p j - 1 2 P1 (q j ) -ν (α) ∞,0 -ν (α) ∞,-1 q j - i =j µ (α) j + µ (α) i q j -q i . (5-1)
The next step is to determine the coefficients (H ∞,j ) 0≤j≤r∞-4 that remain unknown in L 2,1 (λ). To achieve this task, we look at order (λ -q j ) -2 in L α [L 2,1 (λ)] using [START_REF] Bergère | Determinantal formulae and loop equations[END_REF][START_REF] Bertola | Hamiltonian structure of rational isomonodromic deformation systems[END_REF][START_REF] Bertola | The partition function of the two-matrix model as an isomonodromic τ function[END_REF]. We obtain

-pjLα[qj] = -2 µ (α) j   -P2(qj) + r∞-4 k=0 H ∞,k q k j - i =j pi qj -qi   + 2 pj   ν (α) ∞,-1 qj + ν (α) ∞,0 + i =j µ (α) i qj -qi   -µ (α) j pj   P1(qj) + i =j qj -qi   .
(5-2)

Inserting (5-1) provides, for all j ∈ 1, g , r∞-4 k=0 H ∞,k q k j = p 2 j -P1 (q j )p j + P2 (q j ) + i =j p i -p j q j -q i (5-3)
where it is obvious that the r.h.s. is independent of the deformation vector α. The last relation can be rewritten into a matrix form.

Proposition 5.1. We have

(V ∞ ) t    H ∞,0 . . . H ∞,r∞-4    =       p 2 1 -P1 (q 1 )p 1 + P2 (q 1 ) + i =1 p i -p 1 q 1 -q i . . . p 2 g -P1 (q g )p g + P2 (q g ) + i =g p i -pg qg-q i       (5-4)
Finally, in order to obtain the evolution equation for (p j ) 1≤j≤g we look at order (λ

-q j ) -1 of the entry L α [L 2,1 (λ)]. We get, for all j ∈ 1, g , L α [p j ] = i =j (µ (α) i + µ (α) j )(p i -p j ) (q j -q i ) 2 + µ (α) j p j P 1 (q j ) -P 2 (q j ) + r∞-4 k=1 kH ∞,k q k-1 j + ν (α) ∞,-1 p j + r∞-1 k=1 kc (α) ∞,k q k-1 j .
(5-5)

Thus, we have obtained the general evolutions for (p j , q j ) 1≤j≤g through (5-1) and . We may now formulate our first main Theorem showing that these evolutions are Hamiltonian.

(α) (q, p) = r∞-4 k=0 ν (α) ∞,k+1 H ∞,k - g j=1 r∞-1 k=1 c (α) ∞,k q k j -ν (α) ∞,0 g j=1 p j -ν (α) ∞,-1 g j=1 q j p j , (5-6)
the evolutions for j ∈ 1, g ,

L α [q j ] = 2µ (α) j p j - 1 2 P1 (q j ) -ν (α) ∞,0 -ν (α) ∞,-1 q j - i =j µ (α) j + µ (α) i q j -q i , L α [p j ] = i =j (µ (α) i + µ (α) j )(p i -p j ) (q j -q i ) 2 + µ (α) j p j P 1 (q j ) -P 2 (q j ) + r∞-4 k=1 kH ∞,k q k-1 j + ν (α) ∞,-1 p j + r∞-1 k=1 kc (α) ∞,k q k-1 j (5-7)
are Hamiltonian in the sense that

∀ j ∈ 1, g : L α [q j ] = ∂Ham (α) (q, p) ∂p j and L α [p j ] = - ∂Ham (α) (q, p) ∂q j .
(5-8)

Quantities involved in the Hamiltonian evolution are defined by Propositions 4.1, 4.2, 4.3 and 5.1.

Proof. Proof is done in Appendix E.

Remark 5.1. Note that there is an alternative expression for the Hamiltonian (5-6):

Ham (α) (q, p) = - 2 (i,j)∈ 1,g 2 i =j (µ (α) i + µ (α) j )(p i -p j ) q i -q j - g j=1 (ν (α) ∞,0 p j + ν (α) ∞,-1 q j p j ) + g j=1 µ (α) j p 2 j -P1 (q j )p j + P2 (q j ) - g j=1 r∞-1 k=1 c (α) ∞,k q k j (5-9)
Theorem 5.1 shows that the Hamiltonian expression for a general isomonodromic deformation may be split into several contributions

• A linear combination of the (H ∞,k ) 0≤k≤r∞-4 whose coefficients are given by ν (α) ∞,k+1 1≤k≤r∞-3 .
Note that the coefficients (H ∞,k ) 0≤k≤r∞-4 do not depend on the isomonodromic deformations and correspond to the unknown coefficients of L 2,1 (λ).

• A linear combination of

g j=1 q k j 1≤k≤r∞-1
whose coefficients are given by c

(α) ∞,k 1≤k≤r∞-1
.

As we will see in the next sections, these terms will vanish for a suitable choice of non-trivial isomonodromic deformations.

• Two additional terms g j=1 p j and g j=1 q j p j that are respectively proportional to -ν

(α) ∞,0 and -ν (α)
∞,-1 . These terms shall be removed after a suitable symplectic rescaling of (q i , p i ) 1≤i≤g .

Expression of the Hamiltonian and Lax matrices in terms of symmetric Darboux coordinates

In this section, we show that we may use the symmetric polynomials (e i ({q 1 , . . . , qg })) 1≤i≤g to obtain polynomial Hamiltonians and polynomial explicit formulas for the matrices L and Ã(ατ ) .

Notations and identities regarding symmetric polynomials

In the rest of the paper we need to introduce elementary symmetric polynomials and other basis of symmetric polynomials.

Definition 6.1 (Basis of symmetric polynomials). We shall introduce the following basis of symmetric polynomials:

• Elementary symmetric polynomials are denoted by (e i ({x 1 , . . . , x n })) i≥0 with the convention that e 0 ({x 1 , . . . , x n }) = 1 and e k ({x 1 , . . . , x n }) = 0 if k > n. By definition we have:

e k ({x 1 , . . . , x n }) = 1≤i 1 <•••<i k ≤n x i 1 . . . x i k , ∀ k ∈ 1, n (6-1)
• Complete homogeneous symmetric polynomial are denoted by (h i ({x 1 , . . . , x n })) i≥0 with the convention that h 0 ({x 1 , . . . , x n }) = 1. By definition we have:

h k ({x 1 , . . . , x n }) = 1≤i 1 ≤•••≤i k ≤n x i 1 . . . x i k , ∀ k ∈ 1, n (6-2)
• k th symmetric power sum polynomials are denoted by (S k ({x 1 , . . . , x n })) k≥0 . By definition, we have:

S 0 ({x 1 , . . . , x n }) = n S k ({x 1 , . . . , x n }) = n j=1 x k j , ∀ k ≥ 1 (6-3) (e k ({x 1 , . . . , x n })) 0≤k≤n , (h k ({x 1 , . . . , x n })) 0≤k≤n and (S k ({x 1 , . . . , x n }))
0≤k≤n are some basis of symmetric polynomials in the variables {x 1 , . . . , x n }. We also have the relations

n j=1 (λ -x j ) = n k=0 (-1) n-k e n-k ({x 1 , . . . , x n })λ k = n k=0 (-1) k e k ({x 1 , . . . , x n })λ n-k 1 n j=1 (λ -x j ) = ∞ k=0 h k ({x 1 , . . . , x n })λ -n-k (6-4)
The relation between the various sets are given by

h 0 ({x 1 , . . . , x n }) = e 0 ({x 1 , . . . , x n }) h k ({x 1 , . . . , x n }) = k j=1 (-1) j b 1 ,...,b j ∈ 1,k j b 1 +•••+b j =k j m=1 (-1) bm e bm ({x 1 , . . . , x n }) , ∀ k ∈ 1, n (6-5)
and ∀ m ≥ 1:

S m ({x 1 , . . . , x n }) = (-1) m m m k=1 1 k Bm,k (-e 1 ({x 1 , . . . , x n }), . . . , -e m-k+1 ({x 1 , . . . , x n })) = (-1) m m b 1 +2b 2 +•••+mbm=m b 1 ≥0,...,bm≥0 (-1) b 1 +•••+bm (b 1 + • • • + b m ) b 1 + • • • + b m b 1 , . . . , b m m i=1 e i ({x 1 , . . . , x n }) b i (6-6)
where Bm,k m≥k≥0 are the ordinary Bell polynomials. Finally, we also have the identities

(n -k)e k ({x 1 , . . . , x n }) = k i=0 (-1) i e k-i ({x 1 , . . . , x n })S i ({x 1 , . . . , x n }) , ∀ k ∈ 0, n S k ({x 1 , . . . , x n }) = k-1 i=k-n (-1) k-1+i e k-i ({x 1 , . . . , x n })S i ({x 1 , . . . , x n }) , ∀ k ≥ n (6-7)
Elementary symmetric polynomials satisfy some useful relations: Lemma 6.1. For any (i, m) ∈ 1, g 2 :

∂e i ({x 1 , . . . , x g }) ∂x m = i-1 j=0
(-1) j e i-1-j ({x 1 , . . . , x g })x j m (6-8) Proposition 6.1. For any i ∈ 1, g , we have

g k=1 ∂e i ({x 1 , . . . , x g }) ∂x k g j =k λ -x j x k -x j = i-1 j=0 (-1) j e i-j-1 ({x 1 , . . . , x g })λ j (6-9)
These relations allow to express Q(λ).

Corollary 6.1. We have

Q(λ) = g-1 j=0 (-1) j-1   g i=j+1 P i Q i-j-1   λ j (6-10)
Moreover, the elementary symmetric polynomials satisfy:

∀ (i, j) ∈ 1, n 2 : e n-i ({x 1 , . . . , x n } \ {x j }) = n m=i (-1) m-i e n-m ({x 1 , . . . , x n })x m-i j ∀ j ∈ 1, n : 0 = g m=0 (-1) n-m e n-m ({x 1 , . . . , x n })x m j (6-11)
so that we obtain Lemma 6.2. For any i ∈ 1, n and any M ≥ 0 we have:

x M i = n j=1 n m=Max(j,j+n-1-M ) (-1) n-m e n-m ({x 1 , . . . , x n })h M +m-j-n+1 ({x 1 , . . . , x n })x j-1 i (6-12)
In particular, we may invert the Vandermonde matrix with the following Proposition.

Proposition 6.2. For any i ∈ 0, n and M ≥ 0 we have:

n j=1 (-1) n-i e n-i ({x 1 , . . . , x n } \ {x j }) m =j (x j -x m ) x M j = n m=Max(i,i+n-1-M ) (-1) n-m e n-m ({x 1 , . . . , x n })h M +m-i-n+1 ({x 1 , . . . , x n }) (6-13)
In particular, for M ≤ n -1 we get:

∀ M ∈ 0, n -1 : n j=1 (-1) n-i e n-i ({x 1 , . . . , x n } \ {x j }) m =j (x j -x m ) x M j = δ i,M +1 (6-14)
Proof. For completeness, the proofs are presented in Appendix F.

Symmetric Darboux coordinates

Let us first recall the well-known result from symplectic geometry.

Lemma 6.3. If we define new coordinates ( Q1 , . . . , Qg , P1 , . . . , Pg ) from old symplectic coordinates (Q 1 , . . . , Q g , P 1 , . . . , P g ) by

Qi = f i (α -1 Q 1 + β, . . . , α -1 Q g + β) αP i + h(α -1 Q i + β) = α g k=1 Pk ∂f k (α -1 Q 1 + β, . . . , α -1 Q g + β) ∂Q i , ∀ i ∈ 1, g (6-15)
with α ∈ C\{0} and β ∈ C two given constants, h any function of class C 1 and (f m (x 1 , . . . , x g )) 1≤m≤g any functions of class C 2 then the change of coordinates is symplectic.

Proof. For completeness, the proof is done in Appendix G

We may now apply the lemma with the elementary symmetric polynomials (e i (q 1 , . . . , qg )) 1≤i≤g which is a basis of the symmetric polynomials in (q 1 , . . . , qg ). Definition 6.2 (Symmetric Darboux coordinates). We define (Q 1 , . . . , Q g , P 1 , . . . , P g ) using the elementary symmetric polynomials:

Q i = e i (q 1 , . . . , q g ) = e i (T -1 2 q1 -T 1 , . . . , T -1 g qg -T 1 ) p i = g k=1 P k ∂e k (q 1 , . . . , q g ) ∂q i = T 2 pi + 1 2 P1 (T -1 2 qi -T 1 ) , ∀ i ∈ 1, g (6-16)
We shall denote (Q 1 , . . . , Q g , P 1 , . . . , P g ), the symmetric Darboux coordinates.

It is obvious from Lemma 6.3 that the change of coordinates is symplectic. Indeed, for the old variables (q 1 , . . . , q g , p 1 , . . . , p g ) this is nothing but an application of Lemma 6.3 with α = 1, β = 0 and h(x) = 0 while for the other old variables (q 1 , . . . , qg , q1 , . . . , qg ), this corresponds to an application of Lemma 6.3 with α = T 2 , β = -T 1 and h(x) = 1 2 P1 (x).

Polynomial expression of the Hamiltonian in the symmetric Darboux coordinates

Since the change of coordinates (q i , pi ) 1≤i≤g → (Q i , P i ) 1≤i≤g is symplectic, we may compute the Hamiltonian Ham(Q, . . . , P) by just replacing the coordinates (q i , p i ) 1≤i≤g in terms of (Q i , P i ) 1≤i≤g in Theorem 5.1.

Theorem 6.1 (Expression of the general Hamiltonian in terms of symmetric Darboux coordinates). We have:

Ham (α) (Q, P) = - g j=1 (ν (α) ∞,0 p j + ν (α) ∞,-1 q j p j ) - g j=1 r∞-1 k=1 c (α) ∞,k q k j + g i=1 ν (α) ∞,i H ∞,i+1 = -ν (α) ∞,0 g-1 k=0 (g -k)Q k P k+1 -ν (α) ∞,-1 g k=1 kQ k P k - r∞-1 k=1 c (α) ∞,k S k ({q 1 , . . . , q g }) - g i=1 ν (α) ∞,i g k=i+1 (-1) i (g -i)P k Q k-1-i + k-1 m=i+1 (-1) m P k Q k-1-m S m-i ({q 1 , . . . , q g }) + g i=1 ν (α) ∞,i g k1=1 g k2=1 P k1 P k2 (-1) i-1 Min(k1-1,i-1) r1=Max(0,i-k2) Q k1-1-r1 Q k2-i+r1 + 0≤r1≤k1-1 0≤r2≤k2-1 r1+r2≥g (-1) r1+r2 Q k1-1-r1 Q k2-1-r2 g m=i (-1) g-m Q g-m h r1+r2+m-i-g+1 ({q 1 , . . . , q g }) + g i=1 ν (α) ∞,i g k=1 Min(k-1,i-1) r=0 (-1) r t ∞,2i-2r P k Q k-1-r + k-1 r=0 g+1 s=g-r g m=i (-1) g+r-m t ∞,2s+2 P k Q k-1-r Q g-m h r+s+m-i-g+1 ({q 1 , . . . , q g }) + g i=1 ν (α) ∞,i 2r∞-4 r=g g m=i (-1) g-m P (2) ∞,r Q g-m h r+m-i-g+1 ({q 1 , . . . , q g }) (6-17)
where P

∞,k g≤k≤2g+4 are given by Proposition 2.4, ν

(α) ∞,i -1≤i≤g
are given by (4-8), and coefficients (S k ({q 1 , . . . , q g })) k≥0 , (h k ({q 1 , . . . , q g })) k≥0 are given by Definition 6.1.

The main advantage of the explicit expression [START_REF] Bertola | The partition function of the two-matrix model as an isomonodromic τ function[END_REF][START_REF] Biquard | Wild non-abelian hodge theory on curves[END_REF][START_REF] Boalch | Symplectic manifolds and isomonodromic deformations[END_REF][START_REF] Boalch | Simply-laced isomonodromy systems[END_REF][START_REF] Boalch | Twisted local wild mapping class groups: configuration spaces, fission trees and complex braids[END_REF][START_REF] Chekhov | Hermitean matrix model free energy: Feynman graph technique for all genera[END_REF][START_REF] Chekhov | Hermitian matrix model free energy: Feynman graph technique for all genera[END_REF][START_REF] Chekhov | Matrix eigenvalue model: Feynman graph technique for all genera[END_REF][START_REF] Dumitrescu | Quantization of spectral curves for meromorphic Higgs bundles through topological recursion[END_REF][START_REF] Eynard | Quantization of classical spectral curves via topological recursion[END_REF][START_REF] Eynard | Invariants of algebraic curves and topological expansion[END_REF][START_REF] Eynard | Topological recursion in random matrices and enumerative geometry[END_REF] is that it immediately shows that the general Hamiltonian is polynomial in (Q i , P i ) 1≤i≤g , i.e. it has the same kind of singularities as the initial connection. In particular, it is quadratic in (P i ) 1≤i≤g . Note also that the explicit dependence of the Hamiltonian in the irregular times is contained only in ν

(α) ∞,i -1≤i≤g and P (2) ∞,k g≤k≤2g+1 .
Proof. The proof is done in Appendix H.

Expressing the Lax matrices with the symmetric Darboux coordinates

Symmetric Darboux coordinates (Q 1 , . . . , Q g , P 1 , . . . , P g ) are well-suited for the matrix L given by (2-37) as the following proposition shows Proposition 6.3. Entries of the matrix L(λ) are given by

L1,1(λ) = - g-1 j=0 (-1) j-1 g i=j+1 PiQi-j-1 λ j - 1 2 t∞,2r ∞-2 λ + g0 g j=0 (-1) g-j Qg-jλ j L1,2(λ) = g m=0 (-1) g-m Qg-mλ m L2,2(λ) = g-1 j=0 (-1) j-1 g i=j+1 PiQi-j-1 λ j + 1 2 t∞,2r ∞ -2λ + g0 g j=0 (-1) g-j Qg-jλ j - r∞-2 k=0 t ∞,2k+2 λ k L2,1(λ) = - r∞-2 i=0 2r∞-4 j=g+i P (2) ∞,j hj-g-i({q1, . . . , qg}) λ i + g-1 i=0 g-1 j=i g+1 s=g+i-j (-1) j-1 t∞,2s+2 g r=j+1 PrQr-j-1 hs+j-i-g λ i - g-2 i=0 g-1 j 1 =i+1 g-1 j 2 =g+i-j 1 (-1) j 1 +j 2 g i 1 =j 1 +1 Pi 1 Qi 1 -j 1 -1 g i 2 =j 2 +1 Pi 2 Qi 2 -j 2 -1 hj 1 +j 2 -g-i({q1, . . . , qg}) λ i - 1 2 t∞,2r ∞-2 λ + g0 2 g m=0 (-1) g-m Qg-mλ m + 1 2 t∞,2r ∞ -2λ + g0 g+1 s=0 t∞,2s+2λ s - 1 2 t∞,2r ∞-2 λ + g0 g-1 j=0 (-1) j-1 g i=j+1 PiQi-j-1 λ j (6-18)
where

g 0 = 1 2 t ∞,2r∞-4 + 1 2 t ∞,2r∞-2 Q 1 . Proof.
Only the expression of Q(λ) is non-trivial and is given by Corollary 6.1.

We remind the reader that P

∞,j g≤j≤2r∞-4 are given by Definition 2.4 while the complete homogenous symmetric polynomials (h j ({q 1 , . . . , q g })) 0≤j≤g may be expressed in terms of (Q k ) 1≤k≤g using h 0 ({q 1 , . . . , q g }) = 1

h k ({q 1 , . . . , q n }) = k j=1 (-1) j b 1 ,...,b j ∈ 1,k j b 1 +•••+b j =k j m=1 (-1) bm Q bm , ∀ k ∈ 1, g (6-19)
In particular Proposition 6.3 implies that the entries of the Lax matrix L are also polynomial in the symmetric Darboux coordinates and at most quadratic in (P i ) 1≤i≤g .

We may also obtain the entries of Ãα in terms of the symmetric Darboux coordinates.

Proposition 6.4. Entries of the matrix Ãα (λ) are given by

[ Ãα (λ)] 1,1 = r∞-1 i=0 c (α) ∞,i λ i - g i=0 g-1-i m=Max(-1,-i) (-1) i+m-1 ν (α) ∞,m g r=i+m+1 P r Q r-i-m-1 λ i - 1 2 t ∞,2r∞-2 λ + g 0 g+1 i=0 g-i m=Max(-1,-i) (-1) g-i-m Q g-i-m ν (α) ∞,m λ i [ Ãα (λ)] 1,2 = g+1 j=0   g-j m=Max(-1,-j) (-1) g-j-m ν (α) ∞,m Q g-j-m   λ j [ Ãα (λ)] 2,2 = -[ Ã(α) (λ)] 1,1 - r∞-1 s=1 1 s α ∞,2s λ s + 2c (α) ∞,0 + (g + 1)ν (α) ∞,-1 - r∞-2 j=0 t ∞,2j+2 ν (α) ∞,j [ Ãα (λ)] 2,1 = - 2 ν (α) ∞,-1 t ∞,2r∞-2 λ 2 -ν (α) ∞,-1 g 0 λ -ν (α) ∞,-1 (-1) g P g - 2 gt ∞,2r∞-2 ν (α) ∞,-1 λ -gg 0 ν (α) ∞,-1 + 2 t ∞,2r∞-2 ν (α) ∞,-1 Q 1 - 2 gt ∞,2r∞-2 ν (α) ∞,0 + (r ∞ -1)c (α) ∞,r∞-1 λ + (r ∞ -1)c (α) ∞,r∞-1 Q 1 + (r ∞ -2)c (α) ∞,r∞-2 + ν (α) ∞,-1 H ∞,r∞-4 + 2 α ∞,2r∞-2 + 2 α ∞,2r∞-4 + 2 α ∞,2r∞-2 Q 1 - g i=0 g-1 j=Max(i-1,0) g+1 s=g+i-j-1 g r=j+1 s+j-g-i m=-1 (-1) j t ∞,2s+2 ν (α) ∞,m h s+j-m-g-i P r Q r-j-1 λ i - r∞ i=0 2r∞-4 j=Max(g,g+i-1) j-g-i m=-1 ν (α) ∞,m h j-g-m-i P (2) ∞,j λ i - g i=0 g-1 j1=0 g-1 j2=0 j1+j2-g-i m=-1 (-1) j1+j2 ν (α) ∞,m h j1+j2-g-m-i g r1=j1+1 g r2=j2+1 P r1 P r2 Q r1-j1-1 Q r2-j2-1 λ i + 1 2 t ∞,2r∞-2 λ + g 0 r∞-1 i=0 r∞-2 s=Max(i-1,0) t ∞,2s+2 ν (α) ∞,s-i λ i -(t ∞,2r∞-2 λ + 2g 0 ) g i=0 g-1 j=Max(i-1,0) g r=j+1 (-1) j-1 ν (α) ∞,j-i P r Q r-j-1 λ i - 1 2 t ∞,2r∞-2 λ + g 0 2 g+1 i=0 g j=Max(i-1,0) (-1) g-j Q g-j ν (α) ∞,j-i λ i (6-20) with g 0 = 1 2 t ∞,2r∞-4 + 1 2 t ∞,2r∞-2 Q 1 . Coefficients ν (α) ∞,m -1≤m≤g
are given by (4-8). Coeffi-

cients c (α) ∞,i 1≤i≤r∞-1
given by (4-20) and ν

(α) ∞,r∞-2 = g k=1 (-1) g-k ν (α) ∞,k Q g+1-k from (4-17).
Proof. The proof is rather long and done in Appendix I.

Decomposition and reduction of the space of isomonodromic deformations

The second goal of this paper is now to provide a decomposition of the space of isomonodromic deformations (of dimension 2g + 4) into a subspace of trivial deformations associated to trivial times (i.e. for which rescaled Darboux coordinates are independent) and a subspace of non-trivial deformations of dimension g associated to non-trivial times while providing the corresponding Hamiltonian evolutions.

Subspaces of trivial and non-trivial deformations

In the previous section we considered general isomonodromic deformations relatively to all irregular times by considering L α characterized by a general vector α ∈ C 2g+4 . However, as we will see below, there exists a subspace of deformations of dimension g + 4 for which the evolutions of the Darboux coordinates are trivial, thus leaving only a non-trivial subspace of deformations of dimension g. These non-trivial deformations shall later be mapped to g isomonodromic times whose expressions will be explicit in terms of the initial irregular times. Trivial deformations shall correspond to the fact that only odd irregular times (t ∞,2k-1 ) 1≤k≥r∞-1 are relevant whereas even irregular times (t ∞,2k ) 1≤k≥r∞-1 do not appear in the Hamiltonians. In other words, considering meromorphic connections in gl 2 (C) or in sl 2 (C) is essentially the same at the level of the Hamiltonian systems. This remark shall provide a subspace of trivial deformations of dimension g + 2. Finally, the remaining 2 trivial deformations correspond to the remaining two degrees of freedom in the action of the Möbius transformations. As we will see below, this choice encodes the necessity of a symplectic rescaling (translation and dilatation) of the Darboux coordinates. These two degrees of freedom shall be used to fix the first two leading non-trivial coefficients at infinity: t ∞,2r∞-3 (conventionally set to 2) and t ∞,2r∞-5 (conventionally set to 0).

Let us first recall that the space of isomonodromic deformations, denoted T , is given by:

L α = 2r∞-2 k=1 α ∞,k ∂ t ∞,k (7-1)
We make the identification with C 2r∞-2 by identifying an isomonodromic deformation L α with its vector α ∈ C 2r∞-2 :

L α = 2r∞-2 k=1 α ∞,k ∂ t ∞,k ⇔ α = 2r∞-2 k=1 α ∞,k e k (7-2)
where we shall denote (e k ) 1≤k≤2r∞-2 the canonical basis of C 2r∞-2 .

Definition 7.1. We define the following vectors of C 2r∞-2 and their corresponding deformations.

w k = e 2k , ∀ k ∈ 1, r ∞ -1 u k = 1 2 r∞-k-2 m=1 (2m -1)t ∞,2m+1+2k e 2m-1 + 1 2 r∞-k-2 s=1 2s t ∞,2s+2k+2 e 2s = 1 2 2r∞-2k-4 r=1 r t ∞,r+2k+2 e r , ∀ k ∈ -1, r ∞ -3 (7-3)
and we shall denote:

U trivial = Span {w 1 , . . . , w r∞-1 , u -1 , u 0 } U iso = Span {u 1 , . . . , u r∞-3 } (7-4)
Note in particular that U iso is of dimension g = r ∞ -3 and that (w 1 , . . . , w r∞-1 , u -1 , . . . , u r∞-3 ) is a basis of C 2r∞-2 . The choice of basis is such that the following proposition holds. Proposition 7.1. We have for all k ∈ 1, r ∞ -1 :

ν (w k ) j = 0 , ∀ j ∈ -1, r ∞ -3 µ (w k ) j = 0 , ∀ j ∈ 1, g c (w k ) ∞,j = - 1 2k δ j,k , ∀ j ∈ 1, r ∞ -1 (7-5)
and for all k ∈ -1, r ∞ -3 :

ν (u k ) j = δ j,k , ∀ j ∈ -1, r ∞ -3 c (u k ) ∞,j = 0 , ∀ j ∈ 1, r ∞ -1 (7-6)
Proof. The proof is presented in Appendix J.

Note that Proposition 7.1 and (4-11) imply that

µ (u -1 ) j = 0 = µ (u 0 ) j , ∀ j ∈ 1, g (7-7)
so that for all j ∈ 1, g : Proposition 7.1 combined with Theorem 5.1 provides the following theorem.

Theorem 7.1. For any j ∈ 1, g , we have:

L w k [q j ] = 0 , ∀ k ∈ 1, r ∞ -1 L w k [p j ] = - 2 q k-1 j , ∀ k ∈ 1, r ∞ -1 L u -1 [q j ] = -q j L u -1 [p j ] = p j L u 0 [q j ] = - L u 0 [p j ] = 0 (7-8)
Proof. The proof is done in Appendix K.

Note that L u -1 and L u 0 do not act trivially on (q j , p j ) 1≤j≤g . As we will see below, one needs to rescale the Darboux coordinates in order to have a trivial action. The purpose of the next section is to define trivial and isomonodromic times that are dual to the previous deformations. However, it is not possible to define some times (τ 1 , . . . , τ r∞-3 ) such that ∂ τ k = L u k since the system becomes non compatible for g ≥ 4.

Definition of trivial times and isomonodromic times

The split in the tangent space between trivial and non-trivial subspaces may be translated at the level of coordinates. This corresponds to choosing g non-trivial times and g + 4 trivial times for which the evolutions of the shifted Darboux coordinates are trivial. In particular, one may then choose the values of these trivial times to any arbitrary values without changing the Hamiltonian evolutions. However, it s important to notice that the choice of trivial times and isomonodromic times is not unique since, for example, one may use any arbitrary combination of isomonodromic times to provide a new one. We propose the following set of trivial and non-trivial times that are particularly convenient in our context. Definition 7.2 (Trivial and non-trivial deformation times). Let us define the following "trivial times":

T ∞,k = t ∞,2k , ∀ k ∈ 1, r ∞ -1 T 2 = 1 2 t ∞,2r∞-3 2 2r∞-3 T 1 = t ∞,2r∞-5 2r ∞ -5 1 2 t ∞,2r∞-3 -2r∞-5 2r∞-3 (7-9)
We also define the g = r ∞ -3 "isomonodromic" times (τ k ) 1≤k≤g ,for all k ∈ 1, g , by:

τ k = k-1 i=0 (-1) i i s=1 (2r ∞ -2k + 2s -7) 1 2 t ∞,2r∞-5 i 1 2 t ∞,2r∞-3 - (2r∞-3)i+2r∞-5-2k 2r∞ -3 1 2 t ∞,2r∞-5-2k+2i i!(2r ∞ -5) i + (-1) k k s=1 (2r ∞ -2k + 2s -7) 1 2 t ∞,2r∞-5 k+1 1 2 t ∞,2r∞-3 - (k+1)(2r∞-5) 2r∞-3 (k + 1)(k -1)!(2r ∞ -5) k (7-10)
We shall denote T trivial the set of trivial times and T iso the set of isomonodromic times:

N trivial = {T ∞,1 , . . . , T ∞,2r∞-2 , T 1 , T 2 } , N iso = {τ 1 , . . . , τ g } (7-11)
The previous set of trivial and non-trivial times is trivially in one-to-one correspondence with the irregular times (t ∞,k ) 1≤k≤2r∞-3 . Moreover, the inverse change of coordinates is given by the following proposition. Proposition 7.2. One may recover the irregular times (t ∞,k ) 1≤k≤2r∞-3 from T trivial ∪ T iso with the following formulas:

t ∞,2r∞-3 = 2T 2r∞-3 2 2 t ∞,2r∞-5 = (2r ∞ -5)T 1 T 2r∞-5 2 2 t ∞,2i = T ∞,i , ∀ i ∈ 1, r ∞ -1 (7-12)
and for all k ∈ 1, r ∞ -3 :

t ∞,2k-1 = 2T 2k-1 2 2      r∞-k-2 p=1 r∞-k-2 m=p+1 (2r ∞ -2m -5) 2 r∞-k-p-2 (r ∞ -k -p -2)! T r∞-k-p-2 1 τ p + T r∞-1-k 1 r∞-k-2 m=0 (2r ∞ -2m -5) 2 r∞-1-k (r ∞ -1 -k)!      (7-13)
Proof. The proof is computational and is proposed in Appendix L.

The last proposition allows to obtain immediately the expression of derivatives relatively to trivial and non-trivial times using the chain rule. Proposition 7.3. For all k ∈ 1, r ∞ -3 , we have:

∂ τ k = 2 r∞-2-k i=1 r∞-i-2 m=k+1 (2r ∞ -2m -5) 2 r∞-i-k-2 (r ∞ -i -k -2)! T r∞-i-k-2 1 T 2i-1 2 2 ∂ t ∞,2i-1 (7-14)
and

∂T ∞,i = ∂t ∞,2i , ∀ i ∈ 1, r∞ -1 ∂T 1 = (2r∞ -5)T 2r∞-5 2 2 ∂t ∞,2r∞ -5 +2 r∞-3 k=1 T 2k-1 2 2      r∞-k-3 p=1 r∞-k-2 m=p+1 (2r∞ -2m -5) 2 r∞-k-p-2 (r∞ -k -p -3)! T r∞-k-p-3 1 τp + T r∞-k-2 1 r∞-k-1 m=0 (2r∞ -2m -5) 2 r∞-1-k (r∞ -k -2)!      ∂t ∞,2k-1 ∂T 2 = (2r∞ -3)T 2r∞-5 2 2 ∂t ∞,2r∞ -3 + (2r∞ -5) 2 2 T1 T 2r∞ -7 2
2
∂t ∞,2r∞-5

+ r∞-3 k=1 (2k -1)T 2k-3 2 2      r∞-k-2 p=1 T r∞-k-p-2 1 τp r∞-k-2 m=p+1 (2r∞ -2m -5) 2 r∞-k-p-2 (r∞ -k -p -2)! + T r∞-1-k 1 r∞-k-2 m=0 (2r∞ -2m -5) 2 r∞-1-k (r∞ -1 -k)!      ∂t ∞,2k-1 (7-15)
For clarity, we shall denote α τ k the corresponding vector in the tangent space associated to ∂ τ k for any k ∈ 1, r ∞ -3 . Its entries are given by

α τ k ∞,2i = 0 , ∀ i ∈ 1, r ∞ -1 α τ k ∞,2r∞-3 = 0 α τ k ∞,2r∞-5 = 0 α τ k ∞,2i-1 = δ 1≤i≤r∞-k-2 2 r∞-i-2 m=k+1 (2r ∞ -2m -5) 2 r∞-i-k-2 (r ∞ -i -k -2)! T r∞-i-k-2 1 T 2i-1 2 2 , ∀ i ∈ 1, r ∞ -3 (7-16)
Remark 7.1. Note that inserting (7-16) into (4-8) implies that

∀ k ∈ 1, g , ∀ i ∈ -1, k -1 : ν (α τ k ) ∞,i = 0 (7-17)
In particular

∀ k ∈ 1, g : ν (α τ k ) ∞,-1 = ν (α τ k )
∞,0 = 0 (7-18)

Properties of trivial and isomonodromic times

Trivial and non-trivial times are chosen so that they satisfy the following properties.

Proposition 7.4. For all k ∈ 1, r ∞ -1 :

L w k [T ∞,j ] = δ j,k , L u -1 [T ∞,j ] = jt ∞,2j , L u 0 [T ∞,j ] = jt ∞,2j+2 , ∀ j ∈ 1, r ∞ -1 L w k [T 2 ] = 0 , L u -1 [T 2 ] = T 2 , L u 0 [T 2 ] = 0 L w k [T 1 ] = 0 , L u -1 [T 1 ] = 0 , L u 0 [T 1 ] = T 2 (7-19)
Proof. Results on (T ∞,j ) 1≤j≤r∞-1 follow by straightforward computations using the fact that

L w k = ∂ t ∞,2k , ∀ k ∈ 1, r ∞ -1 L u -1 = 2 2r∞-2 r=1 rt ∞,r ∂ t∞,r L u 0 = 2 2r∞-4 r=1 rt ∞,r+2 ∂ t∞,r (7-20)
Results on T 2 are also straightforward using the fact that T 2 only depends on t ∞,2r∞-3 . Finally since T 1 depends only on t ∞,2r∞-3 and t ∞,2r∞-5 , we get that

L u -1 [T 1 ] = 2 (2r ∞ -3)t ∞,2r∞-3 ∂ t ∞,2r∞-3 [T 1 ] + 2 (2r ∞ -5)t ∞,2r∞-5 ∂ t ∞,2r∞-5 [T 1 ] = 2 (2r ∞ -3)t ∞,2r∞-3 - t ∞,2r∞-5 2r ∞ -5 (2r ∞ -5) (2r ∞ -3) 1 2 -2r∞-5 2r∞-3 (t ∞,2r∞-3 ) -2r∞-5 2r∞-3 -1 + 2 (2r ∞ -5)t ∞,2r∞-5 1 2r ∞ -5 1 2 t ∞,2r∞-3 -2r∞-5 2r∞-3 = 0 (7-21)
and

L u 0 [T 1 ] = 2 (2r ∞ -5)t ∞,2r∞-3 ∂ t ∞,2r∞-5 [T 1 ] = 2 (2r ∞ -5)t ∞,2r∞-3 1 2r ∞ -5 1 2 t ∞,2r∞-3 -2r∞-5 2r∞-3 = 1 2 t ∞,2r∞-3 2 2r∞-3 = T 2 (7-22)

Shifted Darboux coordinates

Theorem 7.1 indicates that deformations L u -1 and L u 0 do not act trivially on the Darboux coordinates (q j , p j ) 1≤j≤g . However, since the action is very simple, we may easily perform a symplectic transformation on the Darboux coordinates to obtain "shifted Darboux coordinates" for which the action of L u -1 and L u 0 becomes trivial.

Definition 7.3. The shifted Darboux coordinates (q j , pj ) 1≤j≤g are defined by

qj = T 2 q j + T 1 pj = T -1 2 p j - 1 2 P1 (q j ) , ∀ j ∈ 1, g (7-23)
Using Theorem 7.1 and Proposition 7.4, we get that the shifted Darboux coordinates satisfy the following proposition. Proposition 7.5. For all j ∈ 1, g :

L w k [q j ] = L w k [p j ] = 0 , ∀ k ∈ 1, r ∞ -1 L u -1 [q j ] = L u -1 [p j ] = 0 L u 0 [q j ] = L u 0 [p j ] = 0 (7-24)
In other words, for any α 0 ∈ U trivial : L α 0 [q j ] = L α 0 [p j ] = 0, hence the terminology "trivial deformations" and "trivial subspace".

Proof. The proof directly follows from Theorem 7.1 and Proposition 7.4 but we detail it in Appendix M.

Note also that the change of coordinates (q j , p j ) 1≤j≤g ↔ (q j , pj ) 1≤j≤g is symplectic in the sense that

g j=1 dq j ∧ dp j = g j=1 dq j ∧ dp j .
We finally get to our second main theorem.

Theorem 7.2. [Independence of the shifted Darboux coordinates relatively to trivial times] The shifted Darboux coordinates (q j , pj ) 1≤j≤g are independent of the trivial times (T ∞,1 , . . . , T ∞,r∞-1 , T 1 , T 2 ). They are only functions of isomonodromic times (τ k ) 1≤k≤g . Moreover, any function

f (t ∞,1 , t ∞,2 , . . . , t ∞,2r∞-3 ) that is solution of ∀ k ∈ 1, r ∞ -1 : L w k [f ] = 0 and L u -1 [f ] = 0 and L u 0 [f ] = 0 (7-25)
is an arbitrary function u of the isomonodromic times:

f (t ∞,1 , t ∞,2 , . . . , t ∞,2r∞-3 ) = u(τ 1 , . . . , τ g ).
Proof. The proof is presented in Appendix N.

Finally let us mention the following observation.

Proposition 7.6. For any isomonodromic deformations (τ j ) 1≤j≤g , associated to vectors α τ j , the trace of the corresponding matrices Ǎα τ j and Ãα τ j are independent of λ because of the compatibility equations. Moreover, the matrices ( Ǎα τ j ) 1≤j≤g (resp. ( Ãα τ j ) 1≤j≤g ) can be set traceless simultaneously by the additional gauge transformation Ψn = Ǧ Ψ (resp. Ψn = G Ψ) with

Ǧ = exp   - 1 2 g j=1 τ j
Tr( Ǎα τ j )(s)ds

  I 2 G = exp   - 1 2 g j=1 τ j
Tr( Ãα τ j )(s)ds

  I 2 (7-26)
Note that these additional gauge transformations do not change neither Ľ nor L. implies that ∂ λ Tr Ǎα τ = 0. Moreover, for g ≥ 2, if we denote (τ i ) 1≤i≤g a set of isomonodromic times, the compatibility of the Lax system also gives

∂ τ j [ Ǎα τ i ] = ∂ τ i [ Ǎα τ j ] + Ǎα τ j , Ǎα τ i , ∀ i = j. ( 7
α τ 1 = Ǎ(2)
α τ 1 and the result follows by induction. Finally, it is obvious that a gauge transformation independent of λ and proportional to I 2 does not change neither L nor L.

The last proposition shall be useful when L and Ľ are traceless. In this case, it is interesting to perform this additional gauge transformation in order to obtain a Lax pair that belongs to sl 2 (C) rather than gl 2 (C). In particular, this is always possible for the canonical choice of trivial times that shall be proposed in Section 8.

Canonical choice of trivial times and simplification of the Hamiltonian systems

The purpose of this section is to select some specific values of the trivial times in order to obtain simpler form of the Hamiltonian evolutions of Theorem 5.1. Indeed, the last section indicates (Theorem 7.2) that the shifted Darboux coordinates are independent of the values of the trivial times so that we may choose them without affecting the Hamiltonian evolutions. As it turns out, there exists a natural choice of the trivial times for which the Hamiltonian evolutions drastically simplify.

Canonical choice of the trivial times and main theorem

Definition 8.1 (Canonical choice of the trivial times). We define the "canonical choice of trivial times" by choosing

T ∞,k = 0 , ∀ k ∈ 0, r ∞ -1 , T 1 = 0, T 2 = 1. (8-1)
In the rest of the paper, we shall always set the trivial times to their canonical values. The canonical choice of trivial times implies that • All even irregular times are set to 0: for all k ∈ 1, r ∞ -1 : t ∞,2k = 0.

• t ∞,2r∞-3 = 2 and t ∞,2r∞-5 = 0.

• P1 is identically null. This is equivalent to say that L and Ľ are traceless. Hence, Proposition 7.6 implies that under a potential additional trivial gauge transformation, we may choose a gauge in which L, Ľ, Ãα τ and Ǎα τ are traceless for any isomonodromic time τ ∈ T iso .

• The shifted Darboux coordinates are identical to the initial Darboux coordinates:

∀ j ∈ 1, g : qj = q j and pj = p j (8-2)

• The isomonodromic times τ k identify with an irregular time:

∀ k ∈ 1, g : τ k = 1 2 t ∞,2r∞-2k-5 ⇔ 1 2 t ∞,2k-1 = τ r∞-k-2 (8-3) • P2 reduces to P2 (λ) = -λ if r ∞ = 3 or for r ∞ ≥ 4: P2 (λ) = -λ 2r∞-5 - 2r∞-7 k=r∞-2   2τ 2r∞-k-6 + r∞-3 m=k-r∞+6 τ r∞-m-2 τ r∞-k+m-5   λ k -2τ r∞-3 + r∞-3 m=3 τ r∞-m-2 τ m-2 λ r∞-3 (8-4)
In other words, for r ∞ ≥ 4, we have

P (2) ∞,2r∞-4 = 0 P (2) ∞,2r∞-5 = -1 P (2) ∞,2r∞-6 = 0 P (2) ∞,k = - 2r∞-7 k=r∞-2 2τ 2r∞-k-6 + r∞-3 m=k-r∞+6 τ r∞-m-2 τ r∞-k+m-5 , ∀ k ∈ r ∞ -2, 2r ∞ -7 P (2) ∞,r∞-3 = -2τ r∞-3 + r∞-3 m=3 τ r∞-m-2 τ m-2 (8-5) • Coefficients c (ατ ) ∞,k 1≤k≤r∞-1
are vanishing for any isomonodromic deformation τ ∈ T iso .

• The gauge matrices G 1 (λ) and J(λ) of Proposition 2.2 simplifies to

G 1 (λ) = I 2 , J(λ) =    1 0 - g i=1 pi λ-q i j =i 1 qi -q j 1 g j=1 (λ-q j )    (8-6)
In particular, Ľ(λ) = L(λ) and for all τ ∈ T iso , Ǎ(α τ ) (λ) = Ã(α τ ) (λ).

We also get the explicit expression Proposition 8.1. Under the canonical choice of trivial times given by Definition 8.1, the Lax matrices L is given by

L1,1 (λ, ) = -Q(λ, ), L1,2 (λ, ) = g j=1 (λ -qj ), L2,2 (λ, ) = Q(λ, ), L2,1 (λ, ) = ∂ Q(λ, ) g j=1 (λ-q j ) ∂λ + L 2,1 (λ, ) g j=1 (λ -qj ) - Q(λ, ) 2 g j=1 (λ -qj ) (8-7) with L 2,1 (λ, ) = -P2 (λ)+ r∞-4 k=0 H ∞,k λ k - g j=1 pj λ-q j and Q(λ, ) = - g i=1 pi j =i
λ-q j qi -q j and L 2,2 (λ, ) = g j=1 λ-q j . Similarly, the matrix A α τ (λ, ) is given by

[A α τ (λ, )] 1,1 = - g j=1 µ (α τ ) j pj λ -qj [A α τ (λ, )] 1,2 = g j=1 µ (α τ ) j λ -qj [A α τ (λ, )] 2,1 = - g j=1 i =j µ (α τ ) j (λ -qj )(λ -qi ) +   g j=1 µ (α τ ) j λ -qj   -P2 (λ) + r∞-4 k=0 H ∞,k λ k [A α τ (λ, )] 2,2 = - g j=1 µ (α τ ) j pj λ -qj + g j=1 i =j µ (α τ ) j (λ -qj )(λ -qi ) (8-8)
We may also simplify Propositions 6.3 and 6.4.

Proposition 8.2. Under the canonical choice of trivial times given by Definition 8.1, the Lax matrices L and Ãα may be expressed in terms of symmetric Darboux coordinates as follow. For any τ ∈ T iso :

L1,1(λ) = - g-1 j=0 (-1) j-1 g i=j+1 PiQi-j-1 λ j L1,2(λ) = g m=0 (-1) g-m Qg-mλ m L2,2(λ) = g-1 j=0 (-1) j-1 g i=j+1 PiQi-j-1 λ j L2,1(λ) = - r∞-2 i=0 2r∞-5 j=g+i P (2) ∞,j hj-g-i({q}) λ i - g-2 i=0 g-1 j 1 =i+1 g-1 j 2 =g+i-j 1 (-1) j 1 +j 2 g i 1 =j 1 +1 Pi 1 Qi 1 -j 1 -1 g i 2 =j 2 +1 Pi 2 Qi 2 -j 2 -1 hj 1 +j 2 -g-i({q}) λ i [ Ãα τ (λ)]1,1 = - g-2 i=0 g-1-i m=1 (-1) i+m-1 ν (α) ∞,m g r=i+m+1 PrQr-i-m-1 λ i [ Ãα τ (λ)]1,2 = g-1 j=0 g-j m=1 (-1) g-j-m ν (α) ∞,m Qg-j-m λ j [ Ãα τ (λ)]2,2 = -[ Ã(α) (λ)]1,1 [ Ãα τ (λ)]2,1 = - g i=0 2r∞-5 j=Max(g,g+i-1) j-g-i m=1 ν (α) ∞,m hj-g-m-i({q}) P (2) ∞,j λ i - g i=0 g-1 j 1 =0 g-1 j 2 =0 j 1 +j 2 -g-i m=1 (-1) j 1 +j 2 ν (α) ∞,m hj 1 +j 2 -g-m-i({q}) g r 1 =j 1 +1 g r 2 =j 2 +1 Pr 1 Pr 2 Qr 1 -j 1 -1Qr 2 -j 2 -1λ i (8-9)
where P

∞,k r∞-3≤k≤2r∞-4

are determined by (8-5) and (h k ({q})) k≥0 are expressed in terms of symmetric Darboux coordinates by h 0 ({q}) = 1 and

h k ({q}) = k j=1 (-1) j b 1 ,...,b j ∈ 1,k j b 1 +•••+b j =k j m=1 (-1) bm Q bm , ∀ k ∈ 1, g (8-10)
Coefficients ν

(α τ ) ∞,k 1≤k≤r∞-3
shall be given by Proposition 8.3 depending on the isomonodromic deformation τ ∈ T iso and ν

(α) ∞,r∞-2 = g k=1 (-1) g-k ν (α) ∞,k Q g+1-k .
We shall now apply Theorem 5.1 for the canonical values of the trivial times and obtain our third main theorem. Theorem 8.1 (Hamiltonian representation for the canonical choice of trivial times). The canonical choice of the trivial times given by Definition 8.1 and the definition of trivial times (Definition 7.2) imply that for any isomonodromic time τ ∈ T iso :

Ham (α τ ) (q, p) = r∞-4 k=0 ν (α τ ) ∞,k+1 H ∞,k (8-11)
In other words, the Hamiltonian is a (time-dependent) linear combination of the isospectral Hamiltonians (H ∞,k ) 0≤k≤r∞-4 that are determined by

        1 q1 . . . . . . qg-1 1 1 q2 . . . . . . qg-1 2 . . . . . . . . . . . . 1 qg . . . . . . qg-1 g               H ∞,0 . . . . . . H ∞,r∞-4       =          p2 1 + P2 (q 1 ) + i =1 pi -p 1 q1 -q i . . . . . . p2 g + P2 (q g ) + i =g pi -pg qg-qi          (8-12)
where P2 is given by (8-4). Coefficients ν

(α τ ) ∞,k 1≤k≤r∞-3
shall be given by Proposition 8.3 depending on the isomonodromic deformation τ ∈ T iso . In terms of symmetric Darboux coordinates, the Hamiltonian is given by:

Ham (α τ ) (Q, P) = - g i=1 ν (α τ ) ∞,i g k=i+1 (-1) i (g -i)P k Q k-1-i + k-1 m=i+1 (-1) m P k Q k-1-m S m-i ({q}) + g i=1 ν (α τ ) ∞,i g k 1 =1 g k 2 =1 P k 1 P k 2 (-1) i-1 Min(k 1 -1,i-1)
r 1 =Max(0,i-k 2 ) Q k 1 -1-r 1 Q k 2 -i+r 1 + 0≤r 1 ≤k 1 -1 0≤r 2 ≤k 2 -1 r 1 +r 2 ≥g (-1) r 1 +r 2 Q k 1 -1-r 1 Q k 2 -1-r 2 g m=i (-1) g-m Q g-m h r 1 +r 2 +m-i-g+1 ({q}) + g i=1 ν (α τ ) ∞,i 2r∞-5 r=g g m=i (-1) g-m P (2) ∞,r Q g-m h r+m-i-g+1 ({q}) (8-13)
where (S k ({q})) 0≤k≤g and (h k ({q})) 0≤k≤g are determined by h 0 ({q}) = 1, S 0 ({q}) = g and for all k ∈ 1, g :

h k ({q}) = k j=1 (-1) j b 1 ,...,b j ∈ 1,k j b 1 +•••+b j =k j m=1 (-1) bm Q bm , ∀ k ∈ 1, g S k ({q}) = (-1) k k b 1 +2b 2 +•••+kb k =k b 1 ≥0,...,b k ≥0 (-1) b 1 +•••+b k (b 1 + • • • + b k ) b 1 + • • • + b k b 1 , . . . , b k k i=1 Q b i i (8-14)
Proof. The proof is obvious since the canonical choice of trivial times implies that the coefficients c

(ατ ) ∞,k 1≤k≤r∞-1
are vanishing for any isomonodromic deformation. Moreover, Remark 7.1 implies that ν

(α τ ) ∞,-1 = ν (α τ ) ∞,0 = 0.
Note that only the coefficients ν

(α τ ) ∞,k 1≤k≤r∞-3
of the linear combination depend on the deformation, since the isospectral Hamiltonians are independent of it. We shall now obtain their explicit values from the simplification of (4-8) depending on the choice of isomonodromic time τ j with j ∈ 1, g .

Proposition 8.3 (Expression of ν (α τ j )
∞,k ). For any j ∈ 1, g , the coefficients ν

(α τ j ) ∞,k 1≤k≤r∞-3
are determined under the canonical choice of trivial times of Definition 8.1 by

            1
τ g-2 τ g-3 . . . τ 1 0 1                    ν (α τ j ) ∞,1 . . . . . . ν (α τ j ) ∞,r∞-3        = 2 2r ∞ -2j -5 e j (8-15)
One may also obtain a simplified expression for µ

(α τ j ) i 1≤i≤g
from [START_REF] Bergère | Determinantal formulae and loop equations[END_REF][START_REF] Bertola | Hamiltonian structure of rational isomonodromic deformation systems[END_REF][START_REF] Bertola | The partition function of the two-matrix model as an isomonodromic τ function[END_REF][START_REF] Biquard | Wild non-abelian hodge theory on curves[END_REF][START_REF] Boalch | Symplectic manifolds and isomonodromic deformations[END_REF][START_REF] Boalch | Simply-laced isomonodromy systems[END_REF][START_REF] Boalch | Twisted local wild mapping class groups: configuration spaces, fission trees and complex braids[END_REF][START_REF] Chekhov | Hermitean matrix model free energy: Feynman graph technique for all genera[END_REF]. We find for the canonical choice of trivial times:

             1
                      1
                µ (α τ j ) 1 . . . . . . µ (α τ j ) g        = 2 2r∞ -2j -5 ej (8-16)
for all j ∈ 1, g .

Finally, the evolution equations under this canonical choice reduce to

∂ τ j qm = 2µ (α τ j ) m pm - i =m µ (α τ j ) m + µ (α τ j ) i qm -qi , ∂ τ j pm = i =m (µ (α τ j ) i + µ (α τ j ) m )(p i -pm ) (q m -qi ) 2 + µ (α τ j ) m -P 2 (q m ) + r∞-4 k=1 kH ∞,k qk-1 m (8-17)
for all (j, m) ∈ 1, g 2 .

Let us now underline a few aspects of Theorem 8.1:

• Note that Theorem 8.1 is valid for any choice of the trivial times and not only the canonical choice of Definition 8.1. Indeed (q j , pj ) 1≤j≤g are independent of the choice of trivial times so that we may choose any value to compute the Hamiltonian system. However, for other choices of trivial times, the connection between Darboux coordinates and shifted Darboux coordinates and the relation between irregular times and isomonodromic times may be more complex and is given by Definitions 7.3 and 7.2. This observation was already made in the non-twisted case in [START_REF] Marchal | Hamiltonian representation of isomonodromic deformations of general rational connections on[END_REF].

• Note that the Hamiltonian system for (q j , pj ) 1≤j≤g does not depend on P1 . Since P1 = Tr L, this means that the non-trivial isomonodromic evolutions are the same in the study of isomonodromic deformations of twisted connections on gl 2 (C) or sl 2 (C). However, in the case of gl 2 (C), apparent singularities are no longer the right Darboux coordinates and a shift by -1 2 P1 (q j ) becomes necessary (Cf. Definition 7.3). This observation was already made in the non-twisted case in [START_REF] Marchal | Hamiltonian representation of isomonodromic deformations of general rational connections on[END_REF].

• Hamiltonian evolutions only depend on τ g and τ g-1 through P2 (because the matrix in Proposition 8.3 does not depend neither on τ g-1 nor τ g ) so that the explicit dependence of the Hamiltonians in τ g and τ g-1 is linear. The explicit dependence of the Hamiltonians in (τ i ) 1≤i≤g-2 is polynomial and the corresponding degrees are given by [START_REF] Boalch | Symplectic manifolds and isomonodromic deformations[END_REF][START_REF] Boalch | Simply-laced isomonodromy systems[END_REF][START_REF] Boalch | Twisted local wild mapping class groups: configuration spaces, fission trees and complex braids[END_REF][START_REF] Chekhov | Hermitean matrix model free energy: Feynman graph technique for all genera[END_REF][START_REF] Chekhov | Hermitian matrix model free energy: Feynman graph technique for all genera[END_REF][START_REF] Chekhov | Matrix eigenvalue model: Feynman graph technique for all genera[END_REF][START_REF] Dumitrescu | Quantization of spectral curves for meromorphic Higgs bundles through topological recursion[END_REF][START_REF] Eynard | Quantization of classical spectral curves via topological recursion[END_REF][START_REF] Eynard | Invariants of algebraic curves and topological expansion[END_REF][START_REF] Eynard | Topological recursion in random matrices and enumerative geometry[END_REF][START_REF] Fuchs | Sur quelques équations différentielles linéaires du second ordre[END_REF][START_REF] Gaiur | Isomonodromic deformations: Confluence, reduction & quantisation[END_REF][START_REF] Gambier | Sur les équations différentielles du second ordre et du premier degré dont l'intégrale générale est à points critiques fixes[END_REF][START_REF] Garnier | Solution du problème de riemann pour les systèmes différentiels linéaires du second ordre[END_REF][START_REF] Hurtubise | On the geometry of isomonodromic deformations[END_REF].

Explicit expressions for the inverse of the matrices

One may invert the Vandermonde matrix in Theorem 8.1 in order to have some explicit expressions for (H ∞,k ) 0≤k≤r∞-4 . For all i ∈ 1, r ∞ -3 we find

H ∞,i-1 = r∞-3 m=1 (-1) g-i e g-i ({q 1 , . . . , qg } \ {q m }) r =m (q m -qr )   p2 m + P2 (q m ) + s =m ps -pm qm -qs   (8-18)
where we have defined the elementary symmetric functions by e m ({y 1 , . . .

y k }) = 1≤j 1 <•••<jm≤k y j 1 . . . y jm , ∀ m ≥ 0, k ≥ 1 (8-19)
Similarly, one may invert the lower triangular Toeplitz matrix of Proposition 8.3 in order to have an explicit expression for ν

(α τ j ) ∞,k . We find ν (α τ j ) k = 2 2r ∞ -2j -5 (δ j,k + F j-k-1 (τ 1 , . . . , τ k-j-1 )δ k≤j-2 ) , ∀ (j, k) ∈ 1, g 2 (8-20)
where we have defined:

            1
τ g-2 τ g-3 . . . τ 1 0 1             -1 =             1
F g-2 (τ 1 , . . . , τ g-2 ) F g-3 (τ 1 , . . . , τ g-3 ) . . . F 1 (τ 1 ) 0 1             (8-21)
with

F i (τ 1 , . . . , τ i ) = (b 1 ,...,b i )∈N i i j=1 (j+1)b j =i+1 i j=1 b j b 1 , . . . , b i (-1) i j=1 b j τ b 1 1 . . . τ b i i , ∀ i ≥ 1 (8-22)
For example, the first values of (F i (τ 1 , . . . , τ i )) 1≤i≤5 are

F 1 (τ 1 ) = -τ 1 F 2 (τ 1 , τ 2 ) = -τ 2 F 3 (τ 1 , τ 2 , τ 3 ) = τ 2 1 -τ 3 F 4 (τ 1 , τ 2 , τ 3 , τ 4 ) = 2τ 1 τ 2 -τ 4 F 5 (τ 1 , . . . , τ 5 ) = -τ 3 1 + 2τ 1 τ 3 + τ 2 2 -τ 5 (8-23)
Finally, we may obtain an explicit expression for

µ (α τ k ) i = 2 (2r ∞ -2i -5) g m =i (q i -qm ) (-1) g-i e g-i ({q 1 , . . . , qg } \ {q i }) + r∞-3 r=i+2 (-1) g-r e g-r ({q 1 , . . . , qg } \ {q i }) (b1,...,br-i-1)∈N r-i-1 i s=1 (s+1)bs=r-i r-i-1 s=1 b s b 1 , . . . , b r-i-1 (-1) r-i-1 s=1 bs τ b1 1 . . . τ br-i-1 r-i-1 (8-24) for all (k, i) ∈ 1, r ∞ -3 2 .
9 Topological type property and formal WKB solutions

Starting from twisted meromorphic connections on gl 2 (C) with a pole at infinity, we have obtained some isomonodromic times (τ j ) 1≤j≤g and some Lax pairs L(λ, τ , ), Ãα τ 1 (λ, τ , ), . . . , Ãα τg (λ, τ , ) corresponding to the compatible differential systems

∂ λ Ψ(λ, τ , ) = L(λ, τ , ) , ∂ τ j Ψ(λ, τ , ) = Ãα τ j (λ, τ , ) Ψ(λ, τ , ) , ∀ j ∈ 1, g (9-1)
These matrices are expressed in terms of the isomonodromic times and the Darboux coordinates (q i , pi ) 1≤i≤g satisfying some Hamiltonian systems. This construction is independent of the type of solutions, in particular in [START_REF] Marchal | Isomonodromic deformations of a rational differential system and reconstruction with the topological recursion: the sl 2 case[END_REF], it is argued that the most general formal solutions are expected to be formal -transseries. However, one may look for a simpler form of solutions. Of particular interests are formal power series solutions of the Hamiltonian systems:

qi (τ , ) = ∞ k=0 q (k) i (τ ) k , pi (τ , ) = ∞ k=0 p (k) i (τ ) k , ∀ i ∈ 1, g (9-2)
that equivalently correspond to formal WKB solutions

Ψ(λ, τ , ) = exp ∞ k=-1 Ψ k (λ, τ ) k (9-3)
of the Lax system. In [START_REF] Marchal | Isomonodromic deformations of a rational differential system and reconstruction with the topological recursion: the sl 2 case[END_REF], the authors proved that, in this formal WKB setup, the Lax systems arising from general isomonodromic deformations (twisted or not) always satisfy the socalled "Topological Type property" of [START_REF] Bergère | Rational differential systems, loop equations, and application to the q th reductions of KP[END_REF]. In particular, the central argument (section 4.2 of [START_REF] Marchal | Isomonodromic deformations of a rational differential system and reconstruction with the topological recursion: the sl 2 case[END_REF]) to prove the topological type property is the existence of an isomonodromic time τ (built from isospectral deformations in [START_REF] Marchal | Isomonodromic deformations of a rational differential system and reconstruction with the topological recursion: the sl 2 case[END_REF]) for which the auxiliary matrix Ãα τ (λ, ) is of the form Ãα

τ (λ, ) = B 1 λ+B 0 p(λ)
where B 0 and B 1 are independent of λ and p is a polynomial. Our formalism generates a similar result without using isospectral deformations. Indeed, it is obvious that the isomonodromic time τ ∞,r∞-3 provides a matrix Ãα τ r∞-3 that satisfies the condition presented above. Thus, in the context of formal WKB solutions (or equivalently of formal power series solutions of the Hamiltonian systems), the Lax pair satisfies the topological type property and one may reconstruct the formal correlation functions (W n (λ 1 , . . . , λ n )) n≥1 built from "determinantal formulas" (see [START_REF] Bergère | Determinantal formulae and loop equations[END_REF] for definitions) of the differential system ∂ λ Ψ(λ, ) = L(λ, ) Ψ(λ, ) by the formal -power series of the Eynard-Orantin differentials (ω k,n ) k≥0,n≥0 produced by the topological recursion on the classical spectral curve (that always reduces in this formal WKB setup to a genus 0 curve):

W n (λ 1 , . . . , λ n ; τ , ) = ∞ k=0 ω k,n (λ 1 , . . . , λ n ; τ ) n-2+2k , ∀ n ≥ 1 (9-4)
Moreover, the formal Jimbo-Miwa-Ueno τ -function τ JMU [START_REF] Jimbo | Monodromy preserving deformation of linear ordinary differential equations with rational coefficients: I. general theory and τ -function[END_REF][START_REF] Bertola | The partition function of the two-matrix model as an isomonodromic τ function[END_REF] is reconstructed by the free energies (ω k,0 ) k≥0 produced by the topological recursion

ln τ JMU (τ , ) = ∞ k=0 ω k,0 (τ ) 2k-2 (9-5)
We stress again that the Hamiltonian systems and Lax pairs obtained in this article do not depend on the type of the solutions considered. As explained above, when considering solutions expressed as formal power series or formal WKB series in , the picture simplifies since the genus of the classical spectral curve drops to 0, the topological type property is verified and one may reconstruct the formal correlation functions or the formal tau-function of the Lax system directly from the topological recursion. Nevertheless, it is presently an open question to prove that the same picture remains valid when considering more general solutions of the Lax system like -transseries solutions. Even in the formal WKB setup, the issue of giving some analytic meaning to the formal solutions is currently a widely open question.

Examples

Let us now apply the general theory to the first cases of the Painlevé 1 hierarchy.

The Airy case: r ∞ = 3

The Airy case corresponds to r ∞ = 3 so that g = 0. The canonical choice of trivial times corresponds to t ∞,4 = t ∞,2 = 0, t ∞,3 = 2 and t ∞,1 = 0 so that P2 (λ) = -λ. There is no Darboux coordinates and any Hamiltonian evolutions. However, one may still write the Lax matrices L and L. They are given by

L(λ) = 0 1 -λ 0 = L(λ) (10-1)
giving the Airy spectral curve: y 2 = λ. Since g = 0, the only interesting result provided by the paper is that the wave function Ψ may be reconstructed by topological recursion after the introduction of the formal parameter . This is of course in agreement with known results about the Airy spectral curve [START_REF] Kontsevich | Intersection theory on the moduli space of curves and the matrix Airy function[END_REF][START_REF] Eynard | Topological recursion in random matrices and enumerative geometry[END_REF][START_REF] Dumitrescu | Quantization of spectral curves for meromorphic Higgs bundles through topological recursion[END_REF].

10.2 Painlevé 1 case: r ∞ = 4
Let us consider r ∞ = 4, i.e. g = 1 corresponding to the Painlevé 1 case. In this setup, the canonical choice of trivial times corresponds to t ∞,6 = t ∞,4 = t ∞,2 = 0, t ∞,5 = 2 and t ∞,3 = 0. The only non-trivial isomonodromic time is τ := τ 1 = 1 2 t ∞,1 . Since g = 1, we shall drop the useless index in this case (i.e. q := q 1 , p := p 1 , etc.). Application of the general results to this case is straightforward and give under the choice of trivial times made in Definition 8.1:

P2 (λ) = -λ 3 -2τ λ H ∞,0 = p2 + P2 (q) = p2 -q3 -2τ q Q(λ, ) = -p ν (α τ ) ∞,1 = 2 , µ (α τ ) = 2 (10-2)
Thus, we get that the Hamiltonian is

Ham (α τ ) (q, p) = 2H ∞,0 = 2p 2 -2q 3 -4τ q (10-3)
It corresponds to the ordinary differential equations ∂ τ q = 4p , ∂ τ p = 6q 2 + 4τ so that q(τ ) satisfies a Painlevé 1 like equation:

2 d 2 q dτ 2 = 24q 2 + 16τ (10-5)
The associated Lax pairs are given by

L(λ, ) = 0 1 λ 3 + 2τ λ + p2 -q3 -2τ q -p λ-q λ-q , A α τ (λ, ) = -2p λ-q 2 λ-q 2 λ-q λ 3 + 2τ λ + p2 -q3 -2τ q -2p λ-q (10-6) or equivalently L(λ, ) = p λ - q λ 2 + qλ + q2 + 2τ -p , Ãα τ (λ, ) = 0 2 2(λ + 2q) 0 (10-7)
Remark 10.1. If we perform t = 2 6 5 τ , q = 2 -2 5 q, p = 2 2 5 p we find that q(t) satisfies the normalized Painlevé 1 equation:

2 d 2 q dt 2 = 6q 2 + t (10-8)
Moreover, one may recover the Jimbo-Miwa Lax pair (eq. C.2 of [START_REF] Jimbo | Monodromy preserving deformation of linear ordinary differential equations with rational coefficients[END_REF]):

L JMU (x) = -z(t) x 2 + y(t)x + y(t) 2 + t 2 4(x -y(t)) z(t) , A JMU (x) = 0 x 2 + y(t) 2 0 (10-9)
Ψ JMU (x) = 0 1 2

6 5 0 Ψ(λ) , x = 2 -2 5 λ , t = 2 6 5 τ , y(t) = q(t) , z(t) = p(t) (10-10)
10.3 Second element of the Painlevé 1 hierarchy: r ∞ = 5

Let us consider r ∞ = 5, i.e. g = 2 corresponding to the second element of the Painlevé 1 hierarchy. In this setup, the canonical choice of trivial times corresponds to t ∞,8 = t ∞,6 = t ∞,4 = t ∞,2 = 0, t ∞,7 = 2 and t ∞,5 = 0. The only non-trivial isomonodromic times are

τ 1 = 1 2 t ∞,3 and τ 2 = 1 2 t ∞,1 . We have also P2 (λ) = -λ 5 -2τ 1 λ 3 -2τ 2 λ 2 Q(λ, ) = - p1 (λ -q2 ) q1 -q2 - p2 (λ -q1 ) q2 -q1 = (p 2 -p1 )λ + p1 q2 -p2 q1 q1 -q2 (10-11)
Coefficients (H ∞,0 , H ∞,1 ) are determined by [START_REF] Boalch | Symplectic manifolds and isomonodromic deformations[END_REF][START_REF] Boalch | Simply-laced isomonodromy systems[END_REF][START_REF] Boalch | Twisted local wild mapping class groups: configuration spaces, fission trees and complex braids[END_REF][START_REF] Chekhov | Hermitean matrix model free energy: Feynman graph technique for all genera[END_REF][START_REF] Chekhov | Hermitian matrix model free energy: Feynman graph technique for all genera[END_REF]:

H ∞,0 = (q 1 p2 2 -q2 p2 1 ) (q 1 -q2 ) 2 - (p 1 -p2 ) (q 1 -q2 ) + (q 1 + q2 )q 1 q2 (q 2 1 + q2 2 + 2τ 1 ) + 2τ 2 q1 q2 H ∞,1 = p2 1 -p2 2 q1 -q2 -2τ 1 (q 2 1 + q1 q2 + q2 2 ) -2τ 2 (q 1 + q2 ) -q4 1 -q3 1 q2 -q2 1 q2 2 -q1 q3 2 -q4 2 (10-12) Coefficients ν (α τ 1 ) ∞,1 , ν (α τ 1 ) ∞,2 , ν (α τ 2 ) ∞,1 , ν (α τ 2 ) ∞,2
are determined by Proposition 8.3 whose l.h.s. is trivial for g = 2 so that

ν (α τ 1 ) ∞,1 = 2 3 , ν (α τ 1 ) ∞,2 = 0 ν (α τ 2 ) ∞,1 = 0 , ν (α τ 2 ) ∞,2 = 2 (10-13) Coefficients µ (α τ 1 ) 1 , µ (α τ 1 ) 2 , µ (α τ 2 ) 1 , µ (α τ 2 ) 2
are determined by [START_REF] Boalch | Symplectic manifolds and isomonodromic deformations[END_REF][START_REF] Boalch | Simply-laced isomonodromy systems[END_REF][START_REF] Boalch | Twisted local wild mapping class groups: configuration spaces, fission trees and complex braids[END_REF][START_REF] Chekhov | Hermitean matrix model free energy: Feynman graph technique for all genera[END_REF][START_REF] Chekhov | Hermitian matrix model free energy: Feynman graph technique for all genera[END_REF][START_REF] Chekhov | Matrix eigenvalue model: Feynman graph technique for all genera[END_REF][START_REF] Dumitrescu | Quantization of spectral curves for meromorphic Higgs bundles through topological recursion[END_REF][START_REF] Eynard | Quantization of classical spectral curves via topological recursion[END_REF][START_REF] Eynard | Invariants of algebraic curves and topological expansion[END_REF]:

µ (α τ 1 ) 1 = - 2q 2 3(q 1 -q2 ) , µ (α τ 1 ) 2 = 2q 1 3(q 1 -q2 ) µ (α τ 2 ) 1 = 2 q1 -q2 , µ (α τ 2 ) 2 = - 2 q1 -q2 (10-14)
The Hamiltonians are

Ham (α τ 1 ) (q, p) = 2 3 (q 1 p2 2 -q2 p2 1 ) (q 1 -q2 ) 2 - (p 1 -p2 ) (q 1 -q2 ) + (q 1 + q2 )q 1 q2 (q 2 1 + q2 2 + 2τ 1 ) + 2τ 2 q1 q2 = 2 3 P 2 2 (Q 2 1 -Q 2 ) + P 2 1 + 2P 1 P 2 Q 1 + P 2 Ham (α τ 2 ) (q, p) = 2 p2 1 -p2 2 q1 -q2 -2τ 1 (q 2 1 + q1 q2 + q2 2 ) -2τ 2 (q 1 + q2 ) -q4 1 -q3 1 q2 -q2 1 q2 2 -q1 q3 2 -q4 2 = 2 -P 2 2 Q 1 -2P 1 P 2 -Q 4 1 + 3Q 3 1 Q 2 -Q 2 2 + 2(Q 2 -Q 2 1 )τ 1 -2Q 1 τ 2 (10-15)
where

Q 1 = q1 + q2 , Q 2 = q1 q2 , P 1 = q1 p1 -q2 p2 q1 -q2 , P 2 = - p1 -p2 q1 -q2 (10-16)
The Lax matrices are

L(λ, ) = 0 1 λ 5 + 2τ 1 λ 3 + 2τ 2 λ 2 + H ∞,1 λ + H ∞,0 λ-q1 + λ-q2 [A α τ 1 (λ, )] 1,1 = 2p 1 q2 3(q 1 -q2 )(λ -q1 ) - 2q 1 p2 3(q 1 -q2 )(λ -q2 ) [A α τ 1 (λ, )] 1,2 = - 2q 2 3(q 1 -q2 )(λ -q1 ) + 2q 1 3(q 1 -q2 )(λ -q2 ) [A α τ 1 (λ, )] 2,1 = 2 - q2 3(q 1 -q2 )(λ -q1 ) + q1 3(q 1 -q2 )(λ -q2 ) λ 5 + 2τ 1 λ 3 + 2τ 2 λ 2 + H ∞,1 λ + H ∞,0 - 2 3(λ -q1 )(λ -q2 ) [A α τ 1 (λ, )] 2,2 = 2q 2 p1 3(q 1 -q2 )(λ -q1 ) - 2q 1 p2 3(q 1 -q2 )(λ -q2 ) + 2 3(λ -q1 )(λ -q2 ) [A α τ 2 (λ, )] 1,1 = 2 - p1 (q 1 -q2 )(λ -q1 ) + p2 (q 1 -q2 )(λ -q2 ) [A α τ 2 (λ, )] 1,2 = 2 (q 1 -q2 )(λ -q1 ) - 2 (q 1 -q2 )(λ -q2 ) [A α τ 2 (λ, )] 2,1 = 2 1 (q 1 -q2 )(λ -q1 ) - 1 (q 1 -q2 )(λ -q2 ) λ 5 + 2τ 1 λ 3 + 2τ 2 λ 2 + H ∞,1 λ + H ∞,0 [A α τ 2 (λ, )] 2,2 = 2 - p1 (q 1 -q2 )(λ -q1 ) + p2 (q 1 -q2 )(λ -q2 ) (10-17)
or equivalently

L(λ) = p1 -p2 q1 -q 2 λ + q1 p2 -q 2 p1 q1 -q 2 (λ -q1)(λ -q2) λ 3 + (q1 + q2)λ 2 + (q 2 1 + q1 q2 + q2 2 + 2τ1)λ -( p1 -p2 ) 2 (q 1 -q 2 ) 2 + (q 2 1 + q2 2 + 2τ1)(q1 + q2) + 2τ2 -p1 -p2 q1 -q 2 λ -q1 p2 -q 2 p1 q1 -q 2 = -P2λ -Q1P2 -P1 λ 2 -Q1λ + Q2 λ 3 + Q1λ 2 + (Q 2 1 -Q2 + 2τ1)λ -P 2 2 + Q1(Q 2 1 -2Q2) + 2Q1τ1 + 2τ2 P2λ + Q1P2 + P1 (10-18) Ãα τ 1 (λ) = 2(p 1 -p 2 ) 3(q 1 -q 2 ) 2 3 (λ -q1 -q2 ) 2 3 λ 2 + (q 1 + q2 )λ + q2 1 + q2 2 + 2τ 1 -2(p 1 -p 2 ) 3(q 1 -q 2 ) = -2 3 P 2 2 3 (λ -Q 1 ) 2 3 λ 2 + Q 1 λ + Q 2 1 -2Q 2 + 2τ 1 2 3 P 2 Ãα τ 2 (λ) = 0 2 2(λ + 2(q 1 + q2 )) 0 = 0 2 2(λ + 2Q 1 ) 0 (10-19)
10.4 Third element of the Painlevé 1 hierarchy: r ∞ = 6

Let us consider r ∞ = 6, i.e. g = 3 corresponding to the third element of the Painlevé 1 hierarchy. In particular, this case is the first case where the Hamiltonians are non-trivial linear combinations of the coefficients (H ∞,k ) 0≤k≤r∞-4 , this is due to the fact that the matrix of Proposition 8.3 is no longer diagonal. Indeed we have

Ham (α τ 1 ) = 2 5 H ∞,0 - 2 5 τ 1 H ∞,2 Ham (α τ 2 ) = 2 3 H ∞,1 Ham (α τ 3 ) = 2H ∞,2 (10-20)
In this setup, the canonical choice of trivial times corresponds to t ∞,10 = t ∞,8 = t ∞,6 = t ∞,4 = t ∞,2 = 0, t ∞,9 = 2 and t ∞,7 = 0. The only non-trivial isomonodromic times are

τ 1 = 1 2 t ∞,5 , τ 2 = 1 2 t ∞,3 and τ 3 = 1 2 t ∞,1 .
For compactness, we shall only present results expressed in terms of the symmetric Darboux coordinates (Q 1 , Q 2 , Q 3 , P 1 , P 2 , P 3 ). One may recover the expression in terms of shifted Darboux coordinates using

Q 1 = q1 + q2 + q3 , Q 2 = q1 q2 + q1 q3 + q2 q3 , Q 3 = q1 q2 q3 P 1 = q2 1 (q 2 -q3 )p 1 -q2 2 (q 1 -q3 )p 2 + q2 3 (q 1 -q2 )p 3 (q 1 -q2 )(q 1 -q3 )(q 2 -q3 ) P 2 = - q1 (q 2 -q3 )p 1 -q2 (q 1 -q3 )p 2 + q3 (q 1 -q2 )p 3 (q 1 -q2 )(q 1 -q3 )(q 2 -q3 ) P 3 = (q 2 -q3 )p 1 -(q 1 -q3 )p 2 + (q 1 -q2 )p 3 (q 1 -q2 )(q 1 -q3 )(q 2 -q3 ) (10-21)
We have

P2 (λ) = -λ 7 -2τ 1 λ 5 -2τ 2 λ 4 -(τ 2 1 + 2τ 3 )λ 3 Q(λ) = -P 3 λ 2 + (P 3 Q 1 + P 2 )λ -P 3 Q 2 -P 2 Q 1 -P 1 (10-22)
The Hamiltonians are

Ham (α τ 1 ) = 2 5 (-Q1Q3 + Q 2 2 )P 2 3 + 2(Q1 + Q2)P1P3 + Q 2 1 P 2 2 + 2Q1P1P2 + P 2 1 + 2(Q1Q2 -Q3)P2P3 -2Q1P1P3 -Q3(Q 4 1 -3Q 2 1 Q2 + 2Q1Q3 + Q 2 2 ) + 2(Q 2 1 -Q2)τ1τ2 + 2Q1τ1τ3 + Q1τ 3 1 + (2Q 3 1 -4Q1Q2 + Q3)τ 2 1 -(-Q 5 1 + 4Q 3 1 Q2 + 2Q1P2P3 + Q2P 2 3 -Q 2 1 Q3 -3Q1Q 2 2 + 2P1P3 + P 2 2 )τ1 -2Q1Q3τ2 -2Q3τ3 + (Q1P3 + 2P2) Ham (α τ 2 ) = 2 3 -2Q 2 1 P2P3 -2Q1P1P3 -2P1P2 + (Q3 -Q1Q2)P 2 3 -2Q1P 2 2 + 4Q1Q2Q3 + Q 4 1 Q2 -Q 2 1 Q 2 2 -Q 3 1 Q3 + Q 3 2 -Q 2 3 + Q2τ 2 1 + 2(Q 2 1 Q2 -Q1Q3 -Q 2 2 )τ1 + 2(Q1Q2 -Q3)τ2 + 2Q2τ3 -P3 Ham (α τ 3 ) = 2 2Q1P2P3 + Q2P 2 3 + 2P1P3 + P 2 2 -Q 5 1 + 4Q 3 1 Q2 -3Q 2 1 Q3 -3Q1Q 2 2 + 2Q2Q3 -Q1τ 2 1 + 2(2Q1Q2 -Q 3 1 -Q3)τ1 + 2(Q2 -Q 2 1 )τ2 -2Q1τ3 (10-23)
The Lax matrices are

L1,1 (λ) = P 3 λ 2 -(P 2 + Q 1 P 3 )λ + P 1 + Q 2 P 3 + Q 1 P 2 L1,2 (λ) = λ 3 -Q 1 λ 2 + Q 2 λ -Q 3 L2,1 (λ) = λ 4 + Q 1 λ 3 + (Q 2 1 -Q 2 + 2τ 1 )λ 2 + (-P 2 3 + Q 3 1 + Q 3 -2Q 1 Q 2 + 2Q 1 τ 1 + 2τ 2 )λ 2P 2 P 3 + Q 1 P 2 3 + Q 4 1 -3Q 2 Q 2 1 + 2Q 3 Q 1 + Q 2 2 + τ 2 1 + 2(Q 2 1 -Q 2 )τ 1 + 2Q 1 τ 2 + 2τ 3 L2,2 (λ) = -P 3 λ 2 + (P 2 + Q 1 P 3 )λ -P 1 -Q 2 P 3 -Q 1 P 2 (10-24)
and

Ãα τ 1 (λ) = 2 5 (P3λ -Q1P3 -P2) 2 5 λ 2 -Q1λ + Q2 -τ1 2 5 λ 3 + Q1λ 2 + (Q 2 1 -Q2 + τ1)λ -P3 + Q 3 1 -2Q1Q2 + 2Q3 + 2τ2 -2 5 (P3λ -Q1P3 -P2) Ãα τ 2 (λ) = 2 3 P3 2 3 (λ -Q1) 2 3 λ 2 + Q1λ + Q 2 1 -2Q2 + 2τ1 -2 3 P3 Ãα τ 3 (λ) = 0 2 2(λ + 2Q1) 0 (10-25)

Outlooks

In this article, we complemented the results of [START_REF] Marchal | Hamiltonian representation of isomonodromic deformations of general rational connections on[END_REF] by dealing with the case of twisted meromorphic connections in gl 2 (C) obtaining explicit expressions of the Hamiltonians and Lax pairs in various sets of Darboux coordinates. Moreover, we provided a reduction of the initial space of irregular times (of dimension 2g + 4) to a set of non-trivial isomonodromic times of dimension g.

In particular, we recover in the twisted case the fact that meromorphic connections in gl 2 (C) are equivalent at the level of Hamiltonian systems to meromorphic connections in sl 2 (C), a point that was already raised in [START_REF] Marchal | Hamiltonian representation of isomonodromic deformations of general rational connections on[END_REF] in the non-twisted case. The method used in the present article opens the way to several generalizations:

• This article and [START_REF] Marchal | Hamiltonian representation of isomonodromic deformations of general rational connections on[END_REF] ends the study of meromorphic connections in gl 2 (C). Thus, a natural issue is to know if the present setup extends to gl d (C) with d ≥ 2. In principle, a similar strategy shall be used but it is unclear if all technical issues might be overcome when the dimension is arbitrary, especially in the twisted case where the underlying geometric construction is far less understood. We let this very exciting question for future works.

• This article deals with meromorphic connections in gl 2 (C). However, one may be interested in a more general abstract setup with any Lie algebra and not only a matrix representation of it. In this case, we believe that most of the results shall be generalized upon the adequate quantities and terminology.

• As mentioned in Section 3 and similarly to [START_REF] Marchal | Hamiltonian representation of isomonodromic deformations of general rational connections on[END_REF], this article makes the connection with formal WKB expansions and -transseries via topological recursion. At the level of isomonodromic deformations, the introduction of the formal parameter is done by a simple rescaling of irregular times (Section 2.5). A better understanding of the role of and its limit to 0 in the context of meromorphic connections would be interesting in order to address the issue of resumation and analytical properties associated to the formal (trans)series.

A Proof of Proposition 2.3

As mentioned in Proposition 2.1 there exists a local gauge transformation

G ∞ (z) = G ∞,-1 z + G ∞,0 + G ∞,1 z + . . . with z = λ 1 2 and G ∞,-1 of rank 1 such that Ψ ∞ = G ∞ Ψ is
given by (we recall that we added the extra parameter by rescaling)

Ψ∞(λ) = Ψ (reg) ∞ (z) diag exp - 1 2r∞-2 k=1 t ∞,k k z k + 1 2 ln z , exp - 1 2r∞-2 k=1 (-1) k t ∞,k k z k + 1 2 ln z =     Ψ (reg) ∞ (z) 1,1 exp -1 2r∞-2 k=1 t ∞,k k z k + 1 2 ln z Ψ (reg) ∞ (z) 1,2 exp -1 2r∞-2 k=1 (-1) k t ∞,k k z k + 1 2 ln z Ψ (reg) ∞ (z) 2,1 exp -1 2r∞-2 k=1 t ∞,k k z k + 1 2 ln z Ψ (reg) ∞ (z) 2,2 exp -1 2r∞-2 k=1 (-1) k t ∞,k k z k + 1 2 ln z     (A-1)
Since Ψ is normalized at infinity by (2-13), the verification is straightforward when one considers the highest order in the expansion

G∞,-1 = X 0 X 0 (A-2)
In particular, the expansion of

G -1 ∞ is of the form G -1 ∞ = Ĝ∞,0 + Ĝ∞,1z -1 + Ĝ∞,2z -2 +. . . with Ĝ∞,0 = 0 0 X X .
Thus, multiplying on the left by

G -1 ∞ provides Ψ(z) =     R (reg) ∞ (z) 1,1 exp -1 2r∞-2 k=1 t ∞,k k z k -1 2 ln z R (reg) ∞ (z) 1,2 exp -1 2r∞-2 k=1
(-1) k t ∞,k k z k -1 2 ln z R (reg) ∞ (z) 2,1 z exp -1 2r∞-2 k=1 t ∞,k k z k -1 2 ln z R (reg) ∞ (z) 2,2 z exp -1 2r∞-2 k=1 (-1) k t ∞,k k z k -1 2 ln z     (A-3)
where R

(reg)

∞ (z) i,j (i,j)∈ 1,2 2
are regular at infinity. Let us now prove that these asymptotics are consistent with the one proposed for Ψ. Since Ψ is the solution to a companion-like system, we have Ψ = ψ1 ψ2 ∂ λ ψ1 ∂ λ ψ2 .

Hence, equation (2-38) is equivalent to

Ψ1,1(λ) = exp - 1 2r∞-2 k=1 t ∞,k k z k - 1 2 ln z + O(1) Ψ1,2(λ) = exp - 1 2r∞-2 k=1 (-1) k t ∞,k k z k - 1 2 ln z + O(1) Ψ2,1(λ) = - 1 2 2r∞-2 k=1 t ∞,k z k-2 - 1 4z 2 + o(z -3 ) exp - 1 2r∞-2 k=1 t ∞,k k z k - 1 2 ln z + O(1) , Ψ2,2(λ) = - 1 2 2r∞-2 k=1 (-1) k t ∞,k z k-2 - 1 4z 2 + o(z -3 ) exp - 1 2r∞-2 k=1 (-1) k t ∞,k k z k - 1 2 ln z + O(1) , (A-4)
Thus, Ψ = G 1 JΨ is given by

Ψ(λ) =    1 0 1 2 t ∞,2r∞-2 λ + g 0 + Q(λ) g j=1 (λ-q j ) 1 g j=1 (λ-q j )    Ψ 1,1 (λ) Ψ 1,2 (λ) Ψ 2,1 (λ) Ψ 2,2 (λ) (A-5)
and behaves like

Ψ1,1(λ) = exp - 1 2r∞-2 k=1 t ∞,k k z k - 1 2 ln z + O(1) Ψ1,2(λ) = exp - 1 2r∞-2 k=1 (-1) k t ∞,k k z k - 1 2 ln z + O(1) Ψ2,1(λ) = exp - 1 2r∞-2 k=1 t ∞,k k z k - 1 2 ln z + O(1) 1 2 t∞,2r ∞-2 z 2 + g0 + - 1 2 t∞,2r ∞ -2z 2r∞-4 - 1 2 t∞,2r ∞-3 z 2r∞-5 + O(z 2r∞-6 ) z -2r∞+6 + g j=1 qjz 2r∞+4 + O(z 2r∞+2 ) = - 1 2 t∞,2r ∞-3 z + O(1) exp - 1 2r∞-2 k=1 t ∞,k k z k - 1 2 ln z + O(1) Ψ2,2(λ) = exp - 1 2r∞-2 k=1 (-1) k t ∞,k k z k - 1 2 ln z + O(1) 1 2 t∞,2r ∞-2 z 2 + g0 + - 1 2 t∞,2r ∞ -2z 2r∞-4 + 1 2 t∞,2r ∞-3 z 2r∞-5 + O(z 2r∞-6 ) z -2r∞+6 + g j=1 qjz 2r∞+4 + O(z 2r∞+2 ) = 1 2 t∞,2r ∞ -3 z + O(1) exp - 1 2r∞-2 k=1 (-1) k t ∞,k k z k - 1 2 ln z + O(1) (A-6)
in accordance with (A-3). Moreover, asymptotics (A-6) of Ψ implies by direct computations that ( ∂ λ Ψ) Ψ-1 = L may be set in the form given by (2-13).

B Proof of Proposition 2.4

From Proposition 2.3, the local asymptotics of the wave functions ψ 1 and ψ 2 are

ψ 1 (λ) λ→∞ = exp - 1 2r∞-2 k=1 t ∞,k k λ k 2 - 1 4 ln λ + O(1) ψ 2 (λ) λ→∞ = exp - 1 2r∞-2 k=1 (-1) k t ∞,k k λ k 2 - 1 4 ln λ + O(1) (B-1)
In particular, the Wronskian W (λ) = ψ 1 (λ) ∂ λ ψ 2 (λ) -ψ 2 (λ) ∂ λ ψ 1 (λ) is given by Definition 2.3:

W (λ) = W 0 g j=1 (λ -q j ) exp 1 λ 0 P1 (s)ds (B-2)
The standard relation between Tr L(λ) and the logarithmic derivative of the Wronskian provides the expected result of L 2,2 (λ, ). As an intermediate step we define

Y i (λ) = ψ i (λ) ∂ λ ψ i (λ). Then L 2,1 (λ) = -Y 1 (λ)Y 2 (λ) - Y 2 (λ)∂ λ Y 1 (λ) -Y 1 (λ)∂ λ Y 2 (λ) Y 2 (λ) -Y 1 (λ) (B-3)
One may thus study the asymptotic behavior of L 2,1 at λ → ∞. We have

Y 1 (λ) λ→∞ = - 1 2 2r∞-2 k=1 t ∞,k λ k 2 -1 - 4λ + O λ -3 2 Y 2 (λ) λ→∞ = - 1 2 2r∞-2 k=1 (-1) k t ∞,k λ k 2 -1 - 4λ + O λ -3 2 (B-4) so that - Y 2 (λ)∂ λ Y 1 (λ) -Y 1 (λ)∂ λ Y 2 (λ) Y 2 (λ) -Y 1 (λ) λ→∞ = 4 t ∞,2r∞-2 λ r∞-3 + O λ r∞-4 (B-5)
Hence we obtain:

L 2,1 (λ) λ→∞ = - 1 4 2r∞-2 i=1 2r∞-2 j=1 (-1) j t ∞,i t ∞,j λ i+j 2 -2 - 2 t ∞,2r∞-2 λ r∞-3 + O λ r∞-4 (B-6)
Note that terms with odd values of i + j cancel by symmetry, we end up with

L 2,1 (λ) λ→∞ = - 1 4
(i,j)∈ 1,2r∞-2 2 i+j even

(-1) j t ∞,i t ∞,j λ i+j 2 -2 + O λ r∞-4 λ→∞ = - 1 4 2r∞-2 k=r∞-2 2r∞-2 j=2k-2r∞+6 (-1) j t ∞,j t ∞,2k-j+4 λ k - 1 4 2r∞-3 j=1 (-1) j t ∞,j t ∞,2r∞-j-2 λ r∞-3 + O λ r∞-4 = -P2 (λ) + O λ r∞-4 (B-7)
Since L 2,1 (λ, ) is a rational function of λ with only poles at infinity, we get that it is a polynomial in λ and the previous asymptotics provide its leading coefficients.

C Proof of Proposition 4.1

Let us first observe that the entry [A α (λ)] 1,2 is given by

[A α (λ)] 1,2 = W α (λ) W (λ) = Z α,2 (λ) -Z α,1 (λ) Y 2 (λ) -Y 1 (λ) (C-1)
where we have defined

Z α,i (λ) = L α [ψ i (λ)] ψ i (λ) , ∀ i ∈ 1, 2 W α (λ) = L α [ψ 2 (λ)]ψ 1 (λ) -L α [ψ 1 (λ)]ψ 2 (λ) (C-2)
and from Proposition 2.3 we have:

Z α,1 (λ) = - 2r∞-2 k=1 α ∞,k k λ k 2 + O(1) Z α,2 (λ) = - 2r∞-2 k=1 (-1) k α ∞,k k λ k 2 + O(1) (C-3) so that Z α,2 (λ) -Z α,1 (λ) = 2 r∞-1 j=1 α ∞,2j-1 2j -1 λ j-1 2 + O(1) (C-4) Thus, since [A α (λ)] 1,2 (Y 2 (λ) -Y 1 (λ)) = Z α,2 (λ) -Z α,1 (λ) (C-5) using (B-4) Y 2 (λ) -Y 1 (λ) = r∞-1 j=1 t ∞,2j-1 λ j-3 2 + O λ -3 2 (C-6)
we obtain the form of the entry

A α (λ, ) 1,2 = r∞-3 j=-1 ν (α) ∞,j λ -j + O λ -(r∞-2) (C-7)
where the coefficients ν

(α) ∞,k -1≤k≤r∞-3
are recursively determined by

  r∞-3 j=-1 ν (α) ∞,j λ -j + O λ -r∞+2   r∞-1 k=1 t ∞,2k-1 λ k-1 + O λ -1 = 2 r∞-1 k=1 α ∞,2k-1 2k -1 λ k + O (1)
.

(C-8) These relations may be rewritten in a (r ∞ -1) × (r ∞ -1) lower triangular Toeplitz matrix form:

M ∞     ν (α) ∞,-1 . . . ν (α) ∞,r∞-3     =       2α ∞,2r∞-3 (2r∞-3) 2α ∞,2r∞-5 (2r∞-5) . . . 2α ∞,1 1       with M ∞ =          t ∞,2r∞-3 0 . . . 0 t ∞,2r∞-5 t ∞,2r∞-3 0 . . . . . . . . . . . . . . . . . . t ∞,3 . . . . . . 0 t ∞,1 t ∞,3 . . . t ∞,2r∞-3          (C-9)

D Proof of Proposition 4.3

A straightforward computation shows that

[A α (λ)] 1,1 = Z α,1 (λ)Y 2 (λ) -Z α,2 (λ)Y 1 (λ) Y 2 (λ) -Y 1 (λ) (D-1)
which we rewrite as

[A α (λ)] 1,1 (Y 2 (λ) -Y 1 (λ)) = Z α,1 (λ)Y 2 (λ) -Z α,2 (λ)Y 1 (λ) (D-2)
We proceed using (B-4) and (C-3) that give

Z α,1 (λ)Y 2 (λ) -Z α,2 (λ)Y 1 (λ) = 1 2 2r∞-2 j=1 2r∞-2 i=1 α ∞,j j t ∞,i ((-1) i -(-1) j )λ i+j 2 -1 + O λ r∞-2 = r∞-1 s=1 r∞-1 m=1 α ∞,2s-1 2s -1 t ∞,2m - α ∞,2s 2s t ∞,2m-1 λ m+s-3 2 + O λ r∞-2 (D-3)
We make use (C-6) also in order to obtain

A α (λ, ) 1,1 = r∞-1 j=1 c (α) ∞,j λ j + O (1) (D-4)
where the coefficients c

(α) ∞,k 1≤k≤r∞-1
are recursively determined by

r∞-1 j=1 c (α) ∞,j λ j + O (1) r∞-1 k=1 t ∞,2k-1 λ k-3 2 + O λ -3 2 = r∞-1 s=1 r∞-1 m=1 α∞,2s-1 2s -1 t∞,2m - α∞,2s 2s t∞,2m-1 λ m+s-3 2 +O λ r∞-2 (D-5)
These relations may be rewritten in a (r ∞ -1) × (r ∞ -1) lower triangular Toeplitz matrix form:

M ∞          c (α) ∞,r∞-1 . . . c (α) ∞,k . . . c (α) ∞,1          =             α ∞,2r∞-3 2r∞-3 t ∞,2r∞-2 - α ∞,2r∞-2 2r∞-2 t ∞,2r∞-3 . . . r∞-1 m=k α ∞,2k+2r∞-2m-3 2k+2r∞-2m-3 t ∞,2m - α 2k+2r∞-2m-2 2k+2r∞-2m-2 t ∞,2m-1 . . . r∞-1 m=1 α ∞,2r∞-2m-1 2r∞-2m-1 t ∞,2m - α 2r∞-2m 2r∞-2m t ∞,2m-1             (D-6)
Finally, the coefficients ρ

(α) j 1≤j≤g
are obtained by looking at order (λ-q j ) -3 of L α [L 2,1 (λ)].

E Proof of Theorem 5.1

This appendix is devoted for the proof of Theorem 5.1.

E.1 Preliminary results

We start with the following lemma:

Lemma E.1. For all j ∈ 1, g :

r∞-4 k=0 g i=1 H ∞,k q k i ∂ q j µ (α) i = -µ (α) j r∞-4 k=0 kH ∞,k q k-1 j (E-1)
Proof. The proof follows from the expression relating the coefficients (ν

(α)
p,k ) and (µ

(α)
p,k ) given by [START_REF] Bergère | Determinantal formulae and loop equations[END_REF][START_REF] Bertola | Hamiltonian structure of rational isomonodromic deformation systems[END_REF][START_REF] Bertola | The partition function of the two-matrix model as an isomonodromic τ function[END_REF][START_REF] Biquard | Wild non-abelian hodge theory on curves[END_REF][START_REF] Boalch | Symplectic manifolds and isomonodromic deformations[END_REF][START_REF] Boalch | Simply-laced isomonodromy systems[END_REF][START_REF] Boalch | Twisted local wild mapping class groups: configuration spaces, fission trees and complex braids[END_REF][START_REF] Chekhov | Hermitean matrix model free energy: Feynman graph technique for all genera[END_REF]. Taking the derivative relatively to q j and using the fact that (ν (α) ∞,k ) -1≤k≤r∞-3 are independent of q j gives:

∀ k ∈ 0, r ∞ -4 : g i=1 (∂ q j µ (α) i )q k i = -kµ (α) j q k-1 j (E-2) Thus r∞-4 k=0 g i=1 H ∞,k q k i ∂ q j µ (α) i = -µ (α) j r∞-4 k=0 kH ∞,k q k-1 j (E-3)
so that the lemma is proved.

We may now provide an alternative expression for L α [p j ]:

Proposition E.1. Let j ∈ 1, g , we have an alternative expression for L α [p j ]:

L α [p j ] = i =j (µ (α) i + µ (α) j )(p i -p j ) (q j -q i ) 2 + 2 (r,s)∈ 1,g 2 r =s (p s -p r )(∂ q j µ (α) r + ∂ q j µ (α) s ) q s -q r -µ (α) j P 2 (q j ) -p j P 1 (q j ) - g r=1 (∂ q j µ (α) r ) P2 (q r ) + p 2 r -P1 (q r )p r + ν (α) ∞,-1 p j + r∞-1 k=1 kc (α) ∞,k q k-1 j (E-4)
Proof. Using Lemma E.1, the expression (5-5) for L[p j ] becomes:

L α [p j ] = i =j (µ (α) i + µ (α) j )(p i -p j ) (q j -q i ) 2 + µ (α) j p j P 1 (q j ) -P 2 (q j ) + ν (α) ∞,-1 p j + r∞-1 k=1 kc (α) ∞,k q k-1 j - g i=1 (∂ q j µ (α) i ) r∞-4 k=0 H ∞,k q k i (E-5)
We now use (5-3) to get

L α [p j ] = i =j (µ (α) i + µ (α) j )(p i -p j ) (q j -q i ) 2 + µ (α) j p j P 1 (q j ) -P 2 (q j ) + ν (α) ∞,-1 p j + r∞-1 k=1 kc (α) ∞,k q k-1 j - g i=1 (∂ q j µ (α) i ) p 2 i -P1 (q i )p i + P2 (q i ) + r =i p r -p i q i -q r = i =j (µ (α) i + µ (α) j )(p i -p j ) (q j -q i ) 2 + µ (α) j p j P 1 (q j ) -P 2 (q j ) + ν (α) ∞,-1 p j + r∞-1 k=1 kc (α) ∞,k q k-1 j - g i=1 (∂ q j µ (α) i ) p 2 i -P1 (q i )p i + P2 (q i ) + g i=1 (∂ q j µ (α) i ) r =i p r -p i q r -q i (E-6)
The last sums may be split into a symmetric and anti-symmetric part:

∂ q j µ (α) i = 1 2 (∂ q j µ (α) i - ∂ q j µ (α) i ) + 1 2 (∂ q j µ (α) i + ∂ q j µ (α) i ). The term involving ∂ q j µ (α) i -∂ q j µ (α) i
is trivially zero because the sum is anti-symmetric so that we end up with

L α [p j ] = i =j (µ (α) i + µ (α) j )(p i -p j ) (q j -q i ) 2 -µ (α) j P 2 (q j ) -p j P 1 (q j ) + ν (α) ∞,-1 p j + r∞-1 k=1 kc (α) ∞,k q k-1 j - g i=1 (∂ q j µ (α) i ) p 2 i -P1 (q i )p i + P2 (q i ) + 2 g i=1 r =i (p r -p i )(∂ q j µ (α) i + ∂ q j µ (α) r ) q r -q i (E-7)
proving Proposition E.1.

E.2 Proof of the Theorem 5.1

We may now proceed to the proof of Theorem 5.1. We recall that the Hamiltonian is given by:

Ham (α) (q, p) = - 2 (i,j)∈ 1,g 2 i =j (µ (α) i + µ (α) j )(p i -p j ) q i -q j - g j=1 (ν (α) ∞,0 p j + ν (α) ∞,-1 q j p j ) + g j=1 µ (α) j p 2 j -P1 (q j )p j + P2 (q j ) - g j=1 r∞-1 k=0 c (α) ∞,k q k j (E-8)
A straightforward computation from [START_REF] Bertola | Hamiltonian structure of rational isomonodromic deformation systems[END_REF][START_REF] Bertola | The partition function of the two-matrix model as an isomonodromic τ function[END_REF] and from the fact that the ν

(α) ∞,k -1≤k≤r∞-3 and c (α) ∞,k 1≤k≤r∞-1
are independent of q j provides -∂Ham (α) (q, p)

∂q j = i∈ 1,g i =j (µ (α) i + µ (α) j )(p i -p j ) (q i -q j ) 2 + 2 (r,s)∈ 1,g 2 r =s (∂ q j µ (α) r + ∂ q j µ (α) s )(p r -p s ) q r -q s + ν (α) ∞,-1 p j - g i=1 ∂ q j (µ (α) i ) p 2 i -P1 (q i )p i + P2 (q i ) -µ (α) j P 2 (q j ) -p j P 1 (q j ) + r∞-1 k=1 kc (α) ∞,k q k-1 j Prop. (E.1) = L α [p j ] (E-9)
Similarly a direct computation using the fact that ν

(α) ∞,k -1≤k≤r∞-3 and c (α) ∞,k 1≤k≤r∞-1 and µ (α) i 1≤i≤g
are independent of p j gives:

∂Ham (α) (q, p) ∂p j = - i∈ 1,g i =j µ (α) i + µ (α) j q j -q i -ν (α) ∞,0 -ν (α) ∞,-1 q j + µ (α) j 2p j -P1 (q j ) (E-10)
which is exactly L α [q j ] given by (5-1).

The last step is to verify that from Propositions 4.2 and 5.1:

g j=1 µ (α) j   p 2 j -P1 (q j )p j + P2 (q j ) + i =j p i -p j q j -q i   = µ (α) 1 , . . . , µ (α) g       p 2 1 -P1 (q 1 )p 1 + P2 (q 1 ) + i =1 p i -p 1 q 1 -q i . . . p 2 g -P1 (q g )p g + P2 (q g ) + i =g p i -pg qg-q i       = µ (α) 1 , . . . , µ (α) g (V ∞ ) t    H ∞,0 . . . H ∞,r∞-4    = ν (α) ∞,1 . . . ν (α) ∞,r∞-3    H ∞,0 . . . H ∞,r∞-4    = r∞-4 k=0 ν (α) ∞,k+1 H ∞,k (E-11) so that (E-8) becomes Ham (α) (q, p) = r∞-4 k=0 ν (α) ∞,k+1 H ∞,k - g j=1 r∞-1 k=0 c (α) ∞,k q k j -ν (α) ∞,0 g j=1 p j -ν (α) ∞,-1 g j=1 q j p j (E-12)

F Proofs of the identities involving elementary symmetric polynomials

Let us prove Lemma 6.1. By definition we have

g j=1 (λ -x j ) = g i=0 (-1) i e i ({x 1 , . . . , x g })λ g-i (F-1)
Taking the derivative relatively to x m provides:

j =m (λ -x j ) = g i=0 (-1) i-1 ∂e i ({x 1 , . . . , x g }) ∂x m λ g-i (F-2)
Thus, (-1) i-1 ∂e i ({x 1 ,...,xg}) ∂xm is the coefficient of order λ g-i of j =m (λ -x j ). However, the last quantity is also given by

j =m (λ -x j ) = g j=1 (λ -x j ) λ -x m = g r=0 (-1) r e r ({x 1 , . . . , x g })λ g-r ∞ s=0 x s m λ -s-1 = g r=0 ∞ s=0 (-1) r e r ({x 1 , . . . , x g })x s m λ g-r-s-1 (F-3)
Identifying the coefficient of order λ g-i provides (-1) i-1 ∂e i ({x 1 , . . . , x g }) ∂x m = i-1 r=0 (-1) r e r ({x 1 , . . . , x g })x i-1-ir m (F-4)

proving the lemma.

Let us now prove Proposition 6.1. From, Lemma 6.1, we get that both polynomials take the same value at x m (the value is ∂e i ({x 1 ,...,xg}) ∂xm

) for any m ∈ 1, g . Since both sides are obviously polynomials of order at most g -1, we immediately get that they are equal.

Finally, let us prove Corollary 6.1. By definition

Q(λ) = - g i=1 p i j =i λ -q j q j -q i = - g i=1 P i   g k=1 ∂e i ({q 1 , . . . , q g }) ∂q k j =k λ -q j q k -q j   Prop 6.1 = - g i=1 P i i-1 j=0 (-1) j Q i-j-1 λ j = g-1 j=0 (-1) j-1   g i=j+1 P i Q i-j-1   λ j (F-5)
ending the proof.

Let us now prove [START_REF] Bertola | The partition function of the two-matrix model as an isomonodromic τ function[END_REF][START_REF] Biquard | Wild non-abelian hodge theory on curves[END_REF][START_REF] Boalch | Symplectic manifolds and isomonodromic deformations[END_REF][START_REF] Boalch | Simply-laced isomonodromy systems[END_REF][START_REF] Boalch | Twisted local wild mapping class groups: configuration spaces, fission trees and complex braids[END_REF][START_REF] Chekhov | Hermitean matrix model free energy: Feynman graph technique for all genera[END_REF]. We have for all j ∈ 1, n :

n-1 i=0 (-1) n-1-i e n-1-i ({x 1 , . . . , x n } \ {x j })λ i = k =i (λ -q k ) = n k=1 (λ -x k ) (λ -x j ) = g m=0 (-1) n-m e n-m ({x 1 , . . . , x n })λ m ∞ p=0 x p j λ -1-p = g m=0 ∞ p=0 (-1) n-m e n-m ({x 1 , . . . , x n })x p j λ m-1-p i=m-1-p = g m=0 m-1 i=-∞ (-1) n-m e n-m ({x 1 , . . . , x n })x m-1-i j λ i = g-1 i=-∞ g m=i+1 (-1) n-m e n-m ({x 1 , . . . , x n })x m-1-i j λ i (F-6)
Identifying the coefficient λ i for i ∈ Z leads to

(-1) n-1-i e n-1-i ({x 1 , . . . , x n } \ {x j }) = g m=i+1 (-1) n-m e n-m ({x 1 , . . . , x n })x m-1-i j , ∀ i ∈ 0, n -1 0 = g m=i+1 (-1) n-m e n-m ({x 1 , . . . , x n })x m-1-i j , ∀ i < 0 (F-7)
which is equivalent to

e n-i ({x 1 , . . . , x n } \ {x j }) = g m=i (-1) i-m e n-m ({x 1 , . . . , x n })x m-i j , ∀ i ∈ 1, n 0 = g m=i (-1) n-m e n-m ({x 1 , . . . , x n })x m-i j , ∀ i ≤ 0 (F-8)
Let us now prove Lemma 6.2. We have the trivial identity for large λ:

∞ M =0 x M i λ -M -1 = 1 λ -x i = n j=1 (λ -x j ) (λ -x i ) n j=1 (λ -x j ) = n m=0 ∞ r=0 ∞ j=1 (-1) n-m e n-m ({x 1 , . . . , x n })h r ({x 1 , . . . , x n })x j-1 i λ m-n-r-j (F-9)
Identifying the coefficient of order λ -M -1 gives

x M i = n m=0 ∞ j=1 (-1) n-m e n-m ({x 1 , . . . , x n })h m-n+M -j+1 ({x 1 , . . . , x n })x j-1 i (F-10)
but only terms with m -n + M -j + 1 ≥ 0, i.e. m ≥ Max(j, j + n -1 -M ) contribute. In particular, since j ≤ m ≤ n, we also have j ≤ n ending the proof.

Let us now prove Proposition 6.2. Let M ≥ 0, the system

     1 x 1 . . . x n-1 1 1 x 2 . . . x n-1 2 . . . . . . . . . 1 x n . . . x n-1 n           C 1 C 2 . . . C n      =      x M 1 x M 2 . . . x M n      (F-11)
is equivalent to say, by inverting the Vandermonde matrix, that for all i ∈ 1, n :

C i = n j=1 (-1) n-i e n-i ({x 1 , . . . , x n } \ {x j }) m =j (x j -x m ) x M j (F-12)
The system is also equivalent to, for all i ∈ 1, n :

n j=1 x j-1 i C j = x M i (F-13)
Lemma 6.2 implies that the last identity is equivalent to, for all j ∈ 1, n

C j = n m=Max(j,j+n-1-M ) (-1) n-m e n-m ({x 1 , . . . , x n })h M +m-j-n+1 ({x 1 , . . . , x n }) (F-14)
Identifying (F-12) and (F-13) finally proves Proposition 6.2.

G Proof of Lemma 6.3

Let J = 0 I -I 0 be the canonical 2g symplectic matrix and define

M = A B C D with                A i,j = ∂Q i ∂ Qj B i,j = ∂Q i ∂ Pj = 0 C i,j = ∂P i ∂ Qj = g r=1 Pr ∂ 2 fr(α -1 Q 1 +β,...,α -1 Qg+β) ∂ Qj ∂Q i -α -1 h (α -1 Q i + β) ∂Q i ∂ Qj D i,j = ∂P i ∂ Pj = ∂f j (α -1 Q 1 +β,...,α -1 Qg+β) ∂Q i = ∂ Qj ∂Q i (G-1)
The change of coordinates is symplectic if and only if M t JM = J which is equivalent to prove that A t D = I and (A t C) is a symmetric matrix. We have for all (i, j) ∈ 1, g 2 :

(A t D) i,j = g k=1 A k,i D k,j = g k=1 ∂Q k ∂ Qi ∂ Qj ∂Q k = ∂ Qj ∂ Qi = δ i,j (G-2)
and

(A t C) i,j = g k=1 A k,i C k,j = g r=1 Pr g k=1 ∂Q k ∂ Qi ∂ 2 f r (α -1 Q 1 + β, . . . , α -1 Q g + β) ∂ Qj ∂Q k -α -1 g k=1 ∂Q k ∂ Qi ∂Q k ∂ Qj h (α -1 Q k + β) = g r=1 Pr ∂ 2 f r (α -1 Q 1 + β, . . . , α -1 Q g + β) ∂ Qi ∂ Qj -α -1 g k=1 ∂Q k ∂ Qi ∂Q k ∂ Qj h (α -1 Q k + β) (A t C) j,i = g k=1 A k,j C k,i = g r=1 Pr g k=1 ∂Q k ∂ Qj ∂ 2 f r (α -1 Q 1 + β, . . . , α -1 Q g + β) ∂ Qi ∂Q k -α -1 g k=1 ∂Q k ∂ Qj ∂Q k ∂ Qi h (α -1 Q k + β) = g r=1 Pr ∂ 2 f r (α -1 Q 1 + β, . . . , α -1 Q g + β) ∂ Qj ∂ Qi -α -1 g k=1 ∂Q k ∂ Qj ∂Q k ∂ Qi h (α -1 Q k + β) (G-3)
so that (A t B) i,j = (A t B) j,i proving the lemma.

H Proof of Theorem 6.1

The proof is based on the computation of each term appearing in the Hamiltonian. Let us compute for i ∈ 1, g :

A i := g j=1 r =j
(-1) g-i e g-i ({q 1 , . . . , q g } \ {q j })

a =j (q j -q a ) p r -p j q j -q r = -g j=1 r =j

(-1) g-i e g-i ({q 1 , . . . , q g } \ {q j })

a =j (q j -q a ) p r -p j q r -q j (6.1)

= -g j=1 r =j

(-1) g-i e g-i ({q 1 , . . . , q g } \ {q j })

a =j (q j -q a ) g k=1 P k k-1 m=0 (-1) m e k-1-m ({q 1 , . . . , q g }) q m r -q m j q r -q j = - g j=1 r =j (-1) g-i e g-i ({q 1 , . . . , q g } \ {q j }) a =j (q j -q a ) g k=1 P k k-1 m=0 (-1) m e k-1-m ({q 1 , . . . , q g }) m-1 s=0 q s r q m-1-s j = - g j=1 g r=1
(-1) g-i e g-i ({q 1 , . . . , q g } \ {q j })

a =j (q j -q a ) g k=1 P k k-1 m=0
(-1) m e k-1-m ({q 1 , . . . , q g })

m-1 s=0 q s r q m-1-s j + g j=1 (-1) g-i e g-i ({q 1 , . . . , q g } \ {q j })

a =j (q j -q a ) g k=1 P k k-1 m=0 (-1) m e k-1-m ({q 1 , . . . , q g }) m-1 s=0 q m-1 j = - g k=1 k-1 m=0 m-1 s=0 (-1) m P k e k-1-m ({q 1 , . . . , q g }) g r=1 q s r g j=1 (-1) g-i e g-i ({q 1 , . . . , q g } \ {q j }) a =j (q j -q a ) q m-1-s j + g k=1 k-1 m=0 m-1 s=0 (-1) m P k e k-1-m ({q 1 , . . . , q g }) g j=1 (-1) g-i e g-i ({q 1 , . . . , q g } \ {q j }) a =j (q j -q a ) q m-1 j (6.2) = - g k=1 k-1 m=0 m-1 s=0 (-1) m P k e k-1-m ({q 1 , . . . , q g }) g r=1 q s r δ i,m-s + g k=1 k-1 m=0 m-1 s=0 (-1) m P k e k-1-m ({q 1 , . . . , q g })δ i,m = - g k=1 k-1 m=0 m-1 s=0 (-1) m P k e k-1-m ({q 1 , . . . , q g })S s ({q 1 , . . . , q g })δ s,m-i + g k=1 k-1 m=0 (-1) m mP k e k-1-m ({q 1 , . . . , q g })δ i,m = - g k=i+1 k-1 m=i (-1) m P k e k-1-m ({q 1 , . . . , q g })S m-i ({q 1 , . . . , q g }) + g k=1 δ k-1≥i (-1) i iP k e k-1-i ({q 1 , . . . , q g }) = - g k=i+1 k-1 m=i
(-1) m P k e k-1-m ({q 1 , . . . , q g })S m-i ({q 1 , . . . , q g })

+ g k=i+1 (-1) i iP k e k-1-i ({q 1 , . . . , q g }) = - g k=i+1 k-1 m=i+1 (-1) m P k e k-1-m ({q 1 , . . . , q g })S m-i ({q 1 , . . . , q g }) - g k=i+1 (-1) i gP k e k-1-i ({q 1 , . . . , q g }) + g k=i+1 (-1) i iP k e k-1-i ({q 1 , . . . , q g }) = - g k=i+1 k-1 m=i+1 (-1) m P k Q k-1-m S m-i ({q 1 , . . . , q g }) - g k=i+1 (-1) i (g -i)P k Q k-1-i (H-1) B i := - g j=1 (-1) g-i e g-i ({q 1 , . . . , q g } \ {q j }) a =j (q j -q a ) P1 (q j )p j = g j=1 g+1 s=0 g k=1 P k k-1 r=0
(-1) r e k-1-r ({q 1 , . . . , q g }) (-1) g-i e g-i ({q 1 , . . . , q g } \ {q j })

a =j (q j -q a ) t ∞,2s+2 q s+r j = g k=1 k-1 r=0 g+1 s=0
(-1) r P k e k-1-r ({q 1 , . . . , q g })t ∞,2s+2 g j=1 (-1) g-i e g-i ({q 1 , . . . , q g } \ {q j })

a =j (q j -q a ) q s+r j = g k=1 k-1 r=0 g-1-r s=0 (-1) r P k e k-1-r ({q 1 , . . . , q g })t ∞,2s+2 g j=1 (-1) g-i e g-i ({q 1 , . . . , q g } \ {q j }) a =j (q j -q a ) q s+r j + g k=1 k-1 r=0 g+1 s=g-r (-1) r P k e k-1-r ({q 1 , . . . , q g })t ∞,2s+2 g j=1 (-1) g-i e g-i ({q 1 , . . . , q g } \ {q j }) a =j (q j -q a ) q s+r j (6.2) = g k=1 k-1 r=0 g-1-r s=0 (-1) r P k Q k-1-r ({q 1 , . . . , q g })t ∞,2s+2 δ s+r+1,i + g k=1 k-1 r=0 g+1 s=g-r (-1) r P k Q k-1-r t ∞,2s+2
g m=i (-1) g-m Q g-m h r+s+m-i-g+1 ({q 1 , . . . , q g }) (-1) g+r-m t ∞,2s+2 P k Q k-1-r Q g-m h r+s+m-i-g+1 ({q 1 , . . . , q g }) (H-2)

C i := g j=1
(-1) g-i e g-i ({q 1 , . . . , q g } \ {q j }) a =j (q j -q a ) P2 (q j ) = g j=1 2r∞-4 r=g P

∞,r (-1) g-i e g-i ({q 1 , . . . , q g } \ {q j })

a =j (q j -q a ) q r j (6. P

k 1 P k 2 k 1 -1 r 1 =0 k 2 -1 r 2 =0
(-1) r 1 +r 2 Q k 1 -1-r 1 Q k 2 -1-r 2 (-1) g-i eg-i({q1, . . . , qg} \ {qj})

a =j (qj -qa) q r 1 +r 2 j = g k 1 =1 g k 2 =1 P k 1 P k 2 k 1 -1 r 1 =0 k 2 -1 r 2 =0 (-1) r 1 +r 2 Q k 1 -1-r 1 Q k 2 -1-r 2 g j=1
(-1) g-i eg-i({q1, . . . , qg} \ {qj}) a =j (qj -qa) q r 1 +r 2 j δ r 1 +r 2 ≤g-1

+ g k 1 =1 g k 2 =1 P k 1 P k 2 k 1 -1 r 1 =0 k 2 -1 r 2 =0 (-1) r 1 +r 2 Q k 1 -1-r 1 Q k 2 -1-r 2 g j=1
(-1) g-i eg-i({q1, . . . , qg} \ {qj}) a =j (qj -qa) q r 1 +r 2 j δ r 1 +r 2 ≥g (6.2)

= g k 1 =1 g k 2 =1 P k 1 P k 2 k 1 -1 r 1 =0 k 2 -1 r 2 =0
(-1) r 1 +r 2 Q k 1 -1-r 1 Q k 2 -1-r 2 δi,r 1 +r 2 +1δr 1 +r 2 ≤g-1

+ g k 1 =1 g k 2 =1 P k 1 P k 2 k 1 -1 r 1 =0 k 2 -1 r 2 =0
(-1) r 1 +r 2 Q k 1 -1-r 1 Q k 2 -1-r 2 δ r 1 +r 2 ≥g g m=i (-1) g-m Qg-mhr 1 +r 2 +m-i-g+1({q1, . . . , qg})

= (-1) i-1 g k 1 =1 g k 2 =1
P k 1 P k 2

Min(k 1 -1,i-1)

r 1 =Max(0,i-k 2 ) Q k 1 -1-r 1 Q k 2 -i+r 1 + g k 1 =1 g k 2 =1 P k 1 P k 2 0≤r 1 ≤k 1 -1 0≤r 2 ≤k 2 -1 r 1 +r 2 ≥g (-1) r 1 +r 2 Q k 1 -1-r 1 Q k 2 -1-r 2 g m=i
(-1) g-m Qg-mhr 1 +r 2 +m-i-g+1({q1, . . . , qg}) (H-4)

We also have: (-1) j P m e m-1-j ({q 1 , . . . , q g })q j i = - (-1) j e m-1-j ({q 1 , . . . , q g })S j ({q 1 , . . . , q g })

I Proof of Proposition 6.4 (-1) g-j-m ν (α) ∞,m Q g-j-m   λ j (I-4)

Let us now compute [ Ãα (λ)] 1,1 . We have from Proposition 2.2:

[ Ãα (λ)] 1,1 = [A α (λ)] 1,1 -   Q(λ) + 1 2 t ∞,2r∞-2 λ + g 0 g j=1 (λ -q j )   [A α (λ)] 1,2 (I-5)
Since we know that [ Ãα (λ)] 1,1 is a polynomial in λ and using (4-18), we get Note that the polynomial part of -∂ λ [A α (λ)] 1,2 + [A α (λ)] 1,2 g j=1 λ-q j Q(λ) P (λ) + 1 2 t ∞,2r∞-2 λ + g 0 is given by (From corollary 6.1, we have Q(λ) = (-1) g P g λ g-1 + O(λ g-2 ))

- 2 ν (α) ∞,-1 t ∞,2r∞-2 λ 2 -ν (α) ∞,-1 g 0 λ -ν (α) ∞,-1 (-1) g P g - 2 gt ∞,2r∞-2 ν (α) ∞,-1 λ -gg 0 ν (α) ∞,-1 + 2 t ∞,2r∞-2 ν (α) ∞,-1 Q 1 - 2 gt ∞,2r∞-2 ν (α) ∞,0 (I-19)
Keeping only the polynomial part of (I-18) leads to (-1) g-j Q g-j ν

[ Ãα (λ)] 2,1 = - 2 ν (α) ∞,-1 t ∞,2r∞-2 λ 2 -ν (α) ∞,-1 g 0 λ -ν (α) ∞,-1 (-1) g P g - 2 gt ∞,2r∞-2 ν (α) ∞,-1 λ -gg 0 ν (α) ∞,-1 + 2 t ∞,2r∞-2 ν (α) ∞,-1 Q 1 - 2 gt ∞,2r∞-2 ν (α) ∞,0 -[A α (λ)]
+ 2 α ∞,2r∞-2 + 2 α ∞,2r∞-4 + 2 α ∞,2r∞-2 Q 1 = - 2 ν (α) ∞,-1 t ∞,2r∞-2 λ 2 -ν (α) ∞,-1 g 0 λ -ν (α) ∞,-1 (-1) g P g - 2 gt ∞,2r∞-2 ν (α) ∞,-1 λ -gg 0 ν (α) ∞,-1 + 2 t ∞,2r∞-2 ν (α) ∞,-1 Q 1 - 2 gt ∞,2r∞-2 ν (α) ∞,0 + (r ∞ -1)c (α) ∞,r∞-1 λ + (r ∞ -1)c (α) ∞,r∞-1 Q 1 + (r ∞ -2)c (α) ∞,r∞-2 + ν (α) ∞,-1 H ∞,r∞-4 + 2 α ∞,2r∞-2 + 2 α ∞,2r∞-4 + 2 α ∞,2r∞-2 Q 1 + [A α (λ)] 1,2 -Q(λ) P1 (λ) -P2 (λ) -Q(λ) 2 P (λ) - 1 2 t ∞,2r∞-2 λ + g 0 P1 (λ) -(t ∞,2r∞-2 λ + 2g 0 ) Q(λ) - 1 2 t ∞,2r∞-2 λ + g 0 2 P (λ) = - 2 ν (α) ∞,-1 t ∞,2r∞-2 λ 2 -ν (α) ∞,-1 g 0 λ -ν (α) ∞,-1 (-1) g P g - 2 gt ∞,2r∞-2 ν (α) ∞,-1 λ -gg 0 ν (α) ∞,-1 + 2 t ∞,2r∞-2 ν (α) ∞,-1 Q 1 - 2 gt ∞,2r∞-2 ν (α) ∞,0 + (r ∞ -1)c (α) ∞,r∞-1 λ + (r ∞ -1)c (α) ∞,r∞-1 Q 1 + (r ∞ -2)c
(α) ∞,j-i λ i (I-20)

where we have ν We recognize that the r.h.s. is -1 2k times the (r ∞ -k) th column of the matrix on the l.h.s. so that we immediately get c We recognize the (k + 2) th column of M ∞ . Hence, we conclude that (ν

(u k ) -1 , . . . , ν (u k )
r∞-3 ) t = e k+2 , i.e. ν (u k ) j = δ j,k for all j ∈ -1, r ∞ -3 .

1 4 2 2 . 2 2 . 3 2 . 4 2 . 5

 422232425 Introduction and summary of the results Twisted meromorphic connections at infinity 5 2.1 Twisted meromorphic connections and irregular times . . . . . . . . . . . . . . . 5 Choice of representative normalized at infinity . . . . . . . . . . . . . . . . . . . . 7 Darboux coordinates . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8 Scalar differential equation and companion gauge . . . . . . . . . . . . . . . . . . 9 Introduction of a scaling parameter . . . . . . . . . . . . . . . . . . . . . . . . . 10 2.6 Explicit expressions of the gauge transformation . . . . . . . . . . . . . . . . . . 10 2.7 Wronksians and asymptotics of the wave functions . . . . . . . . . . . . . . . . . 12 2.8 Explicit expression for the Lax matrix L . . . . . . . . . . . . . . . . . . . . . . . 13

Remark 2 . 2 .

 22 Since Tr L ∞ = -r∞-1 k=1 t ∞,2k z 2k-2 + O(1), the relation Tr L = Tr L ∞ + O(1) shall also provide the diagonal coefficients in

Theorem 5 . 1 (

 51 Hamiltonian evolution). Defining Ham

Proof.

  For any isomonodromic deformation τ we have ∂ τ [ P1 ] = 0 because the coefficients of P1 are trivial times. From the expression of the Wronskians in Definition 2.3, we get that Tr Ľ = Tr L = P1 (λ). Thus, we get that ∂ τ [Tr Ľ] = ∂ τ [Tr L] = 0. The compatibility equation (4-4)

- 27 ) 2 τ 1 α τ 2 )

 27212 This leads to ∂ τ j [Tr Ǎα τ i ] = ∂ τ i [Tr Ǎα τ j ]. It is obvious that the additional gauge transformation Ψn,1 = Ǧ(1) Ψ with Ǧ(1) = exp -1 Tr( Ǎα τ 1 )(s)ds I 2 defines a gauge in which the corresponding Ǎ(1) α τ 1 is traceless. In this new gauge, (7-27) implies that ∂ τ 1 [Tr Ǎ(1) α τ i ] = 0 for all i ≥ 2. In particular, a new gauge transformation Ψn,2 = Ǧ(2(s)ds I 2 does not change the value of Ǎ(1)

(- 1 )

 1 r t ∞,2i-2r P k Q k-1-r

  g-m P (2) ∞,r e g-m ({q 1 , . . . , q g })h r+m-i-g+1 ({q 1 , . . . , q g }) g-m P (2) ∞,r Q g-m h r+m-i-g+1 ({q 1 , . . . , q g }) g-i eg-i({q1, . . . , qg} \ {qj})

(- 1 )

 1 j P m e m-1-j ({q 1 , . . . , q g })

From Proposition 2 . 2 1 ν

 221 we have[ Ãα (λ)] 1,2 = [A α (λ)] λ -m + O(λ -g-1 ) (I-2)so that we get from (F-14):[ Ãα (λ)] 1,2 = g m=-(α) ∞,m λ -m + O(λ -m-1 ) g r=0 (-1) g-r Q g-r λ r g-r ν (α) ∞,m Q g-r λ r-m + O(λ -1 ) g-r ν (α) ∞,m Q g-r λ r-m + O(λ -1 ) g-j-m ν (α) ∞,m Q g-j-m   λ j + O(λ -1 ) (I-3)Since we know that [ Ãα (λ)] 1,2 is a polynomial in λ, we end up with [ Ãα (λ)] 1,2 =

2 Cor. 6 . 1 = 1 ν- 1 ν(- 1 )(I- 6 ) 2 s=0L=

 26111162 g-s Qg-sλ s [Aα(λ)]1,(α) ∞,m λ -m + O(λ -g-1 ) -1 2 t∞,2r ∞-2 λ + g0 g s=0 (-1) g-s Qg-sλ s g m=(α) ∞,m λ -m + O(λ -g-1 ) g-s Qg-sν (α) ∞,m λ s-m + O(λ -1 ) g-i-m Qg-i-mν (α) ∞,m λ iLet us now consider [ Ãα (λ)] 2,2 . The compatibility equation readsL α [ L] = ∂ λ Ãα + Ãα , L (I-7)Taking the trace on both sides leads toL α [Tr( L)] = L α [ P1 (λ)] = ∂ λ Tr( Ãα ) (I-8) Since L α [ P1 (λ)] = -r∞α [t ∞,2s+2 ]λ s =λ)] 2,2 = -[ Ãα (λ)] 1,1 -α ∞,2s+2 λ s+1 + c0 = -[ Ãα (λ)] 1,1 -2s λ s + c0(I-10) Let us now compute c0 , we observe from (4-5) thatTrA α (λ) = 2[A α (λ)] 1,1 + ∂ λ [A α (λ)] 1,2 + [A α (λ)] 1,2 L(λ) 2,2 (I-11)Moreover, from the gauge transformation we get:Tr Ãα (λ) = TrA α (λ) + TrA α (λ) + g j=1 L α [q j ] λ -q j (I-12)But since Ãα (λ) is a polynomial in λ, we may discard the last term and thus we end up with Tr Ãα (λ) = TrA α (λ) = 2[A α (λ)]

1 H 2 P 2 P

 122 ∞,r∞-4 -P2 (λ) [A α (λ)] 1,(λ) -[A α (λ)] 1,2 Q(λ) Q(λ) P (λ) + t ∞,2r∞-2 λ + 2g 0 -1 2 t ∞,2r∞-2 λ + g 0 (λ)[A α (λ)] 1,2

∞ 1 - 1 (- 1 ) 1 (- 1 )P

 11111 ,r∞-2 + ν (α) ∞,-1 H ∞,r∞-4 + 2 α ∞,2r∞-2 + 2 α ∞,2r∞-4 + 2 α ∞,2r∞-2 Q j t ∞,2s+2 ν (α) ∞,m h s+j-m-g-i P r Q r-j-1 λ i -j1+j2 ν (α) ∞,m h j1+j2-g-mr1 P r2 Q r1-j1-1 Q r2-j2-1 λ i

J

  Proof of Proposition 7.1Let k ∈ 1, r ∞ -1 and consider w k . It is obvious from (4-8) (whose r.h.s. only implies odd entries of α) that ν(w k ) j = 0 for all j ∈ -1, r ∞ -3 . Consequently, (4-11) implies that µ (w k ) j = 0 for all j ∈ 1, g . For j ∈ 1, r ∞ -1 , the j th line of the r.h.s. of (44r∞-2j-2m-3 4r∞-2j-2m-3 t ∞,2m -α 4r∞-2j-2m-2 4r∞-2j-2m-2 t ∞,2m-1 . For α = w k it reduces only to -12k t ∞,4r∞-3-2j-2k δ j≥r∞-k Thus we get:

(t

  w k ) ∞,j = -1 2k δ j,k for all j ∈ 1, r ∞ -1 . Let us now take k ∈ -1, r ∞ -3 and consider u k . The (k + 2) th column of M ∞ is given by r∞-1 i=k+2 t ∞,2r∞-2i+2k+1 e i .Similarly, the r.h.s. of (4-8) only implies odd indexes of u k and it ∞,2r∞-2i+1+2k e i (J-2)

  ∂ λ [A α (λ)] 1,2 + [A α (λ)] 1,2 + ∂ λ [A α (λ)] 1,1 + [A α (λ)] 1,2 L 2,1(λ)

	+ = -	2  α ∞,2r∞-2 +   2 α ∞,2r∞-4 + α ∞,2r∞-2 Q 1 2  P1 (λ) + g j=1 λ -q j	    Q(λ) P (λ)	+	1 2	t ∞,2r∞-2 λ + g 0
				P (λ)
	-[A α (λ)] 1,2 Q(λ)	Q(λ) P (λ)	+ t ∞,2r∞-2 λ + 2g 0 -	1 2	t ∞,2r∞-2 λ + g 0	2	P (λ)[A α (λ)] 1,2
	+	2	α ∞,2r∞-2 +	2	α ∞,2r∞-4 +	2	α ∞,2r∞-2 Q 1	(I-18)

1,1 + ∂ λ [A α (λ)] 1,2 + [A α (λ)] 1,2 L(λ) 2,2 (I-13)

For clarity in the exposition, results of[START_REF] Marchal | Hamiltonian representation of isomonodromic deformations of general rational connections on[END_REF] are restricted to the stronger assumption that the matrices defining the Lax matrix are assumed to have distinct eigenvalues at all orders, but as remarked in the paper, the results may easily be generalized to the assumption that only leading orders at each pole are assumed to be diagonalizable after a suitable restriction of the deformation space.

Acknowledgements

The authors would like to thank deeply N. Orantin for suggesting many improvements of the present paper. The authors would also like to thank G. Rembado, J. Douçot for fruitful discussions.

(-1) j e k-j ({q 1 , . . . , q g })S j ({q 1 , . . . , q g })

(g -k)Q k P k+1 (H-5) and g i=1

q i P m ∂e m ({q 1 , . . . , q g }) ∂q i Lemma (6.1)

(-1) j P m e m-1-j ({q 1 , . . . , q g })q j+1 i = -

(-1) j P m e m-1-j ({q 1 , . . . , q g })

(-1) j e m-1-j ({q 1 , . . . , q g })S j+1 ({q 1 , . . . , q g })

(-1) j e k-j ({q 1 , . . . , q g })S j+1 ({q 1 , . . . , q g })

(-1) i-1 e k+1-i ({q 1 , . . . , q g })S i ({q 1 , . . . , q g })

(k + 1)P k+1 e k+1 ({q 1 , . . . , q g })

where we have used a modified version of [START_REF] Bertola | The partition function of the two-matrix model as an isomonodromic τ function[END_REF][START_REF] Biquard | Wild non-abelian hodge theory on curves[END_REF]. Indeed, [START_REF] Bertola | The partition function of the two-matrix model as an isomonodromic τ function[END_REF][START_REF] Biquard | Wild non-abelian hodge theory on curves[END_REF] implies that for any k ∈ 0, g -1 :

k+1 i=0 (-1) i e k+1-i ({q 1 , . . . , q g })S i ({q 1 , . . . , q g }) = (n -k -1)e k+1 ({q 1 , . . . , q g }) (H-7)

Isolating the term i = 0 and reminding that S 0 ({q 1 , . . . , q g }) = g gives:

-(g -k -1)e k+1 ({q 1 , . . . , q g }) = k+1 i=0 (-1) i-1 e k+1-i ({q 1 , . . . , q g })S i ({q 1 , . . . , q g })

(-1) i-1 e k+1-i ({q 1 , . . . , q g })S i ({q 1 , . . . , q g }) -ge k+1 ({q 1 , . . . , q g })

(-1) i-1 e k+1-i ({q 1 , . . . , q g })S i ({q 1 , . . . , q g }) = (k + 1)e k+1 ({q 1 , . . . , q g }) (H-9)

The coefficient c0 is thus the coefficients of λ 0 in the r.h.s.

where ν

Let us now turn to [ Ãα (λ)] 2,1 . From the gauge transformation, we have denoting

The last quantity gives

q j . Since we know that [ Ãα (λ)] 2,1 is a polynomial in λ, we get:

Replacing [A α (λ)] 1,1 and [A α (λ)] 2,1 using (4-5) gives:

Let us now take i ∈ 1, r ∞ -1 and compute the i th line of the r.h.s. of [START_REF] Bergère | Determinantal formulae and loop equations[END_REF][START_REF] Bertola | Hamiltonian structure of rational isomonodromic deformation systems[END_REF][START_REF] Bertola | The partition function of the two-matrix model as an isomonodromic τ function[END_REF][START_REF] Biquard | Wild non-abelian hodge theory on curves[END_REF][START_REF] Boalch | Symplectic manifolds and isomonodromic deformations[END_REF][START_REF] Boalch | Simply-laced isomonodromy systems[END_REF][START_REF] Boalch | Twisted local wild mapping class groups: configuration spaces, fission trees and complex braids[END_REF][START_REF] Chekhov | Hermitean matrix model free energy: Feynman graph technique for all genera[END_REF][START_REF] Chekhov | Hermitian matrix model free energy: Feynman graph technique for all genera[END_REF][START_REF] Chekhov | Matrix eigenvalue model: Feynman graph technique for all genera[END_REF][START_REF] Dumitrescu | Quantization of spectral curves for meromorphic Higgs bundles through topological recursion[END_REF][START_REF] Eynard | Quantization of classical spectral curves via topological recursion[END_REF][START_REF] Eynard | Invariants of algebraic curves and topological expansion[END_REF][START_REF] Eynard | Topological recursion in random matrices and enumerative geometry[END_REF][START_REF] Fuchs | Sur quelques équations différentielles linéaires du second ordre[END_REF][START_REF] Gaiur | Isomonodromic deformations: Confluence, reduction & quantisation[END_REF][START_REF] Gambier | Sur les équations différentielles du second ordre et du premier degré dont l'intégrale générale est à points critiques fixes[END_REF]. It is given by:

Thus, the r.h.s. of (4-20) is null for α = u k so that since M ∞ is invertible, we get c

Similarly, we also have from Proposition 7.1 that c

Let us now consider L u -1 . From Proposition 7.1 we have ν

Similarly, since from Proposition 7.1 we have c

Let us now consider L u 0 . From Proposition 7.1 we have ν

Similarly, since from Proposition 7.1 we have c

we first observe that we may rewrite for all k ∈ 1, g :

We now insert the ansatz

which gives:

Let us now take

and prove that for all k ∈ 1, r ∞ -3 :

(L-6) Indeed we have:

Thus, under the choice (L-5), equation (L-4) simplifies into

We now rewrite:

and we take . Thus, the sum from r = 0 to k -p in (L-8) reduces to

. We obtain:

Thus, we conclude that for all j ∈ 1, r ∞ -3 :

M Proof of Proposition 7.5

= 0 from Proposition 7.4, we immediately get L w k [q j ] = 0. Moreover from Proposition 3-3, we observe that:

so that Theorem 7.1 provides for all j ∈ 1, g :

Let us now consider L u -1 . From Proposition 7.4, we have

We also have from Theorem 7.1 L u -1 [q j ] = -q j and L u -1 [p j ] = p j . Finally we observe that

Thus, we get for all j ∈ 1, g :

Let us now consider L u 0 . From Proposition 7.4, we have

We also have from Theorem 7.1 L u 0 [q j ] = -and L u 0 [p j ] = 0. Finally we observe that

Thus, we get for all j ∈ 1, g :

N Proof of Theorem 7.2

Let us first observe that a function

Hence, the function f may only depend on odd irregular times:

Let us now consider L u -1 [g] = 0. It is equivalent to

whose solutions are arbitrary functions of

In other words:

Let us now translate this result to L u 0 [h] = 0. We find

We proceed using the following lemma.

Lemma N.1. The general solutions of the differential equation

are arbitrary functions of

where k ∈ 1, r ∞ -3 . with k ∈ 1, r ∞ -3 .