Generalized N=2 Topological Amplitudes and Holomorphic Anomaly Equation

In arXiv:0905.3629 we described a new class of N=2 topological amplitudes that depends both on vector and hypermultiplet moduli. Here we find that this class is actually a particular case of much more general topological amplitudes which appear at higher loops in heterotic string theory compactified on K3 x T^2. We analyze their effective field theory interpretation and derive particular (first order) differential equations as a consequence of supersymmetry Ward identities and the 1/2-BPS nature of the corresponding effective action terms. In string theory the latter get modified due to anomalous world-sheet boundary contributions, generalizing in a non-trivial way the familiar holomorphic and harmonicity anomalies studied in the past. We prove by direct computation that the subclass of topological amplitudes studied in arXiv:0905.3629 forms a closed set under these anomaly equations and that these equations are integrable.


Introduction
Over the last couple of decades, considerable effort has been put into the study of 1/2-BPS effective couplings in four-dimensional extended supersymmetric theories. Such couplings depend only on half of the superspace, generalizing the well-known chiral F-terms of N = 1 supersymmetry. They usually enjoy particular non-renormalization theorems which is the reason why they are easier to study from a field theoretic point of view. Moreover, these couplings play a predominant role in many very interesting applications, ranging from string phenomenology to entropy corrections of supersymmetric black holes. It is therefore of considerable interest to classify different kinds of 1/2-BPS couplings and study their properties in as much generality as possible.
In the string theory effective action, 1/2-BPS couplings are expected to be captured by topological amplitudes, i.e. amplitudes which only depend on the zero-mode structure of the internal manifold [1,2,3]. The best studied example in this respect is the celebrated series F g W 2g where W is the chiral N = 2 Weyl superfield of the supergravity multiplet. The coefficients F g are computed as g loop amplitudes in type II string theory compactified on a Calabi-Yau threefold [2,3] and are functions of the vector multiplet moduli in the Coulomb phase of the theory. In particular, their independence of N = 2 hypermultiplet moduli (and therefore also of the type II dilaton) implies a non-renormalization theorem that F g is g-loop exact. Moreover, classically their BPS nature dictates that the F g are holomorphic functions. This property is broken at the quantum level, however, in a controlled manner captured by a first order differential equation, termed the holomorphic anomaly equation [3]. The latter takes the form of a recursion relation in g which in fact allows in some cases to find explicit expressions for F g up to considerably high genus [4,5,6].
The heterotic versions of the F g 's have been considered in [7]. They are semi-topological expressions (i.e. topological only in the supersymmetric sector), related to F-terms of the formŴ 2g whereŴ is now the N = 1 gauge superfield. Again the BPS-nature of these couplings classically implies a (first order) holomorphicity condition for the F g which is broken at the quantum level. The difference to the type II side is, however, that the corresponding (integrable) holomorphic anomaly equation, which captures this breaking, does no longer close on the class of functions F g , but introduces new (semi)-topological objects. The latter have been understood to be physically related to BPS couplings of the form Π nŴ 2g where Π's are chiral projections of non-holomorphic functions of chiral superfields.
More recently, after studying a particular world-sheet involution of a certain class of N = 4 topological amplitudes [8,9] (see also [10] as well as the review [11]), we have presented a completely new class of N = 2 topological amplitudes [12]. The latter compute a corresponding class of 1/2-BPS terms in the low energy effective action of the form F g K 2g − , where K − is a particular superdescendant of the vector superfield in N = 2 harmonic superspace [13]. The novel property is that their coupling coefficients depend on both holomorphic vector multiplet as well as particular (Grassmann) analytic hypermultiplet moduli, despite the common wisdom of mutual decoupling. Classically, this particular analytic dependence -which is again a direct consequence of the BPS nature of the couplings -can be captured by differential equations, namely a holomorphicity relation for the vector moduli and a harmonicity relation together with another second order equation for the hypermultiplet moduli. As before, at the quantum level these differential equations get modified by anomalous string world-sheet boundary contributions.
In this work, we show that these new topological objects are in fact a particular subclass of an even further generalized set of N = 2 topological amplitudes. We make a systematic analysis of the latter on the heterotic string side (compactified on K3 × T 2 ) by first computing them as genus-g amplitudes involving an appropriate number of fermions (gauginos and hyperinos) and establish explicitly their connection with the (semi)-topological theory obtained by twisting only the supersymmetric (left-moving) sector. We then show that this generic amplitude corresponds to a new series of higher order couplings in the effective action, involving both vector multiplets and neutral hypermultiplets. They depend on the moduli (both vectors and hypermultiplets) in a quite generic manner, which particularly means that there is no immediate direct generalization of the holomorphicity or harmonicity equation. However, the BPS nature of these couplings manifests itself in another very useful property: The particular (Grassmann)-analytic projection which is necessary for consistency of these couplings leads us to establish relations between different component couplings. On the string theory side, these relations turn out to be highly non-trivial (first order) differential equations which again get modified by anomalous world-sheet boundary contributions. We explicitly compute these anomaly terms for these more general amplitudes (up to possible curvature dependent contact terms) as well as for the subclass considered in [12] and show that the equations close in the sense that no new topological objects appear. This indicates that the amplitude we consider captures the most generic expression in this particular class of BPS-couplings. We study the integrability of these equations and show that for the latter, the equations we have obtained are integrable including also the curvature terms.
The paper is organized as follows. In Section 2, we compute a particular multi-fermion genus-g physical amplitude and we show that it acquires a (semi)-topological expression as a correlation function in the twisted K3×T 2 heterotic σ-model. This amplitude has four types of indices, labeling the two different helicities (i.e. Weyl spinor components) of the antichiral gauginos and hyperinos involved. Contracting these indices (which have particular symmetry properties) with four fermionic variables, one can introduce a generating functional, in terms of which the results take a compact form and are greatly simplified. In Section 3, we describe the interpretation of the above amplitudes as a particular class of 1/2-BPS terms in the string effective action. In Sections 4 and 5, we derive differential equations for respectively the generalized amplitudes and the reduced ones that were obtained in [12]. As we show in Section 6, these relations follow from the particular (Grassmann)-analytic projection which we used in formulating the couplings. These equations are a direct consequence of the 1/2-BPS structure of these couplings and relate the anti-holomorphic dependence on the vector moduli with the dependence on the hypermultiplets of 'wrong' harmonicity in different component couplings. In string theory we find that (as usual) these equations get modified by an anomaly due to world-sheet boundary contributions. As a result, one obtains recursion relations for the non-holomorphic/harmonic moduli dependence of the above amplitudes. In Sections 4 and 5 we also study the integrability conditions of the anomaly equations. Our conclusions are presented in Section 7. Although the notation and conventions of this paper are those of Ref. [12], to make it self-contained, we include several appendices. Appendix A gives a brief review of N = 2 and N = 4 superconformal algebras, Appendix B summarizes the main properties of N = 2 harmonic superspace, Appendix C presents some essential features of the heterotic K3 × T 2 compactification, Appendix D describes the gauge freedom associated to one of the two general series of couplings, and Appendix E contains the direct string derivation of the differential equations, as well as the computation of a particular world-sheet boundary contribution which we deemed too lengthy to be presented in the main body of the paper.

A Further Extension of a Class of Topological Amplitudes
In the paper [12] by studying higher derivative couplings of holomorphic N = 2 vector multiplets with hypermultiplets of a particular analyticity in heterotic string theory compactified on K3 × T 2 , a new class of N = 2 topological amplitudes was discovered. They can be expressed as the following (semi-)topological correlators Here Mg denotes the integral over the moduli space of Riemann surfaces of genus g. Deformations of the latter are parametrized by a total of 3g − 3 Beltrami differentials µ. The operators sewed with these Beltrami differentials in (2.1) are part of a twisted N = 2 superconformal algebra, which we have briefly reviewed in appendix A. The free fermion ψ 3 , which is inserted at an arbitrary position α on the Riemann surface, serves the purpose of soaking up the holomorphic zero mode of the T 2 torus. Finally, Q 1,2 are the zero modes of the corresponding right-moving currentsJ 1,2 associated with the gauge group in the heterotic theory.
The correlators (2.1), however, are not yet the most generic objects one might envisage. Indeed, already in [12] a generalization was worked out by allowing also a dependence on non-holomorphic vector multiplets. Indeed, the following semi-topological expression was found in a g-loop heterotic string amplitude for n ≥ 0 (for details of the notation see [12]) where M (g,n) denotes the integral over the moduli space of a genus g Riemann surface with n punctures v j . Note that deformations of the latter are parametrized by 3g − 3 + n Beltrami differentials µ. Moreover, it was also noted that the hypermultiplets of the 'wrong' analyticity (i.e. different than the one induced by the F-term like measure in harmonic superspace) are associated with operators of the form (Ξ A ,Ξ A ). The latter are labeled by an Sp(2n) index A and transform as doublets under the SU(2) current algebra. Besides they are N = 4 primary operators in the sense that each of them is annihilated by half of the supercharges This suggests that the most generic amplitude, with a dependence on holomorphic and anti-holomorphic vector multiplets and hypermultiplets of all analyticities, should structurally be of the form 1 with g ≥ n ≥ m 1 and n > m 2 arbitrary integers. As we will see in the remainder of this work, this is essentially correct and we will prove it by showing that there is indeed a string amplitude, which (after some manipulations) gives rise to this expression. More specifically, we will calculate in the next subsection a heterotic g-loop amplitude with insertions of 2g chiral gauginos, 2n − m 1 − m 2 anti-chiral gauginos, two chiral and m 1 + m 2 + 2 anti-chiral hyperfermions and further discuss its relation to (2.4).

BPS-Saturated Amplitudes
Following the discussion in [12], we will compute the heterotic amplitude at a generic point in the moduli space of the K3 × T 2 compactification. For a brief review of the notation as well as for explicit expressions of the vertex operators see appendix C. Here, for convenience, we just give table 1 containing the fermion charges with respect to space-time and torus fermions (bosonized in terms of scalars φ 1 , φ 2 and φ 3 respectively) and the H-charge of the vertex operators for the physical fields as well as the picture changing operators (PCO). The last column RM denotes the right moving (bosonic) piece of the vertex operator. We have not displayed the superghost part in the table, however, all matter fields are understood to be inserted in the (−1/2) ghost picture, i.e. their vertices come with e −φ/2 , while the PCO come with a factor of e φ . The PCO's are arranged in such a way that g + m 1 + m 2 − n of them at r a ∈ {s 3 } contribute the torus part while 2g + 2n − m 1 − m 2 of them at r a ∈ {s H } the K3 part. In fact, by charge conservation, this is the only possible grouping of all PCOs -of course we still have to take into account all possible distributions of the total number 3g − 2 + n of the picture changing operators into these two classes. Notice moreover that we have put c-ghosts together with the vertex operators at the points v j and t k which means that the latter will not be integrated over the world-sheet. Instead, however, the genus g world-sheet will have n + 1 punctures.
With this table we can write down the amplitude in a straight-forward manner (for our notation concerning the spin structure dependence see appendix C.2). In order to save field pos.
number writing we define the following world-sheet positions with which the amplitude can be written in the following form For simplicity we will consider in the following differences of vector multiplet gauge groups in which case all of the operatorsJ I andJ J with holomorphic indices can only contribute zero modes. Since therefore their contribution is trivial, we will mostly drop them in the following. Furthermore, we are still free to choose a gauge and by picking the condition essentially both the ϑ-functions cancel in (2.5). For the remaining expression we can directly perform the spin-structure sum using (C.6) with the result that where, after using the gauge condition (2.6), we obtain Using the same manipulations as in [12] we can bring the amplitude into the following form We can further simplify this expression by collapsing two of the PCO-insertion points with z 1 and z 2 respectively. Just as in [12], this converts these insertions to (0)-picture hyperscalar vertex operators which we can write as covariant derivatives in the following manner operator dimension charge statistics Performing now also the integration over x i and y i gives a factor of (det(Imτ )) 2 which cancels the one already present. The left-over expression therefore becomes In a final step we can again follow [12] and commute the positions of the PCO to the Beltrami differentials yielding Here we have introduced numerical normalization factors which will turn out helpful in reducing unnecessary writing in some of the later computations. Since this topological expression has four different types of insertions we have compiled the dimensions and U(1) charges for all of them in table 2 for the reader's convenience. Notice that all insertions with dimension 1 are integrated over the world-sheet, while those with dimension 0 are inserted at fixed points associated to the n + 1 punctures. We also note that this is not quite the object anticipated in (2.4), however, we shall see in section 2.3 that indeed it is very closely related.

Notation and Simplifications
Due to the great number of insertions, (2.10) is unfortunately a rather complicated expression. We therefore would like to introduce certain simplifications (and shorthand notations) in order to save writing and moreover also make the computations more transparent. Let us first introduce the operator Although this operator does not have a well defined harmonic U(1)-charge, there is no actual inconsistency because of the following reason. Upon writing (2.10) with the help of (2.11) it is clear that in order to be able to provide the correct number of contractions withψ 3 (and to soak up the torus zero mode) from the 3g + n − 2 operators G − precisely (g + m 1 + m 2 − n) have to contribute the G − T 2 piece. All the remaining ones have to contribute the G − K3,+ part which makes sure that H has exactly the correct harmonic charge 2g + 2n − m 1 − m 2 − 2. Therefore, (2.11) provides us with a helpful shorthand notation without loosing any information about the amplitude (2.10).
Eq.(2.12) is in fact identical to the one appearing in [7] in the context of N = 1 heterotic topological amplitudes where there are n + 1 charge (−1) operators and n + 1 charge (+2) operators together with the required (3g + n − 2) G − insertions folded with the Beltrami differentials. These charge (−1) and (+2) operators correspond to N = 1 anti-chiral matter multiplets (with two different helicities). When the N = 1 theory is embedded inside an N = 2 theory, the N = 1 anti chiral multiplets can come from N = 2 vector multiplets or the hypermultiplets. In (2.12), the (n + 1) charge (−1) (or charge (+2)) operators are correspondingly split into m 1 + 1 (or m 2 + 1) that come from hypermultiplets and n − m 1 (or n − m 2 ) that come from vector multiplets. In other words, in (2.12), the specific features of the N = 2 theory have not been used apart from splitting the operators in terms of their vector or hyper origins. In the next subsection we will define certain reduced topological amplitudes that are obtained from (2.12) by using the special properties of the N = 2 theory namely the N = 4 world-sheet superconformal structure and as a result they will not have an N = 1 counterpart.
In [7], it was also also pointed out that this class of amplitudes extends to genus zero surfaces, where they correspond to Π n+1 terms. On genus zero, there must be a net charge +3 and three of the dimension zero operators are unintegrated which implies that the number of G − is n − 2. It is clear that also in the N = 2 case under consideration we will have the genus zero counterpart of (2.12) with the number of G − being n − 2.
To proceed further, it will turn out to be very useful to write (2.12) in a totally integrated form. To this end, for g > 1, we will localize the n + 1 operators G − , which are sewed with operator dimension charge statistics we can rewrite the topological amplitude (2.12) as (2.15) For convenience we have again listed dimensions and charges of all insertions of this expression in table 3. Notice that all operators have now dimension 1 and are integrated over the world-sheet as it is appropriate. Let us also remind again that the operator Φ B does not have well defined harmonic U(1) charge. However, just as before, after performing all contractions of the ψ 3 free fermions, only those contributions will survive in (2.15) which have a well defined charge.

Symmetry Properties of the Amplitude
For later use we will consider now the symmetry properties of the topological expression (2.10). First of all, we notice that the latter is anti-symmetric in each set of indices A, B, I andJ separately. This is manifest in the expression (2.10) for the indices A andĪ and follows from the fermionic properties of the operatorsΞ A andψ 3 respectively. As we can see from the fully integrated expression (2.15), in fact the same is also true for the indices B andJ as it is expected from the physical component amplitude which we started from in section 2.2.1.
Although not manifestly visible from the expression (2.12), the amplitudes H g also have symmetries upon exchanging, say,Ī withJ indices. To see this, we consider the expression (2.15) with integrated insertions and rewrite for example Deforming the contour integral the only non-vanishing residues are such that we immediately get We can rewrite the right hand side of this equation to obtain which is the direct generalization of identity (3.7) in [7].

Generating Functional
In view of the symmetry properties derived in the previous section, we can introduce a shorthand notation which will allow us to greatly simplify many of the following computations. Indeed, observing the particular fermionic nature of all insertions in table 3 it seems appropriate to introduce four types of Grassmann valued quantities (ϑĪ , ηJ, χ A , ξ B ) and to define the following generating functional -just like in [7] H g : In this form, every summand in H g can be thought of as a differential form with respect to (ϑĪ , ηJ, χ A , ξ B ) of degree (n − m 1 , n − m 2 , m 1 + 1, m 2 + 1). In order to also capture the particular case m 2 = 0, which will be discussed in more detail in section 2.4.1, let us also introduce a reduced version of H g . By restricting to the term m 2 = 0 in the summation of (2.22) we can define Notice that every summand of this expression can be interpreted as a differential form with respect to (ϑĪ , ηJ, χ A , ξ B ) of degree (n − m 1 , n, m 1 + 1, 1).
In order to demonstrate the usefulness of this notation let us just note that the symmetry property (2.20) can be written more compactly as Here we have used the fact that (2.20) is valid for every single summand in the sum over m 2 in (2.22). Therefore it follows that the equation holds for H g as well as for the reduced generating functional (2.23).

χ-Exactness and PotentialH
An observation which will turn out important for the remainder of this work is the fact that H g is in a certain sense an exact form with respect to χ A . A first hint towards this is to realize that H g is closed, i.e.
To see this, we consider the anti-symmetrized derivative Deforming the contour integral G + K3,+ , however, it produces no pole with any of the other operators and therefore vanishes. This proves relation (2.25), however, we can derive an even stronger statement by writing G − K3,+ = G + K3,+ J −− K3 in (2.10) and deforming the contour integration The insertion of G +

K3,+Ξ
A 1 corresponds to a derivative with respect to D A 1 + which we can pull out. Writing the remaining expression in terms of the G − we find Upon introducing the quantity (which is indeed exactly the object we had anticipated in we obtain the relation Comparing, however, to (2.10) we notice that the prefactor on the left hand side is identical to the harmonic charge of H (i.e. it counts the total number of G + K3,+ operators). We can thus replace it by where D 0 is given in (B.3). Introducing the generating functional , equation (2.31) can be more compactly written in the following form Notice that this relation is in a certain sense equivalent to stating that H g is an exact form as concerning the χ A variables. One should realize that similar to H g in (2.22), also every summand ofH g can be interpreted as a differential form in (ϑĪ, ηJ , χ A , ξ B ). The important difference, however, is that the corresponding degrees are (n − m 1 , n − m 2 , m 1 , m 2 + 1), i.e. the degree in χ A is shifted by one relative to (2.22).

ξ-Exactness and the Particular Case
The discussion of the previous sub-section focused entirely on the properties of H g as a function of the χ A variables. Since in the amplitude computation of section 2.2.1 the indices A and B (an therefore also the variables χ A and ξ B ) enter on the same footing, one would expect similar properties of H g also with respect to the ξ B variables. This is essentially the case, although, as we will discover now, there is a slight subtlety. Let us again begin by proving that H g is closed with respect to ξ B , i.e.
To this end, we consider the anti-symmetrized derivative Here we have inserted the vertex operator in its un-integrated form at an additional puncture of the genus g Riemann surface and added a further Beltrami differential. Writing now Ξ B 1 = ψ 3 (ψ 3 Ξ B 1 ) and deforming theψ 3 contour integral it produces no pole with any of the other insertions. Therefore this expression vanishes, which proves equation (2.33).

The Case m 2 = 0
In order to understand in how far H g is also an exact form with respect to ξ A (in the sense of the previous section), we will first discuss the particular case m 2 = 0, which has some interesting properties on its own. The corresponding H takes the form of In writing this expression we have used the fact that the only left-over ψ 3 (originally situated at t 1 ) is not able to contract with any other operator but is forced to soak the torus zero mode. It can therefore be inserted at an arbitrary position on the world-sheet which we call α. We can now proceed to transform the Ξ B insertion into its integrated form by localizing one of the Beltrami differentials The insertion of G − K3,+ Ξ B is, however, just a zero picture vertex operator, which we can pull out in the form of a moduli derivative We can also find a similar interpretation for equations (2.36) and (2.45) involving the reduced generating functional. Indeed, equation (2.36) can be rewritten in the following form with the newly introduced potential For completeness, let us also mention that H g red andĤ g red are of course still exact with respect to χ A in the following sense with the corresponding potentials which satisfies the following identity The degrees of every single summand in (2.40) are respectively We remark that (2.44) is precisely the topological expression conjectured in (2.4). Hence we have shown that there is also an actual string amplitude related to (2.44). This raises the question what the corresponding term in the string effective action looks like and whether it is possible to extract from it information about the 1/2-BPS nature of H along the lines of [12,8,11]. We will return to this question in section 3.

The Case m 2 > 0 and Gauge Freedom
It is important to point out that all manipulations of the previous section were only possible for m 2 = 0 since only then the string theory amplitude allows to put the ψ 3 insertion at some arbitrary position. In the case of m 2 > 0, there would be additional contractions which invalidate these steps. However, independent of its origin in string theory, even for m 2 > 0 we can consider the following object (2.49) Notice in particular that this quantity depends on the position α (which is a clear indication that it cannot come from a string theory computation). However, as we show in appendix D, changing the position α just amounts to adding an expression which is a derivative with where G g (α, β) is written in (D.4). Therefore, taking an anti-symmetrized derivative ofĤ we in fact obtain Notice that this expression is independent of α: According to (2.50) changing the position of α results in an extra contribution which lies in the kernel of ξ B D B + . Thus we can put α at the position of Ξ B 1 , such that which in the language of the generating functionals is written as Therefore, H g is also an exact form with respect to ξ, however, the corresponding potential has no nice interpretation within the framework of BPS-saturated string amplitudes. We note in passing that the subtle difference between χ A and ξ A is due to the particular choice of setup (and gauge fixing condition) in section 2.2.1, where we have computed the topological amplitude.

Harmonic description and Effective Action Coupling
In the previous Section we have considered particular amplitudes in heterotic string theory which are captured by correlation functions in a twisted two-dimensional theory. In this Section we would like to understand which terms in the heterotic effective action these amplitudes correspond to and whether they have any interesting properties with respect to their moduli dependence. It turns out that the effective action is best formulated in N = 2 harmonic superspace [13], for which we have reviewed our conventions and discussed the relevant Grassmann analytic on-shell superfields in appendix B.
Using the G-analytic superfields (B.7) and (B.16), we now want to construct effective action couplings which correspond to the superstring amplitudes computed in section 2.2.1. As it turns out, just as the amplitudes themselves were generalizations of the expressions in [12], also the corresponding effective action couplings will be a generalization of those discussed in [12]. The coupling discussed in the latter work was of the form where we have introduced the G-analytic superspace measure . . , n is an SO(n) vector index labeling the coordinates of the coset of physical scalars (for a discussion see [12]), while a = 1, 2 is an SU(2) index. For a closer description of these superfields see appendix B.2. Notice finally, the superfield K α − being a fermion, one needs a gauge group of sufficiently high rank, so that (K α − ǫ αβ K β − ) g = 0. We will assume that this requirement is met.
We can generalize (3.1) in various ways. The most general is to allow the coupling function to depend on holomorphic W I and anti-holomorphic vector multipletsWĪ as well as hypermultiplets of both analyticities i.e. (q + ,q − ) ↔ q + a and (q − ,q + ) ↔ q − a , a = 1, 2. To make this coupling manifestly G-analytic, however, we need to project it with four spinor derivatives where we have a great number of possibilities to distribute the derivatives, as we shall encounter in section 6. Note also that the U(1) charge of the function F has changed, to compensate for the extra charges brought along by the derivatives. In addition to (3.3) we can also put a number of fermionic factors such that the most generic coupling takes the form where d = 2g + 2n − m 1 − m 2 − 2 and 2 − d is the total harmonic charge of the function F . As usual, the operator (D − · D − )(D + ·D + ) is the projector on the G-analytic superspace, in order to match the properties of the measure. Moreover, the newly introducedΨ denotes the antichiral hyperino. Since we now have both types of hypermultiplets, G-analytic and anti-analytic, we can obtain the same hyperino in two ways, (both superfields start with the same spinor component). The labels (±) in (3.4) indicate the helicities of the corresponding spinors. In addition, in order to match the string theory computation of the previous section, we make a distinction between the labels A, B ofΨ and I,J ofK, in association with the helicity, to indicate their symmetry properties. The origin of this helicity splitting of the coupling (3.4) can be understood as follows. Let us start by denoting the two antichiral fermions collectively as Then we can write the expression in the square brackets in (3.4) in the form where is the Lorentz invariant contraction of two spinors. Here (±) denotes the two projections of the spinor indexα = ((+), (−)), or two helicities. Note that the expression in (3.7) must have a total harmonic charge −2g + 2 to match the charge of the first line in (3.4).
The coefficient function H in (3.7) has two obvious symmetries. One of them reflects the M ↔ N symmetry of the contraction in (3.8) (recall that Φ are fermions), the other corresponds to swapping two such contractions less obvious symmetry is due to the two-component spinor cyclic identity These symmetries do not make the rank (n − 1) tensor The reduction is achieved by breaking up the contractions (3.8) in (3.7) into pieces, first according to the two ingredientsΨ,K in (3.6), and then according to helicity. At the first step we get three types of contractions, where 2n − 2 is the total number of fermions and d = −2g + 2 − b − 2c is the harmonic charge of the function H. Assuming that a ≤ c, i.e. d ≤ 3 − 2g − n, and using the identity (3.9), we can rewrite all contractions of the first type as contractions of the second type. So, from now on we set a = 0, b + c = n − 1. Further, we split each contraction into helicity states according to (3.8), and expand the resulting binomials. This gives rise to a number of terms of the type in the square brackets in (3.4), In the process, each term of the new type picks a particular irreducible projection of the tensor H M 1 M 2 ...M n−2 M n−1 . As an example, consider the following simple case: The function H is symmetric under the exchange of pairs of indicesĪ A andJ B . This symmetry can be realized in two ways, H (ĪJ) and H [AB] [ĪJ] , each of which is associated with a particular helicity distribution: [ĪJ] . (3.12) 4 Differential Equations from String Theory for m 2 > 0 In this section we will derive differential equations for the amplitudes H with respect to the vector multiplet and hypermultiplet moduli. All the equations in this section will be true only modulo possible contact terms that could give rise to moduli space curvature terms. In fact we have not been able to find the structure of the contact terms in these equations that would satisfy the integrability conditions in the most generic setting. Violation of the integrability conditions therefore turns out to be proportional to curvature terms. Even so the fact that the equations are integrable in the flat limit is highly non-trivial and this will be proven in Section (4.3). However in Section 5 we will obtain differential equations for the reduced topological amplitudesĤ andĤ where we have been able to get the correct contact terms so that the resulting equations are exactly integrable even in the the case of general curved moduli spaces.
Based on the relations found in [12], we start out be applying the 'harmonicity operator' toH g . As we will see, this equation will lead us directly to further interesting equations. In particular, sinceH g is related to H g , application of a projected derivative χ A D A + will yield a (weaker) condition for H g . 3

Differential Equation forH g
We will start by applying the harmonicity operator to the 'χ-potential'H The covariant form of the harmonicity operator was given in [12] to be of the form only depends onū i + , we can also use the more practical expression involving -strictly speaking illegal -partial harmonic derivatives (4.2) Applying this operator toH g , we obtain in the language of the generating functionals the following equation (the explicit calculation is rather lengthy and technical and we have therefore relegated it to appendix E.2) where (as already mentioned in appendix E.2) we have converted the prefactor (2g + 2n − m 1 − m 2 − 2), which accompanies the boundary terms, into the action of the operator (−D 0 + 2). We should mention at this point, that this relation is correct up to possible contact terms which stem from two of the punctures coming close together. Since these terms are notoriously difficult to handle, they have not been considered. We note, however, for completeness, that these terms will correspond to curvature contributions and are a consequence of the fact that the moduli space of the heterotic compactification is not flat.
In fact, it will be much more convenient to slightly rewrite the harmonicity operator by moving the (projected) hyper-derivatives to directly hitH g Using relation (2.32), we may further write Notice that here (−D 0 + 2) acts on the whole equation. Since none of the topological amplitudes lies in the kernel of (−D 0 + 2) we can drop it and formulate the first order differential equation where for later convenience we have introduced the first order differential operator andC bdy as given in Appendix E is

Differential Equation for H g
After the discussion ofH g we will now also derive a differential equation for H g . We will show that combining the equations we have already obtained so far results in a new non-trivial differential equation for H g . In appendix E.3, we have shown that the same equation can also be obtained from a direct computation within string theory.
We will start from equation (4.6) and act with the differential χ A D A + on it. Inserting the explicit form of the boundary contribution (E.17) we get Our strategy is now to commute the χ C D C + through until it hits either H g orH g and then use either relation (2.25) or (2.32). Using the definition (4.7) of the operator D, we find (4.10) The commutator terms can be computed using (B.25) following the outline of quaternionic geometry in appendix B.3 where χ 2 := χ A χ B Ω AB and we have used that the contribution of the Sp(n) curvature is symmetric in (AC). Using moreover the identification Υ 0 = D 0 and Υ ++ = D + − , we obtain the simpler relation 4 Here we have also used that D + − H g = D + −Hg = 0, since H g andH g only depend onū i + . As we can see, D 0 becomes an overall operator which can be dropped, thus implying the relation We would like to point out that this is also a first order differential equation, just as (4.6). Even though we have neglected possible contact terms in (4.6), we have shown the additional curvature dependent term R in (4.13). This is because the derivation of the equation forĤ starting from that forĤ in Section 5, will follow the same steps as above and there we will indeed keep track of all the contact terms correctly.

Integrability
An important consistency check for our differential equations (in particular the boundary terms) is integrability. For simplicity we will consider here only equation (4.13) and check its self-consistency by computing D 2 H g . Since the square of the operator D is given by 14) and since in this Section we are ignoring curvature terms, the relation we have to prove is Ignoring the curvature terms we can further move D across various covariant derivatives which results in the following equation. (4.16) In order to further streamline this expression, we will use equation (4.13) Explicitly expanding this expression (and using anti-symmetry of Ω AB ) we find the following terms Indeed, upon reshuffling the sums over the genera, all the terms on the right hand side of this expression mutually cancel, which proves integrability of (4.13) (up to possible curvature terms).

Differential Equations from String Theory for m 2 = 0
Since the case m 2 = 0 is of particular interest, we will discuss the differential equations for H g red andĤ g red separately and show some interesting relations to the earlier work performed in [12]. Contrary to the previous Section, here we will keep all the curvature terms and show that the resulting equations are integrable.

Differential Equation forĤ g red
Generically, the case m 2 = 0 can be directly read off from the differential equations (E. 24) and (E.16). 5 For completeness, however, we will also work out the equations explicitly for the particular objectĤ g,A 1 ...Am I 1 ...Ī n−m 1 ,J 1 ...Jn defined in (2.44). Here, in order to save writing, we simply relabel m 1 by m. The explicit computation is very similar to the case m 2 > 0 and only some salient features are outlined in appendix E. 4. The final answer can again be written in the language of the generating functionals with the boundary contribution as given in (E.28). Following the same steps as in section 4.1, this equation can in fact be rewritten as which entails the first order differential equation

Differential Equation forĤ g red
As in the case m 2 = 0 the equations (2.41) and (5.3) put together, allow us to derive a further equation forĤ g red . To this end we act with χ C D C + on (5.3). Following the same procedures as in section 4 we obtain Using the fact thatĤ g−gs red is only a function ofū i + , we can simplify this expression which implies the differential equation

Comparison to Earlier Results
From the differential equations derived in the previous subsections we can make direct contact to similar equations worked out in [12]. The first equation derived in the latter work was a holomorphicity relation with respect to the vector-multiplet moduli. Translated into our language, this is the relation 6 ϑĪ∂ĪĤ g red = 0 . (5.5) 6 In fact in [12] it was erroneously stated that there would be additional boundary contributions on the right hand side of this equation, which would come from the exchange of vector multiplets in the limit when the genus g Riemann surface factorizes into a genus g − 1 surface and a torus. However, this contribution should be cancelled by a further exchange of a state of charge 3, which gives the same contribution but must appear with opposite sign. In fact we have checked that this (highly non-trivial) cancellation is indeed necessary to guarantee integrability (and thus consistency) of the differential equations derived in the previous section, as we shall see in the next subsection.

Indeed, this equation is precisely the term in (5.3) which is independent of χ
Here we have used thatĤ g red is at least linear in χ A and the only part of D which is independent of χ A is exactly the vector derivative, as in (5.5).
Moreover, a second order differential equation with respect to the hypermultiplet scalars has been derived in [12]. In preparation of reproducing this equation, we consider the commutator Singling out the term which is independent of ϑĪ we obtain which can equivalently be rewritten in the following form This is precisely the equation derived in [12], including also the boundary terms. Notice particularly the operator (D 0 − 2) in the first term on the right hand side. Translated to the setup in [12], this becomes 2(g − 1), which is exactly the prefactor found there.

Integrability
We would also like to briefly check integrability of the differential equations derived in the previous section. For simplicity we will focus on equation (5.4) for which we have to prove For the first term on the right hand side we can use (5.3), while the remaining terms will also contain commutator terms of covariant derivatives. The latter will result in curvature contributions, in particular (5.10) After some algebra we get Using at this stage thatĤ g red is independent of ξ A , we can streamline this condition in the following manner which is indeed satisfied, thus entailing that equation (5.4) is integrable. Notice that integrability holds also including the curvature contributions, which lends strong evidence to the fact that we have indeed obtained the full boundary contribution including the contact terms. We should point out that for m 2 = 0 using the relation ∂ ∂ξ AH | ξ=0 = D A +Ĥ we can obtain the differential equation for former in terms of the equation for the latter (5.4).The resulting equation has an additional curvature dependent term as compared to the one obtained from (4.6). Including this additional piece we can prove integrability for equations for H andH to the first order in ξ (i.e. for m 2 = 0) for a general curved moduli space. This indeed shows that some contact terms are missing in (4.6). While for m 2 = 0 this extra contact term can be deduced from (5.3), we have not been able to extend this to the general case m 2 = 0.

Field Theory Explanation of Differential Equations
After having derived several differential equations for the string theory amplitudes H g ,H g , H g red andĤ g red we would also like to obtain a better understanding from a field theoretic point of view. In the past [8,9,11,12] such equations could always be explained to root in particular analyticity properties of the effective action couplings with respect to moduli as well as harmonic coordinates. Put differently, these equations were a consequence of the fact that the couplings did not depend on all multiplets in a random manner. As explained in section 3, the situation now is slightly different in so far as, say, a generic H depends on all projections of the hypermultiplets and holomorphic and anti-holomorphic vector multiplets alike. Thus, all equations derived in the previous sections indeed will not reflect the analyticity properties of a single coupling, but rather encode relations between two or more couplings which are induced by supersymmetry. We will explain this idea in the following for equation (4.6) while a similar argumentation will also hold for the remaining equations.

Extracting the Component Coupling
As a starting point let us first translate equation (4.6) back from the language of generating functionals into the language of effective couplings. To this end we recall that the coefficient of (χ A ) m 1 +1 (ξ B ) m 2 +1 (ϑĪ) n−m 1 (ηJ ) n−m 2 of the left hand side of (4.6) reads We note in passing that this expression is consistent in so far as each of the terms has the same harmonic charge −(2g + 2n − m 1 − m 2 − 4). We would now like to prove that this expression is zero. 7 To this end we realize that ( is defined in terms of the 'master' coupling where an appropriate antisymmetrization of the indices is understood. Notice that the coupling H indeed has charge −2g − 2n + m 1 + m 2 + 2 as required. The twoH-couplings in (6.1) now correspond to slightly different couplings with a slightly modified distribution of the spinor derivatives. To be concrete we may for example consider the following distributions These couplings have charges −(2g + 2n − m 1 − m 2 − 4) and −(2g + 2n − m 1 − m 2 − 3) respectively, which is indeed as required. We should note that all expressions are still arbitrary functions with respect to the scalar moduli f A i and ϕ, as well as the harmonic variablesū i ± . Before proceeding, however, we notice an opportunity for simplifying our notation and saving a lot of writing: Since all superfields which appear in (6.2), (6.4) and (6.5) are fermionic in nature, all four sets of indices (A, B,Ī andJ) will be separately antisymmetrized. We can thus, just as in section 2.2.4 introduce Grassmann valued expressions (χ A , ξ B , ϑĪ, ηJ ) and formulate everything in terms of generating functionals H g andH g .

Harmonic Gauge Fixing
Our next step after preparing the component couplings is to perform the harmonic u-integral which is most easily done after fixing the harmonic gauge freedom. In order to do this, we proceed in a similar fashion as in [12] and reduce all couplings to SU(2) irreducible representations. The main novelty, however, is that, due to the dependence on the anti-G-analytic hypermultiplets q − ,q + , the couplings can have an arbitrary scalar dependence. Thus the generic gauge fixed form for the three couplings respectively will look like (as announced, we will be using the language of the generating functionals) where χ a 1 ϑ a 2 ξ a 3 η a 4 means extraction of the appropriate coefficient in the power series expansion in terms the Grassmann variables. The precise power also determines the harmonic charge, i.e. d = 2g + 2n − m 1 − m 2 − 2. Notice that after taking into account the definition in terms of the master couplings (6.6) and (6.7), we may immediately state the relation which is exactly equation (2.32) which we find in string theory. After these preparations we are now in a position to explain equation (4.6). To this end, we consider (the equation is understood to be the coefficient of χ m 1 +1 ϑ n−m 1 ξ m 2 +1 η n−m 2 in the Grassmann expansion) Using at this stage relation (6.11) as well as (6.6) and (6.7) we may rewrite In the last term, however, we recognize again (6.9), such that 14) which is exactly relation (6.1) that we wished to prove. Notice, if we treat the auxiliary fields properly (instead of putting them to zero in the on-shell approach we are using here), the right-hand side will be modified by additional boundary terms, which is indeed what we find from the explicit string computation in equation (4.6).
There is one further comment we would like to make: The crucial step in this derivation was the second line of (6.12). There we have used that a singleū i − combined with set of (totally symmetric)ū i + gives rise to two irreducible SU(2) tensor structures corresponding to the two terms in the square bracket of the second line in (6.12). Since both of them have different isospin, they cannot be related to each other and thus, from a harmonic-SU(2) point of view, equation (6.1) decomposes into two distinct equations. This in fact can also be directly verified in (4.6) by hitting it with the two projection operators D + − and −D 0 + 2 − D − + D + − respectively. In the first case we obtain which is just equation (2.32). Here we have used that D + − annihilates H g ,H g as well as every term contained inC bdy . Applying the operator −D 0 + 2 − D − + D + − , we in turn obtain which is just the original harmonicity equation (4.3). Thus the two SU(2) structures of (6.12) indeed encapsulate all the relevant information about the couplingsH g and H g .

Conclusions
In this work we have studied a generalization of the new class of topological amplitudes introduced in our previous work [12]. These amplitudes compute the coupling-coefficients of a very broad class of 1/2-BPS terms in the string effective action, which mix vector multiplets with neutral hypermultiplets in a non-trivial manner. In fact, we find a pair of (series of) amplitudes, which we combine into the generating functionals (H g ,H g ). These two objects are related through a hyper-moduli derivative andH g plays the role of 'gaugepotential' for H g . In particular the moduli dependence of these coupling functions satisfies a set of differential equations emerging from their 1/2-BPS structure, generalizing in a sense the holomorphic anomaly and harmonicity equations studied in [12,9]. In fact, the analytic projection which we use to obtain a consistent BPS action coupling requires non-trivial relations between different terms at the component level. The latter predict highly nontrivial differential equations between different topological amplitudes which we have directly checked in string theory. The main novelty is that they are first order differential equations (as compared to the second order relations derived in [12,9]), however, mixing H g with H g . As usual, we find that they get violated by world-sheet boundary contributions which we have evaluated explicitly. In addition, however, due to the fact that the hypermultiplet and vector multiplet moduli space of the N = 2 string compactification is not flat, we also obtain contact terms, which turn out to be very difficult to control for the generic amplitudes. However, we were able to identify a particular subset of the series of topological amplitudes (denoted (Ĥ g red ,Ĥ g red )) for which we managed to control the curvature contributions and prove integrability (and thus consistency) of our equations. We should also note that the equations close on the subset (Ĥ g red ,Ĥ g red ) in the sense that no new classes of topological objects are introduced. Thus one might hope, taking into account the first order nature of these equations as well as its iterative structure, that it is possible to iteratively solve for (Ĥ g red ,Ĥ g red ), using similar methods as in the case of pure vector multiplet dependence (see e.g. [3,4,5,6]). In this way one might be able to obtain further insight into the structure of the moduli spaces of N = 2 string compactifications. From a field theoretic point of view, the origin of the anomalous terms can probably be also understood via the process of integrating out auxiliary fields (for a discussion see [12]). It would be very interesting in the future to make this point more precise.

A Superconformal Algebras
For heterotic string theory compactified on K3 × T 2 the world-sheet theory is a product of an N = 2 and an N = 4 superconformal field theory representing T 2 and K3 respectively. In the following we will briefly review both theories mainly in order to fix our notation.

A.1 The N = 2 Superconformal Algebra
The (untwisted) N = 2 superconformal algebra (SCA) contains besides the energy momentum tensor T T 2 two supercurrents G ± T 2 which are positively and negatively charged with respect to a U(1) Kac-Moody current J T 2 . Here we have added the subscript T 2 to all operators in order to indicate that they represent the torus piece of the heterotic world-sheet theory. The conformal weights of these operators are given by order to indicate that these operators represent the K3-piece of the heterotic world-sheet theory. The conformal weights of these operators are The non-trivial OPEs of these objects are given by In view of the harmonic coordinates introduced in appendix B.1 we will rewrite the N = 4 superconformal algebra in an SU(2) covariant manner. To this end, we group the supercurrents into SU(2)-doublets in the following manner We are now also free to project the SU(2)-indices with harmonic variables to define the following operators

B N = 2 Harmonic Superspace
In this section of the appendix we discuss our conventions for the N = 2 harmonic superspace and introduce the Grassmann analytic superfields which will be needed throughout this paper. Our conventions are identical to [12].
The harmonic functions have harmonic expansions homogeneous under the action of the subgroup U(1). The harmonic expansions are organized in irreps of SU(2), keeping the balance of projected indices so that the overall U(1) charge is always the same. An example of a harmonic function is e.g.
The first component in this expansion is a doublet of SU(2). The higher components give rise to higher-dimensional irreps.
The harmonic derivatives can be viewed as the covariant derivatives on the harmonic coset SU(2)/U(1), or equivalently, as the generators of the algebra of SU(2) written in a U(1) basis. This means that they are invariant under the left action of the group SU(2), but covariant under the right action of the subgroup U(1). They can be split into generators of the subalgebra U(1): and of the coset: and The harmonic derivatives are differential operators preserving the defining algebraic constraints (B.1) and (B.2).
The derivative (B.3) acts homogeneously on the harmonic functions. For instance, the function f + (u) above has U(1) charge +1, hence The harmonic expansion of this function defines an infinitely reducible representation of SU (2). It can be made irreducible by requiring that the raising operator D − + annihilates the function: In other words, such a function is a highest-weight state of a doublet of SU(2). The irreducibility condition (B.6) is also called a condition for harmonic (H-) analyticity.

B.2 Grassmann Analytic On-shell Superfields
The introduction of harmonic variables allows us to define '1/2-BPS short' or Grassmann (G-) analytic superfields. 8 They depend only on half of the Grassmann variables which can be chosen to be θ + α = θ i α u + i andθα − =ū i −θα i with α = (±) (α = (±)) an (anti-)chiral spinor index. In this work we will encounter two different examples which we will briefly review now (for more information see e.g. [12]).

B.2.1 Linearized On-shell Hypermultiplet
The first superfield which we will discuss is the linearized on-shell hypermultiplet (N = 2 matter multiplet) Here f i are the two complex scalars, χ α andψα are the two fermions of the on-shell multiplet. To exhibit manifest G-analyticity, one has to choose the appropriate analytic basis in superspace, analogous to the familiar chiral basis. Then q + satisfies the conditions Note that the harmonic dependence in (B.7) is cut down to linear. This is typical for on-shell multiplets which, in addition to the G-analyticity condition (B.9), also satisfy the H-analyticity condition Here the harmonic derivative is supersymmetrized by going to the manifestly G-analytic superspace coordinates (B.8). One can show that the 'ultrashort' on-shell superfield (B.7) is a solution to the simultaneous conditions for G-and H-analyticity [13,15,16].
Note that in the N = 2 G-analytic superspace there exists a special conjugation combining complex conjugation with a reflection on the harmonic coset, such that G-analyticity is preserved. In this sense we can define the conjugate hypermultiplet In what follows it will be convenient to combine the two versions of the hypermultiplet into a doublet of an external SU(2) (not the R symmetry one), (q + ,q − ) ↔ q + a , a = 1, 2.

B.2.2 Linearized On-shell Vector Multiplet
Another example of a G-analytic superfield is the linearized on-shell vector multiplet. It is obtained from the off-shell chiral field strength Here ϕ is the complex physical scalar and F (+) is the self-dual part of the gluon field strengths, while S is a triplet of auxiliary fields. On-shell the latter must vanish, hence the additional constraint Now, define the superfield (a superdescendant of W ) where we have projected the SU(2) index of D α i with the harmonicū i − . This superfield is annihilated by half of the spinor derivatives and hence is 1/2 BPS short. Indeed, this is true for the projectionsD + β since {D + , D − } = 0 andD iβ W = 0 (chirality). Further, hitting (B.14) with D β − we obtain zero as a consequence of the projection of the on-shell constraint (B.13) withū i −ū j − . We conclude that K α − satisfies the same G-analyticity constraints as the on-shell hypermultiplet (B.9) which imply that K − depends only on half of the θ's: In addition, the harmonic dependence of K α − is restricted to be linear. As in (B.10), this follows from the condition for H-analyticity in turn derived from the harmonic independence of W (D − + W = 0) and the commutator [D − + , D α − ] = 0. This is another example of an ultrashort superfield. Note, however, that it is not a primary object but rather a superdescendant of the chiral on-shell vector multiplet.

B.3 Quaternionic Geometry
For completeness, we will review in this appendix some relevant aspects of quaternionic manifolds. We will closely follow the discussion in [17] and use the same notation and conventions as in [12]. Following [18,19,20], a quaternionic manifold M is defined as a 4n dimensional Riemann manifold whose holonomy group is restricted to a subgroup of Sp(n) × Sp(1). Therefore, we can locally choose a coordinate frame {x M k } where M = 1, . . . , 2n and k = 1, 2 are indices of Sp(n) and Sp(1) ∼ SU(2) respectively. In this coordinate frame, the covariant derivatives can be written as where we have introduced the Sp(n) and Sp(1) connections ω Ai(CD) and ω Ai(lk) respectively. Our convention for the Sp(n) and Sp(1) algebra of the generators Γ (CD) and Γ (lk) follows [17] and is given by where Ω AB denotes the Sp(n) invariant tensor. In this paper we will only consider the fundamental spinor representations of Sp(n) and Sp(1), in which case Particularly, we will be interested in the commutator of two covariant derivatives leading to the Sp(n) and Sp(1) components of the curvature tensor According to [17], restricting the holonomy group of M to (a subgroup of) Sp(n) × Sp(1) implies the covariant constraints 9 Hence we have for the commutator where we have introduced as the Sp(1) generators analog to (B.3) and (B.4). Projecting with the corresponding harmonic variables, we particularly obtain where we have introduced the projected quantities In this appendix we recall some important features of heterotic K3 × T 2 compactifications, in particular expressions for various vertex operators and picture changing operators as well as our conventions for the spin-structure sum.

C.1 Vertex Operators
The space-time and T 2 torus piece of all vertex operators is fairly standard and will be written using the free bosons φ 1,2 and φ 3 respectively (for more information see e.g. [12]). For the K3 piece we recall that the SU(2) current algebra inside the N = 4 superconformal algebra can be bosonized in terms of a free boson which we call H: HereĜ K3,i are dimension 5/4 operators which have non-singular OPE with H and have no spin structure dependence. The latter enters through the projections and shifts in the U(1) charge lattice of J K3 which in turn is given by the momentum lattice of H. Therefore in the N = 4 internal theory, only correlation functions and the partition function containing H depend on the spin-structure. The term in the picture changing operator containing a N = 4 superconformal generator is where dots indicate the remaining terms and φ bosonizes the superghost.
Concerning the vertex operators of physical fields, the chiral vector multiplet scalar in the (−1)-ghost picture is simply ψ 3 and in particular does not depend on the K3 fields. The gaugino vertex however involves, besides the space-time and torus spin-field, also e Thus the vertex operator for the gauginos λ ∓ andλ ± in the (−1/2)-picture also carry e Finally the vertex operators (at zero momentum) for the hyperscalar f ± A in the (−1) and (0) ghost pictures are given by The vertex operators in the (−1/2) ghost picture for the hyperfermions χ A andψ A are where S α are the space-time spin fields.
As mentioned above the spin structure dependence enters only through the superghost, spacetime and torus fermions and the charge lattice of H. It does not depend in particular on the rest of the details of the N = 4 superconformal theory. On the other hand, the topological theory (besides shifting the dimensions of the torus fermion) involves precisely twisting by adding an appropriate background charge for the field H and the rest of the internal N = 4 theory is insensitive to this twisting.

C.2 Spin Structure
Let Γ be the U(1) lattice of H charges. The space-time and torus fermions define an SO(2)× SO(2)×SO(2) lattice. If one takes an SO(2)×SO(2) sublattice thereof and combines it with Γ, then (see [22,23,24]), the resulting 3-dimensional lattice is given by the coset E 7 /SO (8). The characters are given by the branching functions F Λ,s (τ ) and satisfy: where χ Λ and χ s are the E 7 and SO(8) level one characters respectively, Λ denotes the two conjugacy classes of E 7 and s represent the four conjugacy classes of SO(8) in the spin structure basis. The characters of the internal N = 4 superconformal field theory times two free complex fermions can therefore be expressed as Λ F Λ,s (τ )Ch Λ (τ ) where Ch Λ (τ ) is the contribution of the rest of the internal theory and most importantly does not depend on the spin-structure. The generalization to higher genus is obtained by assigning an E 7 conjugacy class Λ for each loop and we will denote this collection by {Λ}. We can define a more general character F {Λ},s (a 1 , a 2 , a 3 ) by introducing chemical potentials for the three charges; a 1 and a 2 coupling to the two SO(2) charges and a 3 to H-charge. For a genus g surface, the couplings a 1 , a 2 and a 3 each are g-dimensional vectors and represent the coupling to charges going through each loop. In the calculation of the amplitudes (a 1 , a 2 , a 3 ) are related to the positions of various vertex operators weighted by the corresponding charges via Abel map. The spin structure sum is given by the formula: where we have introduced where λ i (with i = 1, ..., g) are 0 and 1 for the E 7 conjugacy classes 1 and 56 respectively. In fact, apart from the non-zero mode determinant of a scalar, Θ is just the character valued genus g partition function of level one SU(2), with the two classes above corresponding to the two representations of level one SU(2) Kac-Moody algebra based on SU(2) representations 1 and 2.

D Gauge Freedom
To see the gauge freedom of the potential forĤ we consider the (α-dependent) potential functionĤ This indeed shows that with the quantity For further convenience and in order to save writing, we will also introduce a generating functional for G

String Derivation of Differential Equations
Our strategy in deriving differential equations with respect to N = 2 moduli directly within the setup of topological string amplitudes follows closely [3]: Moduli derivatives are directly translated into vertex operator insertions within the twisted correlators. In particular we will consider derivatives with respect to vector multiplet and hypermultiplet moduli, for which the relevant insertions take the form We will then further deform the contour integrals (within the topologically twisted setup this is indeed possible on dimensional grounds) thereby trying to manipulate the correlation functions. A generic feature in this type of computations is the appearance of world-sheet boundary terms. The latter appear when the above mentioned contour integrals act on the measure of the twisted correlator, where they might give rise to insertions of the total energy momentum tensor T = T T 2 + T K3 sewed with some Beltrami differential that parametrizes a particular deformation of the genus g-world-sheet. The full energy-momentum tensor can in fact be written as a total derivative with respect to this particular modulus of the genus g-Riemann surface, thus leading to a boundary contribution in the integral over M g . These boundary contributions can pictorially be quite directly figured as degeneration limits of the genus g Riemann surface and (see e.g. [7,25,26,8,11]) there are two homologically distinct types corresponding to the pinching of either a dividing geodesic or a handle. Examples for genus g = 3 are depicted in Figure 1 and 2 respectively. Assuming that we are at a generic  point in the vector multiplet moduli-space with no charged massless states present there will be no contribution from pinching a handle. It therefore remains to study the pinching of a dividing geodesic where the genus g surface splits into two surfaces of genus g s and g − g s respectively 10 and to determine which states may propagate on the long and thin tube. In principle there are two different types, namely hyper-and vector multiplet states. We will discuss both in the following.
In our particular case, there is in fact a further type of degeneration besides these purely geometric ones. Indeed, since our world-sheets have in addition an arbitrary number of punctures (and the moduli space is thus rather M (g,n) ) we also need to take into account contributions when two of these punctures collide with each other. Such contact terms are notoriously difficult to handle within the string approach. In our differential equations they will manifest themselves as curvature contributions and take into account the fact that the moduli space of the string compactification is in general not flat. We will not be able to say much about these terms in the general case, but we will propose consistent expressions for the case m 2 = 0 as we go along.

E.2 String Differential Equation forH
This appendix is devoted to the derivation of a differential equation forH g,A 1 ...Am 1 ,B 1 ...B m 2 +1 I 1 ...Ī n−m 1 ,J 1 ...J n−m 2 along the lines of section 4.1. The derivation will be divided into first considering the bulk contributions and afterwards the boundary terms, which appear from a particular deformation of the genus g world-sheet.

E.2.1 Bulk Equation
We start by explicitly applying the operator (4.2) to the potentialH given in (2.29) where for convenience we introduce the shorthand notations in order to save writing. Deforming the contour integral, we find two different terms The second term in fact vanishes due to the anti-symmetrization of the Beltrami differentials and the fermionic nature of the G − K3,i operators. For the remaining term we can perform the ū i + differential and furthermore use T = T T 2 + T K3 . Denoting the boundary contribution as C (bdy),A 1 ...Am 1 ,B 1 ...B m 2 +1 I 1 ...Ī n−m 1 ,J 1 ...J n−m 2 we obtain Deforming the contour integration, we obtain The insertion G + T 2ψ3JĪ n−m 1 is in fact a derivative with respect to DĪ n−m 1 . Pulling it out and comparing with the definition (2.29) we find where we have explicitly pulled out the factor of (2g + 2n − m 1 − m 2 − 2) from the boundary contribution. In the language of the generating functional, this latter factor is converted into the action of the operator (−D 0 + 2) such that we can write the more compact relation It still remains to determine the boundary contributionC bdy .

E.2.2 Boundary Contributions
The total contribution from the boundary deformations can be divided into two terms corresponding to either hypermultiplet or vector multiplet states propagating on the long thin tube when pinching a dividing geodesic. We can calculate both contributions in a straight forward manner.

Hypermultiplet Propagation
There is not a unique way to distribute the operator insertions on the two Riemann surfaces but rather we have to consider several rather complicated contributions Fortunately, this expression can be written in terms of topological amplitudes in the following manner  . To this end we start by taking a derivative with respect to an anti-holomorphic vector multiplet scalar where we have introduced the operator Notice that similar to G − also this operator is strictly speaking not well defined regarding its harmonic charge. However, just as before the necessity of soaking up all torus fermionic zero modes ψ 3 takes care of this ambiguity in a well defined way. Deforming now the G + -contour we obtain Here the propagation of the vector multiplet states is in fact cancelled against the contribution stemming from the exchange of a state ρ of charge +3 and the unit operator. As already remarked in footnote 6 , this subtle cancellation was overlooked in the equations derived in [12], but is strictly necessary for the integrability of the final differential equation. For later convenience, we introduce a generating functional for the quantity (E.27)