Parabolic degrees and Lyapunov exponents for hypergeometric local systems

Consider the flat bundle on $\mathrm{CP}^1 - \{0,1,\infty \}$ corresponding to solutions of the hypergeometric differential equation $ \prod_{i=1}^h (\mathrm D - \alpha_i) - z \prod_{j=1}^h (\mathrm D - \beta_j) = 0$ where $\mathrm D = z \frac {d}{dz}$. For $\alpha_i$ and $\beta_j$ distinct real numbers, this bundle is known to underlie a complex polarized variation of Hodge structure. Setting the complete hyperbolic metric on $\mathrm{CP}^1 - \{0,1,\infty \}$, we associate $n$ Lyapunov exponents to this bundle. We compute the parabolic degrees of the holomorphic subbundles induced by the variation of Hodge structure and study the dependence of the Lyapunov exponents in terms of these degrees by means of numerical simulations.


Introduction
Oseledets decomposition of flat bundles over an ergodic dynamical system is often referred to as dynamical variation of Hodge structure. In the case of Teichmüller dynamics, both Oseledets decomposition and a variation of Hodge structure (VHS) appear. Two decades ago it was observed in [Kon97] that these structures were linked, their invariants are related: the sum of the Lyapunov exponents associated to a Teichmüller curve equals the normalized degree of the Hodge bundle. This formula was studied extensively and extended to strata of abelian and quadratic differentials from then (see [FMZ14], [Kri04], [BM10], [EKZ14]). Soon this link was observed in other settings. In [KM16] it was used as a new invariant to classify hyperbolic structures and distinguish Deligne-Mostow's non-arithmetic lattices in SL 2 (C). In [Fil14] a similar formula was observed for higher weight variation of Hodge structures. The leitmotiv in this work is the study of the relationship between theses two structures in a broad class of examples with arbitrary weight. This family of examples will be given by hyperelliptic differential equations which yield a flat bundle endowed with a variation of Hodge structure over the sphere with three punctures. A recent article [EKMZ16] shows that the degrees of holomorphic flags of the Hodge filtration bound by below the sum of Lyapunov exponents. Our investigation will start by computing these degrees and then explore the behaviour of Lyapunov exponents through numerical simulations and their distance to the latter lower bounds. This will enable us to bring out some simple algebraic relations under which there is a conjectural equality.
Hypergeometric equations. Let α 1 , α 2 , . . . , α n and β 1 , β 2 , . . . , β n be two disjoint sequences of n real numbers. We define the hypergeometric differential equation corresponding to those parameters This equation originates from a large class of special functions called generalized hypergeometric functions which satisfies it. For more details about these functions see for example [Yos97].
It is an order n differential equation with three singularities at 0, 1 and ∞ hence the space of solutions is locally a dimension n vector space away from singularities and can be seen in a geometrical way as a flat bundle over P 1 − {0, 1, ∞}. This flat bundle is completely described by its monodromy matrices around singularities. We will be denoting monodromies associated to simple closed loop going counterclockwise around 0, 1 and ∞ by M 0 , M 1 and M ∞ . We get a first relation between these matrices observing that composing the three loops in the same order will give a trivial loop: M ∞ M 1 M 0 = Id. The eigenvalues of M 0 and M ∞ can be expressed with parameters of the hypergeometric equation (1) and M 1 has a very specific form as stated in the following proposition. Proof. See Proposition 2.1 in [Fed15] or alternatively Prop. 3.2 and Theorem 3.5 in [BH89] This proposition determines the conjugacy class of the representation associated to the flat bundle π 1 P 1 − {0, 1, ∞} → GL n (R) thanks to the rigidity of hypergeometric equations (see [BH89]). They will be computed explicitly in section 3.2.
Lyapunov exponents. We now endow the 3 punctured sphere with its hyperbolic metric. As this metric implies an ergodic geodesic flow, for any integrable norm on the flat bundle we associate to it , using Oseledets theorem, a measurable flag decomposition of the vector bundle and n Lyapunov exponents. These exponents correspond to the growth of the norm of a generic vector in each flag while transporting it along with the flat connection.
According to [EKMZ16] there is a canonical family of integrable norms on the flat bundle associated to the hypergeometric equation which will produce the same flag decomposition and Lyapunov exponents. This family contains the harmonic norm induced by the VHS structure and the norm we will use in our algorithm.
Variation of Hodge structure. Hypergeometric equations on the sphere are well known to be physically rigid (see [BH89] or [Kat96]) and this rigidity together with irreducibility is enough to endow the flat bundle with a VHS using its associated Higgs bundle structure (see [Fed15] or directly Cor 8.2 in [Sim90]). Using techniques from [Kat96] and [DS13], Fedorov gives in [Fed15] an explicit way to compute the Hodge numbers for the underlying VHS. We extend this computation and give a combinatorial point of view that will be more convenient in the following to express parabolic degrees of the Hodge flag decomposition.
We introduce a canonical way to describe combinatorics of the intertwining of α's and β's on the circle R/Z. Starting from any eigenvalue, we browse the circle counterclockwise (or in the increasing direction for R) and denote α 1 , . . . , α n , β 1 , . . . , β n by order of appearance η 1 , η 2 , . . . , η 2n and definef : Z ∩ [0, 2n] → Z recursively by the following properties, Let f be defined on the eigenvalues by f (η k ) =f (k). It depends on the choice of starting point up to a shift. For a canonical definition, we shift f such that its minimal value is 0. It is equivalent to starting at the point of minimal value. This defines a unique f which we intertwining diagram of the equation.
Then we have the following theorem, Theorem (Fedorov). The h 1 , h 2 , . . . , h n are the Hodge numbers of the VHS after an appropriate shifting.
Remark. If the α's and β's appear in an alternate order then f (α) ≡ 1 and f (β) ≡ 0 thus their is just one element in the Hodge decomposition and the polarization form is positive definite. In other words the harmonic norm is invariant by the flat connection. This implies that Lyapunov exponents are zero.
In general, this Hodge structure endows the flat bundle with a pseudo-Hermitian form of signature (p, q) where p is the sum of the even Hodge numbers and q the sum of the odd ones. This gives classically the fact that the Lyapunov spectrum is symmetric with respect to 0 and that at least |p − q| exponents are zero (see Appendix A in [FMZ14]).
Pushing further methods of [Fed15] and [DS13], we compute the parabolic degree of the sub Hodge bundles. This computation was done with the help of computer experiments in section 4.2 which yielded a conjectural formula for these degrees. Besides from the intertwining diagram, another quantity emerges to express them; relabel α and β by order of appearance after choosing α 1 such that f (α 1 ) = 0, then take the representatives of α and β in R which are included in [α 1 , α 1 + 1[ and define γ := β − α. The formula will depend on the floor value of γ. As 0 < γ < n we have n possible values 0 ≤ [γ] < n. Theorem 1.2. Let 1 ≤ p ≤ n and E p the p-th graded piece of the Hodge filtration on P 1 − {0, 1, ∞}. We denote by δ p the degree of the Deligne compactification of E p on the sphere. Then, Acknowledgement. I am very grateful to Maxim Kontsevich for sharing this problem and taking time to discuss it. I thank dearly Jeremy Daniel for his curiosity to the subject and his answer to my myriad of questions as well as Bertrand Deroin; Anton Zorich for his flawless support and attention, Martin Möller and Roman Fedorov for taking time to explain their understanding of the parabolic degrees and Hodge invariants at MPIM in Bonn. I am also very thankful to Carlos Simpson for his kind answers and encouragements. A (C-)VHS on a curve C consists of a complex local system V C with a connection ∇ and a decomposition of the Deligne extension V = p∈Z E p into C ∞subbundles, satisfying: are holomorphic (resp. antiholomorphic) subbundles for every p ∈ Z. • The connection shifts the grading by at most one, i.e.
Up to a shift, we can assume that there is a n such that E i = 0 for i < 0 and i > n. We call n the weight of the VHS. We can introduce for convenience the notation E p,n−p := F p /F p−1 . Then we have a C ∞ isomorphism between the bundles: n i=0 E i,n−i ∼ = V 2.2. Decomposition of an extended holomorphic bundle. Let C be a complex curve, we assume that its boundary set ∆ := C\C is an union of points. Consider E an holomorphic bundle on C. We introduce structures which will appear on such holomorphic bundle when they are obtained by canonical extension when we compactify C. The first one will take the form of filtrations on each fibers above points of ∆.
Definition 2.1 (Filtration). A [0, 1)-filtration on a complex vector bundle V is a collection of real weights 0 ≤ w 1 < w 2 < · · · < w n < w n+1 = 1 for some n ≥ 1 together with a filtration of sub-vector spaces The filtration satisfies V ≥ν ⊂ V ≥ω whenever ν ≥ ω and the previous weights satisfy V ≥w i + V ≥w i for any > 0.
We denote the graded vector bundles by gr w i := V ≥w i / V ≥w i + for small. The degree of such a filtration is by definition This leads to the next definition, A parabolic bundle is a holomorphic bundle endowed with a parabolic structure. The parabolic degree of (E, G • ) is defined to be 3. Deligne extension. In the following we consider V a flat bundle on P 1 − {0, 1, ∞} associated to a monodromy representation with norm one eigenvalues. We denote by V C the associated holomorphic vector bundle.
We recall the construction of Deligne's extension of V C which defines a holomorphic bundle on C with a logarithmic flat connection. We describe it on a small pointed disk centered at s ∈ ∆ with coordinate q ∈ D * . Let ρ be a ray going outward of the singularity, then we can speak of flat sections along the ray L(ρ) which has the same rank r as V. As all the L(ρ) are isomorphic, we choose to denote it by V 0 . There is a monodromy transformation T : V 0 → V 0 to itself obtained after continuing the solutions. This corresponds to the monodromy matrix in the given representation. For every α ∈ [0, 1) we define These vector spaces are non trivial for finitely many α i ∈ [0, 1). We define These sections are equivariant under z → z + 1 hence they give global sections of V C (D * ). The Deligne extension of V C is the vector bundle whose space of section over D is the O D -module spanned byṽ 1 , . . . ,ṽ r . This construction naturally gives a filtration on V 0 .
In general, we can define various extensions V a ⊂ V −∞ ⊂ j * V where j is the inclusion j : C → C, V ∞ is the Deligne's meromorphic extension and V a (resp. V >a ) for a ∈ R is the free O C -module on which the residue of ∇ has eigenvalues α in [a, a + 1) (resp. (a, a + 1]). The bundle V • is a filtered vector bundle in the definition of [EKMZ16].
If we have a VHS F • on E over C, it induces a filtration of every V a simply by taking F p V a := j * F p V ∩ V a this is a well defined vector bundle thanks to Nilpotent orbit theorem (see [DS13]).
We define over some singularity s ∈ ∆, for a ∈ (−1, 0] and λ = exp(−2iπa), Sim90]) claims that for any local system with all eigenvalues of the form exp(−2πiα) at the singularities endowed with a trivial filtration we associate a filtered D C -module with residues and jumps both equal to α. Thus the sub D C -module corresponding to the residue α has only one jump of full dimension at α, and where we choose α ∈ [0, 1).
2.4. Acceptable metrics and metric extensions. The above Deligne extension has a geometric interpretation when we endow C with a acceptable metric K. If V is a holomorphic bundle on C, we define the sheaf Ξ(V ) α on C ∪ {s} as follows. The germs of sections of Ξ(E) α at s are the sections s(q) in j * V in the neighborhood of s which satisfy a growth condition; for all ≥ 0 there exists C such that |s(q)| K ≤ C |q| α− . In general this extension is a coherent sheaf on which we do not have much information, but Simpson shows in [Sim90] that under some growth condition on the curvature of the metric, the metric induces the above Deligne extensions. When a curvature satisfies this condition it is called acceptable.

Lemma 2.4 (Lemma 2.4 [EKMZ16]). The local system V with non-expanding cusp monodromies has a metric which is acceptable for its Deligne extension V
Proof. For completeness, we reproduce the construction of [EKMZ16]. The idea is to construct locally a nice metric and to patch the local constructions together with partition of unity. The only delicate choice is for the metric around singularities. We want the basis elementsṽ i of the α-eigenspace of the Deligne extension to be given the norm of order |q| α in the local coordinate q around the cusp and to be pairwise orthogonal. Let M be such that e 2iπM = T , where T is the monodromy transformation. Then the hermitian matrix exp(log|q|M t M ) defines a metric such that the elementṽ i has norm |q| α |ṽ i |.
Corollary 2.5. When the monodromy representation goes to zero, the parabolic degree goes to zero.
Proof. In the proof above, it is clear that when T → Id, M → 0 and thus the metric goes to the standard hermitian metric locally. Thus its curvature goes to zero around singularities and its integral on any subbundle, which by definition is its parabolic degree, goes to zero.
2.5. Local Hodge invariants. Our purpose in this subsection is to show the following relation on local Hodge invariants: Theorem 2.6. The local Hodge invariants for equation 1 are : (1) at z = 0, Remark. Computations of (1) and (2) are done in [Fed15]. We give a similar proof with an alternative combinatoric point of view.
2.6. Computation of local Hodge invariants. In the following, we denote by M the local system defined by the hypergeometric equation 1 in the introduction. The point at infinity plays a particular role in middle convolution, thus we apply a biholomorphism to the sphere which will send the three singularity points 0, 1, ∞ to 0, 1, 2. Hereafter, M will have singularities at 0, 1, 2.
Theorem 2.8 (Dettweiler-Sabbah). Let α 0 ∈ (0, 1), for every singular point in ∆ and any α ∈ [0, 1) we have, and, 2.6.1. Recursive argument. We apply a recursive argument on the dimension of the hypergeometric equation. Let us assume that n ≥ 3 and that Theorem 2.6 is true for n − 1.
For convenience in the demonstration, we change the indices of α and β such that α i (resp. β i ) is the i-th α (resp. β) we come upon while browsing the circle to construct the function f .
We apply Lemma 2.7 with α k and β j such that α k ≤ β j . Let us describe what happens to the combinatorial function f after we remove these two eigenvalues.
We denote by f the function we obtain.
Removing α k will make the function decrease by one for the following eigenvalues until we meet β j , thus for any τ = α k , β j , We apply Theorem 2.8 with α 0 = β j − α k . It yields that we have for all m = k, at singularity zero, In terms of M and M k,j , For any integer i, j we denote by δ(i, j) the function which is 1 when i = j and is zero otherwise.
Similarly for all m = j, at singularity 2, And at 1, we setγ : Now if we choose to pick β n > α n for the computation, we now hodge invariants for all values excepts for α n and β n . We will use the computation with for example α 1 < β 1 . From this one we can deduce the invariants at β n and α n . Yet, we should keep in mind that the previous computations are always modulo shifting of the VHS. That is why we need to have dimension at least 3, since in this case α 2 will appear in both computations and will show there is no shift in our formulas.

Continuity of the parabolic degree.
To compute δ p (V ) in equation 2, we show in the following Lemma a continuity property which implies that it is constant on a given domain. Let α 1 , . . . , α n , β 1 , . . . , β n be all disjoint, not integers and such that γ also is not an integer. Fix the intertwining diagram of α and β and the value of [γ], then for any integer p, δ p (V ) is constant.

Lemma 2.9.
Proof. Let L and L be the local system corresponding to equation 1 for eigenvalues α, β (resp. α , β ) satisfying the above hypothesis. We endow them with a trivial filtration. According to Simpson's theory, L corresponds to some Higgs bundle (E, θ) together with a parabolic structure at singularities. As L has eigenvalues of norm one, θ has no residue, moreover its weight filtration is locally the same as the one for the unipotent part of monodromy matrices of L. Consider now (E , θ ) a Higgs bundle with the same holomorphic structure and Higgs form as (E, θ) but a slightly changed parabolic structure for which we keep the initial filtration but modify the parabolic weights α, β to α , β . It is clear that we keep the same residues res s (θ) = res s (θ ) at every singularity s. We also keep locally the same weight filtration for the unipotent part of the monodromy on any eigenspace. This implies that monodromy matrices of the local system associated to (E , θ ) and those of L are locally isomorphic, and by rigidity of the hypergeometric local systems they are globally isomorphic. The same argument applies for Hodge subbundles. As the considered domain is connected this shows that the parabolic degree of the Hodge subbundles is constant.
Together with Corollary 2.5 this is enough to compute δ p (V ). We fix an intertwining diagram and a floor value for γ and make the first [γ] β go to 1 and the rest of eigenvalues to 0 while staying in the given domain. At the limit, the parabolic degree is zero and we can deduce δ p (V ).

Algorithm
In this section, we describe the algorithm used to compute the Lyapunov exponents. We start simulating a generic hyperbolic geodesic and following how it winds around the surface, namely the evolution of the homology class of the closed path. Finally we compute the corresponding monodromy matrix after each turn around a cusp.
3.1. Hyperbolic geodesics. This first question arising to unravel this computation of Lyapunov exponents is how to simulate a generic hyperbolic geodesic. The answer comes from a beautiful theorem proved by Caroline Series in [Ser85] which relates hyperbolic geodesics on the Poincaré half-plane and continued fraction development of real numbers. We follow here the notations of [Dal07] (see part II.4.1).
Let us consider the Farey tessellation of H (see Figure 3). It is the fundamental The sphere minus three points endowed with its complete hyperbolic metric is a degree two cover of the surface associated to Farey's tessellation. This is why we represent the tessellation with two colors : the fundamental domain for the sphere corresponds to two adjacent triangles of different colors. That is why it will be easy once we understand the geodesics with respect to this tessellation to see them on the sphere.
Let us consider a geodesic going through i. It lands to the real axis at a positive and a negative real number. The positive real number will be called x, this number determines completely the geodesic since we know two distinct points on it.
We associate to this geodesic a sequence of positive integers. Look at the sequence of hyperbolic triangles the geodesic will cross. For each one of those triangles, the geodesic has two ways to cross them (see Figure 3). Once it enters it, it can leave it crossing either the side of the triangle to its left (a) or to its right (b).

Remark. The vertices of hyperbolic triangles are located at rational numbers, so this sequence will be infinite if and only if x is irrational (see [Dal07] Lemme 4.2).
We have now for a generic geodesic an infinite word in two letters L and R associated to a geodesic. For example the word associated to the geodesic in Figure 5, is of the form LLRRLR · · · = L 2 R 2 L 1 R . . . . We can factorize each of those words and get R n 0 L n 1 R n 2 L n 3 . . .
Except for n 0 which can be zero the n i are positive integers.  Remark. This theorem states exactly that to study a generic geodesic on the hyperbolic plane, we can consider a Lebesgue generic number in (0, ∞) and compute its continued fraction development.

Figure 5. Crossings of a given geodesic
To compute Lyapunov exponents of the flat bundle, we need to follow how a generic geodesic winds around the cusps. By the previous theorem we can simulate a generic cutting sequence of an hyperbolic geodesic in H. Our goal now will be to associate to such a sequence a product of monodromy matrices following its homotopy class.
Since we will consider universal cover of the sphere minus three points, for convenience we will denote by A, B, C the cusps corresponding in the surface to ∞, 0, 1 respectively and use this latter notation for points in H. Two adjacent hyperbolic triangles of Farey's tessellation, e.g. 0, 1, ∞ and 1, 2, ∞ will form a fundamental domain for this surface. All the vertices of the hyperbolic triangles for this tessellation are associated to either A, B or C. To follow how the flow turns around these vertices in the surface, we will need to keep track of orientation. To do so, we color the triangles according to the order of its vertices, when we browse the three vertices counterclockwise if we have A → B → C → A we color the triangle in white (this is the case for 0, 1, ∞) otherwise we color it in blue (case of ∞, 2, 1).
Let us now consider a point P inside the blue triangle, which will be used as a base point for an expression of the cycles around cusps. We choose a homology marking of the surface by denoting the paths going around A, B, C counterclockwise starting and ending at P , a, b, c (see Figure 6). When we concatenate these paths we get c · b · a = Id and a −1 = c · b. For monodromy matrices we will have the relation In our algorithm we will always follow the cutting sequence until we end up to a blue triangle. Then we will apply an isometry that take the fundamental domain we are in to the (−1, 0, ∞) triangle and the edge the flow will cut when going out of the triangle to be the (0∞) or (−1, ∞) edge in order to place the cusp we are turning around at ∞. We shall warn the reader here that the corresponding cusp on the surface here at points −1, 0 and ∞ may be any of the points A, B, C but their cyclic order will stay unchanged thanks to the orientation. Thus we just need to keep track of the cusp placed at ∞.
When we start with a cutting sequence extracted from the previous theorem we see that the geodesic start by cutting (0, ∞) at i without being counted in the cutting sequence. The first cutting will always be forgotten in the sequence when applying the isometry. Now remark that when the crossing is a sequence of 2n left, we make n turns counterclockwise around the cusp placed at ∞. When it 2n right, we make n turn clockwise. It is a little trickier if the geodesic makes an odd number of the same crossing; we need to take one step further from the next term in the sequence of crossings to end up at P (see Figure 7). There is a last point to consider, since we want to compare the growth of the harmonic norm with regards to the geodesic flow, need to follow its length. Here the discretized algorithm enables us to follow the type of homotopy it will have, but the length will not correspond a priori to the number of iterations of our algorithm. It is proportional to it by the constant.

Monodromy matrices.
In the introduction Proposition 1.1 gave a set of three properties on the monodromy matrices for the hypergeometric differential equation associated to two distinct sequences of real numbers α 1 , . . . , α n and β 1 , . . . , β n . We claim that those properties are sufficient to recover the monodromy matrices up to conjugacy.
For convenience we always assume that the α 1 , . . . , α n are disjoint, otherwise the computation becomes way more tedious, and in our computations we will explore generic domains. We choose a basis in which M 0 is diagonal. Property (3) tells us that M 1 − Id is of rank 1. We can then find two vectors v and w such that M 1 = Id +vw t .
Since M −1 ∞ = M 0 M 1 knowing the eigenvalues of M ∞ we can derive the following n equations, for all j, det M −1 ∞ − e 2iπβ j Id = 0 We can compute this determinant using the particular form of the matrix and the following lemma.
We can conjugate by diagonal matrices so that M 0 v becomes the vector 1 which is one on every coordinates. And obtain the equations det((M 0 − e 2iπβ j Id) + 1w t ) = 0, ∀j Lemma. Let D a diagonal matrix with d 1 , . . . , d n on its diagonal, and x a vector.
Proof. First consider the case where D is the identity matrix. We know that all the eigenvalues except for one are 1. The determinant will then be the eigenvalue of an eigenvector which image through x t is not zero. This vector will be 1 and its eigenvalue (1 + n i=1 x i ). To finish the proof, just factor each column by d i in the determinant.

We obtain
Corollary. The vector w satisfies for all j, We define a matrix N = Thus according to the same article, we know that 2(µ 1 + µ 2 ) is a lower bound for the sum of Lyapunov exponents. We call good cases the equality cases and bad cases the cases where there is strict inequality.  To see what happens in a similar setting for more general hypergeometric equations, we vary C, d and compute the corresponding eigenvalues µ 1 and µ 2 as well as the Lyapunov exponents. On Figure 9a we drew a blue point at coordinate (µ 1 , µ 2 ) if the sum of positive Lyapunov exponents are as close to the parabolic degree 2(µ 1 + µ 2 ) as the precision we have numerically and we put a red point when this value is outside of the confidence interval. Note that according to Figure 9a it seems that all points below the line of equation 3µ 2 = µ 1 + 1 are bad cases. In Figure 9b, we represent the distance of the sum of the Lyapunov exponents to the expected formula. We see that this gives a function that oscillates above zero. More precisely, it seems that good cases are outside of some lines passing through (1/2, 1/2).
According to [BT14] and [SV14], the 7 good cases correspond to cases where the monodromy group of the hypergeometric local system is of infinite index in Sp(4, Z), which is commonly called thin. In the other cases the group is of finite index and is called thick. The three good cases we found by ways of Lyapunov exponents do not seem to have a representation with integers C and d. A lot of questions arise about these points, for example can we find a number-theoretic interpretation of their equality as in Conjecture 6.5 in [EKMZ16]. 4.2. Examples for n = 2. Has we have seen in the introduction the two Lyapunov exponents are symmetric λ 1 and −λ 1 . The sum of the positive Lyapunov exponents is just λ 1 . The parameter space we have for these 2-dimensional flat bundles are α 1 , α 2 , β 1 , β 2 .
The Lyapunov exponents are invariant through translation of the set of parameters. Indeed, we can consider the bundle with e δ M 0 and e −δ M ∞ monodromies, it will have the same set of Lyapunov exponents since both scalar will appear with the same frequency and its parameters will be α 1 + δ, . . . , α h + δ, β 1 + δ, . . . , β h + δ hence without loss of generality we can assume β 1 = 0. Moreover the parameters are given as a set, the order does not matter.
In the following experiments we will consider a set of parameters where the β's will be equidistributed and the α's will be shifted with respect to them. Here we represent the value of the Lyapunov exponent for α 1 = r, α 2 = 2r, β 1 = 0, β 2 = x and we have by definition γ = x − 3r. This will be true whenever the VHS has weight 0.
Another noticeable fact is that zones correspond exactly to different combinatorics for the order of the α and β, and on [γ] introduced in the introduction. Remark that [γ] is 0 in zones 1, 4, and 1 in zones 2, 5. In the following table, we give a relation binding λ 1 , r, x obtained by linear regression. The other column is the formula for the parabolic degree in the given zone. Zone In this case, the VHS is of weight ≤ 1 and thus is in the setting of [Kon97]. In consequence, we have the equality Where deg par is the parabolic degree of the holomorphic bundle and χ(S) = 1 the Euler characteristic of S.
This is a good test for our algorithm and formula on degree. More generally, for any dimension n, this formula will hold as long as the weight is equal to 1.
Using a Monte-Carlo process, we found some values in this configuration for which there is equality with twice the parabolic degree of E 2 ⊕ E 1 . We remarked that several parameter points where there is equality satisfy θ 1 = θ 2 and θ 4 = θ 5 . This motivated us to consider the 2 dimensional subspace of parameters (θ 1 , θ 2 , θ 3 , θ 4 , θ 5 ) = (x, x, 1/2, y, y) For these parameters we can observe a remarkable phenomenon; the difference between the Lyapunov exponent and the formula with parabolic degrees depends only on x + y. We plot this difference in the Figure below and see that for some values of x + y there is equality.