Index-Like Theorems from Line Defect Vevs

In this paper we investigate the relation between complexified Fenchel-Nielsen coordinates and spectral network coordinates on Seiberg-Witten moduli space. The main technique is the comparison of exact expressions for the expectation value of 't Hooft defects in certain 4D $SU(2)$ $\mathcal{N}=2$ gauge theories. We derive an index-like theorem for a class of Dirac operators on singular monopole moduli spaces. Our expression determines the indices of Dirac operators on singular monopole moduli spaces in terms of characteristic numbers for vector bundles over certain Kronheimer-Nakajima quiver varieties.


Introduction and Technical Summary
Exact results in quantum field theories with four-dimensional N = 2 supersymmetry can lead to nontrivial and interesting mathematical predictions. In this paper we consider exact results for line defect expectation values in theories which are of class S but also possess a Lagrangian formulation.
In general, class S theories are closely related to Hitchin systems and in these theories the Hitchin moduli space M is interpreted as a space of vacua when the 4D theory is compactified on a circle. If a line defect L preserves four supersymmetries and is wrapped on this circle then the associated vev ⟨L⟩ is a holomorphic function on M in a complex structure determined by the supersymmetry preserved by L. The preserved supersymmetry may be characterized by a phase ζ, which may be viewed as an element of the twistor sphere: ζ also determines a complex structure on M. We will denote the space M with complex structure determined by ζ as M ζ . In the theories under consideration the holomorphic function ⟨L(ζ)⟩ on M ζ can be computed, exactly, in three different ways. This paper explores some consequences of comparing the resulting exact expressions.
For theories of class S one of the exact methods expresses ⟨L(ζ)⟩ in terms of "spectral network coordinates" on M ζ [17][18][19]. These coordinates are generalizations of well-known cluster, shear, and Fock-Goncharov coordinates. They are functions on the twistor space and, restricted to a fiber M ζ , are holmorphic Darboux coordinates in complex structure ζ. See Section 4 below for a description of "spectral network coordinates." We refer to the exact result for ⟨L(ζ)⟩ in these coordinates as the "Darboux expansion." See equation (1.3) below for the detailed formula.
On the other hand, thanks to the AGT correspondence, ⟨L(ζ)⟩ can also be expressed in terms of complexified Fenchel-Nielsen coordinates [1,10]. These are also holomorphic Darboux coordinates on M ζ . See Section 4.1.2 for a description of Fenchel-Nielsen coordinates. Thus, by comparing answers for a finite set of line defects one can express one coordinate system in terms of the other. The transformation turns out to be rather nontrivial.
In theories of class S that also have a Lagrangian description one can also evaluate ⟨L(ζ)⟩ exactly using localization techniques. These localization techniques are only understood fully for certain line defects known as pure 't Hooft operators, and are only valid in weak-coupling regimes in the Coulomb branch. 1 The function ⟨L(ζ)⟩ for an 't Hooft line defect, when evaluated using localization, turns out, again, to be expressed naturally in terms of complexified Fenchel-Nielsen coordinates and moreover involves interesting equivariant characteristic numbers of Kronheimer-Nakajima varieties.
At the same time, in the weak coupling region of a class S theory the coefficients of the Darboux expansion have an interpretation in terms of indices of certain Dirac-like operators defined on moduli spaces of (singular) magnetic monopoles. Thus, by first comparing line defect vev's to determine the change of variables from spectral network to Fenchel-Nielsen coordinates and then using the above facts, we find unusual expressions for the L 2 -index of certain Dirac operators in terms of characteristic numbers of Kronheimer-Nakajima varieties. This comparison yields our "index-like" theorem.
We should stress how unusual our "index-like" theorem is. In standard index theorems the index of a Dirac operator on one manifold is expressed in terms of a characteristic number of the relevant bundle on the same manifold. In the present case the (L 2 -) index of a Dirac operator is expressed in terms of sums of characteristic numbers on a different (but related) manifold. (Alternatively, an equivariant characteristic number is expressed as a sum of indices of Dirac operators.) Hence the qualifier "index-like theorem." The physical idea here is essentially trivial, but the mathematical statement seems to be fairly nontrivial and would seem to be challenging to prove from first principles.
The remainder of this paper implements the above idea in the special case of 4D G = SU (2) N = 2 supersymmetric gauge theories with adjoint or fundamental matter. In the remainder of this introduction we explain the idea in a bit more detail.

Technical Summary
In the case of 4D G = SU (2) N = 2 supersymmetric gauge theories, there are two complexified-Fenchel-Nielsen coordinates which we will denote a, b. From general principles, the expectation value of the 't Hooft defects can be expressed in Fenchel-Nielsen coordinates as a Fourier expansion in b. More precisely, this can be written as [22,25] ⟨L p,0 (ζ)⟩ = where v = diag(v, −v) and P = diag(p, −p) where v, p ∈ Z. On the right hand side the ζdependence is captured by the use of complexified Fenchel-Nielsen coordinates on M ζ . 2 Here the expectation value above is expressed as a sum over monopole bubbling configurations where cosh(v, b)F (a) v encodes the contribution of bulk fields and Z mono (a, m, ; P, v) describes the contribution from the SQM that arises on the 't Hooft defect from bubbling [4]. See [4][5][6] for more background and explanation of notation.
In the localization computation of ⟨L p,0 ⟩, Z mono (P, v) is given by a characteristic number of a certain resolved Kronheimer Nakajima space 3 Z mono (a, m, ; P, v) = lim Here M ξ KN (P, v) is a certain resolved Kronheimer-Nakajima space determined by the line defect charge (P ) and core magnetic charge (v), e ω+µ T induces the T -equivariant volume form on M ξ KN (P, v),Â T (T M ξ KN ) is the T -equivariantÂ-genus that describes the contribution from the N = 2 vectormultiplet and C T ×T F (V(R)) is a characteristic class related to the matter hypermultiplets where T is the Cartan torus of the conserved global symmetry group of flavor, R-, and global gauge transformations. The equivariant integral can then be evaluated 2 Usually complexified Fenchel-Nielsen coordinates are introduced as holomorphic coordinates, depending on a cutting system, of the character variety X = Hom(π1(C), G C ) conj. for some complex gauge group G C . In our case, M ζ is isomorphic to X for all ζ ≠ 0, ∞ as a complex manifold, but not canonically. Our Fenchel-Nielsen coordinates will therefore also be functions on the twistor space of M (with the fibers above ζ = 0, ∞ removed). When restricted to a fiber M ζ , they are holomorphic Darboux coordinates. It is in this way that they become comparable to spectral network coordinates. 3 There is an additional subtlety with 4D N = 2 SU (N ) theories with N f = 2N . See footnote 10 for more details.
as a contour integral in an algebraic torus whose poles are enumerated by Young tableaux [37,38,42]. 4 On the other hand, using the class S technology, the expectation value of a supersymmetric line defect can be computed by the trace of the holonomy of a flat SL(N ; C) connection along a corresponding curve γ in an associated Riemann surface C. Spectral networks express the expectation value of such 4D line defects as where Ω(γ, L p,0 ; u) are framed BPS indices, Y γ are Darboux functions on the moduli space of flat SL(N ; C) connections on C associated to the physical charge γ ∈ Γ , and Γ is a torsor of the IR charge lattice [17].
In the semiclassical limit of the theories we are considering, the framed BPS indices of 't Hooft defects can be identified with the index of a twisted Dirac operator on singular monopole moduli space [3,20,21,35,36,39,40,45]. Locally on moduli space we can decompose γ non-canonically into magnetic, electric, and flavor charge γ = γ m ⊕ γ e ⊕ γ f . Then we have: Here the superscript γ e ⊕ γ f denotes the associated eigenspace of the L 2 index of D Y , a Dirac operator modified by adding Clifford multiplication by a hyperholomophic vector field defined by Y ∈ t C . The Dirac operator acts on sections of E matter ⊗ SM(P, γ m , u) where SM(P, γ m , u) is the spinor bundle on the singular monopole moduli space M(P, γ m , u) and E matter → M(P, γ m , u) is a vector bundle over it related to hypermultiplet zero-modes. 5 Thus, by comparing the expectation value of 't Hooft defects computed via localization and spectral network techniques in a weak coupling limit, we can derive a relation between characteristic numbers of Kronheimer-Nakajima spaces and indices of Dirac operators on singular monopole moduli space: Since the formula is valid for an infinite number of line defects, we can use it both to express Y γ in terms of a, b (or vice versa) and to determine relations between Dirac indices and characteristic numbers on certain Kronheimer-Nakajima spaces. The outline of the paper will be as follows. We will begin by reviewing 't Hooft defects in 4D N = 2 G = SU (2) asymptotically free gauge theories with adjoint and fundamental matter . Then we will discuss the localization results for their expectation value and the connection to characteristic numbers of Kronheimer-Nakajima spaces. Then we will go on to discuss spectral networks and the relation to the index of a twisted Dirac operator on singular monopole moduli space in the semiclassical limit. Then we will discuss the comparison of localization and spectral networks and provide an index formula for Dirac operators on singular monopole moduli space and a formula for the characteristic numbers of Kronheimer-Nakajima spaces. We will illustrate the derivation of these formulas fully in the example of the SU (2) N f = 0 theory.
We expect that this technique can be extended to more general 4D N = 2 theories of class S. In the more general setting, the calculations analogous to the ones presented in this paper could be used to compute more general indices of the twisted Dirac operator on singular monopole moduli space with higher rank gauge group and coupled to more complicated vector bundles with hyperholomorphic connection.

't Hooft Defects in 4D
Our setting will be 4D G = SU (2) N = 2 gauge theory on R 3 × S 1 with adjoint or fundamental matter. We will be considering the expectation value of reducible 't Hooft defects which are studied in [4,6,25].
A reducible 't Hooft defect is defined in terms of irreducible defects which are 1 2 -BPS operators specified by the data (⃗ x n , P n , ζ). Here ⃗ x n ∈ R 3 specifies the insertion position, P n ∈ Λ mw specifies the 't Hooft charge, 6 and ζ ∈ U (1) specifies the conserved supersymmetries. Given this data, the associated 't Hooft defect is defined by imposing the local boundary conditions is the magnetic field, and X is a real, adjoint valued Higgs field which is related to the complex Higgs field of the N = 2 vectormultiplet Φ as We will denote an irreducible 't Hooft defect as L [Pn,0] (⃗ x n ). We will often suppress the dependence on ⃗ x n in addition to the phase ζ. A reducible 't Hooft defect is then defined as the product of irreducible 't Hooft defects: where h 1 is a simple magnetic weight. As an operator in a quantum theory, the expectation value of a reducible 't Hooft defect inserted at the origin is defined as the limit Note that L p,0 should not to confused with L [ph 1 ,0] . 6 The magnetic weight lattice is defined as Λmw = {h ∈ Λ cochar ⟨µmatter, h⟩ ∈ Z , ∀µmatter} where µmatter are highest weights specifying the representations of the hypermultiplets in the theory and the cocharacter lattice is defined as Λ cochar = {h ∈ t exp{2πh} = 1 G }.

Monopole Bubbling
In order to compute the expectation value of an 't Hooft defect, it is necessary to understand the BPS field configurations. These are generically described by a collection of dyonic states in the presence of an 't Hooft defect. An important feature of the interaction between magnetically charged states and 't Hooft defects is the phenomenon called monopole bubbling [26].
Monopole bubbling is the process in which a smooth monopole is absorbed by an 't Hooft defect, thus screening the magnetic charge of the defect. Additionally, it gives rise to an effective SQM living on the world volume of the screened 't Hooft defect which we will refer to as the bubbling SQM. This can be understood via the following string theory construction [4,5].
Consider embedding the 4D N = 2 SU (2) gauge theory on the world volume theory of a pair of parallel D3-branes which are localized at x 5,6,7,8,9 = 0 and x 4 = ± 1 2 v. This embedding can be achieved by adding a sufficiently large mass deformation to the adjoint hypermultiplet associated with the excitations along the x 6,7,8,9 -directions. 7 This seperation in the x 4 -direction introduces a Higgs vev X ∞ = vH 1 , where H 1 is a simple coroot.
Here, smooth monopoles are described by D1-branes running between the D3-branes in the x 4 -direction [8]. Additionally, a reducible 't Hooft defect at ⃗ x n is described by a collection of NS5-branes that are localized at x 1,2,3 = ⃗ x n and spread out in the x 4 -direction between the D3-branes. If we combine these two elements so that there are m D1-branes running between the D3-branes and 2p NS5-branes localized at ⃗ x n , then we have an asymptotic, relative magnetic charge and 't Hooft chargẽ whereĥ 1 is a simple cocharacter (not magnetic weight), and [P n ] − is the image of the 't Hooft charge in the totally negative Weyl chamber. See Figure 1.
Here, monopole bubbling can be understood by moving D1-branes to be coincident with the collection of NS5-branes in the x 1,2,3 -directions. We can then understand the bubbling SQM living on the world volume of the 't Hooft defect by performing a sequence of Hanany-Witten transformations so that the bubbled D1-branes only end on NS5-branes [4,14,23]. 8 In the case where (p − k) D1-branes have bubbled, the theory living on the 't Hooft defect is described by a N = (0, 4) quiver SQM with quiver:  Figure 1. In this figure we show the brane configuration describing singular monpole configurations on a stack of D3-branes (a). In (b) we show the singular limit in which we have a reducible monopole given by many spatially coincident NS5-branes in the x 1,2,3 -directions.
where the node k is repeated 2p − 2k + 1 times 9 and where v is the effective 't Hooft charge. We will refer to this quiver as Γ (P, v).
In the case of a 4D theory with adjoint matter, the bubbling SQM is enhanced to a N = (4, 4) quiver SQM, and in the case of N f fundamental hypermultiplets, it is modified by coupling to N f short fundamental Fermi multiplets to the central gauge node:

Localization
We can compute the expectation value of an 't Hooft defect on R 3 ×S 1 using localization. Due to the fact that we are on a non-compact space, we must introduce an IR-regulator which is 9 In the case where p = k, there is a single U (k) node with 2 fundamental hypermultiplets adjoined to it.
accomplished via a 1 2 Ω-deformation with corresponding fugacity + . The expectation value of L p,0 can be defined as a supersymmetric index ⟨L p,0 ⟩ = Tr H L p,0 (−1) F e −βH+ +J++iQ⋅Θ+mf ⋅F , (3.1) where H L p,0 is the Hilbert space of the quantum field theory with an insertion of L p,0 , β is the radius of the thermal circle, H is the Hamiltonian, J + and F are the generators of the 1 2 Ω-deformation and flavor symmetries with fugacities + and m f respectively, Q = (γ e , γ m ) is the vector of electric and magnetic charges, and Θ = (θ e , θ m ) is the vector of electric and magnetic theta angles.
Here we will be fixing the electric and magnetic theta angles θ e , θ m . The electric theta angle can be defined in terms of local field configurations by the holonomy of the gauge connection along the thermal circle at infinity (3. 2) The magnetic theta angle, is similarly naturally defined as the holonomy of the dual magnetic gauge field. However, in order to fix both θ e and θ m , we must define the magnetic theta angle as the Fourier dual of path integral with fixed magnetic charge ⟨L ⃗ p,0 ⟩ γm Implementing localization leads to an expectation value which is of the form of (3.5). The details of the localization computation are contained in [4,22,25].

Expectation Value of Line Defects
In the case of G = SU (2), as we consider here, the expectation value of an 't Hooft defect that is computed by localization is naturally written in terms of two complexified Fenchel-Nielsen coordinates: a, b. These form a maximal set of holomorphic coordinates on the Hitchin moduli space M ζ . Associated to a weak coupling description, these coordinates have a semiclassical expansion where θ m and θ e are the magnetic and electric theta angles, ζ −1 Φ ∞ = Y ∞ + iX ∞ are the real and imaginary parts of the phase rotated vev of the adjoint-valued Higgs field Φ of the N = 2 vectormultiplet, and ϑ is the real part of the complex gauge coupling τ . Note that we will generally take Φ ∞ to be fixed so that a, b have ζ-dependence via X ∞ , Y ∞ . Additionally, here β is the radius of the thermal circle and (...) correspond to non-perturbative corrections, which we will discuss later in Section 5.2. More exact expressions for them are given in Section 4.1.2.
In these coordinates, the expectation value of an 't Hooft defect is written [22,25] ⟨L p,0 ⟩ = Here the sum is over monopole bubbling configurations labeled by the effective 't Hooft charge v = diag(v, −v). In each summand, the contribution cosh(v, b)(F (a)) v can be attributed to the contribution of the bulk fields whereas Z mono (P, v) can be attributed to the bubbling SQM living on the world volume of the 't Hooft defect in the given bubbling configuration. As in [4,6], we will focus on the contribution coming from Z mono (P, v). This contribution can be understood as the Witten index I W of the bubbling SQM as described in Section 2.1. 10 These theories naturally have an action which is Q-exact: where Q is a real supercharge that satisfies: where Q a , J + , and F are charges for the flavor symmetry corresponding to global gauge transformations, the R-symmetry corresponding to the 1 2 Ω-deformation, and other flavor symmetries. These generate a group of conserved symmetries of the bubbling SQMs which we will denote T = T G × U (1) + × T f where T G , T f are the maximal tori of the gauge and flavor group respectively.
Because the bubbling SQM is described by a N = (0, 4) quiver SQM, the path integral naturally localizes to a T -equivariant integral over the vacuum moduli space. In order to use localization, we must introduce an FI-parameter in the SQM. This lifts the Coulomb and mixed branches and causes the path integral to localize to a finite dimensional (T -equivariant) integral over the Higgs branch which is isomorphic to the Kronheimer-Nakajima space/quiver moduli space corresponding to the quiver defining the bubbling SQM: Γ (P, v).
It will be more convenient in the following discussion however to identify the quiver Γ (P, v) by its gauge and hypermultiplet nodes. That is to say we will define Γ ( ⃗ k, ⃗ w) ∶= Γ (P, v) where ⃗ k = (k 1 , ..., k n ) encode the gauge nodes and ⃗ w = (w 1 , ..., w n ) encode the fundamental hypermultiplet nodes. In this notation, we can identify the Higgs branch of the bubbling SQMs with the Kronheimer-Nakajima space M KN ( ⃗ k, ⃗ w) associated to the quiver Γ ( ⃗ k, ⃗ w). 10 There is an additional subtlety in the case of the SU (2) N f = 4 theory (and indeed in any SU (N ) Here the Zmono(P, v) has the interpretation as the ground state index I H 0 = lim β→∞ I W (β).
Additionally, there are some issues that arise with applying localization to the related bubbling SQM which requires adding an additional contribution where Iasymp counts ground states on non-compact directions of field space with finite potential energy. See [2,6] for more details.
M KN ( ⃗ k, ⃗ w) is defined as the hyperkähler quotient with respect to the moment maps 10) and the gauge group G = ∏ i U (k i ) acts as Physically, the B i,i±1 can be interpreted as the bosonic fields of the bifundamental hypermultiplets and the I i , J i can be interpreted as the bosonic fields of the fundamental hypermultiplets.

Z mono in terms of Characteristic Numbers
Thus far we have not taken into account the varying matter content of the 4D theories.
Coupling the 4D theory to matter hypermultiplets gives rise to a vector bundle, which we will call the matter bundle, over the same moduli space M KN ( ⃗ k, ⃗ w). This modifies the Tequivariant integral so that we can identify Z mono (P, v) with the T -equivariant integral over where ω is the real symplectic (1,1)-form in some choice of complex structure (which we will define shortly), µ T is the moment map for the action of T , and char T (T M KN ) is a T -equivariant characteristic class on T M KN depending on the matter content of the theory. This equation is not strictly well defined since M KN ( ⃗ k, ⃗ w) is a singular space. However, the singularities can be resolved by introducing an FI-parameter ξ [28]. We will denote the resolved space M ξ KN ( ⃗ k, ⃗ w). In 4D, the parameter ξ comes with a choice of direction of on R 3 (and consequently on the resolved space M ξ KN ) which corresponds to deforming the product in that direction (we have chosen theẑ-direction). The choice of direction defines a complex structure on M ξ KN ( ⃗ k, ⃗ w) and it is in this complex structure that ω is the (1,1) symplectic form. Note that by comparing to the definition of the reducible 't Hooft defect (2.4), computing the expectation value ⟨L p,0 ⟩ requires taking the limit ξ → 0. See [4] for more details. Now we can use the result of [4,25] for the characteristic class for the 4D asymptotically free SU (2) theories with fundamental and adjoint matter to write the monopole bubbling contribution See [5,25,34,37,38,42,43] for details on computing these integrals.

Spectral Networks
In this section we will review the technology of spectral networks [18,19]. Spectral networks are a method for computing the the holonomy of flat connections on a Riemann surface, which, as discussed previously, can be identified with the expectation value of a line defect L in the class S construction. It provides a set of "Darboux coordinates" Y γ on M f lat (C; G C ), the moduli space of flat G C connections on C associated to γ ∈ H 1 (Σ; Z). Spectral networks naturally computes the expectation value of the 't Hooft defect as a Laurent series in these coordinates whose coefficients are framed BPS indices.

Theories of Class S
Theories of class S are those which are constructed by taking the six-dimesional N = (2, 0) theory and compactifying it along an oriented Riemann surface C with a topological twist [13,16,27,46]. For the type SU (N ) theories of class S, this can be described as the low energy effective theory of a stack of N M5-branes wrapped on C ×M 4 with the same topological twist where M 4 is 4D spacetime. When M 4 = R × S 1 where S 1 has radius R the space of quantum vacua can be identified with the moduli space of solutions to the Hitchin equations on C Given a solution of these equations, we can identify the Seiberg-Witten curve and differential as where (x, z) ↦ xdz are coordinates on T * C. Intuitively, the Coulomb branch vacuum is described by separating multiple M5-branes wrapped on C along a transverse direction to produce a single M 5 brane wrapped on Riemann surface Σ which is an N -branched cover Σ → C.

Labeling Line Defects In Class S
As described in Section 7 of [17] a natural class of line defects in class S theories is obtained by considering a semi-infinite M2-brane that ends on a one-dimensional submanifold P ⊂ C times a line in M 4 . The line defect is labeled by a representation R α of SU (N ) associated to each connected component P α of P and a phase ζ, where ζ determines the unbroken supersymmetry. (Geometrically ζ describes how the M2-brane extends in the extra dimensions.) We can denote the resulting line defect by L(R, P, ζ). We now consider the case where the line defect wraps the circle of radius R in M 4 at a fixed point in R 3 . In this case the vacuum is described by a solution to Hitchin's equation to which we can associate a flat G C connection Because these theories are partially topologically twisted, the expectation value of such a line defect can be expressed in terms of the holonomy of this flat gauge field [16,17]: In this paper we will be focusing on the case of theories of class S with gauge group SU (2). We will take all the representations R α to be the fundamental representation. Therefore, the line defects can be labeled by L(P, ζ) where P is a smooth one-dimensional submanifold of C. Isotopy classes of such submanifolds can be conveniently labeled, given a pants decomposition of C in terms of Dehn-Thurston parameters: 11 Theorem (Dehn-Thurston): [7,44] Let C be an oriented Riemann surface with negative Euler characteristic that has genus g and n punctures.
be a maximal set of non-intersecting curves defining a pants decomposition of C and let {γ i } 3g−3+2n i=3g−3+n+1 be a collection of simple closed curves near the punctures. There is a mapping where I(C) is the set of isotopy classes of closed one dimensional submanifolds, q i is the twisting number with respect to γ i , and ⟨ , ⟩ is the intersection number. Elements in the image of D are denoted (⃗ p, ⃗ q) and are called Dehn-Thurston parameters.
above correspond to a weak coupling decomposition of the UV curve C, and specifies a Lagrangian duality frame with gauge algebra su(2) ⊕h with h = 3g − 3 + n. Each curve corresponds to a weakly coupled SU (2) gauge group in the 4D theory. Now consider the line defect associated to a generic 1D submanifold γ ⃗ p,⃗ q with Dehn-Thurston (DT) parameters (⃗ p, ⃗ q) = (p 1 , ..., p h , q 1 , ..., q h ). This submanifold will have a set of connected components γ ⃗ p,⃗ q = ⊕ k α=1 γ (α) ⃗ p,⃗ q labeled by α, each of which has its own Dehn-Thurston parameters: The line defect L(γ ⃗ p,⃗ q , ζ) then decomposes as a product of line defects In [11] it is conjectured that the line defects L(P, ζ) are the same as the 't Hooft-Wilson line defects of the Lagrangian theory with gauge algebra su(2) ⊕h . Moreover, it is proposed that the Dehn-Thurston parameters should be identified with the 't Hooft-Wilson parameters characterizing the magnetic and electric charges. This cannot be true in general, but it seems highly plausible for those Dehn-Thurston parameters that correspond to one-dimensional submanifolds γ ⃗ p,⃗ q with only one connected component. In this case the proposal of Drukker-Morrison-Okuda is that L(P, ζ) corresponds to the 4D line operator where h I(j) is the simple magnetic weight, λ I(j) is the simple weight of the j th factor of the gauge group, and h = 3g − 3 + n. It should be stressed that some more work is needed to make use of this conjecture: In mathematics it is not known what conditions one should put on the Dehn-Thurston parameters (⃗ p, ⃗ q) in order for γ ⃗ p,⃗ q to have a single connected component! The only case where this is known is the once-punctured torus (corresponding to the G = SU (2) N = 2 * theory) and the four-punctured sphere (corresponding to the G = SU (2) N f = 4 theory) [12]. In that case there are only a pair of DT parameters (p, q) and γ (p,q) has g connected components, where g is the gcd of p and q.
More generally, in the case where the four-dimensional gauge group is G = SU (2) we have only a pair of DT parameters (p, q). Here the minimally charged 't Hooft defect corresponds to the line with DT parameters (1, 0) which can be identified with the highest weight representation R h 1 of SU (2) ∨ . Following the decomposition above, a line defect corresponding to DT parameters (p, 0) is the p th power of the simple 't Hooft defect Thus, we see that the 't Hooft defect corresponding to L γ (p,0) is reducible. By equation (4.5) the vev is the trace of the holonomy in the representation R ⊗p h 1 . This is the origin of our notation from (2.3) By contrast L [ph 1 ,0] corresponds to a trace in the representation R ph 1 . If the vev of L p,0 is expressed as (2 cos θ) p then the vev of L [ph 1 ,0] has vev

Complexified Fenchel-Nielsen Coordinates
Because the expectation values of line defects in theories of class S are given by the trace of the holonomy of a flat connection, they are holomorphic functions on Seiberg-Witten moduli space. This allows them to be expressed via the AGT correspondence in terms of a particular set of holomorphic coordinates called complexified Fenchel-Nielsen coordinates. These can be defined as follows.
Choose a weak coupling region of the Coulomb branch. This defines a complex structure and comes with a maximal set of non-intersecting curves that are not isotopic to punctures on the UV curve C which correspond to weakly coupled gauge groups indexed by i. 12 Associated to each γ i , we can define the holomorphic coordinates {a i } ∈ t C defined by (4.12) The {a i } are Poisson commuting with respect to the standard, symplectic (2, 0)-form Ω J on Seiberg-Witten moduli space [9,26] and form a maximal set of Poisson commuting holomorphic functions. Now we can define a set of symplectically dual coordinates {b i } ∈ t C with respect to Ω J such that (4.14) We can then fix the redundancy b where ∂ a i f j = ∂ a j f i by specifying the semiclassical limit as in (3.4).
In the case of a single SU (2) gauge group, as we consider here, the above discussion reduces to a single pair of Fenchel-Nielsen coordinates a, b. These define the Fenchel-Nielsen coordinates a, b used for localization in Section 3.

Spectral Networks in Theories of Class S
For this paper, we are only considering the case of spectral networks in theories of class S of type G = SU (2), G C = SL(2; C). However, the following discussion generalizes to all ADE-type gauge groups [17-19, 32, 33].
A spectral network W subordinate to the covering Σ → C is a collection of oriented, open paths w on C called walls that begin at branch points and flow to punctures of C or other branch points. We will be interested in a special class of spectral networks which arise naturally in theories of class S called WKB spectral networks. These are defined by a meromorphic, quadratic differential ϕ 2 on the closure C of C and a choice of ϑ ∈ R 2πZ. 13 Locally, the quadratic differential can be written This can be used to define a foliation of C by curves γ which satisfy where where t is an affine parameter for γ. The corresponding spectral network W(ϕ 2 , ϑ) is then defined by the critical graph of the foliation defined by ϕ 2 and ϑ -i.e. the set of limiting curves that divide the foliation into distinct sectors. The spectral network technology developed by [18,19] provides a trivialization of the SL(2; C) vector bundle over the complement C W and gives gluing conditions across the walls. This trivialization is subordinate to the covering Σ → C such that we can locally equate the space of flat SL(2; C) connections on C to the space of flat which allows us to compute holonomies of the non-abelian vector bundle E → C in terms of holonomies of the connection of a flat line bundle on L → Σ.
The moduli space M f lat (Σ, GL(1; C)) has a natural set of coordinates: where ∇ ab is the connection on L. These coordinates follow the multiplication rule where ⟨ , ⟩ is the oriented intersection pairing. Further, the Y γ satisfy where γ b is a small loop around a branch point b and ω ∶ Σ → Σ is the map that exchanges sheets of the covering Σ → C. Let us define Σ ′ = Σ {ramification points}. For generic W, we can fix a basis of {γ i } ∈ H 1 (Σ ′ ; Z) ⟨γ + ω * γ⟩ to form our coordinate system on M(Σ, GL(1)).
These Y γ i are the Darboux coordinates related to the spectral network W.
The holonomy of the flat non-abelian gauge connection along a 1-cycle P can now be given in terms of Y γ associated to the decomposition of P with respect to a basis of open paths on C W. These open paths must then be joined across the walls of the spectral network by transition functions such that the complete expression for the holonomy is given by the product of the holonomies along these open paths connected by transition matrices between the different regions of C W.
The holonomy of an open path in a single connected region of C W can be written as where D P corresponds to the path across a simple open region andD P corresponds to a path that crosses a branch cut. Additionally, by making a convenient choice of trivialization, we can write the holonomy across a generic wall [17][18][19] where S w is some function of Y γ i . The case where a path crosses a double wall is a bit more subtle. To compute the parallel transport across a double wall, one must infinitesimally displace the phase ζ (requires picking a resolution convention) so that the double wall is replaced by a pair of generic walls; then one can compute the holonomy across the "double wall" by using the rules above. This prescription allows one to compute the holonomy of a complexified flat gauge connection along any path in terms of Darboux coordinates defined by the spectral network (spectral coordinates).
In addition to these rules, there are also several consistency conditions that restrict the number of free spectral coordinates. These come from abelian gauge symmetry on open path segments and from imposing monodromy conditions around branch points and punctures. This gauge symmetry acts by rescaling the spectral coordinate Y γ ij by a function corresponding to the end points of the curve γ ij where the beginning and end points of γ ij are labeled by i and j respectively. Note that the trace of the holonomy around closed paths are invariant under such gauge transformations. The consistency conditions we impose for monodromy around a branch point b and puncture p is that which come from the condition in (4.20) and the trivialization of the vector bundle at the punctures given in the data specifying the spectral network. Here ∇ is the connection on E → C which locally can be related to the connection on L as explained in [18,19].
In the class S construction, there is a natural choice of quadratic differential, given by the square of the Seiberg-Witten differential In this setting, we have that the charge lattice of the 4D theory is given by Γ ≅ H 1 (Σ ′ ; Z). Here there is a natural central charge function given by This allows us to express the corresponding Darboux coordinate as a function of γ as [15] log where Θ ⋅ Q γ is the Cartesian product of the vector of electric and magnetic theta angles with the vector of electromagnetic charges associated to γ. It is important to note that these coordinates generically have non-perturbative corrections to the semiclassical expression which, while complicated, are known and given in explicit formulas in [15]. 14

Wall Crossing in Spectral Networks
An important feature of spectral networks is that they give us an excellent tool for understanding wall crossing. In this setting, wall crossing is realized by changes of topology of the spectral network W(ϕ 2 , ϑ) as we scan the phase ζ = e iϑ which can be lifted toζ ∈ C * . The locations of the critical phases ζ = ζ c where the spectral network undergoes topology changes lift to a co-dimension-1 "walls" in C * called and are called K-walls [18]. Physically, each K-wall corresponds to a co-dimension loci where ζ is aligned or anti-aligned with the phase of Z γ k . Here, the change in topology of the spectral network causes the Darboux coordinates to undergo a cluster-like transformation/mutation [17] where is a particular choice of quadratic refinement with respect to a choice of splitting of the charge lattice and γ = γ e ⊕ γ m [17]. 15 14 The leading terms were first worked out in unpublished work by B. Pioline and A. Neitzke and in unpublished work by F. Denef and G. Moore. The main point of [15] was to give the full non-perturbative answer. 15 We will be working in the semiclassical limit so that there is always an almost canonical choice of charge lattice splitting. However, since the expectation value of a line defect L P is defined by a path P ⊂ C which is independent of the topology of the spectral network, the expectation value is wall crossing invariant. This means that the Y γ undergo coordinate transformations which exactly cancel the wall crossing of the framed BPS indices. Thus, by studying the wall crossing properties of the Y γ , one can infer the wall crossing of framed BPS states. A nice feature of generic WKB spectral networks is that the walls provide an ideal triangulation of C. In these networks, the associated Darboux coordinates have a natural identification with the edges of the triangulation. These coordinates are given by the holonomy along the lift under the projection π ∶ Σ → C of a path running between the branch points of different triangles through a given edge of the triangulation. See Figure 2. We will use the notation where the Darboux coordinate associated to the edge E is denoted Y E .
In such spectral networks, the fundamental topology shift that occurs in wall crossing is a flip of the triangulation . See Figure 2. Explicitly, in a generic WKB spectral network and consider a quadrilateral with edges E 1 , E 2 , E 3 , E 4 with diagonal edge E, a flip on the edge E ↦ E ′ acts on the corresponding Darboux coordinates by: where the signed intersection pairing of the edges is ⟨E, In the case of theories with vectormultiplets, spectral networks can also undergo a topology change called a juggle. 16 This can be understood as an infinite sequence of flips involving a puncture that has the effect of twisting a wall that runs to a puncture until it completely encircles it [16]. See Figure 3.
The juggle can be understood as follows [16]. Consider an annulus surrounding a puncture, P (which we replace by a disk with a marked point), with a single vertex V of the triangulation on the outer boundary. Now consider lifting the configuration to the simply connected cover which is a triangulated infinite strip as in If we choose an ordering of the lifted images of the vertices, we can define a winding number of an interior edge by the difference of the image number of the end points. Further, we can iteratively increase (decrease) the winding numbers of the interior edges by performing a sequence of simultaneous flips on all of the preimages of the the interior edge with the lowest (highest) winding number. See Figure 5.
After n such flips, the interior edges run between the 0 th exterior vertex preimage to the n th and (n − 1) th interior preimage. We can now make sense of the corresponding Darboux coordinates in the limit as n → ∞. First note that as n → ∞ the interior edges approach a parallel line to the interior and exterior edges. This corresponds to a spectral network where there is a single, double wall circling the puncture of C under consideration. If we define Y + and Y − to be the edges with higher and lower winding number respectively after n flips, then 16 There is another transformation called a "pop" which has to do with changing the decoration of a given puncture, but this will not be important for our story. See [16] for more details. in the n → ∞ limit we can construct the well defined coordinates:

Fenchel-Nielsen Networks
Now we will discuss a special class of spectral networks called Fenchel-Nielsen networks [24]. These spectral networks have only double walls corresponding to a set of minimal cuts necessary to decompose the Riemann surface C into a disjoint product of punctured discs and annuli. This is a WKB spectral network where ϕ 2 is a Jenkins-Strebel differential -ϕ 2 gives a foliation of C by closed paths. Another way of saying this is that a Fenchel-Nielsen spectral network is given by a pants decomposition of C in which on each pair of pants, the spectral network is one of the two networks in Figure 6.  These spectral networks are referred to as Fenchel-Nielsen-type because the a-type Fenchel-Nielsen coordinate has a straightforward interpretation in terms of the associated spectral network coordinates. associated to these networks have a straightforward interpretation as complexified Fenchel-Nielsen coordinates.
Let us take a maximal set of non-intersecting curves that define a pants decomposition of C. On each pair of pants, there are classes of curves which are homotopic to ζ Figure 7. This figure shows the structure of the K-walls in the ζ-plane. There are accumulation points (red) on the imaginary axis where the associated WKB spectral network becomes a Fenchel-Nielsen spectral network. a subset of the {γ i }. The holonomy around a curve that is homotopic to such a γ i is given in terms of the spectral network coordinates 17 However, we see from before, that this is simply the definition of the Fenchel-Nielsen coordinate a: from [17].) The relation to the above two-term expansion is clarified in equation (6.22) above.
In a large class of theories, such as the ones we study here, Fenchel-Nielsen spectral networks can be obtained from a generic WKB spectral network by performing a juggle. This requires changing ζ such that we cross an infinite number of K-walls. In the theories we consider, there are infinite number of such walls which accumulate along co-dimension 1 "accumulation points" in the ζ-plane. See Figure 7. In our setting, sending ζ to an 17 Note that Yγ are defined for γ ∈ H1(Σ; Z) while Lγ is defined for γ ⊂ C. Here we use the loose notation where Yγ for γ ⊂ C is defined as Y π −1 (γ) i the lift under the projection π ∶ Σ → C onto one of the sheets. Due to (4.20), the two choices of lifting are related by inverses and thus are merely a choice of convention. accumulation point is equivalent to undergoing the infinite number of flips that occur in a juggle, leading to a Fenchel-Nielsen spectral network. See Section 5.2 for further discussion.
Remark Recall from the discussion of Section 4, that the spectral network coordinates Y γ is given in (4.28) and has a semiclassical expansion with an infinite number of nonperturbative corrections. Since, as we showed above, we can identify the complexified Fenchel-Nielsen coordinates with spectral network coordinates, the Fenchel-Nielsen coordinates a, b must similarly have an infinite number of non-perturbative corrections to their semiclassical value. We will demonstrate this in the example of the SU (2) N f = 0 theory in Section 6.2 by computing the leading non-perturbative corrections.

Semiclassical Formulation of BPS States
Because the expectation value of a line defect in a theory of class S can be expanded as a series of Darboux coordinates with coefficients that are framed BPS indices (4.1), the expectation value of the line defect is entirely encoded in the spectrum of framed BPS states. In the semiclassical limit of Lagrangian 4D N = 2 theories, the spectrum of framed BPS states can be described by the index of a Dirac operator on singular monopole moduli space [3,39,40].
The identification of BPS states with the kernel of a Dirac operator arises from the effective description of the dynamics of BPS states in the adiabatic limit via collective coordinates. The resulting theory is a SQM on (bundles over) singular monopole moduli space with potential [3,20,21,35,39,40,45]. Solving for the BPS spectrum is reduced to solving the Dirac equation on singular monopole moduli space coupled to certain bundles over M BP S .

Moduli Space Approximation
Consider the semiclassical limit of a four-dimensional N = 2 supersymmetric gauge theory with hypermulitiplet matter and 't Hooft line defect insertions where X ∞ >> Y ∞ . The BPS equations for this theory are given by where ζ −1 Φ = Y + iX is the decomposition of the N = 2 vectormultiplet scalar into real and imaginary parts and ζ is specified by the 't Hooft defect. In the background defined by the asymptotic boundary conditions We can consistently describe the dynamics of BPS states up to order O(g 2 ) as dynamics on the moduli space, M(P, γ m ; X ∞ ). Let us now introduce coordinates {z m } on M. These coordinates will be elevated to time dependent fields 3) whose variation will describe the dynamics of the BPS states.
Since we are considering the dynamics of BPS states which are 1 2 -SUSY, there are still preserved supercharges. This indicates that the z m (t) have super-partners χ m (t) coming from non-trivial fermionic zero modes.
Similarly, there are also fermionic zero-modes ψ s (t) coming from the hypermultiplet fermions that couple the theory to the spin-bundle associated with a vector bundle E matter → M of rank [3]  and {J r } are the complex structures on T M and G(Y ∞ ) m is a triholomorphic vector field that generates a global gauge transformation 19 where ζ is defined by the line defect.
In the Hamiltonian formulation, the supercharge operator becomes a twisted Dirac op- where D n is the spin coviariant derivative on Spin(E matter ) ⊗ SM → M. Thus, BPS states are those that are in the kernel of the supercharge Dirac operators. Now from the supersymmetry algebra 11) we see that if a BPS state is in the kernel of any of the supercharges, then it is in the kernel of all of the supercharges. Therefore, the spectrum of BPS bound states can be determined by considering the kernel of the 4 th superchargê Additionally, since the electric charge and flavor charge operators commute with the Hamiltonian, and hence the supercharge operator, we can simultaneously diagonal the two. This leads to a decomposition of the Hilbert space of BPS states into γ = γ m ⊕ γ e ⊕ γ f -isotypical where here we mean the γ e ⊕ γ f -isotypical component of the L 2 index of the Dirac operator D Y on M(P, γ m ; X ∞ ).

Fenchel-Nielsen Spectral Networks and The Semiclassical Region
Now we can connect the formalism of semiclassical BPS states and spectral networks. Recall that the the expectation value of line defects is determined by the framed BPS index as in (4.1). Therefore, the index of the Dirac operator from the previous section can be used to determine the expectation value of an 't Hooft defect in the semiclassical limit.
Since we want to compare to the localization computation, which is naturally expressed in terms of Fenchel-Nielsen coordinates, one would hope to use Fenchel-Nielsen spectral networks and the associated Dirac operators. In order to implement this we need to know: 1.) if there exists Fenchel-Nielsen spectral networks in the semiclassical limit and 2.) where in parameter space these spectral networks exist so that we can compare to indices of Dirac operators. In this section we will show that such spectral networks exist in the semiclassical limit, but that they only exist in parameter space where the moduli space approximation breaks down.
The question of whether or not a Fenchel-Nielsen network exists is equivalent to the question of whether or not there exists a Jenkins-Strebel differential on C that encodes the data of the theory in some semiclassical limit. The data of the differential is (u, ζ, m) ∈ The existence of Jenkins-Strebel differentials on a Riemann surface C with punctures are studied by Liu [29,30]. There, Liu shows that given a decomposition of C into a collection  of punctured disks {D m } and annuli {R k }, there exists a uniquely determined real Jenkins-Strebel differential with closed trajectory ϕ 2 with fixed monodromy m i ∈ R around each puncture and height h k ∈ R around each annuli where the height is defined as where the infimum is taken over all paths that run between the boundaries of R k . Note that the Fenchel-Nielsen spectral network is exactly given by the union of the boundaries of these component disks and annuli. Now consider as an example the case of the 4D SU (2) N = 2 * theory. This theory is constructed as a theory of class S by taking C to be a torus with a single puncture. This theory comes with a complex 2-dimensional parameter space defined by u ∈ B ≅ C and the complex mass parameter of the hypermultiplet. C can be decomposed as an annulus R a and a punctured disk D m . See Figure 8 for the example of the 4D SU (2) N = 2 * theory where C = T 2 {0}. 20 Thus, there is a 3 dimensional family (specifying m, γ, and ζ) of Jenkins-Strebel differentials which forms a real co-dimension 1 subspace of parameter space. This suggests that there could exist a Jenkins-Strebel differential in the semiclassical limit ( u → ∞ ) and therefore that there could exist a Fenchel-Nielsen spectral networks in the semiclassical limit. This has been confirmed by numerical computations. 21 Now recall that for a WKB Fenchel-Nielsen spectral network, the real Jenkins-Strebel 20 We would especially like to thank Pietro Longhi for providing these figures. 21 We would like to thank Pietro Longhi for sharing his numerical computation for the SU (2) N = 2 * theory and for making the authors aware of Liu's work on Jenkins-Strebel differentials.
differential is related to the Seiberg-Witten differential as Asking that ϕ 2 as defined by this equation is a Jenkins-Strebel differential defines the Fenchel-Nielsen locus in B * × C * . As usual in Seiberg-Witten theory, the periods of λ SW give the vev's of the Higgs field and mass parameters. In our case the UV curve is given by C = T 2 {0}. This means that if we pick a basis of H 1 (Σ; Z) = span Z {A, B}, where D p is a loop circling the puncture and m f is the mass of the adjoint hypermultiplet. In this notation, the condition that ϕ 2 is a Jenkins-Strebel differential (and hence gives rise to a Fenchel-Nielsen-type WKB spectral network) is that which can be rewritten as This locus in parameter space, which we will call the Fenchel-Nielsen locus, is an accumulation point of K-walls in the ζ-plane and we will denote the associate phase in U (1) as ζ F N . Unfortunately, the Fenchel-Nielsen locus is exactly where the moduli space approximation, which gives the identification between the framed BPS index and the index of a Dirac operator on singular monopole moduli space, breaks down. In the limit X ∞ → 0, the space M(P n , γ m ; X ∞ ) (and M(γ m ; X ∞ )) are not defined. The reason is that the semiclassical expression for the central charge is given by Thus the BPS mass M BP S = Re[ζ −1 Z γ ] for a monopole goes to zero as we scan ζ such that X ∞ → 0. However, we know that monopoles do not become massless in the semiclassical limit. Thus, we can deduce that the non-perturbative quantum effects must become large and therefore the effective SQM description above must break down. However, by taking X ∞ , Y ∞ → ∞ as X ∞ Y ∞ → 0, we can still identify framed BPS indices with the index of a Dirac operator for phases which are arbitrarily close to the Fenchel-Nielsen locus. This will allow us to give an index theorem for the supercharge Dirac operators almost everywhere on the ζ-plane. See Figure 1 of [40] or Figure 4 of [39] for more details.
The above analysis makes it clear that there always exists SU (2) Fenchel-Nielsen networks (and in fact all SU (N )-type Fenchel-Nielsen spectral networks) in the semiclassical limit. These exist on the locus where all of the masses and a i = ∮ A i λ SW have the same phase. Such a spectral network can be constructed by gluing together pairs of pants with semiclassical Fenchel-Nielsen spectral networks on them by the procedure of [24]. The only condition here is that the Fenchel-Nielsen spectral networks all have the associated phase.
Remark Recall that a Fenchel-Nielsen spectral network corresponds to a WKB spectral network with a Jenkins-Strebel differential. This is defined by decomposing the Riemann surface C into a collection of annuli and punctured disks. On each component, the flow lines of ϕ give a foliation of curves that are homotopic to the boundary components. If we consider infinitesimally deforming the phase ζ associated to the quadratic differential, we find that the flow lines on each component are no longer homotopic to the boundary components, but rather spiral into them with a very large winding number. Thus, as we send ζ → ζ F N the flow lines of ϕ 2 twist around the boundary components infinitely many times until they form closed paths, producing a Fenchel-Nielsen spectral network. This infinite spiraling indicates that Fenchel-Nielsen spectral networks can be achieved by performing a juggle on a WKB spectral network where all physical parameters have aligned phases. Using the procedure from Section 8.4 of [24], one can identify the limiting coordinates (4.33) with the Fenchel-Nielsen coordinates Y B defines a choice of e b . 22 Thus, the Darboux coordinates associated to Fenchel-Nielsen spectral networks in the cases we are studying can be obtained by acting on a generic set of spectral network coordinates by an infinite number of cluster coordinate transformations. The resulting spectral network coordinates are those that result from the flip (4.33). This gives a recursion formula for the Darboux coordinates that can be "integrated" to give a relation between the Darboux coordinates of a spectral network in any chamber and the Fenchel-Nielsen coordinates which are used in localization computations. This will be the primary computational tool that we will use to construct an index theorem and give a formula for the characteristic numbers in the next section.

Index Theorem and Characteristic Numbers
In this section we will compare the different methods of computing the expectation value of 't Hooft defects in 4D N = 2 G = SU (2) asymptotically free theories with fundamental and hypermultiplet matter. We will outline how this comparison can be used to give an index theorem for Dirac operators on singular monopole moduli spaces and give the characteristic numbers of certain Kronheimer-Nakajima spaces 23 . We will explicitly show these for the SU (2) N f = 0 theory. 22 Note that we could also approach the Fenchel-Nielsen locus in the opposite direction. The procedure from [24] in conjunction with the relation between the two limiting coordinates (4.34), correspond to two different choices of Fenchel-Nielsen coordinates. Equation (7.52) of [16] shows that a is well-defined and Y ± B define two choices of b that are related by (4.34). 23 These are the transversal slice to the stratum of the bubbling locus of singlar monopole moduli spaces.

General Theory
Consider an N = 2 SU (2) Lagrangian theory of class S with mass parameters of identical phase. Now pick a point in the semiclassical limit of the Coulomb branch away from the Fenchel-Nielsen locus. We are interested in computing the expectation values of a 't Hooft defect which is specified by an integer p and a phase ζ. Now consider comparing the localization and spectral network result for the expectation value of 't Hooft defects. Localization requires introducing an IR regulating 1 2 Ω-deformation and expresses the expectation value in terms of Fenchel-Nielsen coordinates. This coordinate expansion is well defined almost everywhere in a simply connected region on the Coulomb branch and is independent of the phase of ζ there due to trivial monodromy. The spectral networks computation however, is not independent of the phase ζ. Rather, it is different in each chamber c n ⊂ C ζ of the ζ-plane. 24 The reason is that the spectral network undergoes topology change at each K-wall and hence has a different set of associated Darboux coordinates in each chamber c n ⊂ C ζ .
Away from the Fenchel-Nielsen locus, the spectral network coordinates are not Fenchel-Nielsen coordinates, but rather are Darboux coordinates which are related to the localization Fenchel-Nielsen coordinates by an infinite sequence of Kontsevich-Soibelman transformations.
Due to the "simple" transformation properties of the spectral network coordinates, these coordinate transformations can be integrated to determine the mapping between Fenchel-Nielsen coordinates and the Darboux coordinates in every chamber. This can be achieved as follows. First, solve for the expectation value of the minimal Wilson and 't Hooft defects in a generic WKB spectral network of choice. We will assign the chamber in the ζ-plane in which we have computed these as the c 0 chamber. Now by tuning the phase of ζ, we will cross walls of marginal stability which takes us from the c n chamber to the c n±1 chamber depending on the direction we tune ζ. Now we can solve for the expectation values in all chambers by solving the recursive K-wall crossing formulas [17]: where ⟨L⟩ ζ∈cn is the expectation value of L computed using the WKB spectral network acssociated to ζ ∈ c n and the K-wall between the c n and c n−1 chamber is W (γ n ).
After we set the 1 2 Ω deformation parameter + → 0, we can then compare the localization expression of the expectation value of the Wilson and 't Hooft defects to their expression in terms of the spectral network coordinates in a generic chamber chamber c i . Inverting these formulas allows us to solve for (a, b) in terms of the Y γ in some fixed chamber c i .
Then, by combining this with the solution with the KS-wall crossing formulas, we then have an expression for the Fenchel-Nielsen coordinates (a, b) in terms of the Y γ in all chambers c n . Inverting the formulas, one obtains an (admittedly complicated) expression for the a, b in terms of the Y γ . 24 Here ζ changes the decomposition of u into X∞, Y∞.
We can then take these expressions and substitute the expression for a, b in terms of the Y γ into the localization expression above. By identifying the coefficients of the Laurent expansion with that of the spectral network computation we arrive at an expression for the framed BPS indices in every chamber. Then by using the relation of the index of the Dirac operator to the framed BPS indices in the semiclassical limit away from the Fenchel-Nielsen locus we get an index formula for the associated Dirac operator in all chambers arbitrarily close to the Fenchel-Nielsen locus. Similarly, we can substitute the expression for the Darboux coordinates Y γ in the c n chamber in terms of the Fenchel-Nielsen coordinates a, b into the spectral network computation. Then, identifying the coefficients of the Laurent expansion in terms of the exponentiated Fenchel-Nielsen coordinates on both sides allows us to express the characteristic numbers that determine Z mono (a, m; P, v) as a rational function of exponentiated Fenchel-Nielsen a-coordinates, masses, and framed BPS indices.
Remark Note that there is an additional subtlety in the case of the N = 2 SU (2) N f = 4 theory. The reason is that Z mono (a, m; P, v) is not entirely given by a characteristic number but rather has an additional contribution from states on the Coulomb branch of an associated SQM [6].
6.2 Example: SU (2) N f = 0 Theory Now we will apply the above, discussion to determine the framed BPS indices for the SU (2) N f = 0 theory. This will produce an index-like formula for a Dirac operator coupled to a hyperholomorphic vector field G n (Y ∞ ) on singular monopole moduli space. The expectation value of the 't Hooft defect in the N f = 0 theory is given in terms of Fenchel-Nielsen coordinates as where In the case where + = 0, which is necessary for comparing with the Fenchel-Nielsen and Dirac operator expressions, these have the simple form As shown in [17,39], the expectation value of the 't Hooft defect of minimal charge in terms of Darboux coordinates in the chamber c n is given by 25 where and Here we take n ∈ Z + to denote the chamber c n in the ζ-plane and the notation U n to denote the Tchebyshev polynomial of the second kind: In this theory, the Fenchel-Nielsen locus is given by ⟨α, X ∞ ⟩ = 0. 26 Let us pick a u ∈ B such that Φ ∞ = ±iX ∞ when ζ ∈ R and Φ ∞ = ±Y ∞ when ζ ∈ iR. We can now identify the imaginary axis as the Fenchel-Nielsen locus. The spectrum of the vanilla BPS states in the semiclassical region are given by γ = ±α , γ ± n = ±H α ⊕ n α , n ∈ Z , (6.11) with BPS indices Thus, the phase of the central charge corresponding to a state with charge γ = γ m ⊕ γ e is given by ) .
(6.13) 25 Recall that L1,0 is the minimal 't Hooft defect and L0,1 is the minimal Wilson defect. 26 Note, that here we have defined X∞ ∈ t to lie in the positive chamber. We will thus assume that ⟨α, X∞⟩ ≥ 0.
Without loss of generality, we can restrict to the case ζ in the positive real half-plane (the other cases follow analogously). We are now only concerned with the phase of BPS states whose K-walls are in the positive real half-ζ plane. These BPS states have charges γ − n with K-walls along the phases , (6.14) to order O(g 4 ). Note that the phases of the chentral charges are ordered 15) in the positive real half-plane. We can now define the chambers c n ⊂Ĉ * as At the K-wall defined by W (γ − n ), the X i,n−1 mutate as: By comparing with the computation of ⟨L 1,0 ⟩ and ⟨L 0,1 ⟩ using localization, we can determine the coordinate transformation relating the Fenchel-Nielsen coordinates to the spectral network coordinates in the c n chamber: We can now invert the expressions (6.18) to get (6.20) Note that both of these pairs of expressions requires matching the semiclassical expressions for X 1,0 , X 2,0 , e a , e b . This can be used to show explicitly that the spectral network coordinates approach the Fenchel-Nielsen coordinates as we approach the Fenchel-Nielsen locus. Sending the phase ζ → ζ F N can then be achieved by sending n → ∞ or n → −∞. Going to the Fenchel-Nielsen locus rotates the phase of Φ ∞ = ζ(Y ∞ + iX ∞ ) so that, in the limit n → ±∞, ±⟨α, Y ∞ ⟩ > 0. Then from the expression for Y γ (4.28), we see that lim n→±∞ Y γn is exponentially suppressed This allows us to identify Y 1 2 α = e ±a . 27 The coordinate identification above can also be seen from taking limits of the coordinate transformation expressions (6.18). As ζ → ζ F N , we pass through chambers with n → ±∞. Using the identification of the Fenchel-Nielsen coordinates with spectral network coordinates in c n (6.20), it is possible to compute explicitly the non-perturbative corrections to a, b. From the results of [15], we know that the corrections to the semiclassical contribution of Y γ are given by solving the recursive formula log Y γ (u, θ, ζ) = log Y sf γ (u, θ, ζ) where Y sf γ is the semi-flat term which is the semiclassical expression in (4.28) and γ ′ is the ray in the ζ-plane along the K-wall K γ ′ .
To compute the first order non-perturbative corrections to a, b, we need to use the fact that the semi-flat expressions for X e , X m are of order X sf e ∼ O(1) while X sf m ∼ O e − 4π g 2 . Using this, we can expand the expressions for a, b as a Laurent series in X m : where f n (X e ) = X e + n(X e − 1 X e ) . (6.28) 27 As we will see below, the signs are correlated with the different limits n → ±∞.
Then by taking into account the non-perturbative corrections to X e , X m and comparing orders in e − 4π g 2 , we can compute the first order non-perturbative corrections to the leading (semiclassical) expression for a, b as Here a (1) n.p. are the first of an infinite series of non-perturbative corrections to the semiclassical values of a, b are the leading order non-perturbative corrections to the semi-flat expressions for X m , X e which are given by [15]: Using the BPS indices (6.12), these integrals simplify to (6.32) where above we have used the notation where the integral over ζ ′ has been mapped to the integral over the positive reals by the phase rotation e iα γ ′ = phase(Z γ ′ ) and ⟪γ, γ ′ ⟫ is the DSZ pairing of charges. Note that the Y γ are functions of u, θ e , θ m , ζ on the Hitchin moduli space. Because of the relation between the a, b and the X e , X m in (6.18) and (6.20), we clearly see that the a, b must also be functions of u, θ e , θ m , ζ. The explicit dependence of a, b on ζ can be seen first fixing a point in Hitchin moduli space with fixed coordinates (u, θ e , θ m ), and then studying the Fenchel-Nielsen coordinates as functions of the complex structure ζ.

Index Theorem
We can now use these coordinate transformations to determine an index formula as follows: 1. Calculate the localization computation for the expectation value of the given line operator: is the corresponding Kronheimer-Nakajima quiver variety as described in Section 2.1. In this example it evaluates to (6.34) 2. Perform the change of coordinates: where ∆ ne,nm is a delta function that restricts the sum over the {m, j, k, , i, q, d i , i i } such that are fixed. Additionally, This index formula for the case of SU (2) SYM theory is also found in [39].

Characteristic Numbers
Now by expressing Y γ in terms of Fenchel-Nielsen coordinates, we can perform a Laurent expansion with respect to the exponential Fenchel-Nielsen coordinate e b . This will allow us to isolate the characteristic number. By using the equations for the Darboux coordinates in terms of Fenchel-Nielsen coordinates (6.18), we get the expansion 2n m n j 1 × j 2 + 2n m (n + 1) − 1 j 2 (−1) i 2 +j 2 +2nm a 2n(i 2 +j 2 )+(2n+2)(i 1 +j 1 )−2nm−ne−2v (1 − a 2 ) 2(nm+v) .  and the sums are restricted such that ∑ b n=a is identically zero for b < a. Note that in both the formulas for the index of D Y and the characteristic numbers on M KN , there is a clear mixing of framed BPS states of magnetic charge γ m among many characteristic numbers for different. This suggests that there is a very non-trivial relationship between framed BPS states and the geometry of singular monopole moduli space since the M KN are transversal slices to singular strata in M(P, γ m ; X ∞ ). It is an interesting challenge to differential geometers to try to prove such relations. in terms of the spectral network coordinates in the c 1 chamber. Note that this matches the direct computation from spectral networks (6.7) [39]. From this expansion we can read off the indices of the Dirac operator:

Explicit
Ind D where P = 1 2 diag(2, −2) and ζ ∈ c 1 . Now let us perform the inverse coordinate substitution to derive the characteristic number from the spectral network computation. Let us start with the expectation value of L 2,0 from the spectral network associated with ζ ∈ c 1 : (6.53) The coordiante transformation (6.18) now takes the form . (6.54) Plugging this into (6.53) we find

Comments on the N = 2 * Theory
Here we would like to make some clarifying comments on the SU (2) N = 2 * theory and the Fenchel-Nielsen locus in this theory. In the case of the N = 2 * theory the UV curve C is given by the once punctured torus. The algebra of line operators of this theory can be generated by the three simple line operators L γ (1,0) , L γ (0,1) , and L γ (1,1) . Note that there are three generating operators because the homology lattice is generated by a cycles that wrap the A-cycle, B-cycle, and the puncture. A generic spectral network associated to the SU (2) N = 2 * theory is given by an ideal triangulation of C as in Figure 9. In each chamber of the ζ-plane c, the charge lattice is spanned by three simple elements γ i [c] for i = 1, 2, 3 such that In such a chamber, the expectation values of the line operators can be expanded in terms as Here Y γ i is the spectral network coordinate corresponding to the edge E i . In the SU (2) N = 2 * theory, the Fenchel-Nielsen locus is defined bȳ where ζ is the phase defining the line operator (and corresponding WKB spectral network).
As it turns out, this coincides with the exceptional locus from [31]. Here mathematical simplifications arise that allow for the exact computation of the spectrum generator which encodes the entire spectrum of BPS states. Figure 11. This figure shows the behavior of the generic WKB spectral network on the punctured torus from Figure 10 as it approaches the Fenchel-Nielsen spectral network (left in Figure 6). Again, while it is not drawn here due to computational limitations, the walls corresponding to the open paths run to the (upper) puncture.