Finite deformations from a heterotic superpotential: holomorphic Chern-Simons and an L∞ algebra

We consider finite deformations of the Hull-Strominger system. Starting from the heterotic superpotential, we identify complex coordinates on the off-shell parameter space. Expanding the superpotential around a supersymmetric vacuum leads to a thirdorder Maurer-Cartan equation that controls the moduli. The resulting complex effective action generalises that of both Kodaira-Spencer and holomorphic Chern-Simons theory. The supersymmetric locus of this action is described by an L3 algebra.


Introduction
A full understanding of the parameter space of string theory is an outstanding mathematical challenge and would lead to powerful constraints on the landscape of string models.
Of the various limits of string theory, the heterotic string has been the focus of much phenomenology thanks to the relative ease with which one can engineer four-dimensional theories with chiral fermions and the Standard Model gauge group [1][2][3][4][5][6][7][8]. Much of this work has been on models where the internal manifold is Calabi-Yau, mostly because such spaces can be constructed using algebraic geometry and then used for compactifications without knowledge of their explicit metrics. Calabi-Yau compactifications are not the most general way to obtain an N = 1 theory in four dimensions that admits a Minkowski vacuum. The general solution to O(α ) is given by compactifying on a complex three-fold X with H flux and a gauge bundle that satisfy an anomaly cancellation condition. The conditions on the geometry and fluxes for such a solution are known as the Hull-Strominger system [9,10]. Known solutions to this system include Calabi-Yau spaces with bundles and a small number of honestly non-Kähler geometries. Generically, a given solution of the Hull-Strominger system will admit deformations of the geometry, flux and bundle that remain N = 1 solutionsthese deformations are known as moduli. These moduli appear in the massless spectrum of the low-energy theory, so it is important that we understand the moduli space of a given compactification.
The moduli spaces of Calabi-Yau compactifications at zeroth order in α are well understood using the language of special geometry. Until recently the general case had not been tackled -this might come as a surprise. Certainly in type II theories the conditions for an N = 1 Minkowski solution are sufficiently complicated (thanks to branes and other ingredients) that their moduli spaces might not admit a general formulation. In the heterotic case, the underlying geometry is relatively straightforward. One might have expected that the gauge sector and anomaly conditions complicate matters somewhat, but that the moduli space might still be understood. Starting with [5,[11][12][13][14][15], this gap is now being filled. (See also [16][17][18][19] for worldsheet approaches.) Infinitesimally, the moduli space is characterised by the existence of a holomorphic structure D on a bundle Q over the three-fold X. The fact that the moduli space is finite dimensional is intimately connected to this holomorphic structure and the Bianchi identity JHEP10(2018)179 for the flux. The infinitesimal moduli are captured by the cohomology H (0,1) D (Q), where Q is defined by a series of extensions. In this way, the complex structure, hermitian and bundle moduli are combined in a single structure. Furthermore, one can define the analogue of special geometry for these general heterotic compactifications and find the metric on the moduli space [20,21].
A natural question to ask is whether one can understand the moduli spaces to higher order. If we think about deformations of a complex structure, we know the infinitesimal moduli are given by H (0,1) ∂ (T (1,0) X), while the higher-order deformations satisfy a Maurer-Cartan equation. A similar thing happens for bundle deformations [22], and simultaneously deformations of the bundle and complex structure [23,24]. In this way, moduli can be obstructed at higher orders and can give non-zero contributions to the superpotential of the four-dimensional effective theory. The aim of this work is to derive the corresponding conditions on the moduli for the Hull-Strominger system at higher orders. In other words, we want to derive the conditions on the moduli when they describe a small but finite deformation of the original heterotic solution.
There are a number of ways one might go about this. One path would be to start with the equations of the Hull-Strominger system and deform the various fields. The deformed fields should still satisfy the Hull-Strominger system (as it describes the most general solution) so one can rewrite the system of equations as conditions on the deformations themselves. This is similar to the path taken in [20,21]. Our approach will be complimentary. It has been shown that supersymmetry of the heterotic system can be described using a four-dimensional superpotential [21,25,26]. The vanishing of the superpotential and its first derivative imposes the F -term conditions in the four-dimensional theory and leads to an N = 1 Minkowski vacuum. Our plan is to deform the fields that appear in the superpotential and then read off the conditions on the moduli for the superpotential and its derivative to vanish. These two approaches will be shown to be equivalent in a future publication [27]. A particular advantage of proceeding this way is that one can use the knowledge that the superpotential is a holomorphic function of the moduli fields to streamline the problem.
In addition to the usual N = 1 lore that the superpotential is holomorphic, we give an argument that the superpotential is holomorphic on the space of moduli fields without requiring that they give a solution to the Hull-Strominger system. This is equivalent to saying that the off-shell parameter space -the space of SU(3) structures, B fields and gauge bundles -is a complex space and the superpotential is a holomorphic function of these parameters. We outline how this follows from generalised geometry where N = 1 NS-NS compactifications are described by a generalised SU(3) × SU(4) structure [28]. We identify the invariant (holomorphic) object which characterises this structure and find that the complex coordinates on the space of structures match with the usual complex structure, complexified hermitian and bundle deformations.
We show that the conditions on the moduli fields from the vanishing of the superpotential and its first derivative can be written as a pair of third-order Maurer-Cartan equations using the holomorphic structure and a number of brackets. Moreover, we show that the JHEP10(2018)179 superpotential itself can be rewritten using these operators in a Chern-Simons-like form: ∆W = X y, Dy − 1 3 [y, y] − ∂b ∧ Ω, (1.1) where y describes the complex, hermitian and bundle moduli, and b is a (0, 2)-form. ∆W is written using a holomorphic structure D, a pairing ·, · on the moduli fields, and a bracket [·, ·]. The brackets can be understood as coming from an underlying holomorphic Courant algebroid that describes the combined deformations of the complex structure, metric and fluxes [29,30]. We then show that the supersymmetry conditions can be recast in terms of an L 3 algebra. We outline how the L 3 algebra gives a C ∞ resolution of the underlying holomorphic Courant algebroid. The natural L 3 field equation reproduces the supersymmetry conditions, and the L 3 structure gives the gauge symmetries of the moduli space in a compact form.
It is known that generic deformation problems have a description in terms of L ∞ structures, so it is not unexpected that our moduli fields are governed by one. What is unexpected is that the structure truncates at finite order leaving us with an L 3 algebra. Why does the deformation truncate in our case? A generic deformation problem can be parametrised in many equivalent ways -some may truncate at finite order while others do not. Essentially, the structure of the heterotic system and its formulation using a superpotential guides us to pick a "nice" parametrisation. Said another way, we know from supergravity that the superpotential should be a holomorphic function of the parameters. Thus when we express the superpotential in the obvious complex coordinates on the parameter space, we get the most natural way to package the deformation problem.
We begin in section 2 with a review of the Hull-Strominger system and the description of its infinitesimal moduli in terms of a holomorphic structure as in [14]. In section 3 we discuss the off-shell parameter space of the theory and give the complex coordinates on the parameter space. We show how the F -term conditions follow from a heterotic superpotential to set the scene for the higher-order deformations. In section 4 we examine the higher-order deformation problem and find the system of equations that govern the moduli of the Hull-Strominger system. We show how this can be written in terms of the holomorphic structure D and a bracket [·, ·] arising from a holomorphic Courant algebroid. In section 5 we rewrite the equations that govern the moduli in terms of an L 3 structure. We give the various multilinear products k that define the L 3 structure and discuss how various properties, such as the moduli equations and gauge symmetries, are naturally encoded in this L 3 language. In section 6 we discuss how the system simplifies under various assumptions and comment on how the effective field theory is encoded in our language. We finish with a discussion of some open questions and avenues for future work.
In the appendices, we lay out our conventions, include a few comments on how flux quantisation works in the heterotic theory, discuss the off-shell parameter space in terms of generalised geometry, show that the D-term conditions do not affect the moduli problem and review how the massless moduli are captured by the the cohomology of the holomorphic structure.

JHEP10(2018)179
2 The Hull-Strominger system and a heterotic superpotential We begin with a review of the Hull-Strominger system [9,10] and the description of its infinitesimal moduli using a holomorphic structure [14,15,31].

N = 1 heterotic vacua and the Hull-Strominger system
The Hull-Strominger system is a set of equations whose solutions are supersymmetric Minkowski vacua of heterotic string theory to order O(α ). The ten-dimensional solution is a product of four-dimensional Minkowski space with a six-dimensional complex manifold X. X admits a vector bundle V with connection A whose curvature F is valued in End V . The tangent bundle T X of X also admits a connection Θ whose curvature R is valued in End T X. X admits an SU(3) structure defined by a nowhere vanishing spinor η or, equivalently, a non-degenerate two-form ω and a nowhere vanishing three-form Ψ that are compatible The invariant objects are defined by bilinears of the spinor as where we are free to normalise the spinor so that Ψ 2 = 8. In what follows it will be useful to define a three-form Ω which is related to Ψ by a dilaton factor as Supersymmetry of the vacuum follows from the vanishing of the supersymmetry variations of the fermionic fields, given in equations (A.7)-(A.9). To first order in α , these conditions are equivalent to the Hull-Strominger system: where ω CS is the Chern-Simons three-form for the connection, The closure of Ω from (2.4) implies that the manifold is complex with a holomorphically trivial canonical bundle, while condition (2.8) tells us that X is conformally balanced.

JHEP10(2018)179
Conditions (2.6) and (2.7) mean V is a polystable holomorphic bundle 1 so the curvature satisfies the hermitian Yang-Mills equations. Finally, (2.5) defines the flux H in terms of the heterotic B field and the anomaly cancellation condition, and links it with the intrinsic torsion of the SU(3) structure. The corresponding Bianchi identity is 2 This set of conditions defines what one might call a heterotic SU(3) structure. Upon considering the four-dimensional N = 1 theory that would follow from compactifying on such a solution, the Hull-Strominger system naturally splits into F -and D-term conditions. As discussed in [26], the F -term equations are It is these equations that the heterotic superpotential reproduces. The remaining equations of the Hull-Strominger system are the conformally balanced condition and the Yang-Mills equations, referred to as the D-term equations.
Modulo certain mild assumptions on the geometry, the infinitesimal deformations are parametrised by the cohomology H (0,1) D (Q), where D and Q are to be defined below. This cohomology is reviewed in appendix E. Under infinitesimal deformations, the D-term equations fix a representative of a certain cohomology class [14], and so should be thought of as gauge fixing conditions that do not affect the moduli problem. This is of course expected from the four-dimensional N = 1 supergravity point of view [33][34][35]. In appendix D we show that preserving the D-term conditions for finite deformations also amounts to fixing a gauge. One might worry about Fayet-Iliopoulos terms appearing, but these are in fact accounted for by modding out by D-exact terms, as shown in [14].

The Atiyah algebroid and a holomorphic structure
The vector bundle V is hermitian in agreement with (0, 2) supersymmetry on the worldsheet [36]. The curvature F of the bundle is given by where A is a one-form connection valued in End V . The exterior derivative on V twisted by A is 1 More precisely, it is the complex vector bundle V C (defined in appendix C.3) that is a holomorphic bundle. 2 The curvature R in the Bianchi identity is the curvature of a connection on T X, satisfying its own hermitian Yang-Mills conditions in order for the equations of motion to be fulfilled [32]. To O(α ), this connection is ∇ − , given by taking the connection in (A.7) with the opposite sign for H.

JHEP10(2018)179
where the action of the bracket on a p-form β is (2.14) A holomorphic structure on V is fixed by the (0, 1) component of d A , which we denote ∂ A . This operator squares to zero if the bundle is holomorphic, that is F (0,2) = 0. Moreover, the Bianchi identity for the curvature is simply A deformation of the Hull-Strominger system corresponds to simultaneous deformations of the complex structure, hermitian structure and gauge bundle. Taking each of these in isolation is not sufficient. In particular, deformations of the hermitian structure alone lead to an infinite-dimensional moduli space. It is surprising that if one considers the full deformation problem together with the anomaly cancellation condition, one finds a finite-dimensional moduli space. Of course, this is what one would expect from string theory, but the precise way in which this happens is rather remarkable.
As discussed in [14], the infinitesimal moduli of the Hull-Strominger system are captured by deformations of a holomorphic structure. The holomorphic structure D acts on a bundle Q. Locally Q is given by Globally, Q is defined by an extension 3 where the bundle Q 1 is defined by The holomorphic structure D on Q is a derivative 4 D : Ω (0,p) (Q) → Ω (0,p+1) (Q), (2.19) where D 2 = 0 if and only if the Bianchi identities for H, F and R are satisfied. The Hull-Strominger system is then equivalent to the data of the extension bundle Q, the nilpotent holomorphic structure D, polystability of V and T X and the conformally balanced condition on X.
The infinitesimal deformations of the holomorphic structure are simply elements of the D-cohomology of Q-valued (0, 1)-forms -H (0,1) D (Q). As shown in [14], this is also the moduli space of heterotic SU(3) structures. We give a short review of this in appendix E.
For the rest of the paper, we make a field redefinition to absorb the explicit α dependence (2.20) 3 Full details can be found in [14]. 4 A similar operator has appeared in the context of generalised Kähler geometry [29].

JHEP10(2018)179
One can restore the proper factors of α by the inverse transformations. We will also suppress the connection on T X -we can reintroduce it in what follows by treating T X as part of the gauge bundle and defining the bundle metric on the T X subspace to be negative definite so that the Bianchi identity comes with a negative sign for the tr R ∧ R term.
The main aim of this work is to understand what happens for finite deformations. In particular we will see the holomorphic structure is an important ingredient in describing higher-order deformations. First let us discuss the off-shell parameter space and how the heterotic SU(3) structure can be rephrased using a superpotential.

The off-shell parameter space
We now show that a subset of the system corresponding to F -term conditions can be derived from a superpotential. We then discuss how deformations of the geometry, flux and bundle can be parametrised using the observation that the superpotential is a holomorphic function.
The four-dimensional effective theory that one finds after compactifying the heterotic string on an SU(3) structure manifold is controlled by a superpotential [21,25,26]. The superpotential W is given in terms of the flux H and the SU(3) invariant forms by 5 As we will review, given the SU(3) structure relations (2.1) and the definition of H in (2.5), W = δW = 0 reproduces the F -term conditions of the Hull-Strominger system [26]. Notice that δW = 0 requires us to vary the superpotential over some space of field configurations. We need to understand what this space is in order to find how the superpotential behaves when we perform a finite deformation of the background fields. We pause briefly to distinguish between this parameter space and the moduli space of solutions to the Hull-Strominger system.
The parameter space or space of field configurations Z is the space of SU(3) structures, B fields and hermitian gauge bundles on the real manifold X. The SU(3) structure is equivalent to the existence of a nowhere-vanishing spinor so that on this space of field configurations the heterotic theory admits an "off-shell" N = 1 supersymmetry of the kind discussed in [37,38]. These fields do not necessarily solve the Hull-Strominger system and so we often refer to them as off-shell field configurations. This is the space over which the superpotential is varied.
The moduli space M of the Hull-Strominger system is a subspace of Z on which the fields also solve the Hull-Strominger system. This set of fields is what one usually means by moduli, and we will often refer to them as on-shell configurations. Another way of saying this is that the superpotential and its derivatives vanish when evaluated on M.
As W is a superpotential for a four-dimensional N = 1 theory we expect it to be a holomorphic function. The holomorphicity of W is a powerful tool for understanding deformations of the Hull-Strominger system and we will see later how its presence greatly 5 Here we have scaled away an overall factor of α /4 that comes from the field redefinition in (2.20).

JHEP10(2018)179
simplifies the problem. On the supersymmetric locus (on-shell), it is known that the superpotential is a holomorphic function of the moduli fields [21] -this is simply the statement that anti-holomorphic derivatives of the superpotential vanish on imposing the F -term conditions. Physics goes further than this and insists that W is a holomorphic function of the off-shell field configurations -the off-shell field space Z must admit complex coordinates and W must be a holomorphic function of these coordinates. In other words, the three-form Ω and the particular combinations of B + i ω and A that appear in (3.1) must be parametrised by these complex coordinates.
We outline a proof that such complex coordinates exist and that Ω is holomorphic on the parameter space Z in appendix C using the formalism of generalised geometry. We also discuss how the hermitian structure on V survives off-shell. For completeness, one should really show that W itself can be expressed as a holomorphic function of the object χ that we define in appendix C -we leave this for a future work.
Note that on-shell, Ω is also holomorphic as a function of the complex coordinates of X. When we talk of Ω being holomorphic we are instead referring to its dependence on the coordinates of the off-shell parameter space.
Let t and t denote holomorphic and anti-holomorphic coordinates on the parameter space Z. The corresponding holomorphic and anti-holomorphic variations are where D is a covariant derivative on the parameter space [20]. As we discuss in appendix C.3, off-shell the gauge bundle admits a real hermitian connection valued in V . This decomposes into (1, 0)-and (0, 1)-forms, with a corresponding decomposition of the Chern-Simons form. Not all components of the Chern-Simons form appear in W (as it is wedged with Ω); only the (0, 1)-form components of A contribute. It is this component of the connection that is the complex coordinate on the off-shell parameter space. A holomorphic deformation of the connection is then given by a (0, 1)-form valued in End V , which we denote α: We show in appendix C.3 that one can write the Chern-Simons form using connections valued in V , V C , V C or a combination of these. From this it is clear that it is equivalent to work with (0, 1)-forms valued in End V or with (0, 1)-forms valued in End V C . A deformation of the complex structure is parametrised by a (1, 0)-vector valued (0, 1)form, µ ∈ Ω (0,1) (T (1,0) X), also known as a Beltrami differential. The complex coordinates on X that define the complex structure deform to dz a → dz a + ∆dz a = dz a + µ a . (3.4) Infinitesimally, the deformed complex structure is J + 2i µ. For a small but finite deformation, the holomorphic three-form becomes [39,40] Ω → Ω + ∆Ω = Ω + ı µ Ω + 1 2 ı µ ı µ Ω + 1 3! ı µ ı µ ı µ Ω, (3.5)

JHEP10(2018)179
where the variations of Ω in coordinates are (3.6) The fact that a variation of Ω is completely captured by µ without needing µ is an indication that Ω is a holomorphic function of the coordinates of the off-shell parameter space. Note that the variation of Ω can in principle have a (3, 0) component. However, the (3, 0) part should be interpreted as a Kähler transformation and so is not part of the physical moduli. Another way of saying this is that ∆ in (3.2) is built from covariant derivatives on the parameter space [20]. Restricted to variations of the complex structure, the (3, 0) component is attributed to a connection in the usual way. The holomorphic deformations of the hermitian and B field moduli are (B + i ∆ω) (1,1) and where a subscript (p, q) denotes the type with respect to undeformed complex structure.
Here B is a combination of variations of the B field and the exact term in the variation of the Chern-Simons term, given by up to a d-closed two-form. The Green-Schwarz mechanism ensures B is gauge invariant. As we show in appendix B, flux quantisation then implies that B is a globally defined two-form, so it can indeed be a modulus. As we will see, the (0, 2) component of ∆ω is actually fixed in terms of the other moduli by the SU(3) relations, but it will be convenient to package this with B (0,2) into (B + i ∆ω) (0,2) .

F -term conditions from the superpotential
Let us review how one derives the F -term conditions from the superpotential. We take δ to be an infinitesimal deformation of the fields, leading to a corresponding variation of the superpotential δW The F -term conditions come from requiring that both W and δW vanish for generic values of the moduli. For arbitrary B + i δω and δA, the vanishing of δW requires ∂Ω = dΩ = 0, F ∧ Ω = 0. With this we see W = δW = 0 reproduces the F -term equations of the Hull-Strominger system. Our plan is to extend this discussion to understand finite deformations around a supersymmetric solution. First we need to understand how the requirement of an SU(3) structure and holomorphicity of the superpotential constrain the possible deformations.

Constraints from the SU(3) structure and holomorphicity
The existence of an SU(3) structure is part of the data that goes into the superpotential, so the deformed geometry should also define an SU(3) structure. Another way of saying this is that the off-shell parameter space on which the superpotential is varied is the space of SU(3) structures (plus bundles, and so on). This means the SU(3) structure compatibility condition must still hold: (ω + ∆ω) ∧ (Ω + ∆Ω) = 0, (3.12) where ∆ is a finite holomorphic variation. Expanding this out according to complex type, we find (∆ω) (0,2) = ı µ ω + ı µ (∆ω) (1,1) − 1 2 ı µ ı µ (∆ω) (2,0) . (3.14) Now consider an anti-holomorphic variation ∆ under which Ω does not vary as it is a holomorphic on the parameter space. We then note that as an anti-holomorphic deformation does not change the complex structure (as ∆Ω = 0) the SU(3) compatibility condition reduces to ∆ω ∧ Ω = 0. (3.15) From this we see (∆ω) (0,2) = 0 and so, taking a conjugate, (∆ω) (2,0) = 0. Combined with the previous result of a holomorphic variation of the compatibility condition, we have From this we see the (0, 2) component of the variation of ω is fixed by the complex structure and hermitian moduli. We can also play the same trick with the superpotential itself. The superpotential is a holomorphic function of the moduli so an arbitrary anti-holomorphic variation of it must vanish without having to impose the supersymmetry conditions -it must vanish JHEP10(2018)179 "off shell". Let us turn off the gauge sector for now, and consider an anti-holomorphic variation of W : As W is a holomorphic functional, this must vanish for all anti-holomorphic deformations without needing to impose the F -term conditions. Without imposing integrability of the complex structure, generically we have dΩ ∈ Ω (3,1) (X)⊕Ω (2,2) (X). For an anti-holomorphic variation of the superpotential to vanish, it is sufficient that where we have used the first condition from (3.16) to remove (∆ω) (0,2) . This agrees with (3.7) where we stated that the holomorphic combinations are the (1, 1) and (0, 2) components of B + i ω (see also discussion in appendix C.1). Taking a conjugate of these conditions we have Taken together these give For what follows, it is useful to definẽ which are our complex coordinates on the parameter space. One can repeat this exercise with gauge sector turned on. W is a holomorphic function of the moduli so it does not change for arbitrary anti-holomorphic variations ∆A. For ∆W = 0 to hold at a generic off-shell point in field space, we find it is sufficient that (∆A) (0,1) = 0. This implies (∆A) (1,0) = 0, so the holomorphic deformations correspond to ∆A = α ∈ Ω (0,1) (End V ). (3.22) In other words, the holomorphic coordinate on the parameter space is α, in agreement with (3.3). Furthermore, one sees that the holomorphic deformations of the complexified hermitian moduli are the (1,1) and (0,2) components of B + i∆ω.
As an aside, we note that there is a schematic way to see that Ω is a holomorphic function of the parameter space coordinates. Consider a generic anti-holomorphic deformation of the superpotential around a point in moduli space where the holomorphic top-form is closed, dΩ = 0: For an infinitesimal deformation, the second term can be dropped, and a sufficient condition for the first term to vanish for generic H and ω is that

JHEP10(2018)179
infinitesimally. Note that this also kills the second term in (3.23) at second order in perturbation theory. A sufficient condition for ∆W to vanish at this order is hence that ∆Ω = 0 to this level as well. This argument can be continued ad infinitum, and we are left with condition (3.24), at least for finite deformations away from a closed Ω. We will assume that this condition holds true at generic off-shell points in the parameter space. Stronger evidence for this is provided by the generalised geometry formulation of the off-shell parameter space presented in appendix C; we also find a matching between the complex coordinates and the natural parametrisation from generalised geometry.

Higher-order deformations
The main aim of this paper is to derive the conditions on the moduli when we move from infinitesimal to finite deformations of solutions to the Hull-Strominger system. In other words, we consider higher-order deformations of the fields that parametrise the supersymmetric Minkowski solution. As we have mentioned we only need to consider the F -term relations to understand the moduli space. We show in appendix D that under some reasonable assumptions the D-term conditions do not constrain the moduli problem and should be thought of as gauge fixing conditions -we expect this to hold in general.

The superpotential
Let us consider the effect on the superpotential of a finite deformation of the background fields away from a point on the supersymmetric locus. In other words, we start with a supersymmetric vacuum solution described by a superpotential W which is a functional of the SU(3) structure, H and the bundle. Let us denote the superpotential evaluated at this point by W | 0 . The vacuum is supersymmetric if both the superpotential and its first derivative vanish when evaluated on the solution. Now move a finite distance from this solution in parameter space by deforming the background. The superpotential evaluated at this new point is W | ∆ = W | 0 + ∆W . We have a supersymmetric solution if both W | ∆ and its first derivative vanish at that point in parameter space, which is equivalent to the vanishing of ∆W and δ∆W . Let us see how this works out. For clarity of presentation, let us ignore the bundle moduli -we will reinstate these in section 4.3.
A finite holomorphic deformation of the parameters gives As we are deforming about a supersymmetric point we have H = i(∂ − ∂)ω and dΩ = 0, so the first term simplifies and the last term vanishes, giving where the ı µ ı µ ı µ Ω term vanishes due to the type of d(∆B + i ∆ω), and we have written the second line in terms ofb and x, the (0, 2) and (1, 1) parts of the complexified hermitian moduli (3.21). As ı µ satisfies a graded Leibniz identity, we can rewrite the above as Our first condition for the deformed background to be supersymmetric is ∆W = 0 when evaluated on the solution. In other words, the terms in the brackets in (4.3) should be zero up to a ∂-exact term.
We also need to impose the vanishing of the first derivative of ∆W . As ∆W is a functional, this amounts to treating it as an action and finding the resulting equations of motion. Varying ∆W , one finds A few comments are in order. The condition in (4.5) is nothing but the Maurer-Cartan equation for finite deformations of a complex structure. This is somewhat expected as we know solutions to the Hull-Strominger system are manifolds with a complex structure. Notice that we also have a second condition on µ in (4.6) which is not usually seen in discussions on the moduli space of complex structures. This condition comes from requiring that the deformed three-form Ω + ∆Ω is closed and thus holomorphic with respect to the new complex structure -this is stronger than requiring a complex structure alone. Note that this same condition that appears in [40] for Kodaira-Spencer gravity -there it is imposed as a constraint from the outset but one should actually think of it as requiring that the deformed three-form Ω remains d-closed. One could make a change of variables which solves this constraint explicitly by taking ı µ Ω = a + ∂b where a is ∂-harmonic. This may be useful when investigating the quantum theory defined by the superpotential but we do not use it in what follows.

A Maurer-Cartan equation from the holomorphic structure
The idea now is that these equations can be interpreted as a Maurer-Cartan equation for the deformations. We know the infinitesimal moduli of the Hull-Strominger system are captured by the D-cohomology of the holomorphic structure, so we expect D will be the differential that appears in such a Maurer-Cartan equation. The other ingredient is a bracket. We introduce the bundle Q as the sum of the holomorphic cotangent and tangent bundles:

JHEP10(2018)179
This is the bundle Q defined in (2.16) with the gauge sector suppressed for now. We write (0, 1)-forms valued in this bundle as The holomorphic structure for the Hull-Strominger system without the bundles is given by an operator D that acts on sections of Q . We can also introduce a bracket [·, ·] on forms valued in Q and a pairing ·, · that traces over the Q indices: We give explicit expressions for D, the bracket [·, ·] and the pairing ·, · in equations (A.11), (A.13) and (A.16). Using these we can rewrite ∆W , given in (4.3), as Note that we have redefined the (0, 2)-form field as b =b − µ a ∧ x a . The equations of motion that follow from varying ∆W can then be written compactly as Let us make a few comments. Looking at ∆W in equation (4.11), we see it resembles a Chern-Simons action. More specifically, the form of the action is that of a holomorphic Chern-Simons theory for y with a Lagrange multiplier b that enforces a constraint for y. This constraint is the same as the gauge choice that is imposed in Kodaira-Spencer theory [40,41]. Note that the conventional Chern-Simons action has appeared as a superpotential in other work [42]; we expect a similar analysis could be applied here.
Notice also that infinitesimal deformations are captured by (4.14) It follows from this and dH ∝ ∂∂ω = 0 that ∂b is ∂-closed. 6 If the underlying manifold X satisfies the ∂∂-lemma or H (0,2) (X) vanishes, ∂b is ∂-exact and can be absorbed in a redefinition of the complexified Kähler moduli x. We then see that infinitesimally the complexified Kähler moduli are counted by H (0,1) (T * (1,0) X).

Including the bundle
We now want to include the bundle degrees of freedom in the superpotential -we do this by adding a Chern-Simons term ω CS (A) for the gauge connection A:

JHEP10(2018)179
The B field transforms in the usual Green-Schwarz manner so that H is gauge invariant: Consider a shift of the gauge connection by A → A + α. The corresponding change of the Chern-Simons form is The exact term in this variation combines with the variation of the B field to give As we show in appendix B.4, B is a globally defined two-form so it can be a modulus [14].
Expanding about a supersymmetric point, the extra terms in the variation of the superpotential from the gauge sector are where the final terms combine with b in ∆W to give , the extra terms in the variation of the superpotential are The full expression for ∆W is given by the sum of expressions (4.20) and (4.3) (with b → B). The equations of motion that follow from varying the full superpotential are Remarkably, the superpotential can still be written in a Chern-Simons fashion as Together with the vanishing of ∆W , these are the conditions for a supersymmetric Minkowski solution.

Vanishing of the superpotential
We now want to understand the condition ∆W = 0 in more detail. In what follows, we will consider the moduli problem with the gauge sector turned off. Everything we say goes through when we replace the D operator, bracket and pairing with those that include the gauge sector.
The deformed vacuum solution is supersymmetric if the equations of motion are satisfied and the superpotential itself vanishes. As X is assumed to be compact, the superpotential vanishes if the terms wedged with Ω are ∂-exact, that is where β is an arbitrary (0, 2)-form. Upon substituting the first equation of motion (4.24) into this expression, it simplifies to . 7 We can integrate this expression to give where k is constant as it is anti-holomorphic and X is compact. Combining this with the vanishing of the superpotential (4.27) gives k Ω = ∂b + ∂β.

JHEP10(2018)179
As Ω is not ∂ exact, we must have k = 0. We then identify β = −b up to a ∂-closed (0, 2)-form. Putting this all together the full set of equations is These equations are equivalent to the vanishing of the superpotential and its derivative, and so their solutions are a supersymmetric Minkowski vacuum. In other words, solutions (y, b) to these equations are precisely the moduli of the Hull-Strominger system. We pause to make a few comments. First note that these equations contain finitely many powers of y and b; the equations do not give an infinite set of relations. This is somewhat striking -generic deformations of geometric structures do not usually truncate at a given order. In our case the fact that the equations depend on terms up to O(y 3 ) is an indication that there is more structure to the Hull-Strominger system than at first sight. This extra structure is the existence of an underlying holomorphic Courant algebroid describing the subsector of deformations given by simultaneous deformations of the complex structure and the three-form flux. In deriving these equations and finding third-order equations for the moduli, we might be encouraged to think there is some sort of algebroid underlying the full heterotic system. Indeed, the ten-dimensional heterotic theory has a description in terms of generalised geometry [31,43,44] and the Hull-Strominger system can be recast as in terms of holomorphic Courant algebroid [30].
One might wonder if the form of these equations can survive α corrections. We derived the equations for the moduli by starting from the superpotential for the four-dimensional theory that one would get by compactifying on a solution to the Hull-Strominger system. Part of the data of such solutions is a complex manifold. A complex manifold admits an SU(3) structure whose torsion is constrained [45]. A special case of such manifolds are those with vanishing torsion so they have SU(3) holonomy and are Calabi-Yau. If the solution to the Hull-Strominger system admits an α → 0 limit, the α = 0 solution is simply Calabi-Yau. In this case it is known that the superpotential receives no α corrections (to finite order) and so, although the α -corrected geometry is not longer Calabi-Yau, the tree-level superpotential is exact [22,46]. This means equations (4.31)-(4.33) will be correct even after α corrections. It is not known what happens if there is no large-volume Calabi-Yau limit.
Up to this point there has been an asymmetry in the way we have treated the vanishing of ∆W and its derivative. In the next section we will see that we can combine these conditions into a single Maurer-Cartan equation for an L 3 algebra.

Moduli and an L 3 algebra
So far we have derived the equations that determine the moduli for solutions to the Hull-Strominger system for finite deformations. We will show in this section that these equations

JHEP10(2018)179
can be reinterpreted as the Maurer-Cartan equation for an L 3 algebra. At first sight there is no obvious reason why the deformations of a system as complicated as the heterotic string should be described by such a "nice" algebra structure. However it is not as surprising if one remembers that the data of the Hull-Strominger system is equivalent to a holomorphic Courant algebroid with a holomorphic vector bundle [30]. The L ∞ structures that govern deformations of Courant algebroids (or Dirac structures) have been found; in particular it is known that the deformation complex of a Dirac structure is isomorphic to a cubic L ∞ or L 3 structure [47][48][49][50].
As we review in appendix A.3, an L ∞ structure is specified by a choice of graded vector spaces Y n and multilinear products k [51,52]. The idea is that the conditions from the superpotential are most naturally written in terms of an L ∞ structure that combines the action of D on the moduli fields y and ∂ on the (0, 2) moduli b. We take the vector spaces Y n to be where y is a (0, 1)-form valued in Q and b is a (0, 2)-form. Using this notation, we write the multilinear products k as We give expressions for the k where the entries are arbitrary elements of Y n in (A.29). One can check that these products have the correct symmetry properties and obey the L ∞ relations, which we write in (A. 26). Note that this is highly non-trivial and is an indication of the underlying holomorphic Courant algebroid.

Quasi-isomorphism to a natural holomorphic L 3 algebra
We briefly remark that these structures have a nice mathematical interpretation: our L 3 algebra (Y, 1 , 2 , 3 ) is L ∞ equivalent to the underlying holomorphic algebra. Neglecting the gauge bundle for a moment, we have the sheaf E of D-holomorphic sections of Q T * (1,0) X ⊕ T (1,0) X and the sheaf of holomorphic functions O X . These form an L 3 algebra, with underlying two-term complex in precisely the same way that a Courant algebroid E together with the real C ∞ functions form an L 3 algebra [53] with two-term complex

JHEP10(2018)179
One can then consider the Dolbeault resolutions of the sheaves E and O X , and extend ∂ to a morphism between them: Our complex Y is then the total complex of the deleted resolution and the differential 1 is the natural differential on this (see e.g. [54]). Our construction gives higher n brackets on Y, providing an L 3 algebra structure on the total complex. As (5.5) is simply a resolution of (5.3), this construction essential provides us with a local reformulation of the holomorphic L 3 algebra (5.3) in terms of C ∞ objects. Explicitly, one has a map of complexes (of sheaves) as follows: As the cohomology of the total complex is the same as the cohomology of the complex it is resolving, this is a quasi-isomorphism. (One can check this explicitly in our case.) However, the morphism in (5.6) also respects the bracket structure of the L 3 algebras on each complex, thus it is a quasi-isomorphism of L 3 algebras. We conclude that, in the L ∞ sense, our L 3 algebra (Y, 1 , 2 , 3 ) is locally equivalent to the holomorphic algebra (5.3).
Including the gauge bundle in this construction is straightforward; one simply replaces the bundle Q above with the full holomorphic bundle (C.34) (this also recently appeared in [30]). One finds an essentially identical two-term complex to (5.3) (see [55] for the analogue of (5.4) including the gauge bundle), giving an L 3 algebra on the local holomorphic sections. Via the Dolbeault resolution, one sees that this is quasi-isomorphic to our L 3 algebra (Y, 1 , 2 , 3 ) (now including the gauge bundles) exactly as above.

An L ∞ field equation
As explained in [52], there is a natural field equation that one can write down for a given L ∞ structure. The constraint on the form of the field equation is that it is covariant under gauge transformations of the fields Y. In terms of the L ∞ products, the field equation is For us this expression truncates at third order as k≥4 = 0. Remarkably, the L 3 field equation coming from (5.2) reproduces the conditions from the vanishing of the superpotential and its derivative:

JHEP10(2018)179
In other words, the conditions for a supersymmetric Minkowski vacuum are equivalent to A particularly nice property of this rewriting is that the L ∞ structure gives us the gauge transformations of the moduli for free and guarantees that the gauge algebra closes. The gauge transformation of Y by a gauge parameter Λ = (λ, ξ) ∈ Y 0 is where the higher-order brackets vanish. In general the gauge transformations take field equations to combinations of field equations -the field equations are covariant. If one could construct an action that has the L 3 field equation as its equations of motion, one would expect that action to be invariant under the L 3 gauge transformations. In contrast to [52], we have not been able to find such a candidate action nor do we expect one to exist; this is due to the fact that the supersymmetry conditions are the vanishing of the superpotential and its first derivative. Note that the superpotential alone is not expected to be invariant under both λ and ξ transformations -we will see in the next section that the superpotential is actually invariant under the ξ transformations alone which correspond to shifts by ∂ a -exact forms.

A reduced L 3 algebra and an effective action
In this section we discuss some consequences of the L 3 algebra. In particular we comment on how the L 3 algebra can be reduced by quotienting by ∂-exact forms and show that this is equivalent to integrating out b. We also discuss the relation of the moduli (y, b) to the effective theory one would find by compactifying on a solution to the Hull-Strominger system. As we have mentioned, the form of the superpotential (4.23) closely resembles that of holomorphic Chern-Simons theory, and is in fact a generalisation of this theory. Holomorphic Chern-Simons theory has several interesting relations with mathematical disciplines such as open and closed topological string theory, knot theory, Donaldson-Thomas invariants and so on, and it would be interesting to look for heterotic generalisations of these relations. This will be the subject of future work. For now, we will restrict ourselves to making some observations about the (semi-) classical effective action (4.23), and its relation to the lower-dimensional effective physics.

Integrating out b
We want to integrate out the (0, 2)-form field b. Looking back at the form of the L 3 gauge transformations (5.10), taking Λ = (0, ξ) gives where the second term can be written as From this we see that b plays the role of a Lagrange multiplier. We shall see below that given certain assumptions about the Hodge diamond of X, specifically h (2,0) = 0, then the field b has no associated massless modes. We can then integrate out b, resulting in where y now satisfies the constraint ∂ı µ Ω = 0. We want to think of this functional as an effective action. Note that for ∂ı µ Ω = 0 there is a gauge symmetry of this action where ξ ∈ Ω 0,1 (X). We are thus led to defineQ as a reduced sheaf of Q whose sections satisfy the constraint and are defined up to ∂ a -exact terms: One can check that the brackets on Q are well defined onQ, and that the L 3 algebra descends to a differential graded Lie algebra (DGLA). A very similar sheaf is also considered in the discussion of βγ systems in [56]. Note that this DGLA is L ∞ quasi-isomorphic to the algebra on (5.1). The second gauge symmetry of (6.4) is a generalisation of the Chern-Simons symmetry. The superpotential is invariant under where λ ∈ Ω 0 (Q) satisfies ∂ı λ Ω = 0. Note that the gauge algebra generated by (5.10) is reducible; a gauge transformation by Λ = (λ, ξ) is trivial if

JHEP10(2018)179
so it is well defined as an equation on the sheafQ. Recall that we already know the superpotential is invariant under δy a = ∂ a ξ. One can show that equation (6.9), together with the condition that the effective action vanishes, is equivalent to the Maurer-Cartan equations (4.31)-(4.33). Indeed, note that (6.9) is equivalent to for some b ∈ Ω (0,2) (X). For solutions to this equation, the condition that the action vanishes can then be written as y, [y, y] = ∂-exact, which can be rewritten as for some β ∈ Ω (0,2) (X). Here we used that the second term on the left-hand side is ∂-exact -it integrates to zero against Ω due to the constraint ∂ı µ Ω = 0. This then gives the same starting point for our derivation of equations (4.31)-(4.33). It is beyond the scope of the present paper to investigate general solutions to (6.9) and ∆W [y] = 0, i.e. integrable deformations of heterotic geometries. We will however make some comments on the couplings derived from (6.4) in the four-dimensional effective field theory. From this we make a conjecture about the obstructions that can appear in the Maurer-Cartan equations.

Effective field theory and Yukawa couplings
Our starting point to derive the effective physics is the superpotential (4.23), where we have re-introduced the field b. When dimensionally reducing the theory, it is common practice to split our fields (y, b) into "massless" and "heavy" modes We imagine performing a formal dimensional reduction of the theory to a four-dimensional Minkowski background where we keep all the massive Kaluza-Klein modes for the time being. The corresponding mass matrix of the reduced theory reads 8 where {α, β, . . .} denote holomorphic directions in the parameter space and K is the Kähler potential. Full knowledge of the Kähler potential is not necessary at this point, but the curious reader is referred to [20,21] for more details. From the form of the mass matrix, it is easy to see that a field direction α is massless if and only if where the field directions γ can in principle be massive. 8 In principle, there is also a potential coming from D-terms. However, as we show in appendix D, the D-terms can be set to zero by a complexified gauge transformation and so they do not lift any moduli.

JHEP10(2018)179
In the end, we are interested in a reduced field theory of massless modes, where the massive modes have been "integrated out". It is easy to see that the field direction corresponding to b is massive (although need not be an eigenmode of the mass matrix). Indeed, from (6.16) it follows that b 0 must satisfy ∂ a b 0 = 0, (6.17) and so b is an anti-holomorphic section of Ω (0,2) (X). We restrict ourselves to geometries where this bundle has no sections, in other words This is true in particular for Calabi-Yau geometries. It follows that we can integrate out the b field as far as the effective theory is concerned, leaving us with the effective superpotential (6.4), where now the Beltrami differential component µ of y satisfies ∂ı µ Ω = 0 as above. From condition (6.16) it follows that the remaining massless fields y 0 then satisfy Dy 0 = 0, (6.19) where this should be viewed as an equation in the sheafQ. It is also natural to decompose the symmetry transformations (6.7) in terms of the massless and massive modes. A suggestive decomposition, given the condition (6.19), is the following This cohomology is isomorphic to H (0,1) D (Q) for manifolds satisfying either the ∂∂-lemma or H (0,1) (X) = 0 [14,15]. We give a brief review of this cohomology and its decomposition in into more familiar cohomologies by means of long exact sequences in appendix E.
Decomposed in terms of massless and heavy modes, the effective action now reads where we denote We see that the heavy y h modes are the only ones that propagate internally. Note that even though we take the expectation value of y h to vanish, by including internal quantum corrections, we see that the coupling between y 0 and y h can generate higher-order couplings of the massless fields. These new couplings are however of quartic order and higher in y 0 , and are hence non-renormalisable in the effective field theory. The only renormalisable

JHEP10(2018)179
coupling we need to worry about from an effective field theory point of view is therefore the Yukawa coupling ∆W Yuk (y 0 ) = − 1 3 X y 0 , [y 0 , y 0 ] ∧ Ω. (6.25) This argument is similar to and generalises Witten's standard argument for the gauge sector [22].
Note that in addition to the standard Yukawa couplings between bundle moduli, the Yukawa couplings (6.25) also contain couplings of gravity-gravity type (couplings of deformations of the geometry) and gravity-bundle type, often referred to as µ-terms in the literature. It would be interesting to investigate the phenomenological implications of such couplings, but it is beyond the scope of the present paper to do so.
Let α A ∈ H (0,1) (Q) denote a set of inequivalent cohomology classes spanning H (0,1) (Q), and expand where the C A now correspond to the four-dimensional fields, including in principle moduli and matter fields. The Yukawa couplings then read A massless field direction α A is then truly free if and only if In particular, this is true if Note that, starting from the Maurer-Cartan equation (6.9), this is simply the condition for an infinitesimal deformation in the field direction α A to be unobstructed. The effective field theory then prompts us to make the following conjecture: when H (2,0) ∂ (X) = 0, 9 the only non-trivial obstructions coming from the Maurer-Cartan equations are given by the constraints (6.29) on the infinitesimal moduli.

Conclusions
In this paper we have considered finite deformations of the Hull-Strominger system. Starting with the four-dimensional N = 1 superpotential, we showed that integrable deformations corresponding to holomorphic directions on the moduli space can be parametrised by solutions of a Maurer-Cartan equation for an L 3 algebra, which we described in detail.
There are many directions one could follow from this work. Firstly, one might wonder which of the infinitesimal deformations parametrised by H are there some special cases where a generalisation of the Tian-Todorov lemma applies? It would be interesting to apply our formalism to some explicit examples, and in this way work out the spectrum of free fields in the low-energy four-dimensional theory.
The superpotential led to a generalisation of holomorphic Chern-Simons theory [57][58][59] that couples hermitian and complex structure moduli. Following [57] and [40], it seems that one should think of this theory by taking spacetime to be spanned by the anti-holomorphic directions with holomorphic Ω-preserving generalised diffeomorphisms playing the role of a gauge group. 10 It would be interesting to study this further. In particular, by starting from the superpotential as an effective action and investigating its quantisation one might hope it defines a consistent quantum theory (cf. [61,62]) and gives analogues of Donaldson-Thomas or holomorphic Casson invariants for heterotic geometries. Note that the superpotential is complex in general so the path integral will not be convergent. Such complex path integrals have appeared before in the study of complex Chern-Simons theory [42,63,64] where they are understood by analytic continuation. We foresee a similar treatment here.
As a step towards a complete understanding of the quantum heterotic moduli space, one could construct a world-sheet AKSZ topological model [65,66] or a topological string model for the effective theory similar to Witten's open string model for ordinary Chern-Simons theory [67]. As a guide, one might start by comparing the heterotic moduli space with the spectrum of holomorphic βγ systems and the chiral de Rham complex [56,[68][69][70][71][72]. Several other approaches to the (0, 2) world-sheet have appeared over the years (see [46,[72][73][74][75][76][77][78][79] and references therein). It would be interesting to investigate how these methods connect with the approach outlined in the present paper. These are all interesting aspects which we hope to explore in future publications.
One might also consider the moduli space of heterotic compactifications on more exotic geometries, such as G 2 or Spin(7) manifolds [80][81][82]. In the case of G 2 compactifications, the form of the moduli space is remarkably similar: for example, the infinitesimal deformations are again captured by a cohomology. Despite this there are notable differences such as the analogue of the bundle Q not appearing as an extension. It would be interesting to investigate the finite deformation algebras in these cases, and in the process identify the corresponding L ∞ structure. This might give a G 2 generalisation of Chern-Simons theory.
We have been concerned with finding the honest supersymmetric deformations of solutions to the Hull-Strominger system. For this we only needed to consider the superpotential in the four-dimensional theory. Of course, the four-dimensional theory also has a Kähler potential which is important for understanding the physical potential of the effective theory. The Kähler potential and the metric on the moduli space have been worked out in recent publications [20,21]. One might wonder how these objects appear in our formalism. It seems that the cleanest description of these objects would follow from a proper analysis using generalised geometry. As outlined in appendix C, the N = 1 heterotic structure is 10 Generalised diffeomorphisms are transformations generated by sections of Q via the Dorfman bracket, whose antisymmetrisation is the Courant bracket. (In the real case, they are simply the diffeomorphisms together with gauge transformations of the supergravity fields.) The Ω-preserving condition is just (4.6) which ensures the deformed three-form is d-closed. Another Chern-Simons like theory featuring generalised diffeomorphisms as a gauge algebra has recently appeared in [60].

JHEP10(2018)179
described by an invariant objectχ so the Kähler potential should be given by a functional of this object, similar to the Hitchin functional for SU(3) and G 2 structures [83]. Indeed, it seems that a similar story applies to heterotic compactifications in other dimensions. We hope to make progress on this in a future work.
Note that even though we have the invariant objectχ we do not have a natural integrability condition for it -the generalised connection is not torsion-free in the heterotic string [43]. Curiously, it appears that when one looks at deformations of this structure there is a nice integrability condition (given by the superpotential). It would be interesting to see if this pattern persists for other generalised geometries.

Acknowledgments
We would like to thank Bobby Acharya, Dmitri Alekseevsky, Philip Candelas, José

A Conventions
In this appendix we set out our conventions and notation. We will use (m, n, . . .) indices to denote real coordinates and (a, b, . . . , a, b, . . .) to denote complex coordinates on the real six-dimensional manifold X. Using this we can expand, for example, a vector as Our fields are form-valued so, to save space, we often omit the wedge symbol where it will not lead to confusion. Our convention for the contraction of a vector-valued one-form with a p-form is that the vector index is contracted as usual and the form components are wedged. In coordinates, for a vector-valued one-form w and a p-form ρ, we have where ρ m is defined as It follows that ı w satisfies a Leibniz rule: The interior product of a vector with a one-form is extended to p-vectors and p-forms using the operation, defined as where u is a p-vector and ρ is a p-form. We indicate the p-vector obtained by raising the indices of a p-form with the metric g by a superscript -for example

JHEP10(2018)179
A.1 Heterotic supergravity The Hull-Strominger system follows from setting the supersymmetry variations of the tendimensional gravitino ψ, dilatino λ and gaugino χ to zero. In our conventions these are where ε is a ten-dimensional Majorana-Weyl spinor and ∇ LC denotes the Levi-Civita connection.

A.2 Holomorphic structure
Ignoring the gauge sector for the moment, the relevant fields are (0, 1)-forms taking values in Q = T * (1,0) X ⊕ T (1,0) X. We write (0, 1)-forms valued in this bundle as where x a and µ a are (0, 1)-forms. More generally we will write y p = (x p , µ p ) ∈ Ω (0,p) (Q ) -we will often drop the subscript denoting the form degree. The holomorphic structure is defined by a D operator that is nilpotent. The action of D on y ∈ Ω (0,p) (Q ) is :

JHEP10(2018)179
The bracket does not satisfy a graded Jacobi identity. As is the case for a Courant algebroid, the Jacobi identity holds up to a ∂-exact term. Evaluated for y ∈ Ω (0,1) (Q ), one finds where the final component is the gauge field piece and F a = F ab dz b . The bracket is where we have written only the extra terms that appear in the bracket. Again, the bracket obeys a Jacobi identity up to a ∂-exact term. The pairing between sections is given by We follow [52] for the conventions of an L ∞ algebra in the " -picture". 11 We start with a graded vector space Y Y = n Y n , n ∈ Z, (A.22) 11 We have swapped n → −n so that the degree of Yn matches the form degree of y

JHEP10(2018)179
where the Y n are of degree n. The L ∞ algebra admits multilinear products 1 , 2 , . . ., where k has degree 2 − k. This means 1 is degree 1, 2 is degree 0, 3 is degree −1, and so on. The products are graded commutative so that for example where a superscript denotes the degree of Y ∈ Y n . More generally we have The sign has two contributions: (−1) σ gives a plus if the permutation is even and a minus if the permutation is odd; (σ; Y ) is determined by In these conventions, the first few L ∞ identities are The field equations and gauge transformations are where Y ∈ Y 1 and Λ ∈ Y 0 . The multilinear products k for the moduli of the heterotic system are

B Comments on heterotic flux quantisation
In this appendix we comment on flux quantisation in the heterotic string, and its relation to the global nature of the deformation of certain quantities that appear in this paper. Our discussion applies to general ten-dimensional heterotic supergravity and can thus be applied to solutions other than four-dimensional Minkowski compactifications. Note that we are not saying anything new here; understanding flux quantisation in the heterotic string is still an open problem for the simple reason that H is not d-closed in general. We proceed in steps, beginning with a toy example of the quantisation of the flux of an abelian gauge bundle. We then present a similar treatment of the two-form gerbe as a warm-up for the case of the heterotic gauge fields. We follow [84] for much of the early part of the discussion.

B.1 A toy example: the abelian bundle
Consider a vector bundle over a manifold X and let A denote an abelian connection with curvature F = dA. Let {U i } denote an open cover of X. We denote the overlaps by U ij = U i ∩ U j and so on for higher intersections. We assume that the cover is "good" so that the U i and their intersections are contractible. We employ the standard notion for theČech co-boundary operator where appropriate: For the curvature to be well defined, we require that on U ij we have where A i and A j denotes the connections on U i and U j . As U ij is contractible, by the Poincaré lemma, we must therefore have on U ij for some zero-forms γ ij . On triple overlaps U ijk we have which is often referred to as taking the co-cycle of dγ ij . It follows that c ijk = (∂γ) ijk are constants: Clearly from (B.4) we have ∂c = 0, so the c ijk define a class of the sheaf cohomology which represents the two-form flux F = dA. From (B.4) this might look like a trivial co-cycle, but this is deceptive since the class is trivial only if the individual γ ij can be chosen to be constant. In this case, we see from (B.2) that the connection A can be made global so that the flux is trivial. In this language flux quantisation is the statement that the c ijk are in fact integers, c ijk ∈ Z(U ijk ), so that they define an integral class [c ijk ] ∈Ȟ 2 (X; Z) ∼ = H 2 (X; Z)

JHEP10(2018)179
Gauge transformations of the connection do not change the class [c ijk ] ∈Ȟ 2 (X; R). To see this consider a gauge transformation which preserves the curvature F : This transformation induces a deformation of the γ ij as where κ ij ∈ R(U ij ) are constants. Thus c = ∂γ = ∂γ + ∂κ = c + ∂κ is shifted by aČech co-boundary, and so [c ijk ] = [c ijk ] ∈Ȟ 2 (X; R) is unchanged. Note that in the case of a quantised flux, the integrality of c ijk holds only in a preferred set of gauges for γ ij under shifts by real constants κ ij . Now consider a general deformation of the above system. We denote the variations by ∆(. . .): for example A = A + ∆A. We fix the (integral) cohomology class of the quantised flux so that ∆c ijk = 0. This means that we must have (∂∆γ) ijk = 0 and thus as C ∞ (R) is acyclic we can find i ∈ C ∞ (U i ; R) with ∆γ ij = (∂ ) ij = i − j . We see the γ ij can be deformed only by a gauge transformation, as in (B.7). Performing a global gauge transformation A i = A i − d i , we find that in the new gauge we havẽ We have shown there exists a gauge in which the variation of the connection is a global one-form∆A i =∆A j . Note that we could have performed the deformation requiring only that ∆c ijk = (∂κ) ijk for κ ij ∈ R(U ij ) -this fixes the real cohomology class [c ijk ] ∈Ȟ 2 (X; R). We would then have deduced that ∆γ ij = (∂ ) ij + κ ij , leading to the same gauge transformation of A as above. This would correspond to deforming away from the gauge (choice of γ ij ) in which c ijk are explicitly integral, while above we restricted ourselves to the gauge with c ijk ∈ Z(U ijk ).
The story becomes more intricate for non-abelian bundles. Recall that a non-abelian connection A on overlaps U ij transforms as For the purpose of the present paper we will assume without proof that we can take ∆g ij = 0, as in the abelian case.

B.2 The two-form gerbe example
The case of a two-form gerbe is a direct generalisation of the abelian bundle. The gerbe is specified by a set of two-forms B i ∈ Ω 2 (U i ) covering the manifold. The field strength H = dB is globally defined. As before, this means that on overlaps U ij we have

JHEP10(2018)179
for some one-forms λ ij ∈ Ω 1 (U ij ). Again we have (∂B) ij = dλ ij on U ij so that (∂λ) ijk is d-closed d(∂λ) ijk = (∂ 2 B) ijk = 0. (B.12) As U ij is contractible (∂λ) ijk is actually exact for some ζ ijk ∈ C ∞ (U ijk ; R). Finally we take the co-cycle of dζ ijkl on quadruple over- 14) It follows that the functions c ijkl = (∂ζ) ijkl are constants. As before, if the flux is quantised the c ijkl can be made integral by appropriate choice of ζ ijk . Indeed, the c ijkl represent a co-cycle in the sheaf cohomology where the class in H 3 (X; Z) is given by the flux H = dB. Consider deformations of the above system. From the quantisation of the constants c ijkl and the acyclicity of C ∞ (R) we have for some κ ij ∈ C ∞ (U ij ; R). This leads immediately to 12 for some i ∈ Ω 1 (U i ). Putting this all together we have Looking at these formulae, we see that the variations are forced to take precisely the form needed to constitute a gauge transformation of B and its descendants λ and ζ. Thus, in an appropriately chosen gauge, we have that ∆B is a global two-form.

B.3 Heterotic flux quantisation
We now come to the example relevant for the present paper: the B field of the heterotic string. Recall that the anomaly cancellation condition reads

JHEP10(2018)179
where ω CS (A) denotes the Chern-Simons three-form and we set α 4 = 1. A gauge transformation of A of the form induces the following transformation of the Chern-Simons form [43] ω where Let us now consider the patching on U ij . From equation (B.20) we get for some one-forms λ ij . Taking a co-cycle of this on triple intersections U ijk gives the relation (∂ω 2 ) ijk = d(∂λ) ijk .

(B.29)
It can further be shown that for some one-forms ω 1 ijk (g) ∈ Ω 1 (U ijk ). The expression for ω 1 ijk (g) is not relevant in the present context, but it can be taken to be independent of the gauge connection A, depending on the transition functions g ij . To see this, simply vary ω 2 ij (g, A) with respect to A. We find that where we have defined ∆A = α, which transforms appropriately in the adjoint representation when the transition functions g ij are kept constant. From this it is clear that which shows that we may take ω 1 ijk (g) to be independent of A without loss of generality. Putting together (B.29) and (B.30), we then get the relation

JHEP10(2018)179
for some functions ζ ijk ∈ C ∞ (U ijk ; R). We again co-cycle this relation on quadruple intersections U ijkl to get d ω 0 ijkl (g) + (∂ζ) ijkl = 0, (B.34) where we have used the fact that for some functions ω 0 ijkl (g). We can co-cycle this relation again on quintuple overlaps U ijklm to get d(∂ω 0 ) ijklm (g) = 0, (B.36) and thus the numbers k ijklm = (∂ω 0 ) ijklm (g) define a class This represents the Pontryagin class tr F ∧ F in the cohomology of the sheaf R. Note however from (B.34) we have for our particular bundle for some constants c ijkl . Computing the co-cycle (∂ω 0 ) ijklm (g) we thus find and so [k ijklm ] = 0 ∈Ȟ 4 (X; R). This is the sheaf cohomology version of the statement that the bundle in question has a trivial first Pontryagin class as tr F ∧ F = dH is exact.

B.4 Deforming the system and a well-defined global two-form
We now consider variations of the above heterotic story. The goal is to show that the heterotic quantisation condition leads us naturally to a global two-form This is an essential part of the moduli of the main text, where α = ∆A is the global one-form variation of the gauge connection. The story is very similar to that of the two gerbe, with some subtleties. Note that as for the abelian bundle, we will assume that we can choose to keep the transition functions g ij constant under deformations, that is even in the case where the bundle is non-abelian. With this assumption, any deformation of the bundle connection α = ∆A can be assumed to be a section of Ω 1 (End V ). We begin by noting that, imposing (B.41), a deformation of (B.38) gives Hence the constants ∆c ijkl define a sheaf cohomology class inȞ 3 (X; R), even though the original c ijkl did not. We will require that this class vanishes [∆c ijkl ] = 0 ∈Ȟ 3 (X; R), so

JHEP10(2018)179
that the deformation does not produce any new third cohomology. One could think of this as fixing the topological data associated to the original numbers c ijkl , assuming that they have some notion of "integrality" which cannot be continuously deformed, even though they do not explicitly define such a class themselves. Indeed, this condition also seems to agree with the general world-sheet arguments on the heterotic B field and flux quantisation made in [85], where one thinks of the B field as a torsor. The lack of the notion of zero in the space of B fields is reflected in our setup by the lack of an explicit cohomology class associated to the c ijkl . However, the variations of B about a given starting point do define a cohomology class naturally. From this requirement we have for some r ijk ∈ R(U ijk ) and κ ij ∈ C ∞ (U ij ). Next we take the variation of (B.33), again imposing (B.41), and use the acyclicity of Ω 1 (U ijk ) to find exactly as for the simple two gerbe case. Finally, we take the variation of (B.28) to obtain We see we can absorb the d terms via a global gauge transformation of B as before, so that on the overlaps U ij we have Thus B defines a global two-form. Note in particular that B is gauge invariant with respect to gauge transformations of the bundle [20]. This follows from the Green-Schwarz mechanism wherein the B field transforms as C The off-shell N = 1 parameter space and holomorphicity of Ω In this appendix, we present a description of the space of off-shell scalar field configurations that appears when we rewrite the ten-dimensional theory with manifest four-dimensional N = 1 supersymmetry. This rewriting is done in the same spirit as the rewriting of elevendimensional supergravity as a four-dimensional N = 8 theory in [86] and the rewriting of type II supergravities as four-dimensional N = 2 theories in [37,38,[87][88][89]. Here we focus only on the scalars, which are described as an (infinite-dimensional) space of generalised geometric structures of a type we will specify. We then show how to recover the threeform Ω from the generalised geometric structure and outline how one can see that it is holomorphic on the off-shell field space, as we claim in section 3.

C.1 SU(3) × SO(6) structures in the NS-NS sector
As a warm up, we consider the common NS-NS sector of ten-dimensional string theories, which is well known to admit a description in the language of SO(10, 10) × R + generalised geometry [90] (see also [91,92]). When rewriting the theory as a four-dimensional theory, we consider the spacetime to admit a product structure, breaking the ten-dimensional Lorentz symmetry to SO(3, 1) × SO (6). The fields which are naturally SO(3, 1) scalars can be packaged as the generalised metric of an SO(6, 6) × R + generalised geometry on the six internal spatial dimensions (see for example [38,87]).
As discussed in [28], the conditions for an N = 1 supersymmetric vacuum can be phrased as the existence of an integrable SU(3) × SO(6) structure on the internal sixdimensional generalised tangent space E T X ⊕ T * X. Such a structure encodes the generalised metric and thus the physical fields. If we require only the presence of an off-shell N = 1 supersymmetry, then one merely restricts the fields to admit a globally defined spinor, corresponding to a (possibly non-integrable) SU(3) × SO(6) structure. The configuration space of off-shell scalar fields that we require in our four-dimensional N = 1 description of the ten-dimensional theory is thus the (infinite-dimensional) space of SU(3) × SO(6) structures on the generalised tangent space E. This is the N = 1 NS-NS sector analogue of the discussion of the N = 2 vector-and hyper-multiplet structures given for the full type II theories in [37,88].
An SU(3) × SO(6) ⊂ SO(6, 6) × R + structure on E is specified by a particular type of complex sectionχ of where L ∧ 6 T * X is the auxiliary R + bundle transforming in the 1 +1 representation of SO(6, 6) × R + as introduced in [90]. It is convenient to writeχ as where Φ ∈ Γ(L) is the generalised density with Φ = √ ge −2φ in the coordinate frame (see [90]). In order forχ to have the correct stabiliser, χ must lie in a particular orbit under the action of SO(6, 6) on ∧ 3 E C at each point of the manifold X.
One can alternatively describe such a generalised structure in terms of a generalised metric, which defines an SO(6) × SO(6) structure on E, which then splits as E C + ⊕ C − , so that To this, one must add a non-vanishing section of the spin bundle for C + which we denote JHEP10(2018)179 η ∈ Γ(S(C + )). Using local basesÊ + m for C + andÊ − m for C − defined by 13 one can write an explicit formula for the object χ where Ψ is the three-form spinor bilinear of η with itself. Writing the object χ in this way, it is guaranteed that it will lie in the correct orbit and its stabiliser is SU(3) × SO(6). Given an N = 2 structure parametrised by two complex pure spinors Φ ± ∈ S(E)⊗L 1/2 (see [38,93]), one can also build an N = 1 structure via the expressioñ thus providing a third description of the structure.
Having defined the structure, we will now show how to extract the ordinary complex three-form on the manifold from it. Recall that the generalised tangent space is an extension with the classes of such extensions labelled by the cohomology class of the three-form flux [H] ∈ H 3 (M ; R). The map π is referred to as the anchor map of the Courant algebroid E. This anchor map induces further maps on tensor products of E. In particular we obtain an induced map, which we also label π, π : ∧ 3 E −→ ∧ 3 T X.

(C.9)
Acting on the bundle L ⊗ ∧ 3 E and using L det T * X we have (C.10) Thus, applying the anchor map to the generalised structureχ we obtain an ordinary threeform on the manifold. We now calculate this three-form explicitly. We first note that which immediately gives us These are the split frames of [90], but constructed from an arbitrary local frame for the tangent bundle rather than a vielbein for the ordinary metric. However, the ordinary metric g is still used to lower the indices on the one-form frames e ±m , such that the O(6, 6) metric has components in these frames.

JHEP10(2018)179
However, we are interested in the objectχ = Φχ and in our frame we have Φ = √ ge −2φ so that π(χ) = 1 3! √ ge −2φ Ψ mnp (ê m ∧ê n ∧ê p ) ∈ ∧ 6 T * X ⊗ ∧ 3 T X C (C.13) Now we use the natural isomorphism To conclude our discussion of the SU(3) × SO(6) structure in SO(6, 6) × R + generalised geometry, we note that standard group theoretical arguments give us that the homogeneous space is diffeomorphic to the orbit ofχ (at a point in M ) under the action of SO(6, 6) × R + on the complex 220 C representation. The homogeneous space (C.16) has a complex structure, with respect to which the elementχ of the 220 C is holomorphic, that is the embedding of (C.16) into C 220 is a holomorphic map. If we imagine that this complex structure naturally extends to the infinite-dimensional space of SU(3) × SO(6) structures on E, then we expect that our generalised three-formχ will be holomorphic on that space. As the anchor map π is linear and fixed by the topology of E, we have that π(χ) = −i e −2φ Ψ = −i Ω is also holomorphic on the space of such structures. Note that the decomposition of the Lie algebras appearing here is so that we can locally parametrise the space (C.16) by exponentiating the action of the complex Lie algebra elements which do not annihilateχ. Taking out the overall scale we haveχ = e c e α+β ·χ, The parameters c, α and β then become local complex coordinates on the space (C. 16). Note that as SU(3) objects, α and β carry the same degrees of freedom as the parameters (µ, (∆ω + i ∆B) (1,1) , (∆ω + i ∆B) (0,2) ) used to parametrise the space of N = 1 field configurations in sections 3 and 3.2 of the main text, while the parameter c is associated to Kähler transformations.

C.2 SU(3) × SO(6 + n) structures in heterotic supergravity
We can perform a similar analysis to the above in SO(6, 6 + n) × R + generalised geometry for heterotic supergravity [43,44], where n is the dimension of the gauge group G. In that geometry the generalised tangent space E is an extension of the real Atiyah algebroid A by the cotangent bundle with the composition of the maps above still giving an anchor map π : E → T X.
The structure group of the bundle E is then the geometric subgroup (GL(6, R) × G) ((T * X ⊗ g) ∧ 2 T * X) of SO(6, 6 + n) × R + , though as usual in generalised geometry we think of E as an SO(6, 6 + n) × R + vector bundle. The generalised metric defines an SO(6) × SO(6 + n) structure on E, with a local frame (Ê + m ,Ê − m ,Ê − α ) which can be built from the physical fields (g, B, φ, A) (see [43] for details). The presence of a single globally defined spinor on M breaks the SO(6) factor to SU(3), so that an N = 1 structure is an SU(3) × SO(6 + n) structure on E. Using the split frame, one can again write an explicit formula for an objectχ ∈ L ⊗ ∧ 3 E defining such an N = 1 structure: As we still have π(Ê + m ) =ê + m , the argument of appendix C.1 still holds to give us that π(χ) = −i e −2φ Ψ, and similar group theoretical reasoning leads us to conclude that bothχ and e −2φ Ψ = Ω are holomorphic on the coset space SO(6, 6 + n) × R + SU(3) × SO(6 + n) (C. 22) at each point of X. In this case, the corresponding Lie algebra decomposition reads so that we can locally parametrise the orbit ofχ viã χ = e c e α+β ·χ c ∈ C * α ∈ (3, 1) +2 β ∈ (3, 6 + n) +1 . (C.24) As SU(3) × G objects, we have the same degrees of freedom as in appendix C.1, but now augmented by β a α , that is the (0, 1)-form part of the deformation of the gauge field. These then match the parametrisation of the off-shell N = 1 field space of section 3 as employed in section 4.3.

JHEP10(2018)179
C.3 The off-shell hermitian structure on V The off-shell parameter space of the Hull-Strominger system is the space of SU(3) structures, B fields and gauge fields for G. On-shell, we know that the gauge bundle V must be a polystable holomorphic bundle. We outline which parts of this structure on the gauge bundle survive off-shell. The gauge group G is a compact unitary real Lie group, which has complex representations. The vector bundle V C is a complex vector bundle with structure group G C , the complexification of G. On-shell (on a solution of the Hull-Strominger system), this complex vector bundle has two structures on it: a holomorphic structure and a hermitian structure.
The hermitian structure on V C is a reduction of the structure group G C to the compact real form G. This is given by a set of local G C frames for V C on patches of X, which, on the overlaps of patches, are related by transition functions in G ⊂ G C only. The holomorphic structure on V C is given by a set of G C frames with respect to which the transition functions are holomorphic functions into the complexified group G C . Clearly these two sets of frames are different for the simple reason that there are no holomorphic maps into the real group G.
From the complex vector bundle V C , we can define a real vector bundle which also admits an action of the complex group G C . The hermitian structure is a positive definite metric on V which pairs V C and V C , such that in the special frames alluded to above its component matrix takes the canonical form where α, β = 1, . . . , dim C V C are indices for V C and α, β = 1, . . . , dim C V C are indices for V C . Off-shell, the manifold X admits an almost complex structure as part of the SU(3) structure. As X is not honestly complex, we lose the holomorphic structure on V C -we cannot define holomorphic maps without an integrable complex structure. Nevertheless, physics tells us we have a real connection on V R and that the physical gauge group G is compact and unitary. This means we have the hermitian structure, even off-shell. 14 This hermitian structure simply says that V C defines a real vector bundle V with a compact unitary structure group.
The physical connection A is a local section of The almost complex structure defines a split of this into (1, 0)-and (0, 1)-form parts. As the (1, 0)-and (0, 1)-form summands are intrinsically complex, it is natural to see them as living in the complexified Lie algebra, i.e. Ω (1,0) (X; g C ) and Ω (0,1) (X; g C ), as these are JHEP10(2018)179 defined via tensor products over the field C. Naively, the connection has four parts with indices as However, as the connection is real hermitian, any one of them defines the rest via complex conjugation and multiplication by the hermitian metric on V C . Explicitly, the real condition fixes while the hermitian condition fixes Together, these allow us to determine all parts of the connection given only (A a ) α β ; for example we have On-shell, the holomorphic structure means that, in each patch, one can choose a G C gauge where (A a ) α β = 0. Off-shell we cannot do this in general, but we can use the previous identities to write any formula purely in terms of (A a ) α β . For example, the terms appearing in the superpotential involve traces, which simplify using identities like These enable us to write all of the needed expressions using only the objects with indices for V C (eliminating the appearance of objects with indices for V C ); for example so that ω CS (A) ∧ Ω features only the components (A a ) α β . Equivalently we could have written this expression purely in terms of V C indices or some combination of the twothis freedom reflects the fact that it is V that appears in the heterotic system, so a split into V C and V C is somewhat arbitrary. The important point is that it is the (0,1) part of the gauge field that appears. This mirrors the generalised geometry argument from the previous section.
Using the freedom to write the Chern-Simons form using only V C indices, one can take the bundle appearing in the deformation complex to be the holomorphic bundle (C. 34) which also appears in [30].

JHEP10(2018)179 D Comments on D-terms
In addition to satisfying the F -term conditions derived from the superpotential, a solution should also satisfy the D-term conditions: The first condition is referred to as the conformally balanced condition, while the second condition is the Yang-Mills condition. In this appendix we want to show that these conditions impose no extra constraints on the heterotic moduli, given some mild assumptions on the geometry and bundle. This is of course common knowledge from the supergravity point of view [33]. 15

D.1 Massless deformations
We begin by considering compactifications without bundles. The D-term condition of relevance is the conformally balanced condition (D.1). Consider first a massless deformation of this condition so that y 0 = (x 0 , µ 0 ) satisfies Dy 0 = 0. Our plan is to show that the D-term conditions can be solved order-by-order using the gauge symmetries of x 0 so that they do not further constrain the moduli. First we note that deformations of the hermitian form ω and the dilaton φ are linked via the SU(3) normalisation condition Remembering ∆ω (0,2) = 0, ∆ω (1,1) = −ix and ∆ω (0,2) = ı µ ω −i ı µ x, a massless holomorphic variation of this condition gives Now we expand the massless deformation in terms of a small parameter : with a corresponding expansion for x and µ. At first and second order the SU(3) normalisation condition fixes

JHEP10(2018)179
This continues to higher order -the deformation of the dilaton at n th order is fixed by the non-primitive part of x (n) and the lower-order fields x (i) for i < n. Now consider a massless holomorphic deformation of the conformally balanced condition. We expand in and separate into complex type. At O( 1 ) the deformation of the conformally balanced condition reduces to Note now that the Gauduchon metric defined byω = e −φ ω is balanced asω ∧ω is d-closed. We do this as on a hermitian manifold the various Laplace operators for a balanced metric agree on functions [94]:∆ We also use that the Hodge stars 16  Using this we can write the previous equations as where we have used the relation between the trace of x (1) and φ (1) given in (D.7), and we have defined The Hodge decomposition for Aeppli cohomology implies that (D.14) and (D.13) determine the (∂ + ∂)-exact part of X (1) . Indeed an equivalent set of equations is We now want to argue that these conditions are simply gauge conditions and so do not impose extra conditions on the moduli. Recalling the form of Dy 0 from (A.11) we see shifts of x 0 by ∂-exact terms drop out explicitly and that shifts by ∂-exact terms fall out as we are working modulo ∂-exact forms. 17 A gauge choice for x 0 then amounts to a choice of 16 We are using the convention for the Hodge star where α ∧ β = α β vol so that ω = 1 2 ω ∧ ω and Ω = −i Ω. The dual of a primitive (1, 1)-form αp satisfying ω αp is αp = −ω ∧ αp. We also have α20 = α20∧ω where α20 is a (2, 0)-form. This choice satisfies 2 = (−1) p on a p-form. The adjoint Dolbeault operators are defined by ∂ † = − ∂ , and we denote the corresponding operators for the Gauduchon metric with a tilde. 17 One can also do this calculation with y and b so that the shift by ∂-exact forms is explicit too. so that we get a Hodge decomposition of the space of (∂ + ∂)-exact forms as Note that it is important that we include the dilaton degrees of freedom: shifts of x (1) are generically not primitive and so they will change the SU(3) normalisation condition, but we can compensate for this by shifting φ (1) (which does not appear explicitly in the first-order conformally balanced condition). We start by shifting x (1) by −∂∂κ (1) , where κ (1) is a function. A short calculation shows that equation (D.16) becomes One can check that the operator acting on κ (1) is a positive semi-definite self-adjoint elliptic operator whose image is given by non-constant functions. This means that (D.21) can always be solved by an appropriate choice of κ (1) . This fixes the ∂∂-gauge symmetry of x (1) .

JHEP10(2018)179
where We see that X (2) depends only on the second-order correction to x 0 while A (2) and B (2) are fixed by the first-order terms (which should be thought of as gauge transformed to solve the first-order conditions). Again (D.24) and (D.25) are equivalent to ∂∂ †(X (2) + A (2) ) = i ∂∂ †ı µ (2)ω + ∂∂ †B (2) . (D.31) As X (2) is a function of x (2) alone, we can perform gauge transformations of x (2) without affecting A (2) and B (2) . Generically these gauge transformations will break the SU(3) normalisation condition, but we can always shift φ (2) to compensate for this (which is what we have implicitly done by eliminating φ (2) from the equations). An analogous argument to the one we gave previously then shows that we can always solve these conditions using the gauge freedom of x (2) . From this it is simple to see that this process can be continued to all orders. The conformally balanced condition at order n is a set of equations for x (n) with x (i<n) fixed. Again, one can always find a gauge transformation of x (n) that solves these conditions. From this we conclude that the D-term conditions are gauge fixing conditions for the moduli and do not further constrain the moduli problem.

D.2 Including bundle moduli
We now want to show that when we include the bundle moduli we can solve the D-terms and Yang-Mills equations simultaneously. As stated at the start of section 3, we have the freedom to work with either End V or End V C . We choose to work with End V C in this appendix and appendix E, so that the Donaldson-Uhlenbeck-Yau and Li-Yau theorems directly apply [95][96][97].
The Yang-Mills condition is where we have definedρ = 1 2ω ∧ω. A massless holomorphic deformation of this condition reads Note that ∆ 0ρ is completely determined by x 0 .

JHEP10(2018)179
For now we assume H 0 (End V C ) = 0. This is true for irreducible bundles with hermitian connections satisfying the Yang-Mills equation -such bundles are stable. If H 0 (End V C ) = 0 the Yang-Mills condition can give rise to Fayet-Iliopoulos D-terms in the lower-dimensional supergravity. Such terms have been studied extensively in heterotic compactifications (see [8,11] and references therein). Their appearance can be understood in the context of the infinitesimal moduli problem as modding out by D-exact terms [14,26], see also appendix E.
The equations of motion Dy 0 = 0 admit gauge transformations of x 0 and α 0 : Upon expanding x 0 and α 0 in powers of , the first order contributions to (D.33) are We now want to show that we can solve this and the conformally balanced condition by an appropriate gauge choice. We denote these conditions schematically as where D is the usual operator, and C and Y are the operators that act on y (1) to give the relevant equations. Explicitly (to first order) we have Let us start with a solution to the equations of motion The equations of motion have a gauge symmetry under shifts by ∂a + ∂b + tr F γ: Generically a solution y (1) to Dy (1) = 0 will not solve either Cy (1) = 0 or Y y (1) = 0 on the nose. Let us make a gauge transformation of y (1) : This still solves the equations of motion -Dy (1) = 0 -but affects the remaining two conditions: We showed in the previous section that there is always a choice of α and β that solves this condition, that is

JHEP10(2018)179
At this point we have a solution y (1) that solves Dy (1) = Cy (1) = 0 but generically has Y y (1) = 0. Note also that the forms α and β are fixed by y (1) . We now want to show that we can always find a solution to Y y = 0 using gauge transformations. Let us define This solves Dy (1) = 0 -what happens to the other two conditions?
Again, we can solve Cy (1) = 0 for any γ (1) by choosing a and b appropriately. We now want to see if we can solve Y y (1) = 0 using the gauge freedom in γ (1) .
Note that C picks out the (∂ + ∂)-exact terms of the input (and has no kernel acting on this) so that Cy (1) = 0 is equivalent to Using this we have where P projects onto Ω (1,1) (X)\{im ∂ + im ∂}. Explicitly we have 19 One can check that the operator acting on γ (1) is a positive semi-definite elliptic operator with trivial kernel, so one can always find a γ (1) that solves Y y (1) = 0. Again, one can see that this argument will hold to all orders in . Recall that we included T X as part of the gauge bundle V . Strictly speaking one should show you can solve the analogue of (D.49) for T X as well. Naively an extra minus sign will appear in the trace term, leading to a question of whether the operator acting on γ (1) has trivial kernel. Physically we expect this equation will not cause any problemsthere is evidence that the degrees of freedom in the connection for T X are not real moduli as they can be eliminated using field redefinitions (at least to O(α )) [98]. Another way of saying this is that the hermitian Yang-Mills condition for T X is automatically satisfied provided the other supersymmetry conditions are solved and one chooses the appropriate connection, which at O(α ) is the connection ∇ − . Though we do not have a proof, we expect this will hold to higher order in α when one chooses the connection on T X to be an instanton.

D.3 Polystable bundles
Let us now consider the situation where the bundle is polystable, that is we have a reducible bundle where each factor is stable of slope zero. If a bundle factor has non-vanishing first Chern class, which is possible for a sum of irreducible components, then it admits a single section proportional to the identity isomorphism on the bundle factor. This can give rise to a non-trivial Fayet-Iliopoulos term in the lower-dimensional supergravity. We now show that equation (D.49) can still be solved when such sections are present.
A polystable bundle is a sum of irreducible components where each factor V (i) C is a stable slope-zero bundle with non-zero trace. The trace-less part of (D.49) can be solved by the methods outlined in the previous subsection. The non-trivial part comes from the trace, which we isolate by taking the trace to get − i˜ tr(Y y (1) ) (i) = e −φ tr F (i)˜ x (1) + ∂ † A tr α (1) + e −φ tr F (i)˜ P[tr F (i) γ (i) (1) ] +∆ tr γ (i) (1) . (D.51) It is simplest to show that this can be solved by reintroducing the α parameter, which we have neglected until this point. To simplify the notation, we imagine the bundle V C is abelian for the remainder of the section. Reintroducing α , equation (D.51) becomes (1) ] +∆γ (1) . (D.52) We can now expand this equation in α . We expand the deformations x (1) and α (1) , together with the gauge parameter γ (1) , while keeping the background geometry fixed. The potential non-zero contributions to γ (1) are then (1) + α γ (1) as a function of x (1) . The remaining non-constant part can then be solved for modulo a constant term, which determines γ

JHEP10(2018)179
It is easy to see that this can be continued to higher orders in α . At the next order, we use the so-far undetermined constant piece of γ (0) (1) to fix the constant piece of the first-order equation. This determines γ (0) (1) , which is now independent of γ (1) (1) . We then solve for the non-constant piece, determining γ (1) (1) modulo a new constant. It is also easy to check that a similar argument will hold at higher orders in the expansion parameter .
The argument we have given here for a single abelian factor is easily generalised to the general polystable case. This shows that the D-term equations can be solved even when there are Fayet-Iliopoulos terms present. The reader may have expected that the D-term equations would obstruct some hermitian moduli when FI terms are present. This is not quite right; instead the expected obstructed moduli will form part of the gauge transformations and so should not be included in the first place. We give more detail on this in appendix E.

D.4 Full Maurer-Cartan equations
The arguments we have given for the massless deformations can be generalised to solutions to the full set of Maurer-Cartan equations: Exactly as before, x (1) and α (1) have a gauge symmetries that allow us to solve the D-term conditions to first order.

JHEP10(2018)179
At second order we have Here (y (1) , b (1) ) are the gauge-transformed fields that solve the first-order D-term conditions. Furthermore λ 1 is fixed by these first-order conditions. We note that the second order field y (2) is shifted by the first order gauge transformation by the term [λ 1 , y 1 ]. This means it is y (2) − [λ (1) , y (1) ] that appears in the second-order D-term condition, which then admits shifts by terms of the form Dλ (2) + 1 2 ∂ξ (2) which we may use to solve the D-term conditions just as in the first order case, fixing λ (2) and ξ (2) in the process. We can see that this structure will continue to hold to all orders. Assuming we have solved the D-terms at i th order where i < n, y (n) is shifted by gauge transformations involving the y (i) . We can then use the shifts of y (n) by (∂ + ∂)-exact forms to solve the D-terms at n th order. Solving the D-term conditions thus fixes the gauge symmetries (though there are residual gauge-of-gauge transformations).

E Massless moduli
In this appendix we clarify the meaning of the massless moduli This condition can be viewed as an equation of motion for the field b appearing in the superpotential (4.11). It is also a necessary condition for preserving the d-closure of the holomorphic three-form, which in general is a stronger condition then preserving an integrable complex structure alone. However, with the extra assumption that H (T (1,0) X) has a representative satisfying this condition [14].

JHEP10(2018)179
Next we must solve the equation This is then the familiar condition that the complex structure deformations µ a should be in the kernel of the Atiyah map where F µ is given by contraction with F as in (E.5). The fields α 0 now decompose as where α c 0 is closed and α µ 0 0 is uniquely determined by the complex structure moduli through equation (E.5). Modulo gauge symmetries, we have  It is easy to check that the right-hand side of both of these equations is ∂-closed. The condition imposed by (E.9) is that the sum should be exact. The number of remaining massless moduli is then given by | ker H| + | ker F α | + | im H ∩ im F α |. (E.15)

JHEP10(2018)179
We can understand this from considering gauge transformations by a D-exact term. For polystable bundles, part of these transformations are formed by sections of H 0 (End V C ). Again we have V C = i V (i) C for stable slope-zero factors V (i) C . Thus we have im F α = span{c 1 F 1 + . . . + c n F n }, (E. 23) for constants c i such that i n i c i = 0, where n i is the dimension of the i th bundle, so that the total section remains traceless. As we saw in the previous section, these are precisely the gauge transformations needed to solve the D-term conditions in the presence of Fayet-Iliopoulos terms, and so hermitian elements of this form should not be thought of as moduli.
Open Access. This article is distributed under the terms of the Creative Commons Attribution License (CC-BY 4.0), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.