Dawei Chen
email: dawei.chen@bc.edu

Martin Möller
email: moeller@math.uni-frankfurt.de

Adrien Sauvaget
email: adrien.sauvaget1@u-cergy.fr

Don Zagier
email: donzagier@mpim-bonn.mpg.de

Masur-Veech volumes and intersection theory on moduli spaces of Abelian differentials

We show that the Masur-Veech volumes and area Siegel-Veech constants can be obtained using intersection theory on strata of Abelian differentials with prescribed orders of zeros. As applications, we evaluate their large genus limits and compute the saddle connection Siegel-Veech constants for all strata. We also show that the same results hold for the spin and hyperelliptic components of the strata.

Introduction

Computing volumes of moduli spaces has significance in many fields. For instance, the Weil-Petersson volumes of moduli spaces of Riemann surfaces can be written as intersection numbers of tautological classes due to the work of Wolpert [START_REF] Wolpert | On obtaining a positive line bundle from the Weil-Petersson class[END_REF] and of Mirzakhani for hyperbolic bordered surfaces with geodesic boundaries [START_REF] Mirzakhani | Weil-Petersson volumes and intersection theory on the moduli space of curves[END_REF]. In this paper we establish similar results for the Mazur-Veech volumes of moduli spaces of Abelian differentials.

Denote by M g,n (μ) the moduli spaces (or strata) of Abelian differentials (or flat surfaces) with labeled zeros of type μ = (m 1 , . . . , m n ), where m i ≥ 0 and where n i=1 m i = 2g -2. Masur [START_REF] Masur | Interval exchange transformations and measured foliations[END_REF] and Veech [START_REF] Veech | Gauss measures for transformations on the space of interval exchange maps[END_REF] showed that the hypersurface of flat surfaces of area one in M g,n (μ) has finite volume, called the Masur-Veech volume, and we denote it by vol ( M g,n (μ)). The starting point of this paper is the following expression of Masur-Veech volumes in terms of intersection numbers on the incidence variety compactification P M g,n (μ) described in [START_REF] Bainbridge | Compactification of strata of abelian differentials[END_REF]. For 1 ≤ i ≤ n we define

β i = 1 m i + 1 ξ 2g-2 j =i ψ j ∈ H 2(2g-3+n) (P M g,n (μ), Q) ( 1
)
where ξ is the universal line bundle class of the projectivized Hodge bundle and ψ j is the vertical cotangent line bundle class associated to the jth marked point (see Sect. 3 for a more precise definition of these tautological classes).

Theorem 1.1 The Masur-Veech volumes can be computed as intersection numbers

vol ( M g,n (m 1 , . . . ,

m n )) = - 2(2iπ) 2g (2g -3 + n)! P M g,n (μ) ξ 2g-2 • n i=1 ψ i (2) = 2(2iπ) 2g (2g -3 + n)! P M g,n (μ) β i • ξ (3)
for each 1 ≤ i ≤ n.

The equality of the two expressions on the right-hand side is a non-trivial claim about intersection numbers on P M g,n (μ). Note that we follow the volume normalization in [START_REF] Eskin | Moduli spaces of Abelian differentials: the principal boundary, counting problems, and the Siegel-Veech constants[END_REF] that differs slightly from the one in [START_REF] Eskin | Asymptotics of numbers of branched coverings of a torus and volumes of moduli spaces of holomorphic differentials[END_REF] (see [START_REF] Chen | Quasimodularity and large genus limits of Siegel-Veech constants[END_REF]Section 19] for the conversion).

Theorem 1.1 is the interpolation and generalization of [42, Proposition 1.3] (for the minimal strata) and [START_REF] Chen | Quasimodularity and large genus limits of Siegel-Veech constants[END_REF]Theorem 4.3] (for the Hurwitz spaces of torus covers). In order to prove it, we show that both sides of Eq. ( 2) satisfy the same recursion formula. On the volume side, the recursion formula is expressed via an operator acting on Bloch and Okounkov's algebra of shifted symmetric functions (see Sect. 4). The recursion for intersection numbers is first proved at the numerical level using the techniques developed in [START_REF] Sauvaget | Cohomology classes of strata of differentials[END_REF] to compute the classes of P M g,n (μ) (see Sects. 2, 3). Then we formally lift this relation to the algebra of shifted symmetric functions and show that it is equivalent to the previous one (see Sect. 5).

In particular, the recursion arising from intersection calculations provides the following useful formula. We define the rescaled volume

v(μ) = (m 1 + 1) • • • (m n + 1)vol ( M g,n (m 1 , . . . , m n )) . ( 4
)
For a partition μ, we denote by n(μ) the cardinality of μ and by |μ| the sum of its entries.

Theorem 1. [START_REF] Athreya | Right-angled billiards and volumes of moduli spaces of quadratic differentials on CP 1[END_REF] The rescaled volumes of the strata satisfy the recursion

v(μ) = k≥1 g,μ h P 1 ((m 1 , m 2 ), p) • k i=1 (2g i -1 + n(μ i ))! v(μ i , p i -1) 2 k-1 k! (2g -3 + n)! ( 5
)
where g = (g 1 , . . . , g k ) is a partition of g, μ = (μ 1 , . . . , μ k ) is a k-tuple of multisets with (m 3 , . . . , m n ) = μ 1 The relevant Hurwitz spaces of P1 covers will be introduced in Sect. 2. Note that h P 1 ((m 1 , m 2 ), p) = 0 only if k i=1 ( p i + 1) = m 1 + m 2 + 2. This implies that k ≤ min(m 1 + 1, m 2 + 1) in the summation of the theorem.

For special μ, the strata M g,n (μ) can be disconnected. There are up to three connected components altogether, at most one of which is hyperelliptic, classified by Kontsevich and Zorich [START_REF] Kontsevich | Connected components of the moduli spaces of Abelian differentials with prescribed singularities[END_REF]. We show the refinements of Theorems 1.1 and 1.2 for the spin and hyperelliptic components respectively, given as Theorems 6.3 and 6.12 (conditional on a technical Assumption 6.1 in Sect. 6, which can be deduced from [START_REF] Bainbridge | The moduli space of multi-scale differentials[END_REF]). 1 Equation ( 5) has a similar form compared to the recursion formula obtained by Eskin, Masur and Zorich [START_REF] Eskin | Moduli spaces of Abelian differentials: the principal boundary, counting problems, and the Siegel-Veech constants[END_REF] for computing saddle connection Siegel-Veech constants (joining two distinct zeros). Consider a generic flat surface with n labeled zeros of orders μ = (m 1 , . . . , m n ). The growth rate of the number of saddle connections of length at most L joining, say, the first two zeros is quadratic in L and the leading term of the asymptotics (up to a factor of π to ensure rationality) is called the saddle connection Siegel-Veech constant. Intuitively, the saddle connection Siegel-Veech constant should be proportional to the cone angles around the two concerned zeros. For quadratic differentials this is not correct as shown by Athreya, Eskin and Zorich [START_REF] Athreya | Right-angled billiards and volumes of moduli spaces of quadratic differentials on CP 1[END_REF]. Nevertheless as an application of our formulas, we show that for Abelian differentials the intuitive expectation indeed holds, if we use a minor modification c hom 1↔2 (μ) of the Siegel-Veech constant counting homologous saddle connections only once. An overview about the variants of Siegel-Veech constants is given in Sect. 7.

Theorem 1. [START_REF] Aggarwal | Large genus asymptotics for volumes of strata of abelian differentials. With an appendix by A. Zorich[END_REF] The saddle connection Siegel-Veech constant c hom 1↔2 (μ) joining the first and the second zeros on a generic flat surface of type μ is given by c hom 1↔2 (μ) = (m 1 + 1)(m 2 + 1) . [START_REF] Bainbridge | The moduli space of multi-scale differentials[END_REF] When the stratum is disconnected, we also show that the theorem holds for each connected component under Assumption 6.1. We remark that as an asymptotic equality as g tends to infinity, the formula [START_REF] Bainbridge | The moduli space of multi-scale differentials[END_REF] for the entire stratum was previously shown in the appendix by Zorich to [START_REF] Aggarwal | Large genus asymptotics for volumes of strata of abelian differentials. With an appendix by A. Zorich[END_REF] for saddle connections of multiplicity one and by Aggarwal [START_REF] Aggarwal | Large genus asymptotics for Siegel-Veech constants[END_REF] for all multiplicities.

Another important kind of Siegel-Veech constants is the area Siegel-Veech constant, which counts cylinders (weighted by the reciprocal of their areas) on flat surfaces and is related to the sum of Lyapunov exponents [START_REF] Eskin | Sum of Lyapunov exponents of the Hodge bundle with respect to the Teichmüller geodesic flow[END_REF] (see Sect. 7 for the definition of area Siegel-Veech constants). We similarly establish an intersection formula for area Siegel-Veech constants.

Theorem 1. [START_REF] Aggarwal | Large genus asymptotics for Siegel-Veech constants[END_REF] The area Siegel-Veech constants of the strata can be evaluated as c area (μ) = -1 4π 2

P M g,n (μ) β i • δ 0 P M g,n (μ) β i • ξ (8)
for each 1 ≤ i ≤ n, where δ 0 is the divisor class of the locus of curves with a non-separating node.

Theorem 1.4 completes the investigation of area Siegel-Veech constants begun in [START_REF] Chen | Quasimodularity and large genus limits of Siegel-Veech constants[END_REF]Section 4] (for the principal strata) and [START_REF] Sauvaget | Volumes and Siegel-Veech constants of H(2g-2) and Hodge integrals[END_REF]Equation (2)] (for the minimal strata). It also justifies a speculation of Kontsevich [START_REF] Kontsevich | Lyapunov exponents and Hodge theory[END_REF]Section 7] about the existence of such a β-class for computing area Siegel-Veech constants and sums of Lyapunov exponents of the strata.

Another application of the volume recursion is a geometric proof of the large genus limit conjecture by Eskin and Zorich [START_REF] Eskin | Volumes of strata of Abelian differentials and Siegel-Veech constants in large genera[END_REF] for the volumes of the strata and area Siegel-Veech constants. A proof using direct combinatorial arguments was given by Aggarwal [START_REF] Aggarwal | Large genus asymptotics for volumes of strata of abelian differentials. With an appendix by A. Zorich[END_REF][START_REF] Aggarwal | Large genus asymptotics for Siegel-Veech constants[END_REF]. Our proof, in addition, gives a uniform expression for the second order term as conjectured in [START_REF] Sauvaget | Volumes and Siegel-Veech constants of H(2g-2) and Hodge integrals[END_REF] (see Sect. 11).

Theorem 1.5 [START_REF] Eskin | Volumes of strata of Abelian differentials and Siegel-Veech constants in large genera[END_REF]Main Conjectures] Consider the strata M g,n (μ) such that all the entries of μ are positive. Then

v(μ) = 4 - 2π 2 3 • n i=1 (m i + 1) + O(1/g 2 ), c area (μ) = 1 2 - 1 2 n i=1 (m i + 1) + O(1/g 2 ),
where the implied constants are independent of μ and g.

Finally we settle another conjecture of Eskin and Zorich on the asymptotic comparison of spin components. where the implied constants are independent of μ and g.

Further directions

Our work opens an avenue to study a series of related questions. First, we point out an interesting comparison with the proofs by Mirzakhani [START_REF] Mirzakhani | Weil-Petersson volumes and intersection theory on the moduli space of curves[END_REF], by Kontsevich [START_REF] Kontsevich | Intersection theory on the moduli space of curves and the matrix Airy function[END_REF], and by Okounkov-Pandharipande [START_REF] Okounkov | Gromov-Witten theory, Hurwitz numbers, and matrix models[END_REF] of Witten's conjecture: the generating function of ψ-class intersections on moduli spaces of curves is a solution of the KdV hierarchy of partial differential equations. Mirzakhani considered the Weil-Petersson volumes of moduli spaces of hyperbolic surfaces and analyzed geodesics that bound pairs of pants, while we consider the Masur-Veech volumes of moduli spaces of flat surfaces and analyze geodesics that join two zeros (i.e. saddle connections). Kontsevich interpreted ψ-classes as associated to certain polygon bundles, while we have the interpretation of Abelian differentials as polygons. Okounkov and Pandharipande used Hurwitz spaces of P 1 covers, while we rely on Hurwitz numbers of torus covers. Therefore, we speculate that generating functions of Masur-Veech volumes and area Siegel-Veech constants should also satisfy a certain interesting hierarchy as in Witten's conjecture.

In another direction, one can consider saddle connections joining a zero to itself (see [START_REF] Eskin | Moduli spaces of Abelian differentials: the principal boundary, counting problems, and the Siegel-Veech constants[END_REF]Part 2]) or impose other specific configurations to refine the Siegel-Veech counting (see e.g. the appendix by Zorich to [START_REF] Aggarwal | Large genus asymptotics for volumes of strata of abelian differentials. With an appendix by A. Zorich[END_REF]). From the viewpoint of intersection theory, such a refinement should pick up the corresponding part of the principal boundary when flat surfaces degenerate along the configuration, hence we expect that the resulting Siegel-Veech constant can be described similarly by a recursion formula involving intersection numbers.

One can also investigate volumes and Siegel-Veech constants for affine invariant manifolds (i.e. SL 2 (R)-orbit closures in the strata). It is thus natural to seek intersection theoretic interpretations of these invariants for affine invariant manifolds, e.g. the strata of quadratic differentials (see [START_REF] Andersen | Topolocial recursion for Masur-Veech volumes[END_REF][START_REF] Chen | Masur-Veech volumes and intersection theory: the principal strata of quadratic differentials[END_REF][START_REF] Delecroix | Masur-Veech volumes, frequencies of simple closed geodesics and intersection numbers of moduli spaces of curves[END_REF] for interesting related results in the case of the principal strata). We plan to treat these questions in future work.

Organization of the paper

In Sect. 2 we introduce relevant intersection numbers on Hurwitz spaces of P 1 covers that will appear as coefficients in the volume recursion. In Sect. 3 we prove that the expression of volumes by intersection numbers satisfies the recursion in [START_REF] Andersen | Topolocial recursion for Masur-Veech volumes[END_REF], thus showing the equivalence of Theorems 1.1 and 1.2. In Sect. [START_REF] Aggarwal | Large genus asymptotics for Siegel-Veech constants[END_REF] we exhibit another recursion of volumes by using the algebra of shifted symmetric functions and cumulants. In Sect. 5 we show that the two recursions are equivalent by interpreting them as the same summation over certain oriented graphs, thus completing the proof of Theorems 1.1 and 1.2. In Sect. [START_REF] Bainbridge | Compactification of strata of abelian differentials[END_REF] we refine the results for the spin and hyperelliptic components of the strata. In Sects. 7, 8 and 9 we respectively review the definitions of various Siegel-Veech constants, prove Theorem 1.3 regarding saddle connection Siegel-Veech constants and interpret the result from the perspective of Hurwitz spaces of torus covers. In Sect. 10 we establish similar intersection and recursion formulas for area Siegel-Veech constants, thus proving Theorem 1.4. Finally in Sect. 11 we apply our results to evaluate large genus limits of volumes and area Siegel-Veech constants, proving Theorems 1.5 and 1.6.

Hurwitz spaces of P 1 covers

In this section we recall the definition of the moduli space of admissible covers of [START_REF] Harris | On the Kodaira dimension of the moduli space of curves[END_REF] as a compactification of the classical Hurwitz space (see also [START_REF] Harris | Moduli of Curves[END_REF]), and prove formulas to compute recursively intersection numbers of ψ-classes on these moduli spaces. These intersection numbers will appear as coefficients and multiplicities in the volume recursion. Along the way we introduce basic notions on stable graphs and level functions.

Hurwitz spaces and admissible covers

Let d, g, and g be non-negative integers. Let = (μ (1) , • • • , μ (n) ) be a ramification profile consisting of n partitions. We define the Hurwitz space H d,g,g ( ) to be the moduli space parametrizing branched covers of smooth connected curves p : X → Y of degree d with profile and such that the genera of X and Y are given by g and g respectively. That is, p is ramified over n points and over the ith branch point the sheets coming together form the partition μ (i) (completed by singletons if |μ (i) | < deg( p)).

The Hurwitz space H d,g,g ( ) has a natural compactification H d,g,g ( ) parametrizing admissible covers. An admissible cover p : X → Y is a finite morphism of connected nodal curves such that 123 (i) the smooth locus of X maps to the smooth locus of Y and the nodes of X map to the nodes of Y , (ii) at each node of X the two branches have the same ramification order, and (iii) the target curve Y marked with the branch points is stable.

The space H d,g,g ( ) is equipped with two forgetful maps H d,g,g ( )

f S f T M g,m M g ,n
obtained by mapping an admissible cover to the stabilization of the source or the target. Here n denotes the number of branch points or equivalently the length of and m denotes the number of ramification points or equivalently the number of parts (of length > 1) of all the μ (i) . The Hurwitz number N • d,g,g ( ) is the degree of the map f T , or equivalently the number of connected covers p : X → Y of degree d with profile and the location of the branch points fixed in Y . We also denote by N d,g,g ( ) the Hurwitz number of covers without requiring X to be connected. We remark that each cover is counted with weight given by the reciprocal of the order of its automorphism group, as is standard for the Hurwitz counting problem.

Intersection of ψ-classes on Hurwitz spaces

From now on in this section we will consider the special case g = 0.

Let μ[0] = (m 1 , . . . , m n ) be a list of non-negative integers and μ[∞] = ( p 1 , . . . , p k ) a list of positive integers. We consider the Hurwitz space with profile given by μ (i) = (m i + 1) for i ≤ n and μ (n+1) = ( p 1 , . . . , p k ) such that d = k i=1 p i , i.e. we consider

H P 1 (μ[0], μ[∞]) = H d,g,0 ((m 1 + 1), . . . , (m n + 1), ( p 1 , . . . , p k )) .
By the Riemann-Hurwitz formula, the genus g of the covering surfaces satisfies that

2 -2g = k + d - n i=1 m i = k + k i=1 p i - n i=1 m i .
The forgetful map f S goes from

H P 1 (μ[0], μ[∞]) to M g,n+k
, where we assume that the first n marked points are the first n ramification points and the preimages of the last branch point are the k last marked points. Since there are n + 1 branch points in the target surface of genus zero, we conclude that dim

H P 1 (μ[0], μ[∞]) = dim M 0,n+1 = n -2.
For g = 0, define the following intersection numbers on the Hurwitz spaces

h P 1 (μ[0], μ[∞]) = H P 1 (μ[0],μ[∞]) f * S n i=3 ψ i . ( 9
)
The definition of ψ-classes will be recalled in Sect. 3. If n = 2, then

H P 1 (μ[0], μ[∞]
) is of dimension zero, and hence the intersection number on the right is just the number of points of the Hurwitz space, i.e.

h P 1 ((m 1 , m 2 ), μ[∞]) = N • d,0,0 ((m 1 + 1), (m 2 + 1), ( p 1 , . . . , p k )) .
Again we emphasize that the Hurwitz number on the right-hand side is counted with weight 1/|Aut| for each cover. Correspondingly the intersection numbers are computed on the Hurwitz space treated as a stack. Our goal for the rest of the section is to show the following result.

Proposition 2.1 For n = 2, h P 1 ((m 1 , m 2 ), μ[∞]
) can be computed by the coefficient extraction

h P 1 ((m 1 , m 2 ), μ[∞]) = (k -1)![t m 1 +1 ] k i=1 t -t p i +1 1 -t ,
and for n ≥ 3,

h P 1 (μ[0], μ[∞]
) can be computed recursively by the sum

h P 1 (μ[0], μ[∞]) = ∈RT(μ[0],μ[∞]) 1,2 h( ) over rooted trees.
The definitions of rooted trees and the local contributions h( ) are given in Sects. 2.3 and 2.4 respectively. The above formula for the Hurwitz number obviously agrees with [START_REF] Bainbridge | Compactification of strata of abelian differentials[END_REF]. Proof of Proposition 2.1, case n = 2. Let S d be the symmetric group acting on [ [1, d]

] = {1, . . . , d}, where d = p 1 + • • • + p k = m 1 + m 2 + 2 -k. Define the set of Hurwitz tuples A(m 1 , m 2 , μ[∞]) = {(σ 1 , σ 2 , σ ∞ )} ⊂ S d × S d × S d such that

123

• the permutation σ ∞ is in the conjugacy class of ( p 1 , . . . , p k ) and we fix a bijection of its cycles with [ [1, k]] such that the ith cycle has length p i , • the partitions σ 1 and σ 2 are cycles of order m 1 + 1 and m 2 + 1 respectively, • the relation σ 1 • σ 2 = σ ∞ holds, and • the group generated by σ 1 , σ 2 and σ ∞ acts transitively on [ [1, d]].

Then the (weighted) Hurwitz number

h P 1 ((m 1 , m 2 ), μ[∞]) = |A(m 1 , m 2 , μ[∞])|/d!.
The second and third conditions above imply that the union of the supports of the cycles σ 1 and σ 2 is [ [1, d]]. Therefore, σ 1 and σ 2 contain exactly (m 1 + 1) + (m 2 + 1)d = k common elements. We can write

σ 1 = (a 1 , . . . , a i 1 -1 , c 1 ; a i 1 +1 , . . . , a i 2 -1 , c 2 ; . . . ; a i k-1 +1 , . . . , a i k -1 , c k ), where 1 ≤ i 1 < i 2 < • • • < i k = m 1 + 1 and c 1 , . . . , c k are the common elements of σ 1 and σ 2 . Since σ 1 • σ 2 = σ ∞ is of conjugacy type ( p 1 , . . . , p k ),
it is easy to see that σ 2 must be of the form

σ 2 = (b 1 , . . . , b j 1 -1 , c k ; b j 1 +1 , . . . , b j 2 -1 , c k-1 ; . . . ; b j k-1 +1 , . . . , b j k -1 , c 1 ), for certain b i , such that { j 1 + i 1 , ( j 2 -j 1 ) + (i k -i k-1 ), . . . , ( j k -j k-1 ) + (i 2 -i 1 )} = {p 1 + 1, p 2 + 1, . . . , p k + 1} . If τ is a permutation on [[1, k]] such that ( j k+1--j k-) + (i +1 -i ) = p τ ( ) + 1, then we have 1 ≤ i +1 -i ≤ p τ ( ) .

Conversely, such τ and i-indices determine the j-indices.

There are m 1 ! d m 1 +1 choices for σ 1 . Fixing σ 1 , to construct σ 2 we first choose • a permutation τ ∈ S k , and then

• a partition (i 1 , i 2 -i 1 , . . . , i k -i k-1 ) of m 1 + 1 such that 1 ≤ i +1 -i ≤ p τ ( ) for all . This gives k![t m 1 +1 ] k i=1 (t +• • •+t p i ) choices.
Choose c 1 out of the elements in σ 1 , which gives m 1 + 1 choices. Along with the i-indices this determines the elements c 2 , . . . , c k as well as the set of b's as the complement of the union of a's and c's. Finally σ 2 is determined by arranging the b's, which gives (m 2 + 1k)! choices. Note that in this process only the cyclic order of (c 1 , . . . , c k ) matters and we cannot actually determine which one is the first c, hence we need to divide the final count by k.

In summary, we conclude that

|A(m 1 , m 2 , μ[∞])| 123 = (m 1 + 1)!(m 2 + 1 -k)! d m 1 + 1 • (k -1)![t m 1 +1 ] k i=1 (t + • • • + t p i ). Since d = m 1 + m 2 + 2 -k, we obtain h P 1 ((m 1 , m 2 ), μ[∞]) = |A(m 1 , m 2 , μ[∞])|/d! = (k -1)![t m 1 +1 ] k i=1 t (1 -t p i ) 1 -t using that (m 1 + 1)! (m 2 + 1 -k)! d m 1 +1 = d!.
We remark that the above Hurwitz counting problem can also be interpreted by the angular data of the configurations of saddle connections joining two zeros z 1 and z 2 of order m 1 and m 2 respectively in the setting of [START_REF] Eskin | Moduli spaces of Abelian differentials: the principal boundary, counting problems, and the Siegel-Veech constants[END_REF]. Suppose f : P 1 → P 1 is a branched cover parameterized in the Hurwitz space H d,0,0 ((m 1 + 1), (m 2 + 1), ( p 1 , . . . , p k )), where we treat f as a meromorphic function with k poles of order p 1 , . . . , p k . Then the meromorphic differential η = d f has two zeros of order m 1 and m 2 as well as k poles of order p 1 + 1, . . . , p k + 1 with no residue. Conversely given such η, integrating η gives rise to a desired branched cover f . Such η can be constructed using flat geometry as in [START_REF] Chen | Principal boundary of moduli spaces of abelian and quadratic differentials[END_REF]Section 2.4]. In particular, it is determined by the angles 2π(a i + 1) between the saddle connections (clockwise) at z 1 and the angles 2π(a i + 1) (counterclockwise) at z 2 , such that k i=1 (a i + 1) = m 1 + 1, k i=1 (a i + 1) = m 2 + 1 and a i + a i + 2 = p i + 1. We see again that the choices involve a partition (a 1 + 1, . . . , a k + 1) of m 1 + 1 such that a i + 1 ≤ p i for all i.

Level graphs and rooted trees

The boundary of the Deligne-Mumford compactification M g,n is naturally stratified by the topological types of the stable marked surfaces. These boundary strata are in one-to-one correspondence with stable graphs, whose definition we recall below. The boundary strata of Hurwitz spaces and of moduli spaces of Abelian differentials are encoded by adding level structures and twists to stable graphs.

Definition 2.2 A stable graph is the datum of

= (V, H, g : V → N, a : H → V, i : H → H, E, L [[1, n]])
satisfying the following properties:

123

• V is a vertex set with a genus function g; • H is a half-edge set equipped with a vertex assignment a and an involution i (and we let n(v) = |a -1 (v)|); • E, the edge set, is defined as the set of length-2 orbits of i in H (self-edges at vertices are permitted); • (V, E) define a connected graph; • L is the set of fixed points of i, called legs or markings, and is identified with [ [1, n]]; • for each vertex v, the stability condition 2g(v) -2 + n(v) > 0 holds. Let v( ) and e( ) denote the cardinalities of V and E respectively. The genus of is defined by v∈V ( ) g(v) + e( ) -v( ) + 1.

We denote by Stab(g, n) the set of stable graphs of genus g and with n legs. A stable graph is said of compact type if h 1 ( ) = 0, i.e. if the graph has no loops, which is thus a tree.

We will use two extra structures on stable graphs, called level functions and twists. As in [START_REF] Bainbridge | Compactification of strata of abelian differentials[END_REF] we define a level graph to be a stable graph together with a level function : V ( ) → R ≤0 . An edge with the same starting and ending level is called a horizontal edge. A bi-colored graph is a level graph with two levels (in which case we normalize the level function to take values in {0, -1}) that has no horizontal edges. We denote the set of bi-colored graphs by Bic(g, n).

Recall the notation μ[0] = (m 1 , . . . , m n ) and μ[∞] = ( p 1 , . . . , p k ) where m i ≥ 0 and p j > 0 for all i and j.

Definition 2.3 Let be a stable graph in Stab(g, n +k). A twist assignment on of type (μ[0], μ[∞]

) is a function p : H ( ) → Z satisfying the following conditions:

• If (h, h ) is an edge, then p(h) + p(h ) = 0. • For all 1 ≤ i ≤ n, the twist of the ith leg is m i + 1 and for all 1 ≤ i ≤ k the twist of the (n + i)th leg is -p i . • For all vertices v of 2g(v) -2 + n(v) = h∈L , a(h)=v p(h) .
Suppose the graph comes with a level structure . We say that a twist p is compatible with the level structure if for all edges (h, h ) the condition p(h) > 0 implies that (a(h)) > (a(h )), and respectively for the cases < and =. In this case we call the triple ( , , p) a twisted level graph. For the reader familiar with related results of compactifications of strata of Abelian differentials, the above definition characterizes twisted differentials (or canonical divisors) in [START_REF] Bainbridge | Compactification of strata of abelian differentials[END_REF] and [START_REF] Farkas | The moduli space of twisted canonical divisors[END_REF] (regarding the m i as the zero orders and p j + 1 as the pole orders of twisted differentials on irreducible components of the corresponding stable curves). In particular, every level graph has only finitely many compatible twists. For a graph of compact type, there exists a unique twist p if the entries of

(μ[0], μ[∞]) satisfy the condition that n i=1 m i -k j=1 ( p j +1) = 2g -2. Definition 2.4 Let 1 ≤ i < j ≤ n.
A stable rooted tree (or simply a rooted tree) is a twisted level graph ( , , p) of compact type satisfying the following conditions:

(i) One vertex v j carries the ith and jth legs and no other of the first n legs;

(The vertex v j is called the root, a vertex on the path from v to the root is called an ancestor of v, and a vertex whose ancestors contain v is called a descendant of v.) (ii) There are no horizontal edges;

(iii) A vertex v is on level 0 if and only if v is a leaf. If v is not a leaf, then (v) = min{ (v ) | v is a descendant of v} -1;
(iv) All vertices of positive genus are leaves and hence on level 0; (v) Each vertex of genus zero other than v j carries exactly one of the first n legs.

Since the root is an ancestor of any other vertex, by definition it is the unique vertex lying on the bottom level, hence it has genus zero. Moreover, it is easy to see from the definition that any path towards the root is strictly going down (Fig. 1). We denote by RT(g, μ[0], μ[∞]) i, j the set of such rooted trees, and sometimes simply by RT(μ

[0], μ[∞]) i, j if g = 0.

The sum over rooted trees

Now we assume that g = 0. Below we define the local contributions from rooted trees in Proposition 2.1 and complete its proof. Consider a graph ∈ RT(μ[0], μ[∞]) 1,2 . Since by assumption every vertex of has genus zero, condition v) implies that has exactly n -1 vertices and n -2 edges. Denote by v 2 , . . . , v n the vertices of such that v i carries the ith leg h i for 3 ≤ i ≤ n and v 2 carries the first two legs. This convention is consistent with our previous notation for the root. We denote by μ[∞] i the list of negative twists at halfedges adjacent to v i . These half-edges are either part of the whole edges joining v i to its descendants (as adjacent vertices to v i on higher level) or part of the k last legs (corresponding to the k marked poles).

If i = 2, then there is a unique (non-leg) half-edge h i = h i adjacent to v i such that m i := p( h i ) -1 ≥ 0. Namely, this half-edge is part of the whole edge joining v i to its ancestor (as the adjacent vertex to v i on lower level). With this notation we define the contribution of the rooted tree as

h( ) = h P 1 ((m 1 , m 2 ), μ[∞] 2 ) • n i=3 h P 1 ((m i , m i ), μ[∞] i ) . ( 10
) Let J ⊂ [[1, n + k]]
be a subset such that the cardinalities of J and J c are at least two. Denote by δ J the class of the boundary divisor of M 0,n+k parameterizing curves that consist of a component with the markings in J union a component with the markings in J c . We need the following classical result (see e.g. [START_REF] Arbarello | Geometry of Algebraic Curves[END_REF]Lemma 7.4]).

Lemma 2.5 For all 3 ≤ i ≤ n + k, the following relation of divisor classes holds on M 0,n+k :

ψ i = i∈J ⊂[[3,n+k]] δ J .
If J is a subset of [ [3, n]], then we denote by δ J = J ⊂[[n+1,n+k]] δ J ∪J . For 3 ≤ i ≤ n, the above lemma implies that

ψ i = i∈J ⊂[[3,n]] δ J . ( 11
)
We also need the following result about the boundary divisors of the Hurwitz spaces

H P 1 (μ[0], μ[∞]).

Lemma 2.6 There is a bijection between the boundary divisors of H P 1 (μ[0], μ[∞]) and the corresponding bi-colored graphs (i.e. with two levels only). Moreover, a boundary divisor drops dimension under the source map f S if its bi-colored graph has more than two vertices.

Proof The first part of the claim follows from the same argument as in [START_REF] Sauvaget | Cohomology classes of strata of differentials[END_REF]Proposition 7.1]. Here the vertices of level 0 in the bi-colored graphs correspond to the stable components of the admissible covers that contain the marked poles. For the other part, suppose that a generic point of a boundary divisor has at least two vertices on level 0 (or on level -1). Then one can scale one of the two functions that induce the covers on the two vertices such that the domain marked curve is fixed while the admissible covers vary. It implies that f S restricted to this boundary divisor has positive dimensional fibers.

Proof of Proposition 2.1, case n ≥ 3. We will prove the result by induction on n. The initial case n = 2 follows from the definition of h( ) in [START_REF] Chen | Principal boundary of moduli spaces of abelian and quadratic differentials[END_REF] and we have also described it explicitly in Sect. 2.2. The strategy of the induction for higher n is by successively replacing the ψ i in [START_REF] Bloch | The character of the infinite wedge representation[END_REF] with the sum over boundary divisors as in the preceding lemma, starting with i = 3. To simplify notation, we write

H (μ[0], μ[∞]) instead of H P 1 (μ[0], μ[∞]
). We also simply write ψ and δ as classes in the Hurwitz space for their pullbacks via f S .

Consider a boundary divisor δ J of M 0,n+k pulled back to

H (μ[0], μ[∞]), which is a union of certain boundary divisors of H (μ[0], μ[∞]
). We would like to compute the intersection number

δ J • n i=4 ψ i on H (μ[0], μ[∞]
). By Lemma 2.6 and the projection formula, the only possible non-zero contribution is from the loci in δ J whose bicolored graphs have a unique edge e = (h, h ) connecting two vertices v 0 and v -1 on level 0 and level -1 respectively, such that the last k markings (i.e. the k marked poles) are contained in v 0 . In this case we can assume that p(h) > 0 (and hence p(h ) < 0 as p(h) + p(h ) = 0 by definition). The admissible covers restricted to v 0 and to v -1 belong to Hurwitz spaces of similar type, where at the node (i.e. the edge e) the ramification order of the restricted maps is given by p(h) -1 = -p(h ) -1. It implies that the locus of such admissible covers can be identified with

H (μ[0] 0 , μ[∞] 0 ) × H (μ[0] -1 , μ[∞] -1 ), where μ[0] 0 is the part of μ[0] contained in v 0 union with p(h) -1, μ[∞] 0 = μ[∞], μ[0] -1 is the part of μ[0] contained in v -1 ,
and μ[∞] -1 has a single entry p(h).

By Eq. [START_REF] Chen | Masur-Veech volumes and intersection theory: the principal strata of quadratic differentials[END_REF] we have

ψ 3 = 3∈J ⊂[[3,n+k]] δ J . Fix a subset J ⊂ [[3, n + k]] such that 3 ∈ J .
If the third marking belongs to v -1 , we claim that

H (μ[0], μ[∞]) • δ J • n i=4 ψ i = 0 .
To see this, note that the dimension of

H (μ[0] -1 , μ[∞] -1 ) is equal to n(μ[0] -1 )-2, as μ[∞] -1
has only one entry. However there are n(μ[0] -1 )-1 markings with label ≥ 4 on v -1 . Consequently the intersection with the product of those ψ i vanishes on

H (μ[0] -1 , μ[∞] -1 ).
Therefore, we only need to consider the case when v 0 contains the third marking (hence all markings labeled by J ), and consequently v -1 contains the first and second markings (hence all markings labeled by J c ). In this case we obtain that

H (μ[0], μ[∞]) • δ J • n i=4 ψ i = h P 1 (μ[0] 0 , μ[∞] 0 ) • h P 1 (μ[0] -1 , μ[∞] -1 ),
where h P 1 is defined in [START_REF] Bloch | The character of the infinite wedge representation[END_REF], where in the first factor on the right-hand side the ψ-product skips the third marking and the marking from the half-edge of v 0 , and where in the second factor the ψ-product skips the first and second markings. Now we use the induction hypothesis to decompose the factors h P 1 (μ[0] a , μ[∞] a ) for a = 0 and a = -1. It leads to a sum over all possible pairs of rooted trees, where the two rooted trees in each pair generate a new rooted tree. More precisely, one rooted tree in the pair contains the markings of J c ∪ {h } whose root v 2 carries the first and second markings, the other rooted tree contains the markings in J ∪ {h} whose root v 3 carries the third marking and h, and they generate a new rooted tree by gluing the legs h and h as a whole edge and by using v 2 as the new root.

Therefore, if J is a subset of [ [3, n]] such that 3 ∈ J , then we obtain that

H (μ[0], μ[∞]) • δ J • n i=4 ψ i = ∈RT(μ[0],μ[∞]) 1,2 , j∈J ⇔ (v 3 )≤ (v j ) h( ),
where the sum is over all rooted trees such that the descendants of v 3 are exactly the vertices v j for j ∈ J \{3}.

In summary if we write ψ 3 = 3∈J ⊂[ [3,n]] δ J , then by the above analysis we thus conclude that

H (μ[0], μ[∞]) • ψ 3 • n
i=4 ψ i is equal to the sum of the contributions h( ) over all rooted trees .

Volume recursion via intersection theory

In this section we show that the two main theorems of the introduction, Theorems 1.1 and 1.2 are equivalent. This section does not yet provide a direct proof of either of them.

We first show that the intersection numbers in Theorem 1.1 are given by a recursion formula of the same shape as in Theorem 1.2. Together with an agreement on the minimal strata this proves the equivalence of the two theorems. Along the way we introduce special classes of stable graphs that are used for recursions throughout the paper.

Intersection numbers on the projectivized Hodge bundle

Fix g and n such that 2g -2 + n > 0. We denote by f : X → M g,n the universal curve and by ω X /M g,n the relative dualizing line bundle. We will use the following cohomology classes:

• Let 1 ≤ i ≤ n. We denote by σ i : M g,n → X the section of f corresponding to the ith marked point and by L i = σ * i ω X /M g,n the cotangent line at the ith marked point. With this notation, we define

ψ i = c 1 (L i ) ∈ H 2 (M g,n , Q). • For 1 ≤ i ≤ g, we denote by λ i = c i ( M g,n ) ∈ H 2 (M g,n , Q) the ith
Chern class of the Hodge bundle. (We use the same notation for a vector bundle and its total space.) • We denote by δ 0 ∈ H2 (M g,n , Q) the Poincaré-dual class of the divisor parameterizing marked curves with at least one non-separating node. • The projectivized Hodge bundle P M g,n comes with the universal line bundle class ξ = c 1 (O( 1)) ∈ H 2 (P M g,n , Q).

Unless otherwise specified, we denote by the same symbol a class in H * (M g,n , Q) and its pull-back via the projection p : P M g,n → M g,n .

Recall that the splitting principle implies that the structure of the cohomology ring of the projectivized Hodge bundle is given by

H * (P M g,n , Q) H * (M g,n , Q)[ξ ]/(ξ g + λ 1 ξ g-1 + • • • + λ g ) .
Let μ = (m 1 , . . . , m n ) be a partition of 2g -2. We denote by P M g,n (μ) the closure of the projectivized stratum P M g,n (μ) inside the total space of the projectivized Hodge bundle P M g,n . This space is called the (ordered) incidence variety compactification 2 .

In this section we study the intersection numbers

a i (μ) = P M g,n (μ) β i • ξ = 1 m i + 1 P M g,n (μ) ξ 2g-1 • j =i ψ j ( 12
)
for all 1 ≤ i ≤ n. The reader should think of the a i (μ) as certain normalization of volumes. In fact, Theorem 1.1 can be reformulated as

vol( M g,n (μ)) = 2(2π) 2g (-1) g (2g -3 + n)! a i (μ), (13)
implying in particular that a i (μ) is independent of i.

We prove a collection of properties defining recursively the a i (μ) as the coefficients of some formal series. As the base case for n = 1, i.e. μ = (2g-2), define the formal series

A(t) = 1 t + g≥1 (2g -1) 2 a 1 (2g -2)t 2g-1 ∈ 1 t Q[[t]] (14)
and set

B(z) := z/2 sinh(z/2) =: j≥0 b j z j = 1 - z 2 24 + 7z 4 5760 + • • • . ( 15
)
For a partition μ, recall that n(μ) denotes the number of its entries and |μ| denotes the sum of the entries.

Theorem 3.1 The generating function A of the minimal stratum intersection numbers a i (2g -2) is determined by the coefficient extraction identity

[t 0 ] 1 j! A(t) j = b j , ( 16
)
while the intersection numbers a(μ) = a i (μ) with n(μ) ≥ 2 are given recursively by

(m 1 + 1)(m 2 + 1)a(m 1 , . . . , m n ) = min k=1 1 k! g,μ h P 1 ((m 1 , m 2 ), p) k j=1 (2g j -1 + n(μ j )) p j a( p j -1, μ j ), ( 17
)
where min = min(m 1 + 1, m 2 + 1), with the same summation conventions as in Theorem 1.2.

The first identity [START_REF] Eskin | Asymptotic formulas on flat surfaces[END_REF] was proved in [START_REF] Sauvaget | Volumes and Siegel-Veech constants of H(2g-2) and Hodge integrals[END_REF] and gives

A(t) = 1 t - 1 24 t + 3 640 t 3 - 1525 580608 t 5 + 615881 199065600 t 8 -• • • .
By Lagrange inversion, this formula can be written equivalently as

A(t) = 1 Q -1 (t)
, where

Q(u) = u exp k≥1 (k -1)!b k u k
and will in fact be proved in this form in Sect. 4.4. We observe in passing that Q(u) is the asymptotic expansion of ψ(u -1 + 1 2 ) as u → 0, where ψ(x) = (x)/ (x) is the digamma function. The proof of the second identity [START_REF] Eskin | Invariant and stationary measures for the SL(2, R) action on moduli space[END_REF] will be completed by the end of Sect. 3.5.

In the course of proving Theorem 3.1 we will prove the following complementary result, justifying the implicitly used fact that a i (μ) is independent of i. Proposition 3.2 For all 1 ≤ i ≤ n, we have a i (μ) = -

P M g,n (μ) ξ 2g-2 • n j=1 ψ j .

Boundary components of moduli spaces of Abelian differentials

In Sect. 2 we introduced several families of stable graphs to describe the boundary of Hurwitz spaces. Here we show how these graphs encode relevant parts of the boundary of moduli spaces of Abelian differentials.

The recursions in Theorems 1.2 and 3.1 can be phrased as sums over a small subset of twisted level graphs, with only two levels and more constraints, that we call (rational) backbone graphs, inspired by Fig. 2.

Recall that a bi-colored graph is a level graph with two levels {0, -1} that has no horizontal edges.

Definition 3.3

An almost backbone graph is a bi-colored graph with only one vertex at level -1. For such a graph to be a (rational) backbone graph we require moreover that it is of compact type and that the vertex at level -1 has genus zero.

We denote by BB(g, n) ⊂ ABB(g, n) ⊂ Bic(g, n) the sets of backbone, almost backbone and bi-colored graphs. We denote by BB(g, n) 1,2 ⊂ BB(g, n) the set of backbone graphs such that the first and second legs are adjacent to the vertex of level -1. Moreover, let BB(g, n) 1,2 ⊂ BB(g, n) 1,2 be the subset where precisely the first two legs are adjacent to the lower level vertex. Similarly, we define ABB(g, n) 1,2 and ABB(g, n) 1,2 and drop (g, n) if there is no source of confusion.

We remark that the backbone graphs will play an important role here, while the graphs in ABB(g, n) appear only in the Hurwitz space interlude in Sect. 9. We fix some notations for these graphs, used throughout in the sequel. For ∈ BB(g, n) we denote by v -1 the vertex of level -1. Given a partition μ = (m 1 , . . . , m n ) of 2g -2, let p be the unique twist of type (μ, ∅) for (see

X 5 X 4 X 3 X 2 X 1 = -1 = 0 X 5 X 3 X 2 X 1 X 4
Fig. 2 A backbone graph and the corresponding stable curve 123 Definition 2.3). We denote by μ[0] -1 the list of m i for all legs i at level -1 and with a slight abuse of notation we denote by p = ( p 1 , . . . , p k ) the list of p(h) for half-edges h that are adjacent to the k vertices of level 0. Said differently, the restriction of the twist to level -1 provides v -1 with a twist of type (μ[0] -1 , μ[∞] -1 = p). Finally if v is a vertex of level 0, we denote by μ v the list of p(h) -1 for all half-edges adjacent to v.

The goal in the remainder of the section is to introduce the classes α , , p in [START_REF] Eskin | Billiards in rectangles with barriers[END_REF] below that will be used in Proposition 3.11 to compute intersection numbers on P M g,n (μ). A stable graph ∈ Stab(g, n) determines the moduli space

M = v∈V ( ) M g(v),n(v)
and comes with a natural morphism ζ : M → M g,n . Let be a level function on such that ( , ) is a bi-colored graph with two levels {0, -1}. We define the following vector bundle

M , = ⎛ ⎝ v∈V ( ), (v)=0 M g(v),n(v) ⎞ ⎠ × ⎛ ⎝ v∈V ( ), (v)=-1 M g(v),n(v)
⎞ ⎠ over M . This space comes with a natural morphism ζ # , : M , → M g,n , defined by the composition M , → ζ * M g,n → M g,n where the first arrow is the inclusion of a vector sub-bundle and the second is the map on the Hodge bundles induced from ζ by pull-back. The morphism ζ # , determines a morphism (denoted by the same symbol) on the projectivized Hodge bundles ζ # , : P M , → P M g,n . The image of ζ # , is the closure of the locus of differentials supported on curves with dual graph such that the differentials vanish identically on components of level -1. In the sequel we will need the following lemma (see [START_REF] Sauvaget | Cohomology classes of strata of differentials[END_REF]Proposition 5.9]).

Lemma 3.4 The Poincaré-dual class of ζ #

, (P M , ) is divisible by [START_REF] Eskin | Isolation, equidistribution, and orbit closures for the SL(2, R) action on moduli space[END_REF] be the loci defined by the following three conditions:

ξ h 1 ( ) in H * (P M g,n , Q).
M 0 ⊂ v∈V ( ), (v)=0 M g(v),n(v) , 123 M -1 ⊂ v∈V ( ), (v)=-1 M g(v),n(v)
(i) A differential in M 0 has zeros of orders m i at the relevant marked points and of orders p(h) -1 at the relevant branches of the nodes. (ii) For each v of level -1 there exists a non-zero (meromorphic) differential ω v on the component X v corresponding to v that has zeros at the relevant marked points of orders prescribed by μ and poles at the relevant branches of the nodes of orders prescribed by p, i.e. such that the canonical divisor class of X v is given by h∈H,a(h)=v (p(h) -1)x h , where x h ∈ X v is the marked point or the node corresponding to the half-edge h. (iii) There exist complex numbers

k v = 0 for all vertices v of level -1 such that ω = (v)=-1 k v ω v satisfies the global residue condition of [6].
In particular for a backbone graph , since it is of compact type with a unique vertex v -1 of level -1, we have the identification

M 0 = v =v -1 M g(v),n(v) (μ v ).
We define P M p , as the Zariski closure of P M 0 × M -1 in P M , and define

α , ,p = ζ # , * [P M p , ] if dim(P M p , ) = dim(P M g,n (μ)) -1 0 otherwise ( 19
)
as the corresponding class in H * (P M g,n , Q). By [43, Proposition 5.9], we can describe α , ,p with the following lemma. ( ,,p) is a bi-colored graph of compact type, then α , ,p = 0 if and only if there is a unique vertex v -1 of level -1, and in this case α , ,p is divisible by

Lemma 3.5 If

ξ g(v -1 ) + ζ * (λ v -1 ,1 )ξ g(v -1 )-1 + • • • + ζ * (λ v -1 ,g(v -1 ) ),
where

λ v -1 ,i = (λ i , 1, . . . , 1) ∈ H * (M g(v -1 ),n(v -1 ) , Q) v∈V ( ),v =v -1 H * (M g(v),n(v) , Q) H * (M , Q) .
123

A first reduction of the computation

Recall the (marked and projectivized) Hodge bundle projection p : P M g,n → M g,n . As before we usually denote by the same symbol a class in M g,n and its pullback via p. In this section we show that many p-push forwards of intersections of α , , p with tautological classes vanish or can be computed recursively. The starting point is the following important lemma proved by Mumford in [START_REF] Mumford | Towards an enumerative geometry of the moduli space of curves[END_REF]Equation (5.4)].

Lemma 3.6

The Segre class of the Hodge bundle is the Chern class of the dual of the Hodge bundle, i.e.

c * ( M g,n ) • c * ( M ∨ g,n ) = 1 .
In particular, we have

λ 2 g = 0 ∈ H 4g (M g,n , Q).
Together with the definition of Segre class, this lemma implies that

p * (ξ k γ ) = s k-g+1 ( M g,n )γ = (-1) k-g+1 λ k-g+1 γ for all γ ∈ H * (M g,n , Q) and all k ≥ g -1.
Another important lemma is the following (see e.g. [1, Section 13, Equation (4.31)]).

Lemma 3.7 Let 1 ≤ k ≤ g and let be a stable graph. Then

ζ * λ k = (k v ) v∈V ∈N V , |(k v )|=k v∈V λ k v ,
where the sum is over all partitions of k into non-negative integers k v assigned to each vertex v ∈ V = V ( ).

In particular if h 1 ( ) > g -k, then v∈V g(v) = g -h 1 ( ) < k, hence the above lemma implies that ζ * λ k = 0 as there exists some k v > g(v) for any partition (k v ) v∈V of k.
As a consequence of the above discussion, we obtain the following result.

Lemma 3.8 Let α = i≥0 ξ i α i be a class in H * (P M g,n , Q) where the classes α i are pull-backs from H * (M g,n , Q). Then we have p * (ξ 2g-1 α) = (-1) g α 0 λ g , p * (ξ 2g-2 α) = (-1) g α 1 λ g + (-1) g-1 α 0 λ g-1 , p * (ξ 2g-2 δ 0 α) = (-1) g-1 α 0 δ 0 λ g-1 .
123

Recall the expressions of the intersection numbers a i (μ) in ( 12) and in Proposition 3.2. In order to compute a i (μ), by Lemma 3.8 we only need to consider the ξ -degree zero and one parts of the class

[P M g,n (μ)] in H * (P M g,n , Q).
Combining Lemmas 3.7 and 3.8 together with Lemmas 3.4 and 3.5 of the previous section, we can already prove the following vanishing result for classes associated with some bi-colored graphs. Proposition 3.9 If ( , , p) is not a backbone graph, then

p * ξ 2g-2 α , ,p = 0 ∈ H * (M g,n , Q)
where α , ,p is defined in [START_REF] Eskin | Billiards in rectangles with barriers[END_REF].

Proof For simplicity we write α = α , ,p in the proof. We assume first that is not of compact type, i.e. h 1 ( ) > 0. Then by Lemma 3.4, the class α is divisible by ξ . Note that the cohomology ring of a projective bundle is generated by the universal line bundle class with the classes pulled back from the base. Therefore, we can write

ξ 2g-2 α = i≥0 ξ 2g-1+i α i ,
where α i is a pullback from H * (M g,n , Q) that is supported on ζ (M ) for all i ≥ 0. Thus by Lemma 3.8, p * (ξ 2g-2 α) = (-1) g α 0 λ g = 0, because is not of compact type. Now we assume that is of compact type. By Lemma 3.5, we only need to consider the case when there is a unique vertex v 1 of level -1. Since is not a backbone graph, v 1 has positive genus g 1 . Still by Lemma 3.5 and simplifying the notation

ζ * (λ v 1 ,i ) by λ v 1 ,i , the class α is divisible by ξ g 1 + ξ g 1 -1 λ v 1 ,1 + • • • + λ v 1 ,g 1 . Consequently we can write α = (ξ λ v 1 ,g 1 -1 + λ v 1 ,g 1 )γ 0 + ξ λ v 1 ,g 1 γ 1 + O(ξ 2 ),
where γ 0 and γ 1 are pullbacks from H * (M g,n , Q) and the O(ξ 2 ) term stands for a class divisible by ξ 2 . By Lemma 3.8, we obtain that

p * (ξ 2g-2 α) = (-1) g (λ v 1 ,g 1 -1 λ g -λ v 1 ,g 1 λ g-1 )γ 0 + (-1) g λ v 1 ,g 1 λ g γ 1 .
Using Lemma 3.7, we also obtain that

ζ * (λ g ) = v∈V λ g v and ζ * (λ g-1 ) = v∈V λ g v -1 v =v λ g v .

123

From the projection formula we deduce that

λ v 1 ,g 1 • λ g = ζ * (λ v 1 ,g 1 • ζ * (λ g )) = ζ * λ 2 g 1 v =v 1 λ g v = 0, because λ 2 g 1 = 0 ∈ H * (M g(v 1
),n(v 1 ) , Q) by Lemma 3.6. Once again the same lemma implies that

λ v 1 ,g 1 -1 • λ g = ζ * λ g 1 λ g 1 -1 v =v 1 λ g v , λ v 1 ,g 1 • λ g-1 = ζ * λ g 1 λ g 1 -1 v =v 1 λ g v + v =v 1 v =v,v 1 λ 2 g 1 ⊗ λ g v -1 v =v λ g v = λ v 1 ,g 1 -1 • λ g .
Putting everything together, we thus conclude that p * (ξ 2g-2 α) = 0.

We define the multiplicity of a twist p to be

m(p) = (h,h )∈E( ) -p(h)p(h ) . ( 20
)
Proposition 3.10 If ( , , p) is a backbone graph in BB(g, n) 1,2 , then

P M g,n α , ,p • ξ 2g-1 • n i=3 ψ i = m(p) • h P 1 (μ -1 , p) • v∈V ( ) (v)=0 a 1 ( p v -1, μ v ),
where p v is the entry of p corresponding to the twist on the unique edge of each vertex v of level 0 and μ -1 is the list of entries in μ whose corresponding legs are adjacent to the vertex of level -1.

As a preparation for the proof we relate the space M -1 defined in [START_REF] Eskin | Isolation, equidistribution, and orbit closures for the SL(2, R) action on moduli space[END_REF] to the Hurwitz space for backbone graphs. The idea behind this relation was already mentioned in the last paragraph of Sect. 2.2. If ( , ) is a backbone graph, then we claim that

H P 1 (μ -1 , p) ∼ = M -1
, where the isomorphism is provided by the source map f S that marks the critical points of the branched covers. To verify the claim, let ω be the meromorphic differential on the unique vertex of of level -1 as in part iii) of the definition for M -1 . Since is of compact type, the global residue condition in [START_REF] Bainbridge | Compactification of strata of abelian differentials[END_REF] imposed to ω implies that all residues of ω vanish. Therefore, a point in M -1 can be identified with such a meromorphic differential ω (up to scale) on P 1 without residues, such that ω has zeros of order m i for m i ∈ μ -1 at the corresponding markings and has poles of order p j + 1 for p j ∈ p at the corresponding nodes. In particular, ω is an exact differential and integrating it on P 1 provides a meromorphic function that can be regarded as a branched cover f parameterized in H P 1 (μ -1 , p). Conversely given f in H P 1 (μ -1 , p), we can treat f as a meromorphic function and taking d f gives rise to such ω. We thus conclude that

H P 1 (μ -1 , p) ∼ = M -1 . Consequently for ∈ BB(g, n) 1,2 , we have h P 1 (μ -1 , p) = H P 1 (μ -1 ,p) f * S 3≤i≤n i →v -1 ψ i = M -1 3≤i≤n i →v -1 ψ i , ( 21
)
where i → v -1 means that the ith marking belongs to the vertex of level -1. Now we can proceed with the proof of Proposition 3.10.

Proof We write α , ,p = i≥0 α i , ,p ξ i where α i , ,p is a pull-back from H * (M g,n , Q). By Lemma 3.8 we deduce that

p * (ξ 2g-1 α , , p ) = (-1) g λ g α 0 , ,p . (22)
Therefore we only need to consider the ξ -degree zero part of α , ,p , which is given by

ζ * ⎛ ⎝ [M -1 ] v∈V ( ), (v)=0 [P M g v ,n v ( p v -1, μ v )] 0 ⎞ ⎠ ,
where

[P M g v ,n v ( p v -1, μ v )] 0 is the degree zero part of the Poincaré-dual class of P M g v ,n v ( p v -1, μ v ) in P M g v ,n v . Therefore, we have λ g • α 0 , ,p = ζ * ⎛ ⎝ [M -1 ] v∈V ( ), (v)=0 λ g v • [P M g v ,n v ( p v -1, μ v )] 0 ⎞ ⎠ .
Multiplying this expression by n i=3 ψ i , we obtain that

λ g • α 0 , ,p • n i=3 ψ i = ζ * ⎛ ⎝ ⎛ ⎝ [M -1 ] • i →v -1 ,i≥3 ψ i ⎞ ⎠ v∈V ( ), (v)=0 λ g v • [P M g v ,n v ( p v -1, μ v )] 0 • i →v ψ i ⎞ ⎠ .
123

For the first factor on the right-hand side, equality [START_REF] Eskin | Asymptotics of numbers of branched coverings of a torus and volumes of moduli spaces of holomorphic differentials[END_REF] implies that

[M -1 ] • i →v -1 ,i≥3 ψ i = h P 1 (μ -1 , p) .
Moreover for all v of level 0, we have

p v a 1 ( p v -1, μ v ) = P M gv ,nv ( p v -1,μ v ) ξ 2g v -1 i →v ψ i = (-1) g v M gv ,nv λ g v • i →v ψ i • [P M g v ,n v ( p v -1, μ v )] 0
where the second identity follows from Lemma 3.8. Since p v is the (positive) twist value assigned to the edge of v, the product of p v over all vertices of level 0 equals m(p) defined in [START_REF] Eskin | Moduli spaces of Abelian differentials: the principal boundary, counting problems, and the Siegel-Veech constants[END_REF]. In addition, the sum of g v over all vertices of level 0 equals the total genus g, because is of compact type and v -1 has genus zero. Putting everything together we thus obtain that

M g,n λ g • α 0 , ,p • n i=3 ψ i = m(p) • h P 1 (μ -1 , p) • v∈V ( ), (v)=0 a 1 ( p v -1, μ v ),
which is the desired statement.

The induction formula for cohomology classes

The main tool of the section is the induction formula in [43, Theorem 6 (1)] which we recall now.

Proposition 3.11

For all 1 ≤ i ≤ n, the relation that

(ξ + (m i + 1)ψ i ) [P M g,n (μ)] = ( , , p) i →v, (v)=-1 m(p) |Aut( , , p)| α , ,p (23)
holds in H * (P M g,n , Q), where the sum is over all twisted bi-colored graphs such that the ith leg is carried by a vertex of level -1.

There are two ways of using Eq. ( 23). First one can compute the Poincaré-dual class of P M g,n (μ) in H * (P M g,n , Q) in terms of the ψ, λ, ξ classes and boundary classes associated to stable graphs. This strategy is used in [START_REF] Sauvaget | Volumes and Siegel-Veech constants of H(2g-2) and Hodge integrals[END_REF] to deduce the first formula in Theorem 3.1.

Alternatively, one can compute relations in the Picard group of P M g,n (μ) to deduce relations between intersection numbers on P M g,n (μ) and intersection numbers on boundary strata associated to twisted graphs. This is the strategy that we will use here. We will use this proposition with i ∈ {1, 2} and multiply the formula by ξ 2g-1 n i=3 ψ i to obtain a 1 (μ) on the left-hand side. Then we will use Propositions 3.9 and 3.10 to compute the right-hand side. A first application of this strategy gives a proof of the complementary proposition.

Proof of Proposition 3.2 We use Proposition 3.11 with i = 1. Multiplying formula ( 23) by ξ 2g-2 • n i=2 ψ i , we obtain that

(m 1 + 1) a 1 (μ) + P M g,n (μ) ξ 2g-2 • n i=1 ψ i = ( , ,p) 1 →v, (v)=-1 m(p) |Aut( , , p)| P M g,n α , ,p • ξ 2g-2 • n i=2 ψ i .
It suffices to check that each summand in the right-hand side vanishes. Proposition 3.9 implies that if ( , , p) is not a backbone graph, then the corresponding summand vanishes. If ( , , p) is a backbone graph, then we have seen (in the paragraph below Proposition 3.10) that M -1 is birational to a Hurwitz space of admissible covers of dimension n -1 -2, where n -1 is the number of legs adjacent to the vertex of level -1. Since the ψ-product restricted to level -1 contains n -1 -1 terms (i.e. it misses ψ 1 only), which is bigger than dim M -1 , it implies that the intersection of α , ,p with n i=2 ψ i vanishes. Now we know that a i (μ) is independent of the choice of 1 ≤ i ≤ n and hence we can drop the subscript i. The second use of the strategy presented above leads to the following induction formula.

Lemma 3.12 The intersection numbers a(μ) satisfy the recursion

(m 1 + 1)(m 2 + 1)a(μ) = k≥1 g,μ h P 1 ((m 1 , m 2 , μ 0 ), p) • 1 k! • k i=1 p 2 i a( p i -1, μ i ) where g = (g 1 , . . . , g k ) is a partition of g, μ = (μ 0 , μ 1 , . . . , μ k ) is a (k + 1)- tuple of multisets with (m 3 , . . . , m n ) = μ 0 • • • μ k and p = ( p 1 , . . . , p k ) has entries p i = 2g i -1 -|μ i | > 0.
We remark that this induction formula is not quite the same as the induction formula of Theorem 3.1, e.g. the sums in the two formulas do not run over the same set. Theorem 3.1 will follow further from a combination of Lemma 3.12 and Proposition 2.1 of the previous section.

Proof We apply the induction formula of Proposition 3.11 with i = 2:

(ξ + (m 2 + 1)ψ 2 ) [P M g,n (μ)] = ( , ,p) 2 →v, (v)=-1 m(p) |Aut( , , p)| α , ,p .
We multiply this expression by ξ 2g-1 n i=3 ψ i and apply p * . Since Lemma 3.6 gives p * (ξ 2g [P M g,n (μ)]) = 0, the above equality implies that

(m 1 + 1)(m 2 + 1)a(μ) = ( , ,p) 2 →v, (v)=-1 m(p) |Aut( , , p)| p * ξ 2g-1 • n i=3 ψ i • α , ,p .
By Proposition 3.9 a term in the sum of the right-hand side vanishes if ( , , p) is not a backbone graph. Suppose ( , , p) is a backbone graph such that the first leg does not belong to the vertex of level -1 (which contains n -1 legs).

Then on level -1 the product of ψ-classes contains n -1 -1 terms (i.e. this product misses ψ 2 only), which exceeds the dimension of M -1 (being n -1 -2), hence the corresponding term in the sum also vanishes. Now we only need to consider the case when ( , , p) is a backbone graph in BB(g, n) 1,2 , i.e. the vertex of level -1 carries both the first and second legs. Then the intersection number ξ 2g-1 • n i=3 ψ i • α , ,p is given by Proposition 3.10. We thus conclude that

(m 1 + 1)(m 2 + 1)a(μ) = ( , ,p)∈BB 1,2 h P 1 (μ -1 , p) m(p) 2 |Aut( , , p)| • v∈V ( ), (v)=0 a( p v -1, μ v ) = k≥1 g,μ m(p) 2 h P 1 ((m 1 , m 2 , μ 0 ), p) • 1 k! • k i=1 a( p i -1, μ i ).
The last equality comes from the fact that

g 1 + • • • + g k = g where g i ≥ 1 and (m 3 , . . . , m n ) = μ 0 μ 1 • • • μ k determines uniquely a graph ( , , p) in BB(g, n) 1,2
and an automorphism of the backbone graph is determined by a permutation in S k that preserves both the partition of g and the sets μ 1 , . . . , μ k .

Sums over rooted trees

The purpose of this section is to combine the preceding Lemma 3.12 with Proposition 2.1 that describes the computation of intersection numbers on Hurwitz spaces. We will show that the numbers a(μ) can be expressed as sums over rooted trees in a similar way as we did for intersection numbers on Hurwitz spaces in Sect. 2.4. Let 2 ≤ i ≤ n and ( , , p) be a rooted tree in RT(g, μ) 1,i (here μ[∞] is empty). Since there is no marked pole, it implies that any vertex of genus zero has at least one edge with a negative twist, hence it is an internal vertex of and lies on a negative level. Denote by μ[∞] 0 the list obtained by taking the (positive) entries p(h) for all half-edges h adjacent to a vertex of level 0. Denote by μ[0] 0 the list of entries of μ from those legs carried by the internal vertices (of genus zero). With this notation we define the rooted tree ( 0 , 0 , p 0 ) in RT(0, μ[0] 0 , μ[∞] 0 ) 1,i obtained by removing the leaves of (i.e. vertices of positive genus and hence on level 0). We also define the multiplicity m 0 (p) of ( 0 , 0 , p 0 ) to be the product of entries of μ[∞] 0 . Now we define the a-contribution of the rooted tree ( , , p) as

a( , , p) = m 0 (p) 2 h( 0 , 0 , p 0 ) v∈V ( ), (v)=0 a( p v -1, μ v ), (24)
where h( 0 , 0 , p 0 ) is the contribution of the rooted tree defined in [START_REF] Chen | Principal boundary of moduli spaces of abelian and quadratic differentials[END_REF].

Lemma 3. [START_REF] Costantini | The area is a good metric[END_REF] The following equality holds:

(m 1 + 1)(m 2 + 1)a(μ) = ( , ,p)∈RT(g,μ) 1,2 a( , , p) |Aut( , , p)| .
Proof Removing the leaves of a rooted tree induces a bijection between RT(g, μ) 1,2 and the set

( , ,p)∈BB 1,2 RT(0, μ[0] 0 , μ[∞] 0 ) 1,2
which is a partition of RT(g, μ) 1,2 over all possible decorations of the leaves of the rooted trees (i.e. each decoration is induced by a graph in BB 1,2 ). Moreover, an automorphism of a rooted tree in RT(g, μ) 1,2 is determined by an automorphism of the backbone graph in BB(g, n) 1,2 , because all internal vertices of the rooted tree (i.e. those of genus zero and hence on negative levels) have marked legs by Definition 2.4. Then we can first use Lemma 3.12 to write a(μ) as a sum over backbone graphs in BB(g, n) 1,2 and then use Proposition 2.1 to express it as the desired sum over the set

( , ,p)∈BB 1,2 RT(0, μ[0] 0 , μ[∞] 0 ) 1,2 RT(g, μ) 1,2
as claimed in the lemma.

We define

RT(g, μ) 1 = {trivial graph} ∪ n i=2 RT(g, μ) 1,i
and the a-contribution of the trivial graph

• as a(•, , p) = (m 1 + 1) 2 a(μ).
End of the proof of Theorem 3.1 We will prove for n ≥ 2 the equality that

k≥1 g,μ h P 1 ((m 1 , m 2 ), p) 1 k! k j=1 (2g j -1 + n(μ j )) p j a( p j -1, μ j ) = ( , ,p)∈RT(g,μ) 1,2 a( , , p) |Aut( , , p)| . ( 25
)
This formula together with Lemma 3.13 thus implies Theorem 3.1. Since by definition n i=1 (m i + 1) = 2g -2 + n, Lemma 3.13 implies that ,,p)| .

(2g -2 + n)(m 1 + 1)a(μ) = (m 1 + 1) 2 a(μ) + n i=2 (m i + 1)(m 1 + 1)a(μ) = (m 1 + 1) 2 a(μ) + n i=2 ( , ,p)∈ RT(g,μ) 1,i a( , , p) |Aut( , , p)| = ( , ,p)∈ RT(g,μ) 1 a( , , p) |Aut(
Therefore, the left-hand side of (25) can be rewritten as

k≥1 g,μ h P 1 ((m 1 , m 2 ), p) 1 k! • k j=1 ( , ,p)∈ RT(g j ,( p j -1,μ j )) 1 a( , , p) |Aut( , , p)| = k≥1 g,μ h P 1 ((m 1 , m 2 ), p) 1 k! • ( j , j ,p j )∈ RT(g j ,( p j -1,μ j )) 1 k j=1 a( j , j , p j ) |Aut( j , j , p j )| .
We claim that there is a bijection

RT(g, μ) 1,2 ( , ,p )∈BB 1,2 v∈V ( ), (v)=0 RT(g v , (p v -1, μ v )) 1 .
Indeed given a rooted tree ( , , p) in RT(g, μ) 1,2 we can construct ( , , p ) ∈ BB(g, n) 1,2 by contracting all edges except those adjacent to the root, and the rooted trees

( v , v , p v ) ∈ RT(g v , (p v -1, μ v ))
1 are the connected components of the graph obtained from ( , , p) by deleting the root. Moreover for a rooted tree ( , , p), by Eq. ( 10) we have

h( 0 , 0 , p 0 ) = h P 1 ((m 1 , m 2 ), p ) • k j=1 h( j0 , j0 , p j0 ),
where as before 0 is obtained from by removing the leaves and the j0 are the connected components after removing the root of 0 . Together with the definition of the a-contribution in [START_REF] Farkas | The moduli space of twisted canonical divisors[END_REF], it implies that a( , , p)

= h P 1 ((m 1 , m 2 ), p ) • v∈V ( ), (v)=0 a( v , v , p v ) .
Note also that

Aut( , , p) = Aut( , , p ) × v∈V ( ), (v)=0 Aut( v , v , p v ) .
Combining the above we thus conclude that equality (25) holds.

Proof of the equivalence of Theorems 1.1 and 1.2 We first assume that Theorem 1.2 holds. By Theorem 3.1, the quantities

vol( M g,n (m 1 , . . . , m n )) and 2(2πi) 2g (2g -3 + n)! a(m 1 , . . . , m n )
satisfy the same induction relation that determines both collections of these numbers starting from the case n = 1. The base case (i.e. the minimal strata)

that vol( M g,1 (2g -2)) = 2(2πi) 2g (2g -2)! a(2g -2) (26)
was proved in [START_REF] Sauvaget | Volumes and Siegel-Veech constants of H(2g-2) and Hodge integrals[END_REF] under a mild assumption of regularity of a natural Hermitian metric on O(-1), and we will give an alternative (unconditional) proof in Sect. 4.4. Consequently we conclude that Theorem 1.2 implies Theorem 1.1.

The converse implication follows similarly.

Volume recursion via q-brackets

In this section we define recursively polynomials in the ring

R = Q[h 1 , h 2 , . . .]
and show that they compute volumes of the strata after a suitable specialization. The method of proof relies on lifting the E 2 -derivative via the Bloch-Okounkov q-bracket and expressing cumulants in terms of this lift. This recursion looks quite different from the recursion given in Theorem 1.2, since it is only defined on the level of polynomials in the variables h i and requires h i -derivatives.

To define the substitution, first recall the numbers b j introduced in (15). We let [START_REF] Harris | On the Kodaira dimension of the moduli space of curves[END_REF] where the denominator denotes the inverse function of u/P B (u). For the recursion we define for a finite set I = {i 1 , . . . , i n } of positive integers the formal series

P B (u) = exp - j≥1 j! b j+1 u j+1 and α = u 1 (u/P B (u)) -1 ,
H I ∈ R[[z i 1 , . . . , z i n ]] if |I | ≥ 2 and H {i} ∈ 1 z i R[[z i ]] by H {i} = 1 z i + ≥1 h z i , H {i, j} = z i H (z i ) -z j H (z j ) H(z j ) -H(z i ) -1 = 2h 1 z i z j + h 2 (3z 2 i z j + 3z i z 2 j ) + 4h 3 z 3 i z j + (2h 2 1 + 4h 3 )z 2 i z 2 j + 4h 3 z i z 3 j + • • • , H I = 1 2(n -1) I =I I I ,I =∅ D 2 (H I , H I ), ( 28
)
where we abbreviate

H n = H [[1,n]] , H = H 1 and h 1 ,..., n = [z 1 1 . . . z n n ]
H n and where the symmetric bi-differential operator D 2 is defined by

D 2 ( f, g) = 1 , 2 ≥1 h 1 , 2 ∂ f ∂h 1 ∂g ∂h 2 . ( 29
)
Theorem 4.1 The rescaled volume of the stratum with signature μ = (m 1 , . . . , m n ) can be computed as

v(μ) = (2πi) 2g (2g -2 + n)! h m 1 +1,...,m n +1 h →α
using the recursion [START_REF] Harris | Moduli of Curves[END_REF] and the values of the α in [START_REF] Harris | On the Kodaira dimension of the moduli space of curves[END_REF].

Three sets of generators for the algebra of shifted symmetric functions

We let * be the algebra of shifted symmetric functions (see e.g. [START_REF] Eskin | Asymptotics of numbers of branched coverings of a torus and volumes of moduli spaces of holomorphic differentials[END_REF][START_REF] Zagier | Partitions, quasimodular forms, and the Bloch-Okounkov theorem[END_REF] or [START_REF] Chen | Quasimodularity and large genus limits of Siegel-Veech constants[END_REF]) and recall the standard generators

p (λ) = ∞ i=1 (λ i -i + 1 2 ) -(-i + 1 2 ) + (1 -2 -) ζ(-) . ( 30
) Note that (1-2 -) ζ(-) = !b +1 .
The algebra * is provided with a grading where each p has weight + 1. For Hurwitz numbers the geometrically interesting generators are

f (λ) = z χ λ ( )/ dim χ λ , ( 31
)
where z is the size of the conjugacy class of the cycle of length , completed by singletons. The first few of these functions are

f 1 = p 1 + 1 24 , f 2 = 1 2 p 2 , f 3 = 1 3 p 3 - 1 2 p 2 1 + 3 8 p 1 + 9 640
.

The third set of generators, defined implicitly by Eskin and Okounkov, will serve as top term approximations of f . We define h ∈ * by

h = -1 [u +1 ]P(u) where P(u) = exp - s≥1 u s+1 p s . ( 32
)
123

Observe that by definition h has pure weight + 1. The first few of these functions are

h 1 = p 1 , h 2 = p 2 , h 3 = p 3 - 3 2 p 2 1 .
Proposition 4.2 [21, Theorem 5.5] The difference fh / has weight strictly less than + 1.

We abuse the notation h for generators of R and for elements in * . This is intentional and should not lead to confusion, since the map h → h induces an isomorphism of algebras R ∼ = * , by the preceding proposition.

The lift of the evaluation map to the Bloch-Okounkov ring

Let f : P → Q be an arbitrary function on the set P of all partitions. Bloch and Okounkov [START_REF] Bloch | The character of the infinite wedge representation[END_REF] associated to f the formal power series

f q = λ∈P f (λ) q |λ| λ∈P q |λ| ∈ Q[[q]], (33)
which we call the q-bracket, and proved that this q-bracket is a quasimodular form of weight k whenever f belongs to the subspace of * of weight k (see [START_REF] Bloch | The character of the infinite wedge representation[END_REF], and [START_REF] Zagier | Partitions, quasimodular forms, and the Bloch-Okounkov theorem[END_REF] or [START_REF] Goujard | Counting Feynman-like graphs: quasimodularity and Siegel-Veech weight[END_REF] for alternative proofs).

In [START_REF] Chen | Quasimodularity and large genus limits of Siegel-Veech constants[END_REF]Section 8] we studied in detail an evaluation map Ev (implicitly defined in [START_REF] Eskin | Asymptotics of numbers of branched coverings of a torus and volumes of moduli spaces of holomorphic differentials[END_REF]) on the ring of quasimodular forms that measures the growth rate of the coefficients of quasimodular forms, or equivalently, their asymptotics as τ → 0 along the imaginary axis [START_REF] Chen | Quasimodularity and large genus limits of Siegel-Veech constants[END_REF]Proposition 9.3]. The purpose of this section is to lift this evaluation map to the Bloch-Okounkov ring and to express it in terms of the generators h i introduced in the previous section.

The map Ev is the algebra homomorphism from the ring of quasimodular forms

M * = Q[E 2 , E 4 , E 6 ] to Q[X ]
, sending the Eisenstein series E 2 (normalized to have constant coefficient one) to X + 12, E 4 to X 2 , and E 6 to X 3 . In this way, the larger the degree of Ev( f ), the larger the (polynomial) growth of the coefficients of f , see [START_REF] Chen | Quasimodularity and large genus limits of Siegel-Veech constants[END_REF]Proposition 9.4] for the precise statement. It is also convenient to work with the evaluation map

3 ev[F]( ) = 1 k/2 Ev[F] X → 1 ∈ Q[1/ ] for F ∈ M k . ( 34
)
We also use the brackets f X := Ev f q (X ) and f := ev f q ( ) for f ∈ * as abbreviation. Note that * admits a natural ring homomorphism to Q, the evaluation at the emptyset, explicitly given by the map p → !b +1 .

Proposition 4.3

There is a second order differential operator : * → * of degree -2 and a derivation ∂ : * → * of degree -1 such that

f = 1 k (e ( -∂ 2 -∂/∂p 1 )/2 f ) (∅) ( 35
)
for f ∈ * k homogeneous of weight k. The differential operators are given in terms of the generators p by

∂( f ) = i≥2 i p i-1 ∂ ∂ p i and ( f ) = k, ≥ 1 (k + ) p k+ -1 ∂ 2 ∂ p k ∂ p . ( 36
)
Proof From the definition and [12, Proposition 9.2] we deduce that the evaluation map can be computed for any

F ∈ M k as ev[F]( ) = 1 k a 0 e d F ( 37
)
where d = 12∂/∂ E 2 and where a 0 : F → F(∞) is the constant term map from M * to Q. From [12, Proposition 8.3] we deduce (note that differentiation with respect to Q i in loc. cit. gives the extra p 1 -derivative here) that the differential operators defined above have the property that

d f q = 1 2 -∂ 2 -∂/∂ p 1 f q ( f ∈ * ) . ( 38
)
Since the constant term of the q-bracket of f is in * , the claim follows from these two equations.

To motivate the next section, we recall the notion of cumulants. Let R and R be two commutative Q-algebras with unit and : R → R a linear map sending 1 to 1. (Of course the cases of interest to us will be when R is the Bloch-Okounkov ring * and is the q-, X -, or -bracket to

R = M * , Q[X ], or Q[π 2 ][ ],
respectively.) Then we extend to a multi-linear map R ⊗n → R for every n ≥ 1, called connected brackets, the image of

g 1 ⊗ • • • ⊗ g n being denoted by g 1 | • • • |g n , that we define by g 1 | • • • |g n = α∈P(n) (-1) (α)-1 ( (α) -1)! A∈α a∈A g a . ( 39
)
123

The most important property of connected brackets is their appearance in the logarithm of the original bracket applied to an exponential:

log e g 1 +g 2 +g 3 +••• = log 1 + i g i + 1 2! i, j g i g j + 1 3! i, j,k g i g j g k + • • • = i g i + 1 2! i, j g i |g j + 1 3! i, j,k g i |g j |g k + • • • .
We specialize to the Bloch-Okounkov ring * and we want to compute the leading terms of the connected brackets associated with the • X -or •brackets. Recall from [12, Proposition 11.1]:

Proposition 4.4 Let g i ∈ * ≤k i (i = 1, . .

. , n) be elements of weight less than or equal to k i and let g i ∈ *

k i be their top weight components. Let

k = k 1 + • • • + k n be the total weight. Then deg( g 1 | • • • |g n X ) ≤ 1 -n + k/2 and [X 1-n+k/2 ] g 1 | • • • |g n X = [X 1-n+k/2 ] g 1 | • • • |g n X . ( 40
)
The leading terms of the brackets are consequently [START_REF] Okounkov | Gromov-Witten theory, Hurwitz numbers, and matrix models[END_REF] We call them rational cumulants.

g 1 | • • • |g n L = [X 1-n+k/2 ] g 1 | • • • |g n X = lim X →∞ Ev[ g 1 | • • • |g n q ](X ) X 1-n+k/2 = [ -k-1+n ] g 1 | • • • |g n = lim →0 k+1-n ev[ g 1 | • • • |g n q ]( ) .

The cumulant recursion

In this section we prove a formula for computing the connected brackets associated with the • q -or rather the • -brackets. The core mechanism for their computation is summarized in the following purely algebraic property.

Let R be an N-graded commutative Q-algebra with R 0 = Q, complete with respect to the maximal ideal m = R >0 . The following statement gives a general recursion for expressions that appear in cumulants. We will specialize R to the Bloch-Okounkov ring subsequently. Key Lemma 4.5 Suppose that D : R → R is a linear map. Then the following statements are equivalent:

(1) We have D(x 3 ) -3x D(x 2 ) + 3x 2 D(x) = 0 for all x ∈ R.

123

(2) For all x ∈ R and all n ≥ 2

D(x n ) = n 2 D(x 2 )x n-2 -n(n -2)D(x)x n-1 .
( (4) If we denote by D 2 : R 2 → R the symmetric bilinear form

D 2 (x, y) = D(x y) -x D(y) -y D(x) then D(x 1 • • • x n ) = n i=1 D(x i )x 1 • • • x i • • • x n + 1≤i< j≤n D 2 (x i , x j )x 1 • • • x i • • • x j • • • x n for all x 1 , . . . , x n ∈ R. ( 5
) For any fixed x ∈ R, the bilinear form D 2 (x, y) is a derivation in y. [START_REF] Bainbridge | Compactification of strata of abelian differentials[END_REF] The map D ∈ Sym 2 (Der(R)), i.e. D is a second order differential operator without constant term. [START_REF] Bainbridge | The moduli space of multi-scale differentials[END_REF] For all X ∈ m there exists L(X ) ∈ R such that

log(e D (e X/ )) = 1 L(X ) + O(1) ( → 0) . ( 42
)
If any of these statements holds, the leading term of (42) is given by L(X ) = L [START_REF] Arbarello | Geometry of Algebraic Curves[END_REF], where

L(0) = X, L (t) = 1 2 D 2 (L(t), L(t)) . ( 43
)
Proof The implication (1) ⇒ (3) follows by passing to the polarization. The implication (3) ⇒ (4) can be proved by induction (replace x 1 by x 0 x 1 ). The implications (2) ⇒ (1) and (4) ⇒ (1), ( 2), (3) follow by specialization. The equivalence (5) ⇔ (3) follows by direct computation. To show (6) ⇔ (5) think deeply. To prove (7) ⇒ (1) it suffices to consider the cubic term: the coefficient of 1/ 2 is the expression in (1).

To prove (6) ⇒ (7) and the final formula for L(X ) we write e D (e x ) = e y( ) . [START_REF] Veech | Gauss measures for transformations on the space of interval exchange maps[END_REF] 123

Then y(0) = x. Note that (2) implies that

D(e x ) = 1 2 D(x 2 )e x + D(x)(1 -x)e x .
Differentiation of ( 44) with respect to implies that y = D(y) + 1 2 D 2 (y, y). Equivalently, writing y = n≥0 f n (x) n , then the initial condition is that f 0 (x) = x, and

(n + 1) f n+1 (x) = D( f n ) + 1 2 n m=0 D 2 ( f m (x), f n-m (x)) .
Recursively this implies that

f n (X/ ) = L n (X ) -n-1 + O( -n ) with L 0 (X ) = X and (n + 1)L n+1 (X ) = 1 2 n m=0 D 2 (L m (X ), L n-m (X )) .
We now let L(t) = n≥0 L n (X )t n and the claims follow.

Application to volume computations

We now return to the proof of the main theorem of this section. Recall the main idea from [START_REF] Eskin | Asymptotics of numbers of branched coverings of a torus and volumes of moduli spaces of holomorphic differentials[END_REF] that the volume of a stratum is given by the growth rate of the number of connected torus covers and thus to the leading terms of cumulants of the f . More precisely for 2g -2 = n i=1 m i , the same argument as in [START_REF] Chen | Quasimodularity and large genus limits of Siegel-Veech constants[END_REF]Proposition 19.1] gives that

vol ( M g,n (m 1 , . . . , m n )) = (2πi) 2g f m 1 +1 | • • • | f m n +1 L (2g -2 + n)! . ( 45
)
Proof of Theorem 4.1, one variable case First, P B (u

) = P(u)| p → !b = P(u)(∅).
Next, recall that Lagrange inversion for a power series

F ∈ uC[[u]]
with non-zero linear term and inverse

G(z) states that k[z k ]G n = n[u -n ]F -k
for k, n = 0. We apply this to F = u/P B (u) and to k = 2g -1 and n = -1 to obtain that

(2g -1)! (2πi) 2g v(2g -2) = (2g -1) f 2g-1 L = h 2g-1 L = -1 2g -1 [u](P B (u)/u) 2g-1 = u 2g-1 1 (u/P B ) -1 (46)
using Proposition 4.2, (32) and Lagrange inversion.

123

We pause for a moment to check the initial condition of the theorem in the previous section independently of the Hermitian metric extension problem along the boundary of the strata.

Proof of ( 26) using [START_REF] Vorobets | Ergodicity of billiards in polygons[END_REF] We want to show that (2g-1) 2 a(2g-2) = h 2g-1 L . Recall that a version of Lagrange inversion (see e.g. [25, Formula (2.2.8)]), in fact the case k = 0 excluded in the version of the previous proof, states that if

F ∈ z + z 2 C[[z]
] with composition inverse G(u), then for any Laurent series φ(z)

[z 0 ]φ(F) = [u 0 ]φ(u) + [u -1 ]φ (u) log(G/u) . ( 47
)
If we let A(z) = 1/z + g≥1 h 2g-1 L z 2g-1 , then we need to show that A(z) = A(z). We apply Lagrange inversion to φ(z) = z -2g and F = 1/ A(z) to obtain that

1 2g! [z 0 ] A 2g = 1 2g! [z 0 ]φ(1/ A(z)) = -1 (2g -1)! [u 2g ] log(1/P B ) = -1 (2g -1)! [u 2g ] s≥1 s!b s+1 u s+1 = [u 2g ]B(u) = 1 (2g)! [z 0 ]A 2g
using [START_REF] Vorobets | Ergodicity of billiards in polygons[END_REF] and ( 16). This implies the claim.

For the general case of the theorem, we apply Sect. 4.3 to the differential operator

D = 1 2 ( -∂ 2 -∂/∂ p 1 ) . ( 48
)
Proposition 4. [START_REF] Bainbridge | Compactification of strata of abelian differentials[END_REF] The bilinear differential operator D 2 defined in [START_REF] Ivanov | Kerov's central limit theorem for the Plancherel measure on Young diagrams[END_REF]

is the polarization of D, namely, D 2 ( f, g) = D( f g) -f D(g) -g D( f
) for all f and g.

Proof In terms of the p -generators the polarization is given by

D 2 ( f, g) = k, ≥1 (k + ) p k+ -1 -k p k-1 p -1 ∂ f ∂ p k ∂g ∂ p . ( 49
)
The definition [START_REF] Kontsevich | Lyapunov exponents and Hodge theory[END_REF] of h in terms of p implies that 123

∂H(z) ∂ p k = - zH (z) H(z) k+1 and n≥2 np n-1 H -n (z) = -H(z)/zH (z) -1. (50)
We compute that

k, ≥1 (k + ) p k+ -1 -k p k-1 p -1 z 1 H (z 1 ) H(z 1 ) k+1 • z 2 H (z 2 ) H(z 2 ) +1 = z 1 H (z 1 )z 2 H (z 2 ) H(z 1 )H(z 2 ) n≥2 np n-1 H(z 1 ) 1-n -H(z 2 ) 1-n H(z 2 ) -H(z 1 ) -1 + H(z 1 ) z 1 H (z 1 ) 1 + H(z 2 ) z 2 H (z 2 ) = H 2 (z 1 , z 2 ),
and this implies the claim by the chain rule.

We now define the partition function of h-brackets

H (u) q = exp ≥1 h u q = n≥0 h 1 • • • h 1 n 1 h 2 • • • h 2 n 2 • • • q u n n! = ∞ n=0 1 n! 1 ,..., n ≥1 h 1 • • • h n q u 1 • • • u n (51)
in the h -variables. Then the partition function of the rational cumulants for the h -generators

H (u) q = n≥0 h 1 | • • • |h 1 n 1 | h 2 | • • • |h 2 n 2 | • • • q u n n! = log H (u) q (52)
is simply the logarithm of H .

Proof of Theorem 4.1, general case

We first show that the pieces of H sorted by total degree in u can be recursively computed using the D 2 -operator. For this purpose we let h i = -i h i . From the definition of cumulants, Eq. ( 37) and a 0 ( g q ) = g(∅), we obtain that

n≥0 h 1 | • • • | h 1 n 1 | h 2 | • • • | h 2 n 2 | • • • u n |n| n! = log exp( 1 i≥1 h i u i ) = log e D exp 1 i≥1 h i u i (∅) . ( 53
)
By applying the Key Lemma with X = i≥1 h i u i and undoing the rescaling of the h i using ( 41) we obtain that

n≥0 h 1 | • • • |h 1 n 1 | h 2 | • • • |h 2 n 2 | • • • L u n n! = ∞ n=0 L n (∅) ( 54
)
with

L 0 = n≥1 h i u i and L n = 1 2n r +s=n-1 D 2 (L r , L s ) ( 55
)
for n > 0. Now define a linear map n :

Q[u] → Q[z] by n (u 1 • • • u n ) = Symm(z 1 1 • • • z n n )
and zero for monomials of length different from n, where Symm denotes symmetrization with respect to the S n action on the variables z i . In this notation H 1 = n L n-1 and

H n = n n≥0 h 1 | • • • |h 1 n 1 | h 2 | • • • |h 2 n 2 | • • • L u n n! .
Consequently, (55) and ( 28) together with [START_REF] Veech | Siegel measures[END_REF] and Proposition 4.2 imply the claim.

Equivalence of volume recursions

In this section we introduce another "averaged volume" recursion that interpolates between the D 2 -recursion introduced in Sect. 4 and the volume recursion in Theorem 1.2. We will show that the averaged volume recursion and the D 2 -recursion give the same generating functions, and then Theorem 1.1 will follow from it.

Recall from [START_REF] Harris | Moduli of Curves[END_REF] the definition of

H {i, j} ∈ R[[z i , z j ]] for i = j, where R = Q[h 1 , h 2 , . . .].
For any list of positive integers p = ( p 1 , . . . , p k ) we define 123

H p {i, j} = ∂ k ∂h p 1 • • • ∂h p k H {i, j} .
For a finite set I = {i 1 , . . . , i n } of positive integers we define the formal series

A I ∈ R[[z i 1 , z i 2 ,

. . .]] inductively by

A I = H I ∈ 1 z i R[[z i ]]
if n = 1, and otherwise

A I = 1 n -1 1≤r <s≤n k≥0 p=( p 1 ,..., p k ) I ={i r ,i s } I 1 ••• I k 1 k! H p i r ,i s • k j=1 A [ p j ] I j , ( 56
)
where

A [ p] I := [z p i ]A I ∪{i} for any i / ∈ I . We set A n = A [[1,n]] . Note that A I = H I by definition if |I | = 2.
We have chosen to sum in definition of A I over all p rather than partitions g = k i=1 g i as e.g. in Theorem 1.2. The two summations are equivalent, since g i and p i determine each other once I has been partitioned. Our goal here is to show the following result.

Theorem 5.1 For all non-empty sets of positive integers I , we have

A I = H I .
For the proof of this theorem we will show that both A I and H I can be written as a sum which ranges over certain oriented trees (see Sect. 5.2). The two recursions can then be viewed as stemming from cutting the trees at a local maximum (a "top") or a local minimum (a "bottom") respectively.

Proof of Theorems 1.2 and 1.1

We assume in this section that Theorem 5.1 holds and finish the proof of Theorems 1.2 and 1.1 under this assumption. We abbreviate λ = ( 1 , . . . , n ) and recall from Sect. 4 that we denoted the coefficients of A n for n ≥ 2 by

A n = H n = 1 ,..., n ≥1 h λ z 1 i 1 • • • z n i n .

Proposition 5.2 The coefficients h λ are uniquely determined by the recursion

h λ = 1≤r <s≤n r + s n -1 k≥1 1 k! g,μ h P 1 ( r -1, s -1, p) • k i=1 h λ i , p i (57)
for n ≥ 2, where the summation is as in Theorem 1.2, except that

λ 1 • • • λ k = [[1, n]] {r, s}.
Proof We begin by showing the formula for λ with two entries, which in view of ( 6) is equivalent to show that

1 , 2 ≥1 h 1 , 2 1 + 2 z 1 1 z 2 2 = k≥1 1 k ⎛ ⎝ p≥1 h p z 1 z 2 z p 1 -z p 2 z 1 -z 2 ⎞ ⎠ k = -log z 1 z 2 z 1 -z 2 (H 1 (z 2 ) -H 1 (z 1 )) .
Since applying (z 1

∂ ∂z 1 + z 2 ∂ ∂z 2 )
to the left-hand side above gives A 2 , and since neither side has a constant term, this in turn follows from

H 2 = A 2 = -1 + z 1 H 1 (z 1 ) -z 2 H 1 (z 2 ) H 1 (z 2 ) -H 1 (z 1 ) = -z 1 ∂ ∂z 1 + z 2 ∂ ∂z 2 log z 1 z 2 z 1 -z 2 (H 1 (z 2 ) -H 1 (z 1 )) .
For λ with more entries, we deduce (for all r , s ≥ 1 and all p) from the preceding calculation that

[z r r z s s ]H {r,s} = k >0 r + s k ! p h P 1 ( r -1, s -1, p ) k i=1 h p i and [z r r z s s ]H p {r,s} = k >0 r + s k ! p h P 1 ( r -1, s -1, p ∪ p ) k i=1 h p i ,
Besides, the recursion formula (56) defining A I can be translated for λ with n parts into

h λ = 1≤r <s≤n 1 n -1 k≥1,g i ,λ i |λ i |>0 1 k! [z r r z s s ]H p {r,s} k i=1 h λ, p i = 1≤r <s≤n r + s n -1 k≥1,g i ,|λ i |>0 k >0,p 1 k!k ! h P 1 ( r -1, s -1, p ∪ p ) k i=1 h λ, p i k j=1 h p j , ( 58
)
where

λ 1 • • • λ k partitions λ\{ r , s }.
Next we remark that the interior sum of ( 57) is over all backbone graphs with the two markings labelled with r and s at the lower level. In the preceding formula (58) the contribution of vertices 123 with at least one marking is separated from the vertices with no markings. This choice results in a binomial coefficient k+k k and transforms 1 k!k ! into 1 (k+k )! , thus showing that the two recursive formulas (58) and (57) are equivalent.

Proof of Theorems 1.2 and 1.1 For μ = (m 1 , . . . , m n ) consider the intersection numbers a(μ) that satisfy the recursion in Theorem 3.1, and recall that a(μ) = a i (μ) is independent of the index i by Proposition 3.2. In particular the a(μ) satisfy the recursion [START_REF] Eskin | Invariant and stationary measures for the SL(2, R) action on moduli space[END_REF] for any distinguished pair of indices, and hence satisfy every weighted average of these recursions. We use the weighted average where the recursion with (i, j) distinguished is taken with weight k / ∈{i, j} (m k + 1). Conversely, the a(μ) are uniquely determined by this weighted average and the initial values for μ of length one given in [START_REF] Eskin | Asymptotic formulas on flat surfaces[END_REF].

On the other hand, the collection of (2g -2 + n) n i=1 (m i + 1)a(μ) and the collection of h m 1 +1,...,m n +1 both satisfy the recursion (57), by observing that

(2g -2 + n)a(μ) n i=1 (m i + 1) = 1≤r <s≤n (m r + 1 + m s + 1) n -1 a(μ) n i=1 (m i + 1). (59)
Note that A 1 | h →α = A by Theorem 4.1 and since we already checked (see [START_REF] Vorobets | Ergodicity of billiards in polygons[END_REF] and the subsequent proof) that the one-variable rescaled volumes v(2g -2) and a(2g -2) agree (see [START_REF] Goujard | Counting Feynman-like graphs: quasimodularity and Siegel-Veech weight[END_REF]) up to the factor (2πi) 2g /(2g -1)!. This implies that

a(μ) = h m 1 +1,...,m n +1 | h →α (2g -2 + n) n i=1 (m i + 1) . ( 60
)
The claim now follows from Theorems 4.1, 5.1 and the conversion (13) of volumes to the a(μ).

Oriented trees

We now start preparing for the proof of Theorem 5.1. An oriented tree is the datum of a graph G = (V, E ⊂ V × V ) whose underlying graph of (V, E) is a tree. In particular it is required to be connected. If (v, v ) ∈ E, we will denote v > v . Moreover, a vertex v ∈ V is called a bottom (respectively a top) if there exists no v ∈ V such that v > v (respectively v < v ). We will denote by B(G) and T (G) the sets of bottoms and tops of G.

For any oriented tree G with n vertices, we define the rational number

f # (G) := Card σ : V ∼ → [[1, n]], s.t.∀(v, v ) ∈ E, σ (v) > σ (v ) n! , ( 61
)
whose numerator is the number of total orderings on the set of vertices compatible with the orientation of G.

Lemma 5. [START_REF] Aggarwal | Large genus asymptotics for volumes of strata of abelian differentials. With an appendix by A. Zorich[END_REF] The function f # can be expressed as

f # (G) = 1 n • v∈B(G) G f # (G ) = 1 n • v∈T (G) G f # (G ) ,
where in both cases the product is over all connected components G of the oriented graph obtained by deleting the vertex v.

Proof In order to define a total ordering on V compatible with the orientation of G, we begin by choosing a minimal element v ∈ V . This element is necessarily a bottom. Let us fix such a choice and denote by (G 1 , . . . , G k ) the connected components of G\{v}. Let n i be the number of vertices of G i for 1 ≤ i ≤ k. A total ordering on V with minimal element v is equivalent to choosing a total ordering on the vertices of G i for all 1 ≤ i ≤ k and a partition of [[1, n -1]] into k sets of size (n 1 , . . . , n k ). Such an ordering on V is compatible with the orientation of G if and only if each ordering on the vertices of G i is compatible with the orientation of G i for all 1 ≤ i ≤ k. This implies that the number of total orderings on V with minimal element v is equal to

n -1 n 1 • • • n k • k i=1 (n i ! f # (G i )) = (n -1)! • k i=1 f # (G i ) .
Summing over all possible choices of a minimal element, the number of total orderings on vertices of G compatible with the orientation of G is equal to

(n -1)! • v∈B(G) G f # (G ) ,
which completes the proof.

123

Decorations of oriented trees

Let I be a non-empty finite set of positive integers. An I -decoration of an oriented tree = (V, E) is the datum of a function dec : I → V such that for each vertex v the number of outgoing edges plus the number of decorations is equal to two, i.e.

#(dec -1 (v)) + #(E ∩ ({v} × V )) = 2
for all v ∈ V . If I has cardinality greater than one, we denote by OT(I ) the set of I -decorated oriented trees. One can easily check that the following two properties hold:

• if I has cardinality n ≥ 2, then has n -1 vertices;

• a vertex of a decorated tree is a bottom if and only if it has exactly two markings.

We denote by OT(I ) v , OT(I ) b and OT(I ) t the sets of decorated trees with a choice of an arbitrary vertex, a choice of a bottom and a choice of a top, respectively. If I = {i} has only one element, we define OT({i}) v = {i} as a trivial graph decorated by i.

Lemma 5.4 If I has cardinality greater than one, then there is a bijection

ϕ t : OT(I ) t → ⎛ ⎝ I ⊂I OT(I ) v × OT(I \I ) v ⎞ ⎠ (I ∼ I \I ) (62)
given by cutting at a top vertex, were the union is over all non-empty proper subsets. Similarly, there is a bijection

ϕ b : OT(I ) b → {i 1 ,i 2 }⊂I,k>0 ⎛ ⎝ I ={i 1 ,i 2 } I 1 ••• I k k j=1 OT(I j ∪ {e j }) ⎞ ⎠ S k (63)
given by cutting at a bottom vertex, were the union is over all partitions of I into k + 1 non-empty sets such that the first distinguished set has precisely two elements and where the element e j = max(I ) + j for all 1 ≤ j ≤ k. We denote by ψ t and ψ b the inverses of ϕ t and ϕ b , respectively.

Proof Given a decorated tree with a chosen top vertex v, we define its image under ϕ t as follows:

• If there are two markings on v, then v has no outgoing edges, hence the graph has v as a unique vertex and I has only two elements. It follows that OT(I ) v has only one element (and so does the right-hand side of (62)). • If there is only one marking i ∈ I on v and one outgoing edge to a vertex v < v, then I = {i} and the corresponding element in OT(I \I ) v is the graph obtained by deleting v and choosing v as the distinguished vertex. • If there are two outgoing edges to vertices v < v and v < v, then v has no I -markings and the graph obtained by deleting v has two connected components. We define I to be the set of markings on the component containing v and define the corresponding elements of OT(I ) v and OT(I \I ) v to be the connected components containing v and v as chosen vertices, respectively.

The inverse of ϕ t in the first two cases is clear, and in the last case is given by adding a top vertex adjacent to the two chosen vertices. Given a decorated tree with a chosen bottom vertex v, in the same spirit we define the function ϕ b as follows. Since v is a bottom, it has no outgoing edges, hence it has exactly two I -markings i 1 and i 2 , and the corresponding k graphs on the right-hand side of (63) are the k connected components of the graph obtained by removing v. The inverse of ϕ b is given by gluing these k graphs back to v along the vertices marked by e 1 , . . . , e k . Now we fix a ring R and a function g : OT(I ) → R . By slight abuse of notation, we write f # : OT(I ) → Q for the composition of f # defined in (61) with the forgetful map of the decorations. As a consequence of the two preceding lemmas we see that the sum

S(g) = ∈OT(I ) f # ( )g( )
can be rewritten in two different ways, namely

S(g) = 1 2(n -1) I ⊂I ( ,v )∈OT(I ) v ( ,v )∈OT(I \I ) v f # ( ) f # ( ) g ψ t ( , v , , v ) (64)
and

S(g) = 1 n -1 {i 1 ,i 2 }⊂I,k>0 I ={i 1 ,i 2 } I 1 ••• I k 1 k! ( j ) k j=1 k j=1 f # ( j ) • g ψ b k j=1 j , ( 65
)
123
where n is the cardinality of I (i.e. ∈ OT(I ) has n -1 vertices as remarked before).

Explicit expansions over decorated trees

In order to finish the proof of Theorem 5.1, we will show that both A I and H I are equal to a generating series S I that is directly defined as a sum over OT(I ).

Let = (V, E, I → V ) be an oriented tree with decoration by I . A twist assignment on is a function p : E → Z >0 . We work over the ring R[[(z i ) i∈I , (z e ) e∈E ]] of formal series in variables indexed by I ∪ E. Given a twisted decorated oriented tree, we define the contribution of a vertex as

H v = H p v 2 (z v,1 , z v,2 ) ∈ R[[(z i ) i∈I , (z e ) e∈E ]],
where p v is the list of twists associated to all vertices v with v > v and where (z v,1 , z v,2 ) are the variables attached to either the markings of v or the outgoing edges from v to vertices v with v < v. Then we define the contribution of the oriented tree as

S I = ∈OT(I ) f # ( ) cont( ) .
End of the proof of Theorem 5.1 We will show that S I = H I and S I = A I for all sets of positive integers I with n = Card(I ) > 2. The equalities in the case n = 2 are obvious from the definition. We assume now that n ≥ 3 and that S I = H I = A I for all I such that Card(I ) < n.

We first prove that S I = H I . We begin by rewriting the defining Eq. ( 28) with two auxiliary "edge" variables z e and z e for distinct indices e , e ∈ N\I as

H I = 1 2(n -
= (E , V , I → V ), v ) ∈ OT(I ) v and = (E , V , I \I → V ), v ) ∈ OT(I ) v . Let ( = (V, E, I → V ), v) = ψ t (
, ) be the combined graph in OT(I ) t as described in Lemma 5.4.

The datum of two twist assignments p and p on and respectively together with a pair of positive integers ( p e , p e ) is equivalent to the datum of a twist assignment p : V → Z >0 on the graph . Moreover, the contributions of the vertices of with the twist assignment p are given by for the two distinguished vertices of and , and

• H v for all other vertices v.

One checks that this is still true if one of the graphs or has only one marking. In summary we obtain that

H I = 1 2(n -1) I ⊂I ( ,v )∈OT(I ) v ( ,v )∈OT(I \I ) v f # ( ) f # ( ) p : E→Z >0 z p(e) e ⎛ ⎝ v∈V H v ⎞ ⎠ = 1 2(n -1) I ⊂I ( ,v )∈OT(I ) v ( ,v )∈OT(I \I ) v f # ( ) f # ( )cont(ψ t ( , v , , v )) = S I ,
where we use (64) to pass from the above second line to the third. Finally, we prove that S I = A I by a similar argument. It suffices to prove for the case I = [ [1, n]]. Using the induction hypothesis that A I j ∪{n+ j} = S I j ∪{n+ j} , we can rewrite the inductive definition (56) of A I as 123

A n = 1 (n -1) k>0,[[1,n]]={i 1 ,i 2 } I 1 ••• I k j ∈OT(I j ∪{n+ j}) k j=1 f # ( j ) k! • p=( p 1 ,..., p k ) p ( j) : E j →Z >0 [z p 1 n+1 • • • z p k n+k ] H p i 1 ,i 2 • k j=1 e∈E j z p ( j)(e) e v∈V 1 ∪•••∪V k H v .
The datum of p together with the twist assignments p ( j) for the split graphs j for 1 ≤ j ≤ k is equivalent to a twist assignment for the combined graph = ψ b ( 1 , . . . , k ) defined in Lemma 5.4. Moreover, given such a twist assignment the contribution of the vertex carrying i 1 and i 2 is

H p i 1 ,i 2 .

Thus we obtain that

A n = 1 (n -1) k>0, [[1,n]]={i 1 ,i 2 } I 1 ••• I k j ∈OT(I j ∪{n+ j}) k j=1 f # ( j ) k! • p : E( )→Z >0 e∈E( ) z p(e) e v∈V ( ) H v = 1 (n -1) k>0, [[1,n]]={i 1 ,i 2 } I 1 ••• I k j ∈OT(I j ∪{n+ j}) k j=1 f # ( j ) k! • cont(ψ b ( 1 , . . . , k )) = S(I ),
where we use (65) to pass from the above second line to the third.

Spin and hyperelliptic components

In this section we prove a refinement of Theorems 1.2 and 1.1 with spin structures taken into account. We also prove the corresponding refinement for hyperelliptic components in Sect. 6.5. Along the way we revisit the counting problem for torus covers with sign given by the spin parity and complete the proof of Eskin, Okounkov and Pandharipande [START_REF] Eskin | The theta characteristic of a branched covering[END_REF] that the generating function is a quasimodular form of the expected weight. We then show that the D 2 -recursion has a perfect analog when counting with spin parity and use the techniques of Sect. 5 to convert this into the recursion for intersection numbers.

In this section we assume that all entries of μ = (m 1 , . . . , m n ) are even. The spin parity of a flat surface (X, x 1 , . . . , x n , ω) ∈ M g,n (μ) is defined as

φ(X, ω) = h 0 X, n i=1 m i 2 x i mod 2 .
The parity is constant in a connected family of flat surfaces by [START_REF] Mumford | Theta characteristics of an algebraic curve[END_REF]. We will denote by M g,n (μ) • with • ∈ {odd, even} the moduli spaces of flat surfaces with a fixed odd or even spin parity and simply call them spin (sub)spaces of the corresponding strata. Moreover, we will denote by M g,n (μ) • their incidence variety compactification and will similarly use this symbol, e.g. in the form v(μ) • , c 1↔2 (μ, C) • , and c area (μ) • for volumes and Siegel-Veech constants.

We remark that with our choice of notation the spin spaces M g,n (μ) • are not necessarily connected. Indeed, for μ = (2g -2) with g ≥ 4 and μ = (g-1, g-1) with g ≥ 5 odd, one of the two spin spaces is disconnected, since it contains an extra hyperelliptic component (see [START_REF] Kontsevich | Connected components of the moduli spaces of Abelian differentials with prescribed singularities[END_REF]Theorem 2]). Nevertheless, the hyperelliptic components will be treated separately in Sect. 6.5. Taking the difference of the volumes thus gives a formula for the volume of each connected component.

To state the refined version of the volume recursion we need a generalization of the spin parity. Let ( , , p) be a backbone graph. A spin assignment is a function

φ : {v ∈ V ( ), (v) = 0} → {0, 1}.
The parity of the spin assignment is defined as

φ( ) := v∈V ( ), (v)=0 φ(v) mod 2. ( 66
)
Our goal is the following refinement of Theorems 1.2 and 1.1 under a mild assumption. Recall that the tautological line bundle O(-1) over P M g,n (μ) has a natural hermitian metric given by the area form h(X, ω) = i 2 X ω ∧ ω. Assumption 6.1 There exists a desingularization f : Y → P M g,1 (2g -2) such that f * h extends to a good hermitian metric on f * O(-1).

This assumption was already present in [START_REF] Sauvaget | Volumes and Siegel-Veech constants of H(2g-2) and Hodge integrals[END_REF] and can be deduced from the smooth strata compactification constructed in [START_REF] Bainbridge | The moduli space of multi-scale differentials[END_REF]. 4 Note that we do not need this assumption for Theorem 1.2, as it is stated for the entire stratum whose cohomology class was computed recursively in [START_REF] Sauvaget | Cohomology classes of strata of differentials[END_REF]. However, currently we do not know the cohomology class of each individual spin subspace.

Theorem 6.2 If n ≥ 2, then the rescaled volumes satisfy the recursion

v(μ) odd = k≥1 g,μ, φ odd h P 1 ((m 1 , m 2 ), p) • k i=1 (2g i -1 + n(μ i ))!v(μ i , p i -1) φ(i) 2 k-1 k!(2g -3 + n)! ,
where the summation conventions for g, μ and p are as in Theorem 1.2 and the superscript φ(i) indicates the corresponding spin subspace.

We remark that the same formula holds when replacing "odd" by "even" in the theorem, which follows simply by subtracting the formula in Theorem 6.2 from that in Theorem 1.2. Theorem 6.3 Let P M g,n (μ) • with • ∈ {odd, even} be the subspaces of P M g,n (μ) with a fixed spin parity. Then the volume can be computed as an intersection number

vol( M g,n (μ) • ) = - 2(2iπ) 2g (2g -3 + n)! P M g,n (μ) • ξ 2g-2 • n i=1 ψ i .
We first show in Sect. 6.1 that the intersection numbers on the right-hand side of Theorem 6.3 satisfy a recursion as in Theorem 6.2. This is parallel to Sect. 3. We then complete in Sect. 6.2 properties of the strict brackets introduced by [START_REF] Eskin | The theta characteristic of a branched covering[END_REF]. The volume recursion in Sect. 6.3 is parallel to Sect. 4 and allows efficient computations of volume differences of the spin subspaces. We do not need to prove an analog of Sect. 5 but can rather apply the results, since the structures of the two recursions are exactly the same as before. Only in Sect. 6.4 we need Assumption 6.1 to prove the initial case of Theorem 6.2, i.e. the case of the minimal strata.

Intersection theory on spin subspaces and hyperelliptic components

With a view toward Sect. 6.5 for the hyperelliptic components, we allow here also the profile μ = (g -1, g -1) (with g -1 not necessarily even) and • ∈ {odd, even, hyp}, and study the corresponding union of connected components P M g,n (μ) • .

Let ( , , p) be a twisted bi-colored graph and D be an irreducible component of the boundary P M 123 Proposition 6.6 For μ of length bigger than one and with even entries, we have

(m 1 + 1)(m 2 + 1)a(μ) odd = ( , ,p,φ), φ odd h P 1 ((m 1 , m 2 ), p) |Aut( , , p, φ)| • v∈V ( ), (v)=0 p v (2g v -1 + n(μ v ))a(μ v , p v -1) φ(v) ,
where the sum is over all choices of backbone graphs with only the first two marked points in the lower level component.

This proposition is a refined combination of Lemma 3.13 and Eq. ( 25). Again we remark that the same formula holds when replacing "odd" by "even" in the proposition, which simply follows from subtracting the above from the corresponding formula for the entire stratum.

Proof We apply Proposition 6.5 to P M g,n (μ) odd . The proposition then follows from the description of the boundary divisors of connected components of P M g,n (μ).

Let (L → X → ) be a one-parameter family of theta characteristics, i.e. L is a line bundle such that L| ⊗2 X ω X for every fiber curve X parametrized by a complex disk (centered at the origin) such that

• the restriction of X to \0 is a family of smooth curves;

• the central fiber X 0 is of compact type.

The second condition above is due to the fact that the graphs in BB(g, n) 1,2 are of compact type. We assume that L is odd, i.e. L restricted to every smooth fiber is an odd theta characteristic. The restriction of L to each irreducible component of X 0 (minus the nodes) is a theta characteristic of that component. Since X 0 is of compact type, the parity of L| X 0 equals the sum of the parities over all irreducible components of X 0 (see e.g. [10, Proposition 4.1]), which implies that the number of components of X 0 with an odd theta characteristic is odd.

Let ( , , p) ∈ BB(g, n) 1,2 . From the above description we deduce that P M p,odd , can be written as

φ odd M -1 × v∈V ( ) P M g v ,n v ( p v -1, μ v ) φ(v)
123

where M -1 and P M g v ,n v ( p v -1, μ v ) are defined as in Sect. 3.2. The arguments in the proof of Lemma 3.12 imply that

(m 1 + 1)(m 2 + 1)a(μ) odd = ( , ,p)∈BB 1,2 φ odd h P 1 ((m 1 , m 2 ), p) |Aut( , , p, φ)| • v∈V ( ), (v)=0 p 2 v a(μ v , p v -1) φ(v) .
Then by the same line of arguments as in Sect. 3.5 (expansions over rooted trees), we get the desired expression.

Strict brackets and Hurwitz numbers with spin parity

Let f : SP → Q be any function on the set of strict partitions (i.e. partitions with strictly decreasing part lengths). The replacement of the q-bracket in the context of spin-weighted counting is the strict bracket defined by

f str = 1 (q) ∞ λ∈SP (-1) (λ) f (λ)q |λ| , (q) ∞ = n≥1 (1 -q n ) = λ∈SP (-1) (λ) q |λ| .
The analog of the algebra * is the algebra * = Q[p 1 , p 3 , p 5 , . . .] of supersymmetric functions, where for odd the functions p are defined by

p (λ) = ∞ i=1 λ i - ζ(-) 2 .
Note the modification of the constant term and the absence of the shift in comparison to [START_REF] Ivanov | The Gaussian limit for projective characters of large symmetric groups[END_REF]. We provide * with the weight grading by declaring p to have weight + 1. On the other hand, [START_REF] Eskin | The theta characteristic of a branched covering[END_REF]

= (μ 1 , . . . , μ n ) p (-1) φ( p) q deg( p) |Aut( p)| = 2 n i=1 ( (μ i )-|μ i |)/2 f μ 1 f μ 2 • • • f μ n str , ( 67
)
123
where the sum is over all covers p : X → E of a fixed base curve and profile . then the difference fh has weight strictly less than + 1. In particular f belongs to the subspace * ≤ +1 of weight less than or equal to + 1. More precisely,

f ( ) = -1 2 [t +1 ] -1 j=1 (1 -jt) • exp j odd 2p j t j j (1 -(1 -t) -j ) . ( 69
)
This statement was missing in the proof of the following corollary, one of the main theorems of [START_REF] Eskin | The theta characteristic of a branched covering[END_REF].

Corollary 6.8 The strict bracket f 1 f 2 • • • f n str
is a quasimodular form of mixed weight less than or equal to n i=1 ( i + 1).

We now prepare for the proof of Theorem 6.7 and prove the corollary along with more precise statements on strict brackets in the next subsection. From [30, Definition 6.3 and Proposition 6.4] we know that the central characters are given by

f ρ = μ∈SP X ρ μ P ↓ μ ( 70
)
where the objects on the right-hand side are defined as follows. We define for any partition λ the Hall-Littlewood symmetric polynomials

P λ (x 1 , . . . , x m ; t) = σ ∈S n σ n i=1 x λ i i i< j, i< (λ)
x it x j x ix j .

These polynomials have cousins where the powers are replaced by falling factorials. That is, writing

n ↓k = n(n -1)(n -2) • • • (n -k + 1), we define P ↓ λ (x 1 , . . . , x m ; t) = σ ∈S n σ n i=1 x ↓λ i i i< j, i< (λ)
x it x j x ix j .

Next, we define X • • (t) to be the base change matrix from the basis of p ρ to the basis P λ (x 1 , . . . , x m ; t), that is, we define them by

p ρ = λ |ρ| X ρ λ (t)P λ ( • ; t) . ( 71
)
The existence and the fact that the X Proof of Theorem 6. [START_REF] Bainbridge | The moduli space of multi-scale differentials[END_REF] We need to prove (69). From there one can then derive (68) by expanding the exponential function (just as in [START_REF] Ivanov | Kerov's central limit theorem for the Plancherel measure on Young diagrams[END_REF] Proposition 3.5 is derived from Proposition 3.3). We use that for ρ = ( ) a cycle, the coefficients X λ ρ in (70) are supported on λ with at most two parts. More precisely, by [START_REF] Macdonald | Symmetric functions and Hall polynomials[END_REF]Example III.7.2] we know that

f ( ) = P ↓ ( ) + 2 /2 i=1 (-i) i P ↓ ( -i,i) = 1≤a,b≤ (λ) a =b i=0 (-1) i λ ↓ -i a λ ↓i b λ a + λ b λ a -λ b i =a,b λ a + λ i λ a -λ i λ b + λ i λ b -λ i . ( 72
)
Using

- j∈N 1 -(-1) j p j (λ)t j j (1 -t) -j = log (λ) i=1 1 -(λ i + )t 1 + (λ i -)t
and the specialization of this formula for = 0, our goal is to show that

f ( ) (λ) = -1 2 [t 1 ] -1 j=0 (t -1 -j) (λ) i=1 t -1 + λ i t -1 -λ i t -1 -(λ i + ) t -1 + (λ i -) = -1 2 [z -1 ] -1 j=0 (z -j) (λ) i=1 z + λ i z -λ i z -(λ i + ) z + (λ i -) = (λ) a=1 -1 j=0 (λ a -j) i =a λ a + λ i λ a -λ i λ a -(λ i + ) λ a + (λ i -) .
Using that for odd

i=0 (-1) i λ ↓ -i a λ ↓i b -λ ↓ -i b λ ↓i a = λ ↓ a (λ a -λ b -) -λ ↓ b (λ b -λ a -) λ a + λ b - ,
we see that our goal and the known (72) agree.

Volume computations via cumulants for strict brackets

We denote by an upper index the difference of the even and odd spin related quantities, e.g. v(μ) = v(μ) even -v(μ) odd . Cumulants for strict brackets are defined by the same formula (39) as for q-brackets. We are interested in cumulants for the same reason as we were for the case of the strata as in [START_REF] Veech | Siegel measures[END_REF].

Proposition 6.9 The difference of the volumes of the even and odd spin subspaces of M g,n (μ) can be computed in terms of cumulants by

vol( M g,n (μ)) = (2πi) 2g (2g -2 + n)! f (m 1 +1) | • • • |f (m n +1) str,L ,
and thus, in combination with Theorem 6.7 we have

v(μ) = (2πi) 2g (2g -2 + n)! h (m 1 +1) | • • • |h (m n +1) str,L .
Here the subscript L refers to the leading term

g 1 | • • • |g n str,L = [ -k-1+n ] g 1 | • • • |g n str, = lim h→0 k+1-n ev[ g 1 | • • • |g n str,q ]( )
for g i homogeneous of weight k i and k = n i=1 k i .

This proposition was certainly the motivation of [START_REF] Eskin | The theta characteristic of a branched covering[END_REF], which stops short of this step. To derive the proposition from (67), we need one more tool, the analog of the degree drop in Proposition 4.4. We use the fact that for strict brackets we have a closed formula (rather than only a recursion as for q-brackets), proved in [22, Section 3.2.2] and in more detail in [START_REF] Bloch | The character of the infinite wedge representation[END_REF]Section 13]. First,

(-1) • p str = G +1 := ζ(-) 2 + n≥1 σ (n)q n (73)
and for the more general statement we define the "oddification" of the Eisenstein series to be

G odd (z 1 , . . . , z n ) = - ∞ r =1 D (n-1) q G 2r s 1 +•••+s n =r +n-1 z 2s 1 -1 1 • • • z 2s n -1 n (2s 1 -1)! • • • (2s n -1)! ,
where D q = q∂/∂q. Then by Proposition 13.3 in loc. cit. the n-point function is given by i ≥1, i odd

p 1 p 2 • • • p n str z 1 1 • • • z n n 1 ! • • • n ! = α∈P(n) A∈α G odd {z a } a∈A . ( 74
)
Consequently, the cumulants are simply given by

p 1 |p 2 | • • • |p n str = z 1 1 • • • z n n 1 ! • • • n ! G odd (z 1 , . . . , z n ) . ( 75
)
Proof of Proposition 6.9 Recall the evaluation map Ev used in [START_REF] Chen | Quasimodularity and large genus limits of Siegel-Veech constants[END_REF]Section 8] and in Sect. While (75) provides an easy and efficient way to compute cumulants of strict brackets, we show that the more complicated way via lifting of differential operators to * and the Key Lemma 4.5 also works here. The analog of Proposition 4.3 is the following result.

Proposition 6.10 With

( f ) = 1 , 2 ≥ 1 ( 1 + 2 ) p 1 + 2 -1 ∂ 2 ∂p 1 ∂p 2 we have f str, = 1 k (e ( -∂/∂p 1 )/2 f ) (∅), ( 76
)
where the evaluation at the empty set is explicitly given by p → -ζ(-) 2 . Note that we can regard the differential operator ( -∂/∂p 1 )/2 appearing in the exponent the same as the operator D defined in [START_REF] Wright | Cylinder deformations in orbit closures of translation surfaces[END_REF] when viewing * as a quotient algebra of * with all the even p set to zero, since the differential operator ∂ sending p to a multiple of p -1 is zero on this quotient.

Proof Using the description (37) of the -evaluation we need to show that d f str, = ( -∂/∂p 1 ) f str, . Contrary to the case of q-brackets we will actually show the stronger statement that d f str = ( -∂/∂p 1 ) f str . It suffices to check this for all the n-point functions. For n = 1 this can be checked directly from (73). For general n, we write W (z) = s≥1 z 2s-1 /(2s -1)!. Using (74) and that the commutator [d, D q ] is multiplication by the weight, we compute that

d n i=1 W (z i ) str = α∈P(n) A 1 ∈α, |A 1 |≥2 dG odd {z a } a∈A 1 • A∈α\{A 1 } G odd {z a } a∈A - 1 2 n i=1 z i • α∈P({1,...,n}\{i}) A∈α G odd {z a } a∈A ( 77
)
where for the factor in the summand with

|A 1 | ≥ 2 dG odd {z a } a∈A 1 = - 1 2 r ≥1 2r D |A 1 |-2 q G 2r • s a ≥1, s a =r +|A 1 |-1 a∈A 1 z 2s a -1 a (2s a -1)! , ( 78
)
and where the summation is over all tuples (s a ) a∈A 1 . Since 1 2 ∂/∂p 1 ( n i=1 W (z i )) = 1 2 n i=1 z i j =i W (z j ), the strict bracket of this expression is precisely the second line on the right-hand side of (77). Since

1 2 n i=1 W (z i ) = 1≤i = j≤n (z i + z j )W (z i + z j ) k∈{1,...,n}\{i, j} W (z k ),
its strict bracket matches the first line on the right-hand side of (77), and the part containing the variable for W (z i + z j ) produces of course the special factor (78).

This proposition provides an efficient algorithm to compute the differences of volumes of the spin subspaces. The definitions below are completely analogous to the beginning of Sect. 4, except that objects with even indices have disappeared and they are written in boldface letters for distinction. For the substitution, we define

P Z (u) = exp j≥1 1 2 j+1 2 ζ(-j) u j+1 and α = u 1 (u/P Z (u)) -1 .
We let R = Q[h 1 , h 3 , . . .] and define for a finite set I = {i 1 , . . . , i n } of positive integers the formal series

H I ∈ R[[z i 1 , . . . , z i n ]] by H {i} = 1 z i + ≥1 h z i , H {i, j} = z i H (z i ) -z j H (z j ) H(z j ) -H(z i ) -1, H I = 1 2(n -1)
I =I I D 2 (H I , H I ), ( 79
)
with

D 2 ( f, g) = 1 , 2 ≥1, odd [z 1 1 z 2 2 ]H {1,2} ∂ f ∂h 1 ∂g ∂h 2 .
We still set

H n = H [[1,n]] and h 1 ,..., n = [z 1 1 • • • z n n ]H n .

Corollary 6.11

The even-odd volume differences of the stratum with signature μ = (m 1 , . . . , m n ) can be computed as

v(μ) = (2πi) 2g (2g -2 + n)! h m 1 +1,...,m n +1 h →α
using the recursion (79).

Proof Thanks to Proposition 6.10 and the subsequent remark, the proof of Theorem 4.1 can be copied verbatim here. The extra factor 2 -/2 in the definition of P Z in comparison to the constant term -ζ(-)/2 of the evaluation of p compensates for the fact that the strict bracket of the f gives the counting function in (67) up to a power of two.

Conclusion of the proofs for spin subspaces

Proof of Theorems 6.3 and 6.2 Theorem 6.2 is a consequence of Corrolary 6.11. Indeed the arguments of Sect. 5 adapted to the series H n show that the recursion in Theorem 6.2 is a consequence of the recursion in (79).

To prove Theorem 6.3, we consider first the case n = 1, i.e. μ = (2g-2). Let ν μ be the push-forward of the Masur-Veech volume ν μ form to P M g,n (μ). Assumption 6.1 implies by the same argument as in [42, Lemma 2.1] that 2(2iπ) 2g (2g-1)! ξ 2g-1 can be represented by a meromorphic differential form (of Poincaré growth at the boundary), whose restriction to

P M g,1 (2g -2) is equal to ν μ . This implies that vol(2g -2) • = 2(2iπ) 2g (2g -1)! P M g,1 (2g-2) • ξ 2g-1 .
Now for the case n ≥ 2 Theorem 6.2 and Proposition 6.6 determine a(μ) • and v(μ) • by the equivalent recursive formulas, and hence they coincide up to the obvious normalizing factors.

Volume recursion for hyperelliptic components

In this subsection we prove the volume recursion for hyperelliptic components, which is analogous to but not quite the same as the recursion in Theorem 1.2.

It is a consequence of the work of Athreya, Eskin and Zorich [START_REF] Athreya | Right-angled billiards and volumes of moduli spaces of quadratic differentials on CP 1[END_REF] on volumes of the strata of quadratic differentials in genus zero.

Recall that only the strata M g (g -1, g -1) and M g (2g -2) have hyperelliptic components. For the hyperelliptic components we still have an interpretation of their volumes as intersection numbers as well as a volume recursion as follows. Theorem 6.12 For the hyperelliptic components we have

vol( M g,1 (2g -2) hyp ) = 2(2iπ) 2g (2g -1)! P M g,1 (2g-2) hyp ξ 2g-1 and vol( M g,2 (g -1, g -1) hyp ) = 2(2iπ) 2g g(2g -1)! P M g,2 (g-1,g-1) hyp ξ 2g-1 ψ 2 ,
provided that Assumption 6.1 holds.

As before we set

v(2g -2) hyp = (2g -1)vol ( M g,1 (2g -2) hyp ), v(g -1, g -1) hyp = g 2 vol ( M g,2 (g -1, g -1) hyp ) .
Proposition 6. [START_REF] Costantini | The area is a good metric[END_REF] The Masur-Veech volumes of the hyperelliptic components of M g,2 (g -1, g -1) satisfy the recursion

v(g-1, g-1) hyp = v(2g-2) hyp + g-1 =1 (2 -1)!v(2 -2) hyp (2(g-)-1)!v(2g-2 -2) hyp 4(2g -1)! .
Note in comparison to Theorem 1.2 that only the terms k = 1 and k = 2 appear and that the Hurwitz number h P is identically one here. As a preparation for the proof recall that the canonical double cover construction provides isomorphisms

Q g (2g -3, (-1) 2g-3 ) ∼ = M g,1 (2g -2) hyp , Q g (2g -2, (-1) 2g-2 ) ∼ = M g,2 (g -1, g -1) hyp
that preserve the Masur-Veech volume and the SL 2 (R)-action. Taking into account the factorials for labeling zeros and poles the main result of [START_REF] Athreya | Right-angled billiards and volumes of moduli spaces of quadratic differentials on CP 1[END_REF] can be translated as

vol ( M g,1 (2g -2) hyp ) = 2 (2g + 1)! (2g -3)!! (2g -2)!! π 2g , vol ( M g,2 (g -1, g -1) hyp ) = 8 (2g + 2)! (2g -2)!! (2g -1)!! π 2g
where the double factorial notation means (2k

)!! = 2 k k! and (2k -1)!! = (2k)!/2 k k!.
Proof Expanding the definition of the double factorials and including the summand v(2g-2) hyp as the two boundary terms of the sum (i.e. = 0 and = g), we need to show that

g =0 1 2 + 1 2 1 2g -2 + 1 2(g -) g - = 2 16 g (g + 1) 2 2g + 2 g + 1 -1 . ( 80
)
For this purpose it suffices to prove the following two identities of generating series

≥0 1 2 + 1 2 x 2 = 1 2x arctan 2x √ 1 -4x 2 (81)
and

2 g≥0 16 g g 2 2g g -1 x 2g = 1 4x 2 arctan 2x √ 1 -4x 2 2 , ( 82
)
so that we can take the square of the first series and compare the x 2g -terms. To prove (81) we multiply it by x, differentiate, and are then left with showing that

123 ≥0 2 x 2 = 1/ √ 1 -4x 2
, which follows from the binomial theorem. To prove (82) we differentiate and are then left with the identity which is already proved in [START_REF] Lehmer | Interesting series involving the central binomial coefficient[END_REF]p. 452,Equation (9)].

The last ingredient is the following straightforward consequence of Proposition 6.5 (analogous to the case of spin subspaces in Proposition 6.6). Proposition 6.14 For μ = (g -1, g -1), we have

g 2 a(g -1, g -1) hyp = (2g -1) 2 a(2g -2) hyp + 1 2 g-1 g 1 =1 h P 1 (g -1, g -1), (2g 1 -1, 2g -2g 1 -1) • (2g 1 -1) 2 a(2g 1 -2) hyp (2g -2g 1 -1) 2 a(2g -2g 1 -2) hyp . ( 83
)
Proof of Theorem 6.12 Since the proof of Theorem 6.12 for the case μ = (2g -2) was already given along with the proof of Theorem 6.3, it remains to show that v(g -1, g -1) hyp and a(g -1, g -1)

hyp = P M g,2 (g-1,g-1) hyp β i •ξ
satisfy the same recursion. It is elementary to check that h P 1 (g -1, g -1), (2g 1 -1, 2g -2g 1 -1) = 1. Then the desired conclusion thus follows from Propositions 6.13 and 6.14.

An overview of Siegel-Veech constants

Let (X, ω) be a flat surface, consisting of a Riemann surface X and an Abelian differential ω on X . Siegel-Veech constants measure the asymptotic growth rate of the number of saddle connections (abbreviated s.c.) or cylinders with bounded length (of the waist curve) in (X, ω). There are many variants that we now introduce and compare.

Saddle connection and area Siegel-Veech constants

For each pair of zeros (z 1 , z 2 ) of ω we let

A phy 1↔2 (
where a collection of homologous saddle connections just counts for one. Quadratic upper and lower bounds for such counting functions were established by Masur [START_REF] Masur | The growth rate of trajectories of a quadratic differential[END_REF]. Fundamental works of Veech [START_REF] Veech | Siegel measures[END_REF] and Eskin-Masur [START_REF] Eskin | Asymptotic formulas on flat surfaces[END_REF] showed that for almost every flat surface (X, ω) in the sense of the Masur-Veech measure (see [START_REF] Masur | Interval exchange transformations and measured foliations[END_REF] and [START_REF] Veech | Gauss measures for transformations on the space of interval exchange maps[END_REF]) there is a quadratic asymptotic, i.e. that

A phy 1↔2 (T ) ∼ c phy 1↔2 (X, ω) π T 2 , A hom 1↔2 (T ) ∼ c hom 1↔2 (X, ω) π T 2 . ( 86
)
The constants c phy 1↔2 (X, ω) and c hom 1↔2 (X, ω) are the first type of Siegel-Veech constants we study here, called the saddle connection Siegel-Veech constants. The difference between these two Siegel-Veech constants becomes negligible as the genus of X tends to infinity, which follows from the results of Aggarwal and Zorich (see [START_REF] Aggarwal | Large genus asymptotics for Siegel-Veech constants[END_REF]Remark 1.1]).

The second type of Siegel-Veech constants counts homotopy classes of closed geodesics, or equivalently flat cylinders. Again, there are two variants, the naive count and the count where each cylinder is weighted by its relative area. As above, the most important counting function with good properties (see e.g. [START_REF] Chen | Quasimodularity and large genus limits of Siegel-Veech constants[END_REF]) and connection to Lyapunov exponents [START_REF] Eskin | Sum of Lyapunov exponents of the Hodge bundle with respect to the Teichmüller geodesic flow[END_REF] is the second variant. For the precise definition we consider

A cyl (T ) = Z ⊂X cylinder w(Z )≤T 1, A area (T ) = Z ⊂X cylinder w(Z )≤T area(Z ) area(X ) , (87)
where w(Z ) denotes the width of Z , i.e. the length of its core curve. We then define the cylinder Siegel-Veech constant and the area Siegel-Veech constant by the asymptotic equalities

A cyl (T ) ∼ c cyl (X, ω)π T 2 , A area (T ) ∼ c area (X, ω)π T 2 . ( 88
)
There is a natural action of GL 2 (R) on the moduli space of flat surfaces M g and the orbit closures are nice submanifolds, in fact linear in period coordinates by the fundamental work of Eskin, Mirzakhani and Mohammadi [START_REF] Eskin | Invariant and stationary measures for the SL(2, R) action on moduli space[END_REF][START_REF] Eskin | Isolation, equidistribution, and orbit closures for the SL(2, R) action on moduli space[END_REF]. We refer to them as affine invariant manifolds, using typically the letter M. The intersection with the hypersurface of area one flat surfaces (denoted by the same symbol M) comes with a finite SL 2 (R)-invariant ergodic measure ν M with support M. This measure is unique up to scale and for affine invariant manifolds defined over Q there are natural choices of the scaling.

It follows from the Siegel-Veech axioms (see [START_REF] Eskin | Asymptotic formulas on flat surfaces[END_REF]) that Siegel-Veech constants for almost all flat surfaces (X, ω) in an SL 2 (R)-orbit closure M agree. We call these surfaces generic (for M) and write e.g. c 1↔2 (M) = c 1↔2 (X, ω) for (X, ω) generic.

The relevant orbit closures in this paper are the connected components of the strata of Abelian differentials and certain Hurwitz spaces inside the strata. We usually abbreviate by c 1↔2 (μ) = c 1↔2 ( M g,n (μ)) the Siegel-Veech constants for the strata with signature μ.

Configurations and the principal boundary

One of the main insights of [START_REF] Eskin | Moduli spaces of Abelian differentials: the principal boundary, counting problems, and the Siegel-Veech constants[END_REF] is that Siegel-Veech constants can be computed separately according to topological types, called configurations. We formalize their notion of configurations briefly so that it also applies to Hurwitz spaces, and in fact to all SL 2 (R)-orbit closures M provided with the generalization of the Masur-Veech measure ν M . The concept of configurations will be used for showing the equivalence between Theorems 1.3 and 1.2 in Sect. 8.

Let (X, z 1 , . . . , z n ) be a pointed topological surface. A configuration C of saddle connections joining z 1 and z 2 for M is a set of simple non-intersecting arcs from z 1 to z 2 up to homotopy preserving the cyclic ordering of the arcs both at z 1 and z 2 . The last condition implies that the tubular neighborhood of the configuration is a well-defined subsurface of X , in fact a ribbon graph R(C) associated with the configuration. (Here the subsurfaces of X after removing a ribbon graph are allowed to have positive genus.) The number of arcs in the configuration is called the multiplicity of the configuration.

We say that a subset of the saddle connections joining z 1 and z 2 on a flat surface (X, ω) belongs to the configuration C, if this subset is homotopic to C. Each configuration gives rise to a counting function A 1↔2 (T, C) for sets of saddle connections belonging to the configuration and having individual length ≤ T . From the counting function we deduce the corresponding Siegel-Veech constant c 1↔2 (M, C), where ∈ {phy, hom} respectively.

A configuration C is relevant if c 1↔2 (M, C) > 0.
A full set of saddle connection configurations for an affine invariant manifold M is a finite set of saddle connection configurations C i , with i ∈ I such that the contributions of the configurations C i sum up to the full Siegel-Veech constant, i.e. such that

i∈I c 1↔2 (M, C i ) = c 1↔2 (M) (89)
for ∈ {phy, hom} respectively. Note that [START_REF] Eskin | Moduli spaces of Abelian differentials: the principal boundary, counting problems, and the Siegel-Veech constants[END_REF]Section 3.2] in their definition of configurations made a further subdivision of the notion by adding metric data, i.e. specifying angles between saddle connections. In that context, Eskin, Masur and Zorich determined a full set of saddle connection configurations for the strata and used the Siegel-Veech transform to connect the computation to volume computations. The following statement summarizes Proposition 3.3, Corollary 7.2 and Lemma 8.1 of [START_REF] Eskin | Moduli spaces of Abelian differentials: the principal boundary, counting problems, and the Siegel-Veech constants[END_REF]. Proposition 7.1 For any stratum M g,n (μ) a full set of saddle connection configurations is the set of collections of pairwise homologous simple disjoint arcs joining z 1 and z 2 (up to homotopy).

In this proposition, several configurations are irrelevant, for example those with a connected component of genus zero after removing the saddle connections in the configuration.

The general strategy to compute Siegel-Veech constants is the following relation to volumes, where the submanifold M is in a stratum with labeled zeros.

Proposition 7.2 The saddle connection Siegel-Veech constants of an affine invariant manifold M can be computed as

c 1↔2 (M) = lim ε→0 1 πε 2 C m (C) ν M (M ε (C)) ν M (M) , ( 90
)
where the sum runs over the full set of saddle connection configurations and where m hom (C) = 1 for all C while m phy (C) is equal to the number of arcs in C.

Proof This is a direct consequence of the Siegel-Veech transform applied to the characteristic function of a disc of radius ε, see [START_REF] Eskin | Moduli spaces of Abelian differentials: the principal boundary, counting problems, and the Siegel-Veech constants[END_REF]Lemma 7.3] together with the Eskin-Masur bound on the number of short saddle connections (Theorem on p. 84 of [START_REF] Eskin | Moduli spaces of Abelian differentials: the principal boundary, counting problems, and the Siegel-Veech constants[END_REF]).

We conclude with remarks on Siegel-Veech constants for general affine invariant manifolds to put the digression on Hurwitz spaces (Sect. 9) in context. There is another variant, besides (84) and (85), of counting saddle connections. Given an affine invariant manifold M we say that two saddle connections on (X, ω) ∈ M are M-parallel if they are parallel and stay parallel in a neighborhood of (X, ω) in M. (The terminology is completely analogous to the notion of M-parallel cylinders introduced in [START_REF] Wright | Cylinder deformations in orbit closures of translation surfaces[END_REF].) We thus define the counting function A

1↔2 (M) = c hom 1↔2 (M) if M is
a connected component of a stratum. For Hurwitz spaces the two values can be different, but we will see (Proposition 9.3) that their difference becomes negligible as the degree of the covers tends to infinity.

In the first part of [START_REF] Eskin | Moduli spaces of Abelian differentials: the principal boundary, counting problems, and the Siegel-Veech constants[END_REF] on recursive computations of Siegel-Veech constants, Eskin, Masur and Zorich called the locus of degenerate surfaces that contribute to the Siegel-Veech counting the principal boundary. At that time the notion of principal boundary was used only as a partial topological compactification. Presently, we dispose of a complete and geometric compactification for the strata [START_REF] Bainbridge | Compactification of strata of abelian differentials[END_REF] and for Hurwitz spaces (by admissible covers), and we can then identify the principal boundary as part of the compactification (see [START_REF] Chen | Principal boundary of moduli spaces of abelian and quadratic differentials[END_REF] for the case of the strata and Sect. 9 for the case of Hurwitz spaces). The reader should keep in mind that the locus "principal boundary" depends on the type of saddle connections under consideration.

Finally we remark that there is a zoo of possibilities of associating weights with saddle connections and cylinders and to define Siegel-Veech constants accordingly. This started with [START_REF] Vorobets | Ergodicity of billiards in polygons[END_REF], and see also [START_REF] Bauer | Geometry of periodic regions on flat surfaces and associated Siegel-Veech constants[END_REF] for computations and conversions of Siegel-Veech constants.

Saddle connection Siegel-Veech constants

In this section we deduce from the volume recursion and its refinement for spin and hyperelliptic components a proof of Theorem 1.3. Almost all we need here has been proven already in [START_REF] Eskin | Moduli spaces of Abelian differentials: the principal boundary, counting problems, and the Siegel-Veech constants[END_REF]. We start with two more auxiliary statements.

Proposition 8.1

The full set of saddle connection configurations for the strata given in Proposition 7.1 is in bijection with (possibly unstable) backbone graphs and a cyclic ordering of its vertices at level zero. The subset of relevant configurations is in bijection with stable backbone graphs.

Although not needed in the sequel, we relate for the convenience of the reader our notion of twists and the angle information that [START_REF] Eskin | Moduli spaces of Abelian differentials: the principal boundary, counting problems, and the Siegel-Veech constants[END_REF] recorded. Let {γ i } k i=1 be the arcs of a configuration realized by (X, ω) labeled cyclically and let a i and a i be the angles between γ i and γ i+1 at z 1 and z 2 respectively. If an edge e is separated by the loop formed by γ i and γ i+1 , then the twist is

p(e) = 1 2π (a i + a i ) -1 . ( 91
)
The above proposition can be seen as follows. Given a collection of k homologous short saddle connections there is a sphere (with z 1 as its south pole and z 2 as its north pole, see [START_REF] Eskin | Moduli spaces of Abelian differentials: the principal boundary, counting problems, and the Siegel-Veech constants[END_REF]Figure 5]) supporting the k saddle connections. This sphere is the source of the backbone of the graph. The components at level zero are bounded by the arcs γ i and γ i+1 . The converse is obvious, given that all the edges of a (stable) backbone graph are separating by definition. Finally formula (91) is just a restatement of the Gauss-Bonnet theorem.

Recall that a backbone graph (being of compact type) is compatible with a unique twist p(•). If the vertices at level zero are labeled as 1, . . . , k as usual and if (h j , i(h j )) is the edge connecting the jth vertex to level -1, we write p j = |p(h j )| as we did in Sect. 3. Proposition 8.2 For each configuration C corresponding to a backbone graph , the volume of the subspace M ε g,n (μ, C) satisfies k≥1 g,μ

πε 2 h P 1 ((m 1 , m 2 ), p) k! (2g -3 + n)! • k i=1 (2g i -1 + n i )! p i vol( M g i ,n i +1 (μ i , p i -1)) = 2 1-k vol( M ε g,n (μ, C)) + o(ε 2 )
as ε → 0, where the summation conventions and p are as in Theorem 1.2 and where n i = n(μ i ).

Proof This is mainly contained in [20, Corollary 7.2, Formulas 8.1 and 8.2], stating that the volume of the locus with an ε-short configuration is πε 2 times the volume of the corresponding boundary. We now explain the combinatorial factors that appear. First, the factor of two and the factorials result from the passage of the boundary volume element in the ambient stratum to the product of the volume elements of the components at the boundary, as explained in detail in [START_REF] Eskin | Moduli spaces of Abelian differentials: the principal boundary, counting problems, and the Siegel-Veech constants[END_REF]Section 6]. The 1/k! stems from labelling the level zero vertices. Second, we need to count the ways to obtain a surface in M ε g,n (μ, C) by gluing a collection of surfaces (X i , ω i ) in M g i ,n i +1 (μ i , p i -1).

Suppose we are given a branched cover b : P 1 → P 1 that has ramification profile = ((m 1 + 1), (m 2 + 1), ( p 1 , . . . , p k )) over the points 0, 1 and ∞. We provide the domain with the differential ω -1 = b * dz. Since this differential has no residues we can glue t • ω -1 with the surfaces (X i , ω i ) by cutting the pole of order p i + 1 and gluing it to an annular neighborhood of the zero of order p i -1 of ω i . For t ≤ ε this provides a surface in M ε g,n (μ, C), see e.g. [START_REF] Bainbridge | Compactification of strata of abelian differentials[END_REF]Section 4] for details of the construction. The plumbing construction also depends on the choice of a p i th root of unity at each pole (from the choice of a horizontal slit at a zero of order p i -1). In total there are h P 1 ((m 1 , m 2 ), p) • k i=1 p i possibilities involved in the construction, thus justifying the remaining combinatorial factors in the formula.

We claim that this construction provides a collection of maps to M ε g,n (μ, C) that are almost everywhere injections if none of the surfaces (X i , ω i ) has a period of length smaller than ε. In fact, if two such plumbed surfaces are isomorphic, this isomorphism restricts to an automorphism of (P 1 , ω -1 ) (see [START_REF] Eskin | Moduli spaces of Abelian differentials: the principal boundary, counting problems, and the Siegel-Veech constants[END_REF]Lemma 8.1] for more details) and this happens only on a measure zero set. The locus where one of the (X i , ω i ) has a short period is subsumed in the 123 o(ε 2 ) [20, Lemma 7.1]. Conversely, for each surface (X, ω) in M ε g,n (μ, C) we can cut a ribbon graph around the configuration C. The restriction of ω has no periods since the boundary curves are homologous by definition of C. It can thus be integrated and completed to a map b : P 1 → P 1 with ramification profile as above.

Proof of Theorem 1.3 We first focus on the case that M g,n (μ) is connected. A decomposition g = k i=1 g i and (m 3 , . . . , m n ) = μ 1 • • • μ k (as in Eq. ( 5) we proved) determines uniquely a configuration and the converse is true up to the labeling of the k vertices at level zero by Proposition 8.1. The configuration is relevant if and only if the volumes on the right-hand side of ( 5) are non-zero. Since the saddle connection Siegel-Veech constant is the sum of ratios of the boundary volumes over the total volume (by Proposition (7.2)), comparing the formula in Proposition 8.2 with Eq. ( 5) thus implies Theorem 1.3. More precisely, note that the rescaled volume v(μ) defined in (4) involves a product of all (m i + 1), while on the right-hand side (m 1 + 1)(m 2 + 1) is missing and this factor gives the right hand side of the desired formula in Theorem 1.3. The factor of π in Proposition 8.2 cancels with the one in Proposition 8.1.

For disconnected strata with components parameterized by S ⊆ {odd, even, hyp} the same proof gives the averaged version that

1 v(μ) •∈S v(μ) • c hom 1↔2 (μ) • = (m 1 + 1)(m 2 + 1)
by rewriting Eq. (90) as a sum over the connected components. It remains to prove the statement for the odd spin subspaces (that may be disconnected for some strata in our notation, in which case there is an extra hyperelliptic component with the same spin parity) as well as for the hyperelliptic components.

Then one obtains the desired result for each connected component by taking suitable differences.

For the odd spin subspaces we now focus on Theorem 6.2 instead of (5) and note that each decomposition of g and μ still determines uniquely a configuration. Moreover, since the Arf-invariant on stable curves of compact type is additive (see [START_REF] Chen | Principal boundary of moduli spaces of abelian and quadratic differentials[END_REF]Proposition 4.6] and [START_REF] Eskin | Moduli spaces of Abelian differentials: the principal boundary, counting problems, and the Siegel-Veech constants[END_REF]Lemma 10.1]) and since by definition the spin assignment is additive on the vertices, the configurations appearing on the right-hand side in Theorem 6.2 are precisely the configurations that contribute (in the sense of Proposition 7.2) to the Siegel-Veech constant of the odd spin subspaces. We thus obtain the analogous statement of Theorem 1.3 for the odd spin subspaces by comparing Proposition 8.2 to Theorem 6.2.

For the hyperelliptic components we use Proposition 6.13 and note that the summands on the right-hand side there correspond exactly to the configurations on hyperelliptic components, as explained in [START_REF] Chen | Principal boundary of moduli spaces of abelian and quadratic differentials[END_REF]Proposition 4.3] and [START_REF] Eskin | Moduli spaces of Abelian differentials: the principal boundary, counting problems, and the Siegel-Veech constants[END_REF]Lemma 10.3].

The viewpoint of Hurwitz spaces

This section is a digression on how to interpret the volume recursion and the saddle connection Siegel-Veech constant from the viewpoint of Hurwitz spaces. The results in this section are not needed for proving any of the theorems stated in the introduction. We will rather explain and motivate

• why the homologous count of saddle connections is more natural than the physical count from the viewpoint of intersection theory, • how to heuristically deduce the value of the saddle connection Siegel-Veech constant in Theorem 1.3 from an equidistribution of cycles in Hurwitz tuples, and • why backbone graphs correspond to configurations.

As usual μ = (m 1 , . . . , m n ) is a partition of 2g -2 and the ramification profile consists of n cycles μ (i) of length m i + 1 unless stated otherwise. Let r (μ) = 2g + n(μ) -1 be the dimension of the (unprojectivized) stratum M g (μ). Theorem 9.1 There exists a constant M(μ) such that the Hurwitz numbers N • d ( ) for connected torus covers of profile can be approximated as

M(μ) N • d ( ) = ∈BB 1,2 d k! h P 1 ((m 1 , m 2 ), p) d 1 +•••+d k =d k i=1 p i N • d i ( i , (p i )) + o(d r (μ)-1 ), ( 92
)
where \{μ (1) , μ (2) } = 1 • • • k is the decomposition of the profile according to the leg assignment in and where p = ( p 1 , . . . p k ) is the unique twist compatible with .

At the end of the section we will show by combining Theorem 1.2 together with Theorem 9.1 the following result.

Proposition 9.2

The constant M(μ) in Theorem 9.1 is (m 1 + 1)(m 2 + 1). In particular, M(μ) depends only on the first two zero orders of μ.

Indeed an independent proof of Proposition 9.2 would provide an alternative proof of Theorem 1.2 (and hence Theorem 1.1) that would bypass the complicated combinatorics in Sects. 4 and 5.

123

The strategy to prove Theorem 9.1 consists of comparing the Hurwitz number N • d ( ), that is the fiber cardinality of the forgetful map f T : H d ( ) → M 1,n to the target curve with the fiber cardinality of the extension of f T to the space of admissible covers H d ( ) = H d,g,1 ( ) over degenerate targets of the following type.

Admissible torus covers

Let E 0,{1,2} be the stable curve of genus one consisting of a P 1 -component carrying precisely the first two marked points and of an elliptic curve E carrying the remaining marked points, joint at a node q E . If p : X → E 0,{1,2} is an admissible cover, we denote by X 0 and X -1 the (possibly reducible) curves mapping to E and to P 1 respectively, both deprived of their unramified P 1components. See Fig. 3 for examples of such admissible covers.

The admissible covers of E 0,{1,2} come in two types. One possibility is that the first two branch points are in the same (hence the unique) component of X -1 . The stable dual graph of the cover is thus a graph ∈ ABB 1,2 (recall the definition of ABB graphs in Sect. 3.2) and we denote by N • d ( , E 0,{1,2} , ) the number of such covers. The second possibility is that each of the two branch points is on its own component

X (i) -1 for i = 1, 2. Consequently, the map p| X (i) -1
is a cyclic cover of degree (m i +1). By contracting the components over P 1 we see that such covers are in bijective correspondence (up to the automorphism group of size |Aut(X -1 /P 1 )| = (m 1 + 1)(m 2 + 1)) with covers of E with the profile [START_REF] Chen | Quasimodularity and large genus limits of Siegel-Veech constants[END_REF] = ((μ (1) , μ (2) ), μ (3) , . . . , μ (n) ), where the first two ramification points are piled over the same branch point.

N • d ( , E 0,{1,2} , ) = 1 k! h P 1 ((m 1 , m 2 ), p) d=d 1 +•••+d k k i=1 N • d i ( i , (p i )), ( 93
)
where i and p are associated with as in Theorem 9.1.

Proof Recall from [START_REF] Eskin | Asymptotics of numbers of branched coverings of a torus and volumes of moduli spaces of holomorphic differentials[END_REF] or the proof of [START_REF] Chen | Quasimodularity and large genus limits of Siegel-Veech constants[END_REF]Proposition 9.4] that if is the profile for a cover π with π * ω ∈ M g,n (μ), then there exist

C 1 , C 2 = 0 such that C 1 • d r (μ)-1 ≤ N • d ( ) ≤ C 2 • d r (μ)-1 as d → ∞.
Suppose that has k components on level 0, each of genus g i and with n i marked points or nodes. Let g 0 be the genus of the component on level -1, and let b = h 1 ( ).

Then b + k i=0 g i = g and k i=1 n i = n -2 + k + b .
The cover of the P 1 -component has finitely many choices independent of d. Over the elliptic component of E 0,{1,2} , the number of choices of covers has asymptotic growth given by

B d 1 +•••+d k =d d 2g i -2+n i i
for some constant B independent of d. This quantity is a polynomial of degree k i=1

(2g i -2 + n i ) + (k -1) = 2g -2g 0 -b + n -3 .
We thus conclude that the total number of admissible covers N 0 d ( , E 0,{1,2} , ) has asymptotic growth given by a polynomial of degree r (μ) -2b -2g 0 ≤ r (μ) -2, with equality attained if and only if b = g 0 = 0, i.e. if and only if ∈ BB 1,2 .

To justify Eq. ( 93) we refer to the computation of the Hurwitz numbers in Proposition 2.1 and divide by k! to account for our auxiliary labeling of the k components of X 0 .

We remark that Proposition 9.3 is not used in any other proofs in the paper, but it reveals the geometric reason for the homologous count of saddle connections behind the recursions in Sects. 3-5. The factor h P 1 ((m 1 , m 2 ), p) (possibly with 1/k! if all branches are labeled) appears in the direct count of admissible covers and in the count of configurations. There is no extra factor k in (93), which corresponds to our setting of the coefficient m hom (C) = 1 (instead of k) in (90) for homologous count of saddle connections (instead of physical count).

Proof of Theorem 9.1 We first show that

N • d ( ) -N • d ( (12) ) = ∈BB 1,2 1 k! h P 1 ((m 1 , m 2 ), p) d 1 +•••+d k =d k i=1 p i N • d i ( i , (p i )) + O(d r (μ)-2 ) . ( 94
)
To see this, note that the ramification order of f T over E 0,{1,2} at the branch through a cover π : X → E is equal to the product of ramification orders at the nodes of X . This results in the factors p i inside the product of the right-hand side and cancels the factor 1/|Aut(X -1 /P 1 )| when counting E 0,{1,2} -covers instead of counting N • d ( ( 12) ). On the other hand, since the volume of the stratum can be approximated by counting covers of profile [START_REF] Chen | Quasimodularity and large genus limits of Siegel-Veech constants[END_REF] and since the generating function of counting these covers is a quasimodular form, arguments as in [START_REF] Chen | Quasimodularity and large genus limits of Siegel-Veech constants[END_REF]Proposition 9.4] imply the existence of a constant M(μ) such that

N • d ( (12) ) N • d ( ) = d -M(μ) d + o(1/d) . ( 95
)
The combination of Eqs. ( 94) and (95) thus implies the desired formula (92).

We now address the equidistribution heuristics for saddle connection Siegel-Veech constants. Recall that N

(m i + 1)-cycles in S 2 d as d → ∞, then M(μ) = (m 1 + 1)(m 2 + 1).
Proof If the non-trivial cycles in γ 1 and γ 2 have no letter in common, then taking (α, β, γ 1 • γ 2 , γ 3 , . . . , γ n ) is a Hurwitz tuple of profile [START_REF] Chen | Quasimodularity and large genus limits of Siegel-Veech constants[END_REF] . Assuming equidistribution and comparing to the total number of Hurwitz tuples, the number of Hurwitz tuples with γ 1 and γ 2 having two letters in common is negligible, while the ratio of those cycles having one letter in common is

(m 1 + 1)(m 2 + 1)/d + o(1/d).
Example 9.5 For the reader's convenience we illustrate the contributions to the right-hand side of (94) for the stratum M 2 (1, 1) explicitly in Fig. 3.

E P 1 E P 1 E P 1

123

The picture on the left gives stable graphs in BB 1,2 , while the pictures in the middle and on the right give graphs in ABB 1,2 \BB 1,2 . The preimages of E in the middle and on the right are unramified and thus again are elliptic curves, while on the left the preimage of E is a curve of genus two.

The saddle connection Siegel-Veech counting in this case was carried out in [START_REF] Eskin | Billiards in rectangles with barriers[END_REF] in a similar way as summarized in Theorem 9.1, despite that only primitive torus covers were considered.

The principal boundary of Hurwitz spaces

We focus on saddle connections joining the first two marked zeros and determine a full set of configurations and the corresponding principal boundary of the Hurwitz spaces. We say that ∈ ABB 1,2 is realizable in H d ( ) if there is an admissible cover p : X → E 0,{1,2} whose stable graph is and such that the vertices with (v) = 0 correspond bijectively to the components of X 0 . Recall also the definition of ribbon graphs associated to configurations in Sect. 7.2.

Proposition 9.6 Associating with a configuration C the boundary curves of the ribbon graph R(C) induces a map ϕ : C → (C) from a full set of saddle connection configurations onto the subset of ABB 1,2 that is realizable in H d ( ). The image of ϕ is independent of d for d large enough. The fibers of ϕ are finite with cardinality bounded independently of d.

Moreover if a graph ∈ BB 1,2 is realizable, then the configurations in ϕ -1 ( ) are in bijection with cyclic orderings of the components at level 0.

Proof To define ϕ, we pinch the boundary curves of R(C) to obtain a pointed nodal curve. The configuration C of saddle connections remains in one component of the curve that contains z 1 and z 2 . We provide the dual graph of the curve with the level structure such that the component containing z 1 and z 2 is the unique one at level -1 and all the other components are on level 0. This way we thus obtain a graph ϕ(C) ∈ ABB 1,2 . We leave the straightforward verification of the other statements to the reader.

An application of the Riemann-Hurwitz formula shows that any configuration in ϕ -1 ( ) has multiplicity |E( )| + 2g(X -1 ). In the special case ∈ BB 1,2 (i.e. if g(X -1 ) = 0 and is of compact type), the configuration consists of k = |E( )| pairwise homologous arcs. However, the cover on the right-hand side of Fig. 3 shows that graphs in ∈ ABB 1,2 \BB 1,2 also contribute. It is not hard to give an example that the fiber cardinality of ϕ over a target graph with g(X -1 ) > 0 can indeed be larger than one, and we leave it to the reader since it is irrelevant to our applications.

Finally we address that Theorem 9.1 and an a priori knowledge that M(μ) = (m 1 + 1)(m 2 + 1) would give an alternative proof of Theorem 1.2. In terms 123 of our volume normalization, [START_REF] Chen | Quasimodularity and large genus limits of Siegel-Veech constants[END_REF]Proposition 9.4]

says that D d=1 N • d ( ) = v(μ) 2r n i=1 (m i + 1) D r + O(D r -1 log D) (96)
as D → ∞, where r = 2g + n -1. To sum the right-hand side of (92) we let

S D ( i ) = D d=1 N • d ( i ). The Euler integral 1 0 t a-1 (1 -t) b-1 dt = (a-1)!(b-1)!
(a+b-1)! used recursively implies the following result.

Lemma 9.7 Suppose that S D (

i ) = v i D r i + O(D r i -1 log D) as D → ∞
and that there exists a constant C depending on i only such that N

• d ( i ) < Cd r i -1 for i = 1, . . . , k. Then lim D→∞ d 1 +•••+d k ≤D (d 1 + • • • + d k ) N • d 1 ( 1 ) • • • N • d k ( k ) D r 1 +•••+r k +1 = k i=1 (r i !v k ) k i=1 r i (r 1 + • • • + r k + 1)! . ( 97
)
Alternative proof of Theorem 9.1 (assuming M(μ) = (m 1 + 1)(m 2 + 1)). With the abbreviation r i = 2g i + n(μ i ) we obtain from (96) that

D d i =1 N • d i ( i , (p i )) = v(μ i , p i -1) 2r i p i m i ∈μ i (m i + 1) D r i + O(D r i -1 log D) . Since r 1 + • • • + r k = 2g + n -2, Lemma 9.7 implies that d 1 +•••+d n ≤D (d 1 + • • • + d n ) k i=1 p i N • d i ( i , (p i )) = D 2g+n-1 k i=1 (2g i + n(μ i ) -1)!v(μ i , p i -1) 2 k (2g + n -1)(2g + n -3)! n i=3 (m i + 1) . ( 98
)
Summing over all backbone graphs and taking the limit after dividing by D 2g+n-1 thus implies the desired formula [START_REF] Andersen | Topolocial recursion for Masur-Veech volumes[END_REF].

Proof of Proposition 9.2 Conversely, the above argument shows that the mere knowledge of Theorem 9.1 gives the recursion in Theorem 1.2 with M(μ) on the left-hand side that replaces (m 1 + 1)(m 2 + 1), and hence the two theorems taken together thus determine the value M(μ) = (m 1 + 1)(m 2 + 1) as claimed in Proposition 9.2.

Area Siegel-Veech constants

The goal of this section is to show that area Siegel-Veech constants are ratios of intersection numbers, i.e. to prove Theorem 1.4. For this purpose we introduce

d i (μ) = P M g,n (μ) β i • δ 0 = 1 m i + 1 P M g,n (μ) ξ 2g-2 • δ 0 • j =i ψ j , ( 99
)
and then Theorem 1.4 can be reformulated as

c area (μ) = -1 4π 2 d i (μ) a i (μ) . ( 100
)
The proof proceeds similarly to the proof of Theorem 1.1 by showing a recursive formula for both the intersection numbers and the area Siegel-Veech constants. The difference in the formulas is that one vertex at level zero of the backbone graphs is distinguished by carrying the Siegel-Veech weight. We remark that in this section area Siegel-Veech constants for disconnected strata are volume-weighted averages of the constants for the individual components.

The intersection number recursion leads to the remarkable formula

(m 1 + 1)(m 2 + 1)c area (μ) = C saddle connection configuration c hom 1↔2 (C)c area (C), (101)
where c area (C) is defined to be the sum of the area Siegel-Veech constants of the splitting pieces induced by the saddle connection configuration. The other recursion (via q-brackets as in Sect. 4) leads to a very efficient way to compute area Siegel-Veech constants, given in Theorem 10.6.

A recursion for the d i (μ) via intersection theory

We have seen that the values of a i (μ) do not depend on the index i. Similarly for d i (μ) it suffices to focus on the case i = 1. To state the recursion, we introduce the generating series

(t) = g≥1 (2g -1) 2 d 1 (2g -2)t 2g = 1
(m 1 + 1)(m 2 + 1)d 1 (μ) = k≥1 g,μ h P 1 ((m 1 , m 2 ), p) (k -1)! d 1 ( p 1 -1, μ 1 ) a 1 ( p 1 -1, μ 1 ) k i=1 (2g i -1 + n(μ i )) p i a 1 ( p i -1, μ i ) ( 103
)
for n = n(μ) ≥ 2, with the usual summation conventions as in Theorems 1.2 and 3.1.

The first identity (102) was proved in [START_REF] Sauvaget | Volumes and Siegel-Veech constants of H(2g-2) and Hodge integrals[END_REF]. The proof of the second identity (103) will be completed by the end of this subsection. This identity together with the conversions in Sect. 8 implies (101) immediately. We start the proof with the following analog of Proposition 3.10. ( ,,p) is a backbone graph in BB(g, n) 1,2 , then

Proposition 10.2 If

P Mg,n α , ,p • ξ 2g-2 • δ 0 • n i=3 ψ i = m(p) • h P1 (μ -1 , (p v ) v∈V ( ), (v)=0 ) • v∈V ( ), (v)=0 d 1 ( p v -1, μ v ) v ∈V ( )\{v}, (v )=0 a 1 ( p v -1, μ v )
with the conventions for p v as in Proposition 3.10.

Proof

We have the equality that

ζ * (δ 0 ) = v∈V ( ), (v)=0 δ v 0 , where δ v 0 = δ 0 ⊗ 1 ∈ H * (M g v ,n v , Q) ⊗ v =v H * (M g v ,n v , Q) H * (M , Q). Combining with the fact that δ 0 λ g = 0, it implies that ζ * (δ 0 λ g-1 ) = v∈V ( ), (v)=0 ⎛ ⎝ δ v 0 λ v,g v -1 v =v, (v )=0 λ v ,g v ⎞ ⎠ .
123

Therefore, we obtain that

λ g-1 • δ 0 • α 0 , ,p = v∈V ( ), (V )=0 ζ * [M -1 ] ⊗ λ v,g v -1 • δ 0 • [P M g v ,n v ( p v -1, μ v )] 0 v =v, (v )=0 λ v ,g v [P M g v ,n v ( p v -1, μ v )] 0 .
Using the last formula in Lemma 3.8, the rest of the proof then follows from the same argument as in Proposition 3.10.

We also need the following analog of Lemma 3.12.

Lemma 10.3 The values of d 1 (μ) satisfy the recursion

(m 1 + 1)(m 2 + 1)d 1 (μ) = k≥1 g,μ h P 1 ((m 1 , m 2 , μ 0 ), p) • d 1 ( p 1 -1, μ 1 ) (k -1)! • k i=2 p 2 i a( p i -1, μ i ) ( 104
)
with summation conventions as in Lemma 3.12.

Proof We use the formula in Proposition 3.11 for i = 2, multiply by ξ 2g-2 δ 0 n i=3 ψ i and apply p * . The left-hand side then evaluates (by Lemma 3.8 and the fact that δ 0 λ g = 0) to the left-hand side of (104). The right-hand side evaluates (by Proposition 3.9) to the weighted sum over all ( , , p) ∈ BB(g, n) 1,2 of the expression in Proposition 10.2. To prove the lemma we interpret as usual a backbone graph as a decomposition of g and the marked points. The factor (k -1)! (instead of k! in Lemma 3.12) comes from the fact that one of the top level vertices of the backbone graph is distinguished.

With the same notation as in Sect. 3.5 we now define the d-contribution of a rooted tree to be

d( , , p) = m 0 (p) 2 h( 0 , 0 , p 0 ) v∈V ( ), (v)=0 d 1 ( p v -1, μ v ) v ∈V ( )\{v}, (v )=0 a( p v -1, μ v ) .
123

Then we can rewrite Lemma 10.3 as

(m 1 + 1)(m 2 + 1)d 1 (μ) = ( , ,p)∈RT(g,μ) 1,2 d( , , p) |Aut( , , p)| ,
which is the analog of Lemma 3.13.

End of the proof of Proposition 10.1 Now the proof of the proposition can be completed similarly to the end of the proof of Theorem 3.1 at the end of Sect. 3.

As a consequence, d i (μ) does not depend on i and we simply write d(μ) from now on.

A recursion for area Siegel-Veech numerators via weighted counting of covers

We recall the main steps of [START_REF] Chen | Quasimodularity and large genus limits of Siegel-Veech constants[END_REF] that reduce the computation of area Siegel-Veech constants to a statement about cumulants. An application of the Siegel-Veech formula [START_REF] Chen | Quasimodularity and large genus limits of Siegel-Veech constants[END_REF]Proposition 17.1] gives the quantity we want to compute by

c area ( M g (m 1 , . . . , m n )) = lim D→∞ 3 π 2 D d=1 c • -1 (d, ) D d=1 N • d ( ) , ( 105
)
where N

c area ( M g (m 1 , . . . , m n )) = 3 π 2 T -1 | f m 1 +1 | • • • | f m n +1 L f m 1 +1 | • • • | f m n +1 L = 3 π 2 T -1 |h m 1 +1 | • • • |h m n +1 L h m 1 +1 | • • • |h m n +1 L , ( 106
)
where T -1 is a hook-length moment function on partitions (but not an element of the ring * , see [START_REF] Chen | Quasimodularity and large genus limits of Siegel-Veech constants[END_REF]Section 13]). Here we use Proposition 4.2 and (40) in the second equality for the denominator. Before passing to the leading term, the numerator

T -1 | f m 1 +1 | • • • | f m n +1
is a linear combination of differences of the form T -1 f q -T -1 q f q . By [12, Theorem 6.1] these differences can be written as explicit linear combinations of derivatives of Eisenstein series and homogeneous differential operators. We can perform the same rewriting with T -1 |h m 1 +1 | • • • |h m n +1 . Again thanks to [START_REF] Mumford | Towards an enumerative geometry of the moduli space of curves[END_REF] the difference is a cumulant with each entry being of strictly smaller weight. Consequently, it does not contribute to the leading order term. We have seen in Sect. 4 how to compute the denominator and related it in Sect. 5 to a(μ). Now we take care of the numerator. Recall that we defined H (u) q in (51) in Sect. 4. Define now

C -1 (u) q = n>0 T -1 • ≥1 h n q u n n! .
By definition of cumulants (or by [START_REF] Chen | Quasimodularity and large genus limits of Siegel-Veech constants[END_REF]Proposition 6.2]) we are interested in the leading term of the quotient

C • -1 (u) q := C -1 (u) q H (u) q = n≥0 T -1 |h 1 | • • • |h 1 n 1 | h 2 | • • • |h 2 n 2 | • • • q u n n! .
To evaluate the numerator of this fraction, recall from [START_REF] Chen | Quasimodularity and large genus limits of Siegel-Veech constants[END_REF]Section 16] the definition of the modified q-bracket f q = T -1 f q -T -1 q f q -1 24

∂ 2 ( f ) q , ( 107
)
where ∂ 2 is the differential operator

∂ 2 : ∂ ∂ p 1 + ≥2 ( -1) p -2 ∂ ∂ p .
This bracket is useful, since its effect can be computed by differential operators acting (contrary to T -1 ) within the Bloch-Okounkov algebra. In fact, [START_REF] Chen | Quasimodularity and large genus limits of Siegel-Veech constants[END_REF]Theorem 16.1] states that

f q = j≥1 G ( j-1) 2 ρ 0, j ( f ) q + i≥2, j≥0 G ( j) i ρ i, j ( f ) q ,
where ρ i, j are differential operators of degree j that shift the weight by -i -2 j, whose definition we recall in (110) below. Motivated by the action of these operators we define

C -1 (u) = n>0 j≥1 j! ρ 0, j ≥1 h n u n n! , ( 108
)
and we let H (u) = exp( ≥1 h u ) such that H (u) q = H (u) q .

Lemma 10. [START_REF] Aggarwal | Large genus asymptotics for Siegel-Veech constants[END_REF] The leading terms of C • -1 (u) q and C • -1 (u) q := -1 24

C -1 (u) q H (u) q agree.

Proof First note that Ev(G ( j-1) 2

)(X ) = -1 24 ( j!X + ( j -1)!) by the defining formulas in [START_REF] Chen | Quasimodularity and large genus limits of Siegel-Veech constants[END_REF]Section 9]. This is the reason for the factor j! in (108). From the non-vanishing of the area Siegel-Veech constant, we know that the leading degree contribution is as in [START_REF] Okounkov | Gromov-Witten theory, Hurwitz numbers, and matrix models[END_REF]. Lower weight terms before passing to the cumulant quotient will contribute to lower order in the growth polynomial. Since ∂ 2 is of degree -2, its contribution in ( 107) is negligible and we can work with the star-brackets. For the same reason, the terms with i > 0 in the definition f q are dominated by the corresponding term with i = 0 and can be neglected.

Our goal is to compute the h-evalutaion of C • -1 (u) and its leading term using Proposition 4.3. Using these identities we can evaluate both sides of (109) to be of the form

e -z 2 N /2 α∈P(n) A∈α W (z A )R({z a } a∈A )
where R({z a }) are polynomials that are visibly different on the two sides, but in fact agree by using the identity

z n+1 A = ∅ =J ⊂A |J |! z J ν∈J z ν z n-|J | A .
To verify this expression, let e i = (-1) i [x n-i ] a∈A (x -z a ) be the elementary symmetric functions in the z a . Then the contribution with |J | = j to the righthand side is e n-j 1 (e 1 e j -( j + 1)e j+1 ). This means that the right-hand side is a telescoping sum where only the first term remains after summation.

The preceding Lemma 10.5, Proposition 4.3 for the computation of the h-brackets, Lemma 10.4 and (106) now imply immediately our goal: Theorem 10. [START_REF] Bainbridge | Compactification of strata of abelian differentials[END_REF] The area Siegel-Veech constants can be computed as c area (m 1 , . . . ,

m n ) = -1 8π 2 [z m 1 +1 1 • • • z m n +1 n ] ∂ 2 (H n ) [z m 1 +1 1 • • • z m n +1 n ] H n h →α ,
where H n (z 1 , . . . , z n ) is recursively defined as in Sect. 4.

Proof of Theorem 1.4

We start with an explicit formula for the ∂ 2 -derivative used in Theorem 10.6 in the case of the minimal strata. for all 1 ≤ i ≤ n, where • ∈ {even, odd}. Using Assumption 6.1 to deal with the case of the minimal strata, the validity of ( 112) is equivalent to the validity of

c area (μ) odd = -1 16π 2 [z m 1 +1 1 • • • z m n +1 n ] (2πi) 2g (∂ 2 (H n ) -∂ 2 (H n )) (2g -2 + n)!v(μ) odd h →α h →α
where H n and H n are recursively defined as in Sects. 4 and 6, and where

∂ 2 = ∂ ∂p 1 + ∞ ≥1 2 (2 + 1)p 2 -1 ∂ ∂p 2 +1
is the analog of the differential operator ∂ 2 on the algebra of super-symmetric functions. There is a clear strategy towards this goal:

• Show that there are operators like the T p for p ≥ -1 odd as in [12, Section 12], whose strict brackets compute Siegel-Veech weighted and spin-weighted Hurwitz numbers. • Show that the action of T p inside strict brackets can be encoded by differential operators like the ρ i j in [START_REF] Chen | Quasimodularity and large genus limits of Siegel-Veech constants[END_REF]Section 14]. • Show that these operators satisfy a commutation relation as in Lemma 10.5, with ∂ replaced by ∂ 2 . • Conclude by comparing the recursions as in the preceding proofs.

Given the length of this paper, we do not attempt to provide details here.

Large genus asymptotics

In this section we study the large genus asymptotics of Masur-Veech volumes and area Siegel-Veech constants and prove the conjectures of Eskin and Zorich in [START_REF] Eskin | Volumes of strata of Abelian differentials and Siegel-Veech constants in large genera[END_REF] by using our previous results.

Volume asymptotics

We recall from [12, Theorems 12.1 and 19.1] and [42, Theorem 1.9] that the asymptotic expansions of v(2g -2) and v(1, . . . , 1) can be computed using the mechanisms of very rapidly divergent series [START_REF] Chen | Quasimodularity and large genus limits of Siegel-Veech constants[END_REF]Appendix] as

v(1 2g-2 ) ∼ 4 1 - π 2 24g - 60π 2 -π 4 1152g 2 -• • • , v(2g -2) ∼ 4 1 - π 2 12g - 24π 2 -π 4 288g 2 -• • • .
123

Let μ = (m 1 , . . . , m n ) be a partition of 2g -2 into n positive integers with n ≥ 2. We write μ = (m 1 + m 2 , m 3 , . . . , m n ) and μ = (m 1 + m 2 -2, m 3 , . . . , m n ). We use Theorem 1.2 and the two obvious backbone graphs, the one with a single top level component of genus g (i.e. k = 1) and the one with two top level components (i.e. k = 2) of genus 1 and g -1 respectively, to deduce the inequality

v(μ) ≥ v(μ ) + π 2 (2g -5 + n)! 6(2g -3 + n)! v(μ ), ( 113
)
where we use h Then the inequality (113) implies that

v(μ) -v(μ ) ≥ - 2π 2 3(2g -3 + n)(2g -4 + n) + π 2 (2g -5 + n)! 6(2g -3 + n)! v(μ ) ≥ - Cπ 2 6g(2g -3 + n)(2g -4 + n) .
In particular for all μ we have v(2g -2) -Cπ 2 (n -1) 6g(2g -1) 2 ≤ v(μ) ≤ v(1, . . . , 1)

+ Cπ 2 (2g -2 -n) 6g(2g -1) 2 .
Since n is bounded by 2g -2, there exists a constant C such that | v(μ)| ≤ C /g 2 for all μ. Thus the first part of Theorem 1.5 holds.

Asymptotics of Siegel-Veech constants

We apply the same strategy to control the asymptotic behavior of area Siegel-Veech constants. We denote by Therefore, using the fact that c area (μ) = d(μ)/v(μ) we thus deduce the second part of Theorem 1.5.

Spin asymptotics

The goal here is to prove that the volumes of the odd and even spin components are asymptotically equal. This is the content of Theorem 1.6 in the introduction that refines the conjecture of Eskin and Zorich [START_REF] Eskin | Volumes of strata of Abelian differentials and Siegel-Veech constants in large genera[END_REF]Conjecture 2]. Recall that we defined in Sect. 6.3 (115)

P Z (u) = exp
Moreover, it has the asymptotics

v(2g -2) ∼ -1 2 g-2 1 + 2π 2 3g + 12π 2 + π 4 18g 2 + • • • (116)
as g → ∞.

For the reader's convenience we give a table for low genus values: For the comparison of area Siegel-Veech constants, since we have shown that the two spin components have the same volume asymptotic, combining with [START_REF] Aggarwal | Large genus asymptotics for Siegel-Veech constants[END_REF]Propositions 4.3 and 4.4] it suffices to compare the odd and even one-cylinder contributions in the part of the principal boundary where the corresponding configurations have either a figure-eight boundary or a pair of holes boundary. Using again the same volume asymptotics of the spin components in genus g -1 and combining with [20, Lemmas 14.2 and 14.4], we see that the odd and even contributions are comparable for those two parts of the principal boundary, thus implying the claim. discussions on related topics. We thank the Mathematisches Forschungsinstitut Oberwolfach (MFO), the Max Planck Institute for Mathematics (MPIM Bonn), the l'Institut Fourier (Grenoble), the American Institute of Mathematics (AIM), the Banff International Research Station (BIRS) and the Yau Mathematical Sciences Center (YMSC) for their hospitality during the preparation of this work. Finally we thank the referee for carefully reading the paper and many helpful suggestions.

Funding Open Access funding enabled and organized by Projekt DEAL.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Theorem 1 . 6 (

 16 [23, Conjecture 2]) The volumes and area Siegel-Veech constants of odd and even spin components are comparable for large values of g.More precisely,v(μ) odd v(μ) even = 1 + O(1/g), c area (μ) odd c area (μ) even = 1 + O(1/g),

Fig. 1 A

 1 Fig.1A rooted tree of genus eleven with three vertices of genus zero (black) and seven legs

  Take a bi-colored graph ( , ) and a partition μ of 2g -2. Now we consider a twist p of type (μ[0] = μ, μ[∞] = ∅) compatible with and construct a subspace P M p , ⊂ P M , such that ζ # , (P M p , ) lies in the boundary of P M g,n (μ). Let

)

  For all x, y, z ∈ R D(x yz) = x D(yz) + y D(xz) + z D(x y)x y D(z)xz D(y)yz D(x) .

  -trivial, and define cont( ) = H v for the trivial graph with a unique vertex v and no edges. Finally, we set S {i} = H {i} = A {i} and for |I | ≥ 2

p,.

  We recall from Sect. 3 that ζ # , (D) is a divisor of P M g,n (μ) if and only if dim(P M p , ) = dim(P M g,n (μ))-1. Hence in this case we define α(D) = ζ # , * (D) ∈ H * (P M g,n ), and define α(D) = 0 otherwise. We will denote by P M p,• , the union of the irreducible components of P M p , that are mapped to P M g,n (μ) • .

  used characters of the modified Sergeev group C(d) = S(d) Cliff(d) to produce elements in * . Here Cliff(d) is generated by involutions ξ 1 , . . . , ξ d and a central involution ε with the relation ξ i ξ j = εξ j ξ i . Irreducible representations of C(d) are V λ indexed by λ ∈ SP. We denote by f μ (λ) the central character of the action of a permutation g μ ∈ S(d) ⊂ C(d) of cycle type μ on V λ by conjugation. The analog of the Burnside formula is [22, Theorem 2] stating that for a fixed profile

Theorem 6 . 7

 67 If we define h = -1 [u +1 ]P(u) where P(u) = exp -s≥1,s odd u s+1 p s , (68)

1 )

 1 and similarly P λ = P λ ( • ; -1)

4 . 2 . 2 n 2 n

 4222 We have deg Ev p 1 p 2 • • • p n str = 1 i=1 i , the highest term being contributed by the partition into singletons. From (75) we deduce that deg Ev p 1 |p 2 | • • • |p n str = 1 i=1 i -(n -1), and thus obtain the expected degree drop. The claim now follows from the usual approximation of Masur-Veech volumes by counting torus covers [21] and [12, Proposition 19.1].

  T ) = |{γ ⊂ X a saddle connection joining z 1 and z 2 , γ ω ≤ T }| (84) 123 be the counting function. The upper index emphasizes that we count all physically distinct saddle connections. It should be distinguished from the version A hom 1↔2 (T ) = |{γ ⊂ X a homology class of s.c. joining z 1 and z 2 , γ ω ≤ T }|,

  M-p 1↔2 (T ) and the Siegel-Veech constant c M-p 1↔2 (M) in analogy to (85) and (86), counting once every M-parallel class of cylinders. [20, Proposition 3.1] can now be restated as c M-p

Fig. 3

 3 Fig. 3 Configurations for Hurwitz spaces in M 2 (1, 1)

Lemma 10 . 5 W 2 N / 2 α∈P

 10522 The commutation relation∂ 2 • e D ( f ) = e D j≥1 j! ρ 0, j ( f ) (109)holds for every f ∈ * .Proof We will check the relation on the n-point function for every n. Since we will recall formulas from[START_REF] Chen | Quasimodularity and large genus limits of Siegel-Veech constants[END_REF] we use the rescalingsQ k = p k-1 /(k -1)! of the generators of * , where Q 0 = 1 and Q 1 = 0.The following identities even hold on the polynomial ringR = Q[Q 0 , Q 1 , Q 2 , . . .] mapping to * . We set W (z) = k≥0 Q k z k-1 .We recall from [12, Theorem 14.2] the action of the operators ρ 0, j , namelyρ 0, j W (z 1 ) • • • W (z n ) = (z ν ) (110)where z J = j∈J z j and N = {1, . . . , n}. On the other hand, in terms of the Q i , the operator D defined in Sect. 4 is just D = 1 2 ( -∂ 2 ), where ∂ is the differential operator sending Q i to Q i-1 . From [12, Proposition 10.5] we know thate D W z 1 ) • • • W z n = e -z

Proposition 10 . 7 8π 2 [u 2g- 1 ]

 10721 The area Siegel-Veech constants for the minimal strata arec area ( M g (2g -2)) = -1 D(u) [u 2g-1 ] A(u) , (111)whereD(u) = (A (u) + uA (u))/uA (u) 2 = t -1 18 t 3 + 91 2304 t 5 -• • • .

P 1 (

 1 (m 1 , m 2 ), (m 1 +m 2 +2)) = h P 1 ((m 1 , m 2 ), (m 1 +m 2 , 1)) = 1 for m 1 , m 2 > 0 and v(0) = π 2 /6. In particular this inequality implies that v(μ) ≥ v(μ ) and thus v(2g -2) ≤ v(μ) ≤ v(1, . . . , 1)for all μ. Consequently, there exists a constant C > 0 such that for all μ we have the inequality |v(μ) -4| < C/g. Now we introduce the notation v(μ) := v(μ) -4 + 2π 2 3(2g -3 + n) .

(m i + 1 )

 1 d(μ) := c area (μ)v(μ) • |d(μ)| where d(μ) is defined in (99) and where the second equality stems from (100). For μ = (m 1 , . . . , m n ) with n ≥ 2, we write μ = (m 1 + m 2 , m 3 , . . . , m n ) and μ = (m 1 + m 2 -2, m 3 , . . . , m n ) as before. Then from Proposition 10.1 we have d(0) = 1 and we obtain the inequalityd(μ) ≥ d(μ ) + π 2 (2g -5 + n)! 6(2g -3 + n)! d(μ ) + (2g -5 + n)! 2(2g -3 + n)! v(μ ) .In particular this inequality implies that d(μ) ≥ d(μ ). Moreover, we know the asymptotic expansionsd(2g -2) ∼ 2 -3 + π 2 6g + • • • and d(1, . . . , 1) ∼ 2 -3 + π 2 12g + • • • from [12,Theorem 19.4] and [42, Theorem 1.9]. Consequently, there exists a constant C such that | d(μ) -2| < C/g for all μ. Then by the same argument as above we can show that there exists a constant C such that d(μ) -2 + 3 + π 2 3(2g -3 + n) < C /g 2 .

j≥1 1 2 jProposition 11 . 1

 2111 /2 ζ(j)2 u j+1 . (114)123 The difference v(2g -2) = v(2g -2) even -v(2g -2) odd can be computed as the coefficient extraction v(2g -2) = 2(2πi) 2g (2g -1)! u 2g-1 1 (u/P Z ) -1 .

1 2 ( 2 - 2 -

 1222 The first statement is a reformulation of a special case of Corollary 6.11.The power series P Z (u) is a very rapidly divergent series, just as P B (u) is, since the coefficients !b and ζ(-)/2 = !b • √ ) differ by a factor that grows only geometrically. The asymptotic statement thus follows from the method of very rapidly divergent series.Proof of Theorem 1.6 Proposition 11.1 together with Theorem 1.5 implies that there exists a constant C such that for all g ≥ 1|v(2g -2) odd -v(2g -2) even | ≤ C /g .Repeated application of Theorem 1.5 implies that v(μ) • ≥ v(2g -2)• for all μ with |μ| = 2g -2. Theorem 1.5 moreover implies that there exists a constant C such that for all μ with |μ| = 2g -2 the inequality |v(μ) -4| ≤ C /g holds. Thus for all μ with |μ| = 2g -2 we havev(2g -2) odd ≤ v(μ) odd = v(μ) -v(μ) even ≤ v(μ) -v(2g -2) even ≤ 2 + (C + 3C )/g .It follows that |v(μ) odd -2| ≤ (C + 3C )/g and the same holds for v(μ) even . This implies the claim for the volume comparison.

  • • • μ k , and p = ( p 1 , . . . , p k ) is defined by p i = 2g i -1 -|μ i | and required to satisfy p i > 0. Here the Hurwitz number h P 1 ((m 1 , m 2 ), p) is defined for any p by h P 1 ((m 1 , m 2 ), p) = (k -1)![t m 1 +1

	1	t m 2 +1 2	]	k i=1	t 1 t 2	(t 2 -t p i 1 ) p i t 2 -t 1	. (6)

  To evaluate the derivative of H I , there are two cases to consider, depending on the cardinality of I . If I = {i}, then

									∂H I ∂h p = z	p i for all p > 0. Otherwise,
	we use the induction hypothesis to compute that
	∂H I ∂h p	=	∂S I ∂h p	=	∈OT(I )	f # ( )	∂ cont( ) ∂h p		⎛
				=	∈OT(I )	f # ( )	p:E→Z >0 e∈E	z	p(e) e	v∈V	∂H v ∂h p	⎝ ⎛	v =v	H v⎞ ⎠
				=	( ,v)∈OT(I ) v	f # ( )	p:E→Z >0 e∈E	z	p(e) e	∂H v ∂h p	⎝	v =v	H v⎞ ⎠ .
									p e e z	p e e ]H {e ,e }	∂H I ∂h p e	∂H I \I ∂h p e	.

1) I ⊂I p e , p e >0 [z Now we assume that both I and I \I have at least two elements. Take two oriented trees (

  • H {e ,e } for the top vertex v, •

	∂H v ∂h p e	and	∂H v ∂h p e

  • ( ) is the number (weighted by |Aut( p)|) of transitive Hurwitz tuples (α, β, (γ i ) If the pairs (γ 1 , γ 2 ) appearing in the Hurwitz tuples of profile equidistribute among pairs of

	γ i of type μ (i) .	n i=1 ) ∈ S n+2 d	with [α, β] = n i=1 γ i and
	Proposition 9.4

  is the sum over those covers with -1st Siegel-Veech weight (see[START_REF] Chen | Quasimodularity and large genus limits of Siegel-Veech constants[END_REF] Section 3]). The relation of the sum of Fourier coefficients to the growth polynomial[START_REF] Chen | Quasimodularity and large genus limits of Siegel-Veech constants[END_REF] Proposition 9.4] and a rewriting of the Siegel-Veech weighted counting[START_REF] Chen | Quasimodularity and large genus limits of Siegel-Veech constants[END_REF] Corollary 13.2] translate this into

• d ( ) is the number of connected torus covers of degree d with ramification profile = (m 1 + 1, . . . , m n + 1) and where c • -1 (d, )

This assumption was later verified in[START_REF] Costantini | The area is a good metric[END_REF].

In[START_REF] Bainbridge | Compactification of strata of abelian differentials[END_REF] the notation P M inc g,n (μ) is used. Here we drop the superscript "inc" for simplicity. In[START_REF] Sauvaget | Cohomology classes of strata of differentials[END_REF] this space is denoted by PH g,n (μ).

This would be (2πi) k ev[F](-4π 2 ) in the notation of[START_REF] Chen | Quasimodularity and large genus limits of Siegel-Veech constants[END_REF] Equation (85)].

Assumption 6.1 was recently verified in[START_REF] Costantini | The area is a good metric[END_REF], to which we refer the reader for more details about the area form being a good metric. See also[START_REF] Chen | Masur-Veech volumes and intersection theory: the principal strata of quadratic differentials[END_REF] for a similar application to the volumes of strata of quadratic differentials with odd zeros.

Acknowledgements

We are grateful to Anton Zorich for valuable suggestions on Siegel-Veech constants and numerical cross-checks. Moreover, we would like to thank Alex Eskin and Andrei Okounkov for their help about spin components. We also thank Amol Aggarwal, Ben Bakker, Qile Chen, Yitwah Cheung, Eleny Ionel, Felix Janda, Xin Jin, Edward Looijenga, Duc-Manh Nguyen, Aaron Pixton, Song Sun, Alex Wright and Dimitri Zvonkine for inspiring

D. Research of the first author is partially supported by the NSF CAREER Grant DMS-1350396 and a von Neumann Fellowship at IAS in Spring 2019. M. Möller: Research of the second author is partially supported by the DFG-project MO 1884/1-1 and by the LOEWE-Schwerpunkt "Uniformisierte Strukturen in Arithmetik und Geometrie". A. Sauvaget: Research of the third author is partially supported by the Dutch Research Council (NWO) Grant 613.001.651.

, is contained in the boundary of exactly one connected component of the stratum (see e.g. [START_REF] Chen | Principal boundary of moduli spaces of abelian and quadratic differentials[END_REF]Corollary 4.4]). Thus the neighborhood of a generic point of D in P M g,n (μ) • is given by [43, Lemma 5.6 and the subsequent formula]. In particular the multiplicity of D in the vanishing locus of s is the same as that of the entire boundary stratum, which implies that m

By the same arguments as in Sect. 3.3 (see Proposition 3.9) one can show that

Proposition 6.5 If n ≥ 2, then the values of a i (μ) • are the same for all 1 ≤ i ≤ n, denoted by a(μ) • , and can be computed as

Proof This follows from the same argument as in the proof of Lemma 3.12.

123

Proof Differentiating (50) gives

Combining these two equalities gives ∂ 2 (H 1 (u)) = (H (u)+uH (u))/uH (u) 2 and the claim by substituting h → α .

Proof of Theorem 1. [START_REF] Aggarwal | Large genus asymptotics for Siegel-Veech constants[END_REF] We start with the case of a single zero. Comparing ( 111) and ( 102) we need to show that D(u) = 2 (u)/u, i.e. in view of ( 16) we need to show that

This equality can be implied by showing that

which in turn follows since the derivative

u has no (-1)-term. Finally, to deal with the case of multiple zeros, we recall from (60) that

H n and hence we need to show that

after knowing that this is true for the case of n(μ) = 1. This follows immediately from differentiating (57), since after substituting h → α this is exactly the sum of the recursions (103) (known to hold for the d(μ)), averaging over all pairs (m r , m s ) of the entries of μ, as in (59).

Given Theorem 1.4 for the area Siegel-Veech computation of the strata on one hand, and the refined Theorem 6.3 for the volume computation of the spin components on the other hand, it is natural to suspect that area Siegel-Veech constants for the spin components can also be computed as ratios of intersection numbers c area (μ) • = -1 4π 2

Publisher's Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.