The locus of the representation of logarithmic connections by Fuchsian equations

The generic element of the moduli space of logarithmic connections with parabolic points on holomorphic vector bundle over the Riemann sphere can be represented by a Fuchsian equation with some singularities and some apparent singularities. We analyze the case of rank $3$ vector bundle which leads to third order Fuchsian equation. We find coordinates on an open subset of the moduli space and we construct a non-trivial part of the moduli space by blowing up along a variety in a special case.


Introduction
Let Σ be a compact Riemann surface with n + 1 distinct fixed points P = {t 0 , . . . , t n }. We will call these points parabolic points. Let I be the set of indices of t i ∈ P . Fix a rank m smooth vector bundle V over Σ and consider a holomorphic vector bundle E with underlying topological vector bundle V. Furthermore consider the logarithmic connections on E. Definition 1.1. Let Σ be a compact Riemann surface with a rank m holomorphic vector bundle E. Let U Σ be an open chart in Σ, let f : U Σ → C be any holomorphic function and let σ ∈ Γ(U Σ , E). A logarithmic connection D over E with parabolic points t i ∈ P is a C-linear map D : E → E ⊗ Ω 1 Σ (log(P )) where Ω 1 Σ (log(P )) is the sheaf of meromorphic 1-forms with at most first order poles at P over Σ which satisfies the Leibniz rule D(σf ) = (Dσ)f + σ ⊗ df.
We introduce a local chart U in Σ around a point t / ∈ P and let z be a coordinate on U . A logarithmic connection on U has the shape D = d + A(z)dz, where A(z) is a matrix-valued holomorphic function. Let U i be a local chart in Σ around a point t i ∈ P and let z i be a coordinate on U i such that t i = {z i = 0}. A logarithmic connection on U i has the shape D = d +Ã (zi) zi dz i , whereÃ(z i ) is a matrix-valued holomorphic function. We refer to theÃ(t i )'s as residues and we fix their adjoint-orbits. We denote the eigenvalues ofÃ(t i ) by ρ i,k , where k = 1, . . . , m. The eigenvalues ρ i,k are called generic if we choose az integer 1 ≤ l ≤ m and any K i ⊆ {1, . . . , m} with |K i | = l ∀i ∈ I then i∈I k∈Ki ρ i,k / ∈ Z.
Consider the pairs (E, D) where D is a logarithmic connection with fixed {ρ i,k } n,m i=0,k=1 . Assumption 1. 3. We suppose that the eigenvalues ρ i,k are generic.
This implies by the residue theorem all connections are irreducible. For the construction of moduli spaces one needs a technical condition called stability. It follows immediately from definitions that irreducible connections are stable. Hence under the Assumption 1.3 all pairs (E, D) with the given eigenvalues ρ i,k are stable. Thus we refrain ourselves from spelling out the definition of stability. Remark 1.4. If we fix the eigenvalues ρ i,k then the rank and degree of underlying vector bundle V will be fixed by the residue theorem. (n − 1) − m + 1 (Remark 6.33. in [8]).
Moreover, there exists a holomorphic symplectic structure on M [6]. The topology on the moduli space M refers to the Zariski topology. It is known that the above connections can be represented by a Fuchsian equation with singular locus in P and with some apparent singularities outside P [8].
Definition 1.5. A Fuchsian equation is a linear homogeneous ordinary differential equation of order m in a complex domain U Σ . Its local form is the following where the B k (z)'s are functions with a pole of order k at some points u l ∈ U Σ (k = 1, . . . , m, l ∈ L with set of indices L).
Definition 1.6. A point q ∈ CP 1 is an apparent singularity of a differential equation (2) if any of the (meromorphic) coefficient of Equation (2) has a pole but the fundamental system of solutions is analytic at this point.
We will label parabolic points by i and apparent singularities by j throughout the paper. Moreover, we will label by i or j the expressions related to t i or q j . Denote the set of apparent singularities by Q. Dubrovin and Mazzocco dealt with Fuchsian equations of arbitrary order in [2]. They introduced apparent singularities and some auxiliary parameters as Darboux coordinates on the symplectic space of Fuchsian system. Szabó specified the number of apparent singularities in [9] and he completely described the case of rank-2 vector bundle, which leads to a second order Fuchsian equation [10]. The present paper provides an alternative approach to determine Darboux coordinates. Szabó also describes the cases q j1 = q j2 and t i = q j for some indices by a blow up method.
The third order case is significantly more involved because the coefficient matrix of system of equations is not block diagonal, hence new techniques are needed for the analysis. The case of 4 parabolic points (Theorem 1.7) is the smallest nontrivial example, where there exists a non-empty set over which the natural Darboux coordinates of the moduli space do not extend. This case differs from the previously studied cases (i.e. q j1 = q j2 and t i = q j cases) in [10]. Conjecture 1.9 describes how to construct a non-trivial part of the moduli space by blowing up and how to choose coordinates.
Iena and Leytem [5] point out a similar phenomenon in the Simpson moduli space of semi-stable sheaves. Namely, they blow up the moduli space along a singular subvariety and they get a closed subset of the moduli space of non-singular sheaves.
1.1. The main results. Before the main result we introduce some notations: let ∆ be the subset where q j1 = q j2 or t i = q j3 for any q j1 , q j2 , q j3 ∈ Q and t i ∈ P , let N be the value in Equation (1) and let Sym(N ) be the N -th symmetric group.
Corollary 1.8. We get a canonical coordinate system on M 0 with coordinates q j on the base of (C \ P ) N / Sym(N ) and coordinates p j on the fiber over q j .
The following conjecture is an extension of the previous theorem.
Conjecture 1.9. Let Σ be the Riemann sphere with n parabolic points, fix generic eigenvalues {ρ i,k } n,3 i=0,k=1 and consider the moduli space of logarithmic connections M with eigenvalues ρ i,k . Let V be the subvariety from Theorem 1.7. Then there exists a Sym(N ) invariant affine subvarietyV = ∅ in T * (C \ P ) N , which intersects V generically transversely and there exists an open subset of V ∩V (denoted by (V ∩V ) 0 ) such that the following blow-up is a subset of the moduli space: We will show a general numerical example in the case n = 3 where this conjecture holds. Another natural question is what is the situation with higher rank (m > 3) vector bundle. The conjecture is that there exists two discriminant varieties V and V such that a result similar to the above holds. The computations in this article were made using Wolfram Mathematica 11.0.

1.2.
Outline of the paper. The next section contains the initial statement and equation which is the starting point of our analysis of moduli space. The section also describes some well-known tools which are useful for the computations. In Section 3, we write the equations explicitly in the case of rank-3 vector bundles and we collect the conditions which guarantee the appearance of parabolic points and apparent singularities. In Section 4, we analyze the case n = 3. We construct the discriminant variety V and make a numerical example. After that we show a blow up procedure in support of Conjecture 1.9. In Section 5, we compute the variety V in the case of arbitrary many parabolic points and we prove Theorem 1.7.
Acknowledgments: The author acknowledges support of NKFIH through thé Elvonal (Frontier) program, with grant KKP126683. The author wants to thank Szilárd Szabó who offered the research topic, gave some useful suggestions and paid attention to the writing of the paper. Thanks to András Stipsicz for many useful comments and encouragement.

Preparatory material
The starting point of this paper is a theorem of Katz [7] extended by Dubrovin and Mazzocco in [2]. Here we summarize and apply their results which are relevant to us.
Let Σ be a compact Riemann surface with t 0 , . . . , t n ∈ Σ, let E be a holomorphic bundle on Σ, let U Σ be an open chart in Σ and let D be a logarithmic connection over E with logarithmic points t 0 , . . . , t n . Consider the sheaf S on Σ \ {t 0 , . . . , t n } given by L(U Σ ) := {e ∈ Γ(U Σ , E)|D(e) = 0}. This gives a functor Ψ from the category of regular connections on Σ to the category of local systems on Σ, i. e. the locally constant sheaves of C-vector spaces.
Let τ : [0, 1] → Σ be a path in Σ. Let L be a local system on Σ. Then τ * L is a local system on [0, 1]. The triviality of this local system yields an isomorphism (τ * L) 0 → (τ * L) 1 . The two stalks (τ * L) 0 , (τ * L) 1 are canonically identified with L τ (0) and L τ (1) . Thus we find an isomorphism L τ (0) → L τ (1) induced by τ . Let b be a base point for Σ and let π 1 denote the fundamental group of Σ \ {t 0 , . . . , t n } with respect to this base point. Let F denote the stalk L b . Then for any closed path τ through b we find an automorphism of F . In this way we have associated to L a representation ρ L : π 1 → GL(F ) of the fundamental group. More generally, we get a functor ρ which associates a representation of the fundamental group of Σ \ {t 0 , . . . , t n } to a local system on Σ.
gives a functor which associates a so-called monodromy representation to a regular connection. [8] Definition 2.2. Let a (Fuchsian) differential equation over an open, connected set U in the complex plane. Let π 1 (U ) the fundamental group of U . The monodromy representation of the differential equation is the linear representation of π 1 (U ). Theorem 2.4. Fix a Riemann surface Σ, a rank m holomorphic vector bundle E and set of points P with fixed eigenvalues ρ i,k of residues of logarithmic connection at t i ∈ P (i = 0, . . . , n, k = 1, . . . , m). Then every element of the above moduli space M can be represented by an order m Fuchsian equation (2) which has singularities at the points t i ∈ P with exponents ρ i,k (mod 2πZ) and has apparent singularities at points of some set {q 1 , . . . , q g } = Q for some value g.
Moreover, this correspondence gives Darboux coordinates (q j , p j ) g j=1 on the moduli space M where p j 's derive from the coefficients of the representing Fuchsian equation.
Katz does not specify the value of g, Dubrovin and Mazzocco specified g = N from 1 but does not describe the variety when the Fuchsian representation is not canonical and depends on some choices. We only use the fact that there exists a Fuchsian representation and Darboux coordinates. We will analyze the moduli space with different techniques.
In the case Σ = CP 1 Szabó specified the value g as the number g = N = 1 2 dim C M in [9]. Moreover, he proved that the eigenvalues ρ i,k of residues and exponents ρ i,k of singularities are equal in generic case (i. e. when q j1 = q j2 ). In the case of second order Fuchsian equation [10] he assigned the auxiliary parameters p j (j = 1, . . . , N ) from the coefficient B 2 (z). We will follow this method in our article and it will turn out that the auxiliary parameters of the third order system also come from B 2 (z).
2.1. The Frobenius method. We will examine Equation (2) in the case Σ = CP 1 and m = 3 with the Frobenius method. For this reason, we provide a short summary about this method. Write the coefficients of Equation (2) as quotients of two polynomials. The roots of the denominator are t i ∈ P and q j ∈ Q for all i and j. Let z be the standard coordinate on the chart of CP 1 away from [1 : 0].
Similar theorem concerns the apparent singularities q j , except that the exponents belonging to q j must be natural numbers. Therefore the exponents belonging to q j are not generic in the sense than ρ i,k 's and we need tools to handle them. The next facts concern the case m = 3 but analogous statements hold for any order of differential equation at any (apparent) singularity.
(2) If ρ i,k1 − ρ i,k2 ∈ Z, ρ i,k1 − ρ i,k3 / ∈ Z and ρ i,k1 ≤ ρ i,k2 for any k 1 , k 2 , k 3 then there exists a fundamental system of solutions of the form for any k 1 , k 2 , k 3 then there exists a fundamental system of solutions of the form (4) If ρ i,k1 = ρ i,k2 for any k 1 , k 2 ∈ {1, 2, 3} then the fundamental system of solutions is the same as in the case ρ i,k1 − ρ i,k2 ∈ Z except that κ i k1,k2 = 0. It is well known that above conditions are not linearly independent, the Fuchs' relation holds for all exponents of all singularities due to residue theorem [1], namely:

Confluent Vandermonde matrices.
A confluent Vandermonde matrix is a certain generalization of Vandermonde matrix. Let n be a natural number, let X = {x 1 , x 2 , . . . , x m } be a set of parameters and let λ be a length m integer partition of n such that λ = (n 1 , n 2 , . . . , n m ) ⊢ n.
Definition 2.9. Define the n × n confluent Vandermonde matrix M corresponding to X and λ in the following way. For all indices i ∈ {1, . . . , m} there is a row of M which contains the consecutive powers of the parameter x i from exponent zero to n − 1. We say this is a Vandermonde-like row: If any of n i > 1 then M contains rows which are 1 st , 2 nd , . . . , (n i − 1) th derivative elementwise of above Vandermonde-like row with respect to x i and divided by 1!, 2!, . . . , (n i − 1)! respectively.
Definition 2.10. In this article we require that the rows follow each other by the index of parameters, within that by order of derivatives. (But we do not require that if i < j then n i ≥ n j .) We will refer to this order as standard sequence.
Proposition 2.11. [4] Let M be the n × n confluent Vandermonde matrix with the set of parameters X (with cardinality m) and λ ⊢ n. Let n i be the cardinality of rows which contain parameter x i (n i − 1 is the maximum order of derivatives with respect to x i 's). The determinant of such a confluent Vandermonde matrix is the following: where we introduce a function ConfVand which gives the determinant of a confluent Vandermonde matrix with parameters x i and order of derivatives n i − 1 [3].

The third order system
First of all, fix the notation in the case of rank-3 vector bundles (i. e. m = 3). The Riemann sphere CP 1 has n + 1 fixed points t i (the parabolic points) with fixed exponents ρ i,k (i = 0, . . . n, k = 1, 2, 3) and N = 3n − 5 points (the apparent singularities of Fuchsian equation) according to Equation (1). They are: For simplicity, choose t 0 at infinity. Now, the differential equation (2) has the following shape: Write the coefficients B i (z) as a quotient of two polynomials as in Equations (3). The denominators of these fractions are powers of the following polynomial: The polynomials in the numerators are: Namely, the differential equation takes the following shape: Equation (9) contains 24n − 33 unknown coefficients.
3.0.1. Conditions from parabolic points. Frobenius method provides a solution of Equation (9) near z = t i as the form w(z) = ∞ l=0 a l (z − t i ) l+ρi and get the indicial equation (i = 0, . . . , n): Similar limits at t 0 = ∞ are given by The indicial equation (10) reads as Denote the (fixed) roots by ρ i,k (k = 1, 2, 3). Vieta's formulas say: Express G ti 0 , H ti 0 , I ti 0 and substitute to Equations (11) and (12). We get the following system of linear equations in coefficients of G, H and I as variables (i = 1, . . . , n): 3.1. Conditions from apparent singularities. We treat the case of apparent singularities similarly. Compute Otherwise, we get a system of linear equations from indicial equations and Vieta's formulas, similarly as in Equation (13) (switch index i to j and t i to q j ). The exponents of apparent singularities are not fixed, but they are natural numbers since Definition 1.6. Hence case of (3) or case of (4) in Fact 2.7 satisfied. We will can choose the values of ρ j,k (j = 1, . . . , 3n − 5; k = 1, 2, 3).
The roots of the indicial polynomial cannot be equal because it leads to a nonzero logarithmic term according to the case of (4) in Fact 2.7.
Definition 3.1. The defect of an apparent singularity q j is The defect measures the complexity of an apparent singularity q j . The defect vanishes at a smooth point.
Indeed the set of exponents with defect zero is ρ j,1 = 0, ρ j,2 = 1 and ρ j,3 = 2 for all j. But then G qj 0 = 0, H qj 0 = 0 and I qj 0 = 0, hence the differential equation has no apparent singularity at q j . Thus we exclude this possibility.
In view of Proposition 2.8 the next simplest possibility is δ = 1.
Assumption 3.2. The defect at all apparent singularities q j is equal to 1 which means: Then G qj 0 = −1, H qj 0 = 0 and I qj 0 = 0 and we get the following system of linear equations for the coefficients (j = 1, . . . , 3n − 5): Conditions for the vanishing of logarithmic terms. The above chosen ρ j,k exponents may cause the appearance of logarithmic terms with coefficients κ j k1,k2 in the solution of the differential equation in (9) according to the case (3) in Fact 2.7. The indices (k 1 , k 2 ) refer to the exponents ρ j,k1 and ρ j,k2 .
For the existence of apparent singularities we need to exclude logarithmic terms; therefore, we must construct conditions for the vanishing of these terms. Note that the method in [Lemma 2.] [10] does not work here.
A fundamental system of solutions of Equation (9) near q j based on the Frobenius method with the above values of exponents are the following: where a 0 = 0, b 1 = 0 and c 3 = 0. Lemma 3.3. The vanishing of logarithmic terms is equivalent to the following equations.
Proof. Write the Wronskian from the three functions of Equations (18) and a function w(z) as a general solution of Equation (9). Convert the Wronskian to the shape of (9) with dividing by coefficient of w ′′′ (z). Note that it is enough to compute terms with non-positive exponents in (z − q j ) in the Wronskian, because the higher order terms do not give contributions to coefficients of lower (z − q j )-powers.
Consider the Laurent series expansions of the coefficients about q j .
In the converse direction we choose κ j k1,k2 for 0 in Equations (18) and after that we make a Wronskian in the previous way. After the conversion and series expansions we get Equations (19).
Take the series expansions by definition of the coefficients of the differential equation (9) and use the known Compare the coefficients from the series expansions of the coefficients of the Wronskian and the coefficients from the series expansions of coefficients of Equation (9), namely Equations (19) and Equations (20). Note that the value H qj 1 depends on b 1 and b 2 only, hence H qj 1 will be the free parameter (possible up to an affine transformation) which is denoted by p j [2]. These will be the remaining N Darboux coordinates on the moduli space next to the q j 's. Elimination of Taylor coefficients from Conditions (19) yield: We neglect terms G qj 1 now, because it will turn out in Subsection 3.2 these are zeros.

3.2.
Computation of G qj 1 . The term G qj 1 appears in some previous equations as (24) and (25a), thus we need to compute this value.
Consider the right hand side of Equation (17a) and introduce the notation The term G qj 1 is the zeroth coefficient of the Laurent series of G(z) ψ(z) about q j . Use Expression (7) and expression η j : Substitute Equation (26) and conclude that we have to differentiate −1 which is zero, hence G The total number of equations for the coefficients of H and I is 24n − 33. This is the same number as the number of coefficients contained in H and I.
We note that the right hand side of this linear system of equations is known, since the parabolic points t i and their exponents ρ i,k are given; the apparent singularities q j and parameters p j are freely chosen.
We need to prove that the linear system of equations in the variables of coefficients of G, H and I has a unique solution, i. e. the coefficient matrix has full rank. First, we study the cases n = 2 and n = 3 because the general case will use the method which will be written in the case n = 3.
Proposition 3.6. The coefficient matrix in question is a block diagonal matrix and the first block which refers to the variables of coefficients of G has full rank.
Proof. The first 4n − 4 rows of the coefficient matrix are independent of other rows, because Equations (14a), (14b) and (17a) do not contain coefficients of H or I. Moreover, the Fuchs relation (Equation (5)) reduces the rank by one, hence we may ignore the first row. The remaining block is a Vandermonde matrix which has nonzero determinant according to Proposition 2.11. According to this proposition the determinant of the block of coefficient matrix which determines the coefficients of G does not vanish and it provides the following (trivial) discriminant variety in T * (C \ P ) 4 : The variety V 0 coincides with the subset ∆ hence we do not deal with the block of coefficient matrix which determines the coefficients of G.

The case n = 3
In this section we will give an argument supporting for Conjecture 1.9. Consider System (T ) from Subsection 3.3 and its coefficient matrix M 1 . The determinant of this matrix will define a discriminant variety V 1 on T * (C \ P ) 4 . The augmented matrix of System (T ) will be M b which will define another discriminant varietyV . Our purpose is to analyze the special case when rk(M 1 ) = rk(M b ) < 20n − 28 = 32 that is the System (T ) is under-determined. We will analyze the case rk(M 1 ) = 31 only. For this purpose, we pick a minor of M 1 which defines a variety W . Finally, we will specify a one-parameter family of solutions of System (T ) in the subvariety (V 1 ∩V ) \ W . 4.1. The case n = 2. Before we analyze the case n = 3, for the sake of completeness, we review the case n = 2 shortly. The parabolic points: P = {t 0 , t 1 , t 2 }, the apparent singularity: Q = {q 1 }, the coefficients of Equation (9): I(z) = I 0 + I 1 z + I 2 z 2 + I 3 z 3 + I 4 z 4 + I 5 z 5 + I 6 z 6 .
The constants, which appear in the coefficient matrix are the following: Use System (T ) and get a coefficient matrix for the variables of coefficients of H and I.
where Λ 1 = 6ω 1 q 1 2 + 4ν 1 q 1 , Λ 2 = 10ω 1 q 1 3 + 10ν 1 q 1 2 , Λ 3 = 15ω 1 q 1 4 + 20ν. Compute its determinant with Laplace expansion along the sixth row. The determinant is −(t 1 − t 2 ) 2 which does not vanish since t i = t j if i = j. Thus the coefficients of G, H and I are determined, the differential equation (9) is also determined and Theorem 1.7 is valid in the case n = 2. The set V is the empty set.

4.2.1.
Searching for confluent Vandermonde matrices. The determinant of M 1 cannot be computed directly even with computers. M 1 is almost a block diagonal matrix with two blocks. It is just almost block diagonal because the Equation (25a) provides for the matrix four rows (the so-called common rows) which have nonzero element in both blocks. These quasi blocks are almost confluent Vandermonde matrices in the sense that they have several rows which together form confluent Vandermonde matrices.
Our aim is to transform M 1 to some block lower/upper triangular matrices where the upper 13 × 13 and lower 19 × 19 blocks are confluent Vandermonde matrices and the remaining 13 × 19 or 19 × 13 blocks are zero matrices.
The first matrix (denoted by M 1 1 ) only differs from M 1 in that the 14th-32nd elements of the 13th row are zeros. The second matrix (say M 2 1 ) is the so-called auxiliary matrix, M 2 1 differs from M 1 in that the 1st-13th elements of the 13th row are zeros. Obviously det(M 1 ) = det(M 1 1 ) + det(M 2 1 ) and the matrix M 1 1 is a block lower triangular with the desired shape, since it has a 19 × 13 zero block. (Note, the proof that M 1 1 is confluent Vandermonde still remains.) R 1 and S 1 denote the diagonal blocks of M 1 1 . In the second step interchange the 13th and the 14th rows in M 2 1 . The determinant will be multiplied by −1. The resulting matrix (denoted by M 2 1 ′ ) has the shape of M 1 except some zero elements in the 14th row. (Note that the 13th row contains powers of q 2 and the 14th row contains powers of q 1 .) Apply the first step to M 2 1 ′ and get two new matrices. One of them is a block lower triangular matrix with diagonal blocks R 2 and S 2 .
Continue the process (second step and first step) until the 13 × 19 or the 19 × 13 blocks are zero matrices. We need two more iterations and denote the diagonal blocks of the created matrices by R 3 , S 3 , R 4 and S 4 . Finally, we get another matrix where the 1st-13th elements of 13th-32nd rows are zeros, thus the determinant of this matrix vanishes.
Hence the determinant of M 1 is the following: where the j also refers to q j which is contained in the 13th row in R j and to the q j which is not contained in the first three rows in S j . The last row of R j has a factor µ j (j = 1, . . . , 4). Denote the matrix without this factor by R ′ j and notice that R ′ j is confluent Vandermonde apart from its first row.
The determinant of S j (j = 1, . . . , 4) does not change if we add a scalar multiple of a row to another row. Based on Equations (25a) and (25b) we can change the matrix S j : where we abuse the notation and also denote the new matrix by S j . The term [k, .] denotes the k-th row of S j . The indexk ∈ {16, . . . , 19} denotes the row which contains the same q l powers as the k-th row. For example, in S 3 we leave k = 1, 2, 3 andk = 16, 17, 19 respectively. Divide S j to two other matrices along the (15 + j)-th row. The first matrix contains the second derivative of I(q j ), and if we factor ω j ν 2 1 ν 2 2 ν 2 3 ν 2 4 /ν 2 j from the rows, then the remaining matrix will be a confluent Vandermonde and will be denoted by S ′ j . The second matrix contains the third derivative of I(q j ). Factor ν 2 1 ν 2 2 ν 2 3 ν 2 4 /ν j and denote the remaining matrix byŜ ′ j . The matrixŜ ′ j is an almost Vandermonde, it misses the second derivative of I(q j ). However, det(Ŝ ′ j ) can be derived from det(S ′ j ) as follows: Hence we showed that M 1 can be decomposed as some block lower triangular matrices with two confluent Vandermonde matrices in its diagonals.

4.2.2.
Matrices and their determinants. The parameters in the rows of R ′ 1 are 1, 1 respectively, where the upper index refers to the order of derivative. This sequence differs from the standard sequence (see Definition 2.10).
Definition 4.1. We introduce a total order (≺) on the set of t i 's and q (k) j 's. We define the order such that the ascending order is just the standard sequence. Let L be a sequence with t i 's and q The number of inversions equals the number of necessary row interchanges which converts the sequence of parameters to the standard sequence. The row interchanges change the determinant by a sign hence it is necessary to calculate the parity of number of inversions.
The number of inversions is 12 in the matrix R ′ 1 , thus apply Equation (6) and get the determinant of R ′ 1 : The determinants of the other R ′ j 's (denoted by r j , j = 2, 3, 4) can be computed similarly but there is a simpler way if we change the indices. Change q j to q 1 and vice versa in Equation (29). We only need to compute how many row interchanges are needed to the conversion between R ′ j to R ′ 1 . The numbers of inversions are even in all R ′ j 's, thus the determinants r 2 , r 3 and r 4 can be calculated with index change from r 1 only. We do not write these terms due to limited space.
The parameters in S ′ 1 are: q 3 , q 4 . The number of inversions is 44. Use Equation (6) again for the determinant of S ′ 1 : For the determinants of other S ′ j 's (j = 2, 3, 4) use the similar method as for R ′ j 's. Namely, interchange parameters q j and q 1 . It is easy to compute the parity of the number of inversions: (−1) 3+(j−2) . Note that the sign (−1) j−1 in Equation (27) cancels the previous sign, hence we can ignore these signs in the determinant of M 1 . Get s j := det S ′ j that is, interchange q 1 and q j in Equation (30). We do not write these terms due to limited space.
Next, produceŝ j with Equation (28). Finally, sum up all r j 's, s j 's andŝ j 's, substitute µ j , ν j and ω j and get the determinant of M 1 .
After long simplification a symmetric polynomial χ 1 can be factored out (namely σ 1 = χ 1 φ 1 ), where The polynomial φ 1 is also symmetric for permutations from Sym(4): where the denominators of the fractions cancel with multiplier factors, hence φ 1 is a polynomial.
The matrix M 1 has full rank if and only if its determinant does not vanish. The polynomial χ 1 is not zero because q j = t i for any i and j. In order to show that the polynomial φ 1 does not vanish, it is enough to compute the coefficient of highest power in q 1 . The highest power is q 6 1 and its coefficient is what is nonzero if the q j 's differ from each other. Hence we proved that M 1 has full rank in the generic case, moreover σ 1 specifies a discriminant variety in T * (C \ P ) 4 (denoted by V 1 ). Consequently, the moduli space M of logarithmic connections of a rank-3 vector bundle over the Riemann sphere with 3 parabolic points can be equipped with Darboux coordinates q j , p j outside of V 1 . This proves Theorem 1.7 in case n = 3.

Further discriminant varieties.
It is well known from the theory of system of linear equations that if the coefficient matrix does not have full rank but this rank equals the rank of augmented matrix then the system of equations is underdetermined and the solution has some free parameters. We want to reproduce this phenomenon related to the moduli space M.
In this subsection we create varieties (denoted by V k , k = 2, . . . , 33) which corresponds to the augmented matrix and we show a numerical example where V 1 = V k and V 1 ∩ V k = ∅ for some index k.
Denote by M b that 32 × 33 matrix which has first column b and the others are the same as M 1 . Denote by M k the 32 × 32 matrix which comes from M b by deleting its k-th column (k = 2, . . . , 33). Let σ k := det M k and let V k be the variety in T * (C \ P ) 4 which is defined by σ k . (Note that the above definition gives back M 1 and σ 1 for k = 1.) It is clear that if one of M k (k ≥ 1) has full rank, then the augmented matrix of System (T ) has full rank. We proceed by method of contraposition. Namely, if the augmented matrix does not have full rank then M k for all 1 ≤ k ≤ 33 does not have full rank. Hence the subvariety ofV := 33 k=2 V k describes the locus where the the moduli space does not admit unique coordinates. The moduli space M can be equipped coordinates on V k \V , but it remains to understand what is happening in the subvariety V 1 ∩V . 4.3.1. Numerical example. We do not compute varieties V k directly due to huge matrices and rational fractions. Thus we make a general numerical example which does not prove the Conjecture 1.9 but supports it. We choose random rational numbers for the parameters (t i , α i , β i , γ i ) and coordinates (q j , p j ). Technically we choose rationals from the interval [1/40, 70] in order for the numbers in calculations be manageable easier.
We can convince ourselves about V 1 = V k by substitution of random numbers to σ 1 and σ k .
First, we will find a point P in V 1 ∩V . Fix all parameters in σ 1 and σ k except for two of them. For the sake of simplicity let the two non-fixed parameters be p 1 and p 2 since deg pj σ 1 = 1 and deg pj σ k = 2 for j = 1, . . . , 4. Let the fixed parameters be random rational numbers such that the factors (q i1 − t i2 ), (q j − t i ) and (q j1 − q j2 ) do not vanish. We note β j = 0 is a possible choice because α j and β j appear only in the factor (α j − β j − 1).
Substitute these numbers to σ 1 and σ k , and solve the equation σ 1 = 0 for the variable p 1 . Substitute the solution to the equation σ k = 0 and solve the resulting quadric equation for the variable p 2 . The random numbers specify some points in V 1 ∩ V k for all k. It turns out there exist three points in V 1 ∩V under the parameters fixed above. One is chosen as P . Geometrically, this means that the p 1 − p 2 coordinate plane intersects V 1 ∩V . In the conjecture we suppose that the point P is a smooth point in the varieties V 1 and V k . We have the freedom to choose a minor, the appropriate choice is an easy computable one. If all varieties, defined by these minors, contain an everywhere dense subset of V 1 ∩V , then the rank of M 1 decreases by at least 2 in V 1 ∩V .
Choose the minor associated to M 1 [13, 13] and denote it by M f . Compute det M f with the method used so far. The minor has two diagonal blocks, denote the upper left by R f and the lower right by S f . The upper right block is the zero matrix thus The associated lower right block is S 1 , hence: The minor σ f defines a variety W with codim W = 1. The polynomial σ f factors as χ f φ f and the two polynomials are: The coefficient of the highest power in q 1 of φ f is −p 1 , hence φ f = 0. Moreover, σ f is generically irreducible.
Finally, substitute the computed point P from Subsection 4.3 to σ f . We get a nonzero number, hence the rank of M 1 decreases by 1 at the point P . This is also true in an open neighborhood of P due to conjecture that the varieties are smooth at P . Hence there exists an open subset of V 1 ∩V , where the rank of M 1 decreases by exactly 1:

4.5.
Blowing up along (V 1 ∩V ) 0 . The (irreducible) varieties V 1 andV specify a 1-dimensional linear system D λ,k = (σ 1 + λσ k = 0) in M which defines a base locus It is just the variety V 1 ∩V . We want to blow up along the base locus, but D λ,k is a pencil, hence we can fix only one parameter in the system of the linear equations (T ). Consequently, we do not get a unique solution where the rank of the system decreases by more than 1.
The variety V 1 ∩V ∩ W \ ∆ defines a multi-dimensional linear system but we present the blow up procedure along a 1-dimensional linear system only, hence we blow up along (V 1 ∩V ) 0 .
Let Ω ⊂ M be a 2N = 8 dimensional ball with holomorphic coordinates q 1 , . . . , q 4 , p 1 , . . . , p 4 . The base locus is defined by equations σ 1 = σ k = 0. (Again k ∈ {2, . . . , 33} \ {2, 14, 15, 33}.) Denote the coordinates in the base locus the following way: Q 1 = σ 1 , P 1 = σ k , z 3 = q 2 , z 4 = q 3 , z 5 = q 4 , z 6 = p 2 , z 7 = p 3 , Let [l 1 , l 2 ] be homogeneous coordinates in CP 1 . Our aim is to determine these coordinates. Denote the blow up by The projection to the first coordinate is π : ∼ Ω → Ω and the exceptional divisor is Introduce two local charts U 1 and U 2 on ∼ Ω. The charts U 1 and U 2 are given by l 1 = 0 and l 2 = 0 respectively. The coordinates on chart U 1 are The coordinates on chart U 2 are Obviously, the coordinate systems on the charts U 1 and U 2 are compatible and the transition functions are easy to work out. Note that, the above ratio P1 Q1 determines the solution of our System (T ) in the variable H k−2 or I k−15 due to Cramer's rule as long as Q 1 = 0 (i. e. outside the variety V 1 ): However, the ratio P1 Q1 is undefined in (V 1 ∩V ) 0 . In other words, H k−2 or I k−15 is not determined by the coordinates (Q 1 = 0, P 1 = 0, z 3 , . . . , z 8 ), hence these will be the free parameter in the solution of System (T ). If we fix the value of the free parameter on the exceptional divisor E, then we have a solution on (V 1 ∩V ) 0 . The coordinates on the charts U 1 and U 2 are: if k=3, . . . ,13. The case k = 16, . . . , 32 is similarly computable. The coordinate y 1 is not equal zero in the chart U 2 due to H k−2 , I k−15 = ∞, because H k−2 and I k−15 is a solution of the system of linear equations. The rank of the augmented matrix is greater than the rank of the coefficient matrix in the set V 1 \V . Consequently, we need to subtract the proper transform of V 1 from the blowing up. If the varieties V 1 andV are generically transverse, then the proper transform of V 1 will be the following:  5. The third order system with arbitrary number of parabolic points 5.1. The existence of the variety. We return to the description of the general case where we interrupted in Section 3. Hence we fix n + 1 parabolic points on the Riemann surface Σ and we fix exponents ρ i,k i = 0, . . . , n, k = 1, 2, 3, we have N = 3n − 5 apparent singularities and the same number of free parameters. In this section we will prove Theorem 1.7 with the generalization of the method in Subsection 4.2. We will not compute the discriminant variety exactly, we will only prove its existence. The aim is to show that the coefficient which belongs to the highest order exponent of q 1 in the polynomial associated to the variety is not zero. As we saw in Proposition 3.6, the first 4n − 4 rows of the coefficient matrix are independent of the other rows. Denote the coefficient matrix of System (T ) in the case of general n by M .
Compute the determinant of M with splits and row interchanges along the 3n − 5 common rows. Make confluent Vandermonde matrices as in Subsection 4.2. Generalize the notation for R j and S j the following way: J will be a multiindex which contains an index j if q j appears in the last n − 2 rows in the upper left (8n − 11) × (8n − 11) block. Hence R J is an (8n − 11) × (8n − 11) matrix where the last n − 2 rows contain parameters q j1 , q j2 , . . . , q jn−2 , where j 1 < j 2 < · · · < j n−2 and j l ∈ J for all l = 1, . . . , n − 2. Similarly, J in S J refers to the set J := {1, 2, . . . , 3n − 5} \ J. Through out this section we will use letter j for index from J and letter for index from J, furthermore i = 1, . . . , n.
Two splittings of M are different if their indices J 1 and J 2 are not the same. On the other hand, all possible multi-indices J appear in the expansion of the determinant of M . This shows that the expansion has 3n−5 n−2 pairs of (R J , S J ). One of these pairs is the diagonal block of a matrix which has a zero block in the upper right or the lower left corner.
It is enough to consider the pairs which contribute to the coefficient of the highest exponent of q 1 .

5.2.
The highest exponent. First, we want to decide when we get higher q 1 power: if the common row with q 1 is in R J or is in S J . The quantity and the order of parameters in R J : n type t k . The first 2n−3 rows of the matrices S J come from Equation (25a) and last 3n−5 rows come from Equation (25b). The determinant does not change if we eliminate the coefficients related to the second derivatives of the I(q k )'s from 2n − 3 rows in the last 3n − 5 rows. The remaining n − 2 rows from the last 3n − 5 rows contain q j 's with j ∈ J. Split S J along these n − 2 rows such that one of them contains only the coefficients related to the second derivatives of the I(q j )'s (denoted by S ′ J ), others contain one more coefficient related to the third derivative of I(q j ) thus deg (q1−q k ) (S ′ J ) is greater than degree of any other matrix which comes from the splitting of S J .
If j = 1 then deg q1 (ν j ) = 9 − 3n > deg q1 (ν j ) = 8 − 3n and if j = 1 then deg q1 (ν j ) = 18 − 12n > deg q1 (ν j ) = 17 − 12n according to Equation (22). Hence the coefficient of the highest q 1 exponent contains the factor j∈J p j . Each polynomials in the first sum in the second row of the Equation (31) have different j∈J p j factors, hence the first sum does not vanish.
Proof of Theorem 1.7. The above discussion shows that there exists a nonzero Sym(N ) invariant polynomial with coordinates {q j , p j } j=N j=1 . The polynomial defines an affine subvariety V in T * (C \ P ) N and a part of T * (C \ P ) N outside V is a dense open set of the moduli space M.