Higgs Bundles for M-theory on $G_2$-Manifolds

M-theory compactified on $G_2$-holonomy manifolds results in 4d $\mathcal{N}=1$ supersymmetric gauge theories coupled to gravity. In this paper we focus on the gauge sector of such compactifications by studying the Higgs bundle obtained from a partially twisted 7d super Yang-Mills theory on a supersymmetric three-cycle $M_3$. We derive the BPS equations and find the massless spectrum for both abelian and non-abelian gauge groups in 4d. The mathematical tool that allows us to determine the spectrum is Morse theory, and more generally Morse-Bott theory. The latter generalization allows us to make contact with twisted connected sum (TCS) $G_2$-manifolds, which form the largest class of examples of compact $G_2$-manifolds. M-theory on TCS $G_2$-manifolds is known to result in a non-chiral 4d spectrum. We determine the Higgs bundle for this class of $G_2$-manifolds and provide a prescription for how to engineer singular transitions to models that have chiral matter in 4d.


Introduction
Geometric engineering is at the heart of many applications of string theory, starting with model building for particle physics, the study of superconformal field theories, or sharpening the boundaries of the string theory landscape. For many of these applications F-theory has been the framework of choice in recent years, culminating for instance in the classification of 6d superconformal field theories [1,2]. Clearly a similarly robust and comprehensive analysis would be desirable for four-dimensional models with minimal supersymmetry. In F-theory, the defining data of 4d N = 1 theories is however not purely geometric, unlike the 6d setup, but includes a choice of G 4 -flux. In particular, G 4 -flux is crucial in order to get chiral 4d N = 1 theories from F-theory [3][4][5][6][7].
An alternative framework that yields minimal supersymmetry in 4d is obtained from Mtheory on G 2 -holonomy manifolds (for a review see [8]). As is well known, the main challenge in this setup is the construction of compact G 2 -holonomy manifolds 1 with singularities, which 1 We shall often simply refer to these as G2-manifolds. yield both gauge (codimension 4) and chiral matter (codimension 7) degrees of freedom in 4d [9][10][11][12][13][14]. To this moment this is an open question.
Until recently, the number of known compact G 2 -manifolds was rather limited: the only concrete examples were the Joyce orbifolds given by resolutions of T 7 /Γ [15] and constructions based on orbifolds of a Calabi-Yau three-fold times S 1 . Recently, a comparatively large class of examples (order millions) of compact G 2 -manifolds was described in [16][17][18] as twisted connected sums (TCS).
The physics of M-theory and string theory on TCS G 2 -manifolds has been investigated in [19][20][21][22][23][24][25][26][27][28]. One key property common to all TCS manifolds, which is a direct consequence of this particular construction, is that singularities will occur (if at all) in codimension 4 and 6, but not 7. From the standard geometric engineering dictionary for G 2 -manifolds it then follows that the resulting models in 4d will not have chiral matter. An obvious question is then which type of deformations or singular transitions are required to remedy this limitation.
The present paper will provide a setting which gives some answers to this question and explores how such transitions would be characterized in TCS geometries, by providing a local model description in terms of a Higgs bundle. To achieve this, we first refine and extend the local model framework of [29], to incorporate the local limit of TCS G 2 -manifolds, and then determine the type of deformations that are required.
The approach of using local Higgs bundle models and their spectral covers in F-theory [3,5,[30][31][32][33][34][35][36][37][38] has proven very successful in model building, and more importantly as a precursor to the study of compact F-theory models. The Higgs bundles characterize the gauge sector of a compactification in terms of the local geometry in the vicinity of an ADE-singularity. In F-theory, this was not only useful in making the geometric engineering dictionary precise, but also subsequently in the constructions of compact geometries with favorable 4d effective field theories. For G 2 -manifolds the local structure close to conical singularities has been studied in [12,14]. Here we will take a slightly different approach, starting much like in F-theory with the statement that a local geometry that realizes in M-theory an ADE gauge group in 4d, will necessarily have a description in terms of an ALE-fibration over a compact supersymmetric cycle. In F-theory the local Calabi-Yau is an ALE-fibration C 2 /Γ ADE → M 4 , with M 4 a Kähler surface in the base of the elliptic Calabi-Yau four-fold. For a G 2 -manifold, the local model is analogously given by a fibration where M 3 is a supersymmetric three-cycle, i.e. an associative cycle, in the G 2 -manifold. This approach was advocated in [29], however much of the details of their paper remained somewhat ad hoc and more importantly, does not e.g. include the case of TCS G 2 -manifolds as we shall explain. We will both provide an in depth exploration of the Higgs bundle associated to this model, that will in particular lend itself to generalizations.
As M-theory compactified on an ALE space gives a 7d super Yang-Mills (SYM) theory with ADE gauge group, the effective 4d N = 1 theory of an ALE-fibration can be found by studying a topologically twisted 7d SYM-theory on a three-manifold M 3 . The BPS equations then determine the field configurations along M 3 that ensure that N = 1 supersymmetry is preserved in 4d. They are given in terms of a Higgs bundle specified by an adjoint valued one-form Higgs field φ and a gauge connection W along M 3 . We will focus entirely on diagonalizable Higgs fields, which implies that the connection W furthermore has to be flat.
The diagonalizability implies that we can equivalently describe the Higgs bundle in terms of its eigenvalues or spectral data.
The BPS equations imply that dφ = d † φ = 0 and so φ = df , where f is a harmonic function. This in turn implies that f is constant as long as we require M 3 to be compact and f to be regular. To obtain interesting solutions we introduce 'sources' or equivalently singularities for f , ∆f = ρ. Alternatively, we may excise the loci where sources are located and study the corresponding f -twisted Laplace equation on the resulting three-manifold with boundary M 3 .
In general the solutions to this zero-mode counting are difficult to determine. However, if we assume a fully factored spectral cover, the problem of finding the zero mode spectrum and interactions maps to Morse-Bott cohomology on M 3 . In this case the resulting 4d gauge theory has U (1)-gauge symmetries, which are determined by the number of factors of the spectral cover. The zero modes can then be computed in terms of relative cohomology of M 3 with respect to its boundary. The Higgs bundle spectral cover provides a construction of the three-cycles in the ALE-fibration, and determines the matter fields and couplings in 4d.
If the spectral cover is not fully factored we only have a formal description of the spectrum in terms of the cohomology of the f -twisted complex. This may be somewhat surprising for the reader more familiar with the F-theory spectral cover description, see e.g. [5,34,35,37], where the factorization is usually achieved with some amount of tuning (in order to have extra U (1) gauge symmetries) and the theories without this are usually simpler to describe.
In the G 2 -setting the factorization is paramount for even computing the 4d spectrum.
There is an alternative description -again in the case of fully factored spectral coversin terms of supersymmetric quantum mechanics (SQM), whose grounds states can be com- Extra two-form Factored C Charges ρ in Cartan U (1) symmetry Table 1: Dictionary between ALE-fibration/G 2 geometric data of the fibration over M 3 , Higgs bundle, Morse-Bott theory on M 3 (alternatively SQM), and the 4d N = 1 low energy effective theory.
puted using Morse (more generally Morse-Bott) theory as in Witten's classic work [39]. This characterization in terms of SQM identifies matter and couplings in terms of gradient flow trees in M 3 . A summary of the dictionary between ALE-geometry, i.e. local G 2 geometry, the Higgs bundle, Morse-Bott theory on M 3 or SQM, and the data of the 4d effective theory is provided in table 1.
This setup in particular allows modelling the local geometry of M-theory compactifications on TCS G 2 -manifolds, which have an ALE-fibration over S 3 (e.g. as in [24]). Moreover it will allow us -in the framework of the local Higgs bundle description of the geometryto make a concrete proposal for the types of deformations and transitions that the geometry needs to undergo. Although we necessarily lose the concrete description of the geometry offered in terms of a twisted connected sum 2 we may nevertheless track what happens to our model in the language of the local geometry, which may be useful in modifying/improving 2 Studying such transitions in a compact setting seems to go beyond the current tools available in geometry, as it can no longer be a TCS. However, see also the recent paper by Chen [40].
the TCS construction.
The plan of this paper is as follows: section 2 starts with a careful derivation of the partially topologically twisted 7d Super-Yang-Mills (SYM) theory on M 3 , which in turn determines the BPS equations. We then discuss solutions in terms of Higgs bundles, which characterize the local geometry, and discuss the spectrum of gauge and bulk matter. The spectral cover approach for these Higgs bundles is set up in section 3 and localized matter is studied in Section 4. A description of abelian Higgs field backgrounds in terms of supersymmetric quantum mechanics and its connection with Morse and more generally Morse-Bott theory is given in section 5. This setup is then applied to the study of matter couplings in section 6. Finally, in section 7 we apply this framework to describe the local models for TCS G 2 -manifolds and study the deformations of the associated local models. A summary of results useful for model building applications together with some concrete models is given in section 8. The reader predominantly interested in the rules of how local G 2 Higgs bundles are constructed for practical purposes can focus almost entirely on this section. We conclude with section 9, which furthermore contains a list of future research directions. A glossary of our notation and further technical details are relegated to the appendices.

The Gauge Theory Sector of M-theory on G -manifolds
M-theory compactified on a G 2 -manifold gives rise to a 4d N = 1 supersymmetric gauge theory with matter fields, coupled to supergravity. In this paper we will be interested in the gauge theories obtained from such compactifications and therefore will decouple gravity.
Gauge degrees of freedom in an M-theory compactification on a holonomy G 2 -manifold are localized on codimension 4 subspaces, which are associative (i.e. calibrated) three-cycles M 3 .
Locally the geometry takes the form of an ALE-fibration over M 3 as in (1.1). A useful way to characterize the gauge sector is to think in terms of the 7d SYM-theory obtained from Mtheory on the ALE-fiber: the gauge bosons in the Cartan subalgebra of the gauge group arise from dimensional reduction of the M-theory three-form C 3 on the two-forms in the ALE-fiber, and the remaining non-abelian gauge bosons arise from wrapped M2-branes. In an adiabatic approximation, where the ALE-fibration varies slowly over M 3 , the 4d effective action can be obtained by dimensionally reducing this 7d SYM-theory on the three-cycle M 3 , with a partial topological twist. In this section we carry out this reduction and determine the spectrum of gauge and matter fields, which are determined by solutions of BPS equations along M 3 (see (2.18)). The solutions are given in terms of a Higgs bundle over M 3 , that is specified by a one-form Higgs field φ and an internal gauge field W .

Partial Topological Twist and BPS Equations
We start with 7d SYM with ADE gauge group G. This theory can be obtained by dimensional reduction of the maximally supersymmetric 10d SYM on R 1,6 × T 3  The 10d Majorana-condition descends to a 7d symplectic Majorana-condition 3 . Denoting the gauge coupling in 7d by g 7 the action becomes whereγ denotes the 7d gamma matrices.
This 7d SYM theory is the starting point for the analysis of gauge degrees of freedom in a local G 2 -holonomy compactification of M-theory. For a given ALE-fiber, the singularity determines the 7d gauge group G. We now reduce this theory further on an associative threecycle M 3 . Since this will be generically curved with holonomy group SO(3), the 4d theory will in turn only retain supersymmetry if we partially topologically twist the local Lorentz group SO(3) M with the R-symmetry SU (2) R . Upon compactification on M 3 the local Lorentz 3 We refer to appendix A.2 for our conventions with regards to spinors and supersymmetry.
where P is a G-principal bundle. We denote the field strengths associated to the gauge fields A µ and the Wilson lines W i by F µν and (F W ) ij , respectively, and their associated covariant derivatives as D µ and D i . The latter can be combined with the scalars φ i , which both transform as a 3 of SU (2) twist , into a complex 1-form Note that ϕ,φ and D,D are related by conjugation in the gauge algebra. We further introduce and its conjugate Fφ = F † ϕ . We assume that the 4d gauge fields A µ are independent of the internal coordinates along M 3 , so that the latter two expressions become standard space-time derivatives of the complex scalars ϕ,φ Define the interaction term The partially twisted 7d SYM action is then (2.13) The supersymmetry variations for the bosonic fields are and for the fermionic ones we find To obtain a 4d supersymmetric theory upon twisted dimensional reduction, the field configuration along M 3 needs to preserve supersymmetry. We further require the background to enjoy 4d Poincaré-invariance and therefore require it to be independent of the coordinates The BPS equations are then obtained by setting δλ = 0 and are where the first equation is obtained by setting the real and imaginary parts of δχ to zero separately. 4d Poincaré invariance requires F µν = 0. Rewriting (2.17) with respect to the notation in (2.6) the BPS equations become the F-and D-term equations Background values for the Higgs field φ and gauge field W along M 3 that solve these equations will determine the effective field theory in 4d 5 . In components the BPS equations are On the left hand side the three-cycle M 3 is shown, with the charge distribution ρ that is located along Γ. On the right hand side, a tubular neighborhood T (Γ) is excised and the resulting manifold is M 3 .
We will generally assume this to be the case. The remaining equations are D W φ = D † W φ = 0. If M 3 is a compact three-manifold without boundaries and φ is regular, there are two cases to consider: In the first case φ has to be a harmonic 1-form and thus must be trivial, in the second case it can be non-trivial.
We will be interested in simply-connected three-manifolds (i.e. a homology three-spheres) in the following. To nevertheless have non-trivial solutions we relax the assumption that φ is regular, which can be achieved by including sources into the D-term equations. Writing φ = df , the function f is then required to satisfy Poisson's equation where ρ models the sources supported on a closed subset Γ. This maps the solution of the BPS equations to an electrostatics problem with the identification f = electrostatic potential ρ = charge density, supported on Γ ⊂ M 3 . (2.23) Alternatively this system can be described by excising a tubular neighborhood T (Γ) of the charge support Γ, and studying the problem of finding solutions on M 3 = M 3 \ T (Γ) -see figure 1. In this case φ needs to be regular, with suitable boundary conditions along ∂M 3 .
In summary we are going to consider the following setup which will be used in section 4 to determine physically interesting solutions to the BPS equations including localized matter. Localized matter is characterized by the vanishing of φ. When f is Morse (i.e. it has no degenerate critical points) these are isolated points and we will discuss this setup in section 4. By relaxing the constraint of f only having isolated critical points this can be generalized to situations where f is a Morse-Bott function and higher-dimensional matter loci can be included as well. We will discuss this in section 5.4 and apply it to TCS G 2 -manifolds in section 7.
In the remainder of this section, we assume that φ is non-trivial and regular, but make no further assumptions on the details of the loci φ = 0.

Higgs Bundles
Before studying the low energy effective theory, let us briefly recall the relation between the Higgs bundle and the local ALE-fibration. The BPS equations (in the absence of sources) (2.18) are in fact precisely the odd dimensional analogue of the Hitchin equations for the Higgs field φ giving rise to the data of a Higgs bundle. In the case [φ, φ] = 0 there is an elegant geometric description of the Higgs field φ in terms of an ALE-fibration over M 3 , which we now summarize [29]. This construction is analogous to the one in F-theory, where the Higgs field specifies the unfolding (a complex structure deformation) of the ALE singularity and is closely connected to the compact Calabi-Yau underlying the F-theory compactification [5,30,34]. Recently this was developed also for Spin(7) manifolds [41]. In our case the Higgs field describes the deformations of the full hyper-Kähler structure of an ALE fiber.
Recall that φ is an adjoint valued 1-form Ω 1 (M 3 ) or a section of T * (M 3 ), and we take it to be non-trivial along the commutant G ⊥ of the 4d gauge group G in The Higgs field is i.e. φ lives in a local geometry in the vicinity of M 3 which is the total space of the cotangent bundle T * (M 3 ). This is a local Calabi-Yau threefold. Since [φ, φ] = 0, we can diagonalize the Higgs field to obtain n 1-forms φ j , where n is the rank of the Lie algebra g ⊥ of G ⊥ . To locally recover the ALE-fibration over M 3 associated to this Higgs field, we use the Kronheimer construction [14,42]. Every ALE-space is of the form C 2 /Γ ADE , where Γ ADE is a finite subgroup of SU (2), which are classified by the corresponding ADE Dynkin diagrams. The second homology of the resolution of singularities of C 2 /Γ ADE is isomorphic to g and we can think of the components φ j as measuring the periods of the hyper-Kähler structure forms.
More explicitly, over a local patch of M 3 we can write the fibration as R 3 × C 2 /Γ ADE . We chose a basis σ j of H 2 (C 2 /Γ ADE , Z) and fix a hyper-Kähler triple (ω I , ω J , ω K ). The 1-form φ j can be written as where we identify This uniquely defines the hyper-Kähler structure on each fiber. Observe that the Higgs field has an SO(3) symmetry arising from the SO(3) acting on ω I , ω J and ω K .
In geometric terms we can describe our situation as follows. For simplicity, assume that we have a G ⊥ = U (1)-valued Higgs field φ. We are considering a local model for a G 2 -manifold with ADE-singularities located along an associative submanifold M 3 , which physically means that gauge degrees of freedom are localized along M 3 and the gauge group is given by the ADE type of the singularity. Consider the gauge group G, which by turning on a non-trivial background vev for φ generically higgses to G → G × U (1). This means that the ALE fiber over a generic point of M 3 will have the singularity corresponding to G via the ADE correspondence and there will be a two-cycle in the U (1) direction with non-zero volume, given by φ. Over the points where φ = 0, the two-cycle collapses and the ALE singularity worsens; equivalently the gauge group enhances from G to G. We will elaborate this point in section 6.
We can in fact make the local geometry of the gauge enhancement fairly explicit. For the moment let us restrict our attention to the case where G = SU (2) which corresponds to a C 2 /Z 2 singularity over M 3 . Giving a non-trivial background vev for φ corresponds to deforming the generic fiber to a smooth Eguchi-Hanson space. More precisely, consider the generator σ of H 2 ( C 2 /Z 2 , Z) of the resolved geometry. Recall that σ is topologically a twosphere. From (2.27) and (2.28) we see that at a generic point x ∈ M 3 (which for this purpose is approxated locally by R 3 ) we have by which we mean the volume of σ in the Eguchi-Hanson space above x. Consider now a neighborhood of a non-degenerate zero of φ, which we can assume to be at 0 ∈ R 3 . We can The signs depend on the eigenvalue of the Hessian at 0. The Higgs field φ now has an isolated zero at the origin. The explicit local description of the ALE-fibration is given by Viewing X as a fibration over R 3 all of the fibers are smooth apart from the fiber over (0, 0, 0) i.e. the zero of φ. Moreover, X is a cone in C 3 × R 3 with the apex at the origin. The link of the cone can be found by intersecting X with the unit sphere in C 3 × R 3 and is in fact P 3 realized as the twistor bundle over S 4 . The approximate G 2 -metric on X is given by where η is the 2-form dual to the two-cycle σ.
This can be generalized to arbitrary ALE-fibrations. The local geometry is of the form , with a C 2 /Γ G fiber over the origin. We again work with G = SU (n + 1). For the deformations of other ADE singularities see [43]. The topology in a neighborhood of an isolated zero is This describes a family of SU (n) singularities, with enhancement to SU (n + 1) at the origin (note that we again write φ = df as above). There are also explicit deformations for other ADE groups. Topologically X is now a cone over the weighted projective space P 3 n,n,1,1 with coordinates (y 1 , y 2 , y 3 , y 4 ) [14]. In the link, there is a family of SU (n) singularities along an S 2 given by y 3 = y 4 = 0. In the ambient space, the location of the singularities is a cone R + × S 2 = R 3 , which is identified in our context with a local patch of the base R 3 of X. As before, the apex of the cone is where the cycle σ collapses to zero volume. This therefore establishes a key piece of the dictionary between properties of φ and the ambient G 2 -geometry. The isolated zeroes of φ give rise to conical singularities of the ALE fibered G 2 -manifold. As we show in section 5, this fits together nicely with the physics side as zeroes of φ which occur at codimension 7 are precisely the loci where chiral fermions are localized.

Massless Spectrum
Given a solution to the BPS equations (2.18) with regular Higgs field we can ask what the spectrum of the 4d gauge theory is. The equations of motion of the fermions follow from 34) which are equivalent to the decoupled equations where the adjoint is taken with respect to the Hermitian inner product Tr (ᾱ ∧ * β) . (2.37) Acting on functions g ∈ Ω 0 (M 3 , ad P ) and written in coordinates, e.g. the operator∆ becomes where we pick up a conjugation due to the inner product. We find that (2.35) may be rewritten where by (2.8), χ and ψ are 0-and 1-forms, respectively. Massless modes are therefore described by the kernels of the Laplacians ∆,∆ or equivalently by closed and co-closed forms with respect to the operators in (2.36) value of φ or W oriented along a subgroup G ⊥ of G breaks the gauge group to its commutant G ⊂ G. The adjoint fermions ψ, χ will decompose accordingly to give matter valued in irreducible representation. In this higgsed theory the fermions are sections of the associated gauge bundles, E. The action of D restricts to each of these subbundles allowing us to make the identification (2.41) We next rewrite these cohomology groups with respect to the same co-boundary operator by dualising H 0 D , H 1 D with the Hodge star. Note that by (2.37) we have D † = * D * and D † = * D * so that taking χ α ∈ H 0 D (M 3 , E) for example we find that * χ α is annihilated by the operators D, D † This precisely states that * χ α ∈ H 3 D (M 3 , E), i.e. we have mapped fromD-cohomology to D-cohomology using the Hodge star. The same observations hold true forψ α . The Hodge star relates This allows us to make the following identifications where now all cohomologies are with respect to D and forms of all degrees are employed.
Note that the Z 2 -grading of the exterior algebra aligns with the 4d chirality of the fermionic zero-modes. The Hodge star depends on the metric of M 3 which itself is induced from the metric of G 2 -holonomy of the ambient 7d manifold.
Since M 3 is associative and so calibrated with respect to Φ ijk we equivalently could have used the G 2 3-form Φ ijk to dualize since it restricts to a volume form of M 3 . Contracting elements of H 0 D and H 1 D with the 3-form Φ ijk is then exactly the same as taking their Hodge dual.

Bulk Matter
The first type of matter we will discuss arises from a background Higgs bundle, where φ = 0, which solves the BPS equations, but W = 0 with F W = 0. This will be referred to as bulk matter, as the modes will not be localized. We will see that for π 1 (M 3 ) = 0 there is no chiral index for this matter type. It may be interesting to extend this to non-trivial π 1 setups, which we relegate to future work, and also has been discussed in earlier works from a different point of view (see e.g. [44]).
Turning on a flat gauge field along a subgroup G ⊥ ⊂ G the gauge group G is Higgsed to the commutant G of G ⊥ in G and the adjoint representation of G decomposes as (2.45) For the fields of the theory this decomposition is lifted to the bundle level, where Ad(P ) decomposes into the vector bundles R n ⊗ S n in the representations R n and S n of G and G ⊥ , respectively. The chiral and conjugate-chiral zero modes transforming in R n are then counted by the cohomology groups Chiral : Their CPT-conjugate zero modes in R n are obtained by Hermitian conjugation in the gauge algebra or equivalently from (2.41) with E =S. In order to rewrite these cohomology groups with respect to the same boundary operator D we again dualise H 0 D , H 1 D using the Hodge star and obtain Chiral : These cohomology groups completely determine the chiral and conjugate-chiral spectrum in 4d transforming in R n of the remnant gauge symmetry G Chiral fermion zero-modes : (2.48) The chiral index of the representation R n is which is nothing other than the Euler characteristic of the D-complex. In the case of trivial fundamental group π 1 (M 3 ), there is no flat bundle to break the gauge group, and reduce to the Betti numbers of the de Rham complex on M 3 . The chiral index is then given by the usual Euler characteristic, which vanishes for odd dimensional closed manifolds This concludes our discussion of 'bulk' matter. In the following we will focus our attention on localized matter modes, which arise from non-trivial φ background values. Since these are best characterized in terms of spectral covers we will first develop the framework for that.
We will briefly discuss interactions between bulk matter and localized matter fields later on.

Defect Description of Matter
Thus far our discussion was based on starting with a 7d SYM theory on M 3 with a gauge group G which is generically broken to a smaller subgroup by the Higgs background. Over the zero locus of the Higgs field some of the gauge symmetry is restored. An equivalent description starts with a bulk 7d SYM with gauge group G × U (1), and an unhiggsing to G by inserting defects at points in M 3 . We discuss this mechanism for rank 1 enhancements.
The starting point is the off-shell formulation of the 7d SYM given in appendix B.1, now with gauge group G × U (1). We take the Higgs field to have a background turned on along the abelian directions as φ = df t where t denotes the U (1) generator. The permitted configurations for f are again determined by the BPS equations, including sources. Whenever the Higgs field vanishes the gauge symmetry could potentially enhance. For this we need to extend the field content by the required degrees of freedom. We therefore add defects, coupled to the bulk fields as where we have introduce a chiral multiplet valued in the representation R q or R −q of G × U (1) depending on the choice of sign +, − in (2.51) respectively. Here K is the auxiliary field of the chiral multiplet. The multiplet Λ has a fixed profile along M 3 which is determined by the condition that Λ descends to a massless 4d field upon reduction and will be subject of section 5. This also fixes the choice of sign in  Each λ i is a one-form and their vanishing thus implies that they intersect the 'zero section', i.e. the locus of the ADE singularity of type G present in every ALE fiber, over points p k on M 3 . Those points p i are precisely the loci where matter is localized. In the generic case of non-factored spectral covers, the sheets are furthermore connected by branch-cuts (orange). As we will see later on: In case the point p is connected by a flow line in M 3 to another critical point, there is a corresponding associative three-cycle which is built by fibering the collapsing S 2 (blue) over the flow line. The resulting contribution to the superpotential gives a mass term for the localized states.

Spectral Cover for the Higgs Field
For the case when a higher rank Higgs bundle is turned on but the Higgs field commutes, it is useful to describe the solution to the BPS equations in terms of the spectral data of the Higgs field. This framework is of course very familiar from F-theory spectral covers, see e.g. [5,34,35,37], and for the Lagrangians in Calabi-Yau threefolds and the associated G 2 -manifolds with pointlike singularities was touched upon in [29]. Here we will prepare the setup to also account for more general Higgs field configurations, with the goal to apply it to the TCS-manifolds.
Recall that φ is a section of Ω 1 (M 3 ) ⊗ Ad(G ⊥ ). For concreteness let G ⊥ = SU (n). For a commuting Higgs field we can choose to diagonalize it and study the resulting spectral cover where b i are symmetric polynomials in the eigenvalues of φ and for SU (n), b 1 = 0. The eigenvalues λ i are one-forms which give rise to an n-sheeted cover of M 3 and where b 0 is the zero-section. A cartoon of this is given in figure 2. Using the spectral cover the associated ALE-fibration is simply Each λ i parametrizes the volumes of a corresponding two-sphere in the G ⊥ -ALE-fiber. The gauge symmetry G is generically higgsed to G, except at the loci when the gauge symmetry enhances. Since λ i is a one-form, this condition implies that this happens generically over points in M 3 , though we will encounter other situations as well.
The relation between the eigenvalues λ i and coefficients in the ALE-fibration b i is not linear, and generically the sheets of the spectral cover will be connected by branch-cuts. This effect implies in particular that the U (1)-symmetries associated to the Higgs bundle are not actually present in the low energy effective theory.
The classic example for spectral cover models starts with an E 8 → SU (5) × SU (5) ⊥ . The spectral cover is five-sheeted and λ i = 0 characterizes the location of 10 matter representations (we refer the reader to the F-theory literature where these models have been studied in depth [5,34,35,37]). We will construct an example of this kind in detail at the end of this paper.

U (1) Symmetries
Generically the sheets of the cover are connected by branch-cuts and therefore, although locally it may appear otherwise, the independent gauge group is G, the commutant of G ⊥ = SU (n) in G. If however the coefficients of the spectral cover are tuned such that it globally factors over M 3 then this corresponds to N independent U (1) factors in the 4d effective theory [37]. The possibilities of factorization depend on the monodromy group that acts on the spectral cover, which for SU (n) covers is S n . If the group acts transitively on the n sheets then there is no additional U (1) symmetry. If it has N + 1 orbits then there are N globally defined two-forms, which define U (1) symmetries. To see this, we consider the difference between the factored cover components C (k) − C (l) . Fibered over M 3 , the associated two-cycles define a non-trivial five-cycle in both the local model and in the compact G 2 manifold J. The Poincaré dual two-form to this five-cycle then gives a U (1) gauge boson in the Kaluza-Klein reduction of the three-form C 3 in compactification of M-theory. This can be also be seen concretely in the context of TCS G 2 , see section 7.
Perhaps the only important difference to the F-theory spectral cover story is that here it will be paramount that the spectral cover factors, in order to determine the spectrum of the 4d theory. Although the general twisted cohomology will continue to compute the zero mode spectrum, we do not have a computational tool to determine these cohomologies, unless the spectral cover factors completely.

Localised Matter
We will now study the more interesting and richer class of matter fields, localized on points or one-dimensional loci of M 3 . So far in section 2.4 we considered only flat gauge fields on along M 3 , which corresponds to bulk matter. Turning on vevs for the Higgs fields φ will enlarge the possible matter structure and will allow us to engineer spectra with non-trivial chiral index. The simplest case is an abelian Higgs field configuration where G is the 4d gauge group and U (1) ⊥ the commutant, along which the Higgs fields is turned on. The expectation is that since φ is in the 3 of SO(3) twist , the condition for local gauge enhancement to G occurs at codimension 3 in the base M 3 , i.e. codimension 7 in the G 2 -manifold J. This is also suggested by the earlier spectral cover discussion. We will discuss this case of codimension 7 localized matter first. In less generic situations, such as the twisted connected sums, however, enhancement occurs at codimension 6 loci.
The solution to the BPS equations (2.24) on M 3 will be constructed by excising a tubular neighborhood T (Γ) of a graph Γ, with boundary conditions, which we will discuss in detail.
The central question is how the zero modes in R +q and R −q are counted. In this section we provide the cohomological answer to this question, which applies to both codimension 6 and 7 gauge enhancements. In the next section we will provide specific solutions to the BPSequations, to which the general analysis in this section can be applied, thereby computing the zero mode spectrum.

Zero Modes from Relative Cohomology
We now turn on a background value for the Higgs field φ, which to begin with is U (1)-valued.
As explained in section 2, we now set out to solve the D-term equation (2.18) for φ = df with sources, i.e. the Poisson equation where the charge density ρ satisfies charge conservation We take ρ to be localized on links Γ i in M 3 of definite signs of the charges, Γ ± , Both the Higgs field φ and f diverge along Γ. We again excise a tubular neighborhood as in section 2. The boundary ∂M 3 splits into connected components Σ i , which correspond to the connected components of the underlying links Γ i and correspondingly the boundary splits as The normal derivatives of f , which are computed with respect to the outward pointing unit normal vector fields, have to be positive (resp. negative) restricted to Σ + (resp. Σ − ).
The zero modes of the fields in the representation R q and R −q in the presence of a background Higgs vev φ = df are obtained from the twisted Laplacian where and H f is the Hessian of f . Furthermore we defined the raising/lowering operators Note that D † is not necessarily adjoint to D on manifolds with boundary Σ as Requiring appropriate boundary conditions fixes this problem. Consider a form α split into its tangential and normal component to the boundary The tangential part α t is defined as the pullback of α to the boundary and the normal part as α n = α − α t . The boundary contribution is sensitive only to the tangential components i.e.
where we have used the fact that ( * α) t = * α n . The two types of boundary conditions are Dirichlet : where we set again D matter (∆ f ) = D matter (∆). Note that the d-complex and D-complex are isomorphic (cf. appendix C), so they have isomorphic cohomology groups. In this case, the d-complex is restricted to forms which vanish on Σ − . This computes the relative cohomology of the pair (M 3 , Σ − ) [46] so we get The sign of the U (1)-charge q is important. Changing it amounts to changing the sign of f , which inverts the signs of normal derivatives and consequently exchanges the boundary conditions imposed on the positive and negative boundaries, and we obtain the cohomology groups with respect to the positive boundary. In terms of the operators defined in section 2.3 changing the sign of q corresponds to computing the cohomology with respect to the operator D, which is isomorphic to the D-cohomology but this time with respect to the conjugate representation R.
Returning to the analysis of the spectrum above, we have seen that it is computed by the relative cohomology with respect to the negatively charged boundary components. (4.16) The net amount of chiral matter transforming in the representation R is therefore given by the relative Euler characteristic We have seen that for an M 3 without boundary there is a 4d vector multiplet in the spectrum. Once we introduce sources along Γ and excise a tubular neighborhood around them, we need to check that the vector multiplets remain in the spectrum. Since these adjoint fields are uncharged under the U (1), the associated forms cannot have any tangential boundary conditions, and we impose purely normal boundary conditions. In this case the domain of the relevant Laplace operator becomes The kernel is then isomorphic to the de Rham cohomology groups [45] and we obtain the required zero modes for the vector multiplets in 4d.

Higher Rank Higgs bundles
Next we generalize to higher rank Higgs bundles in G ⊥ . We still assume that [φ, φ] = 0. If we cannot diagonalize the Higgs bundle globally, (i.e. in the spectral cover language the spectral cover does not fully factor) then we still have a local description in terms of the Higgs field along the CSA, but not globally: i.e. we can only diagonalise locally in a patch U of M 3 . Here f i , ρ i : U → R are functions, n = rk G ⊥ and t i the generators of the CSA. Locally this background breaks the gauge where Q = (q 1 , . . . , q n ) denotes a vector of U (1)-charges. If the spectral cover has N + 1 irreducible components (as in (3.5)), N of these n U (1) factors descend to the gauge group of the 4d effective theory. The operator D defined in (2.36) acts on R Q by (4.23) Let us introduce The zero-modes are counted by (2.44) where E = Ad G ⊥ . If the spectral cover does not factor, i.e. the sheets mix under monodromy, the cohomologies of the operator D cannot be rewritten in terms of e.g. de Rham cohomologies. For the case of rank 1 Higgs bundles the isomorphism given between the corresponding complexes was given by conjugation with e qf (this is explained in greater detail in appendix C). This required a globally defined function f whose role for fully reducible Higgs bundles is played by f Q as we will explain in the next section. This isomorphism cannot be adapted in a straightforward manner to general Higgs bundles.
Restricting D to Ad G or Ad(U (1) n ), it is reduced to the exterior derivative Table 2: The 4d N = 1 matter content for a background given by a U (1) n valued Higgs bundle whose spectral cover is fully factored. Here Σ − Q denotes the negative boundary of M 3 with respect to the function Vector and chiral multiplets transforming in these representations are thus simply counted by the zeroth and first Betti numbers of M 3 , respectively.
However, if the Higgs bundle diagonalizes globally, i.e. if we have rank G ⊥ many U (1) symmetries, then a simple generalization of the rank one case applies. The zero modes are counted with respect to where f Q is globally well-defined and a function. As a consequence the results of section 4.1 carry over upon making the replacement qf → f Q . M 3 is obtained by excising the singularities of all the f i and the boundary decomposes again into positive and negative parts depending on whether f Q → ±∞ when approaching the excised charge. By i.e. for differently charged representation R Q each zero mode counting requires an alternate decomposition of the boundary. We therefore find the fermionic zero-mode spectrum in the representation R Q to be enumerated by the relative Betti numbers   Figure 3. Examples of charged graphs in M 3 . Positive, negative charges are coloured red, blue respectively. Both charge distribution give rise to the same chiral index but a different number of zero-modes.

Example 1: Wires in S 3
We now turn to describing concrete charge configurations -these configurations were studied in [29] and we revisit them here. Let M 3 = S 3 and embed charges in S 3 which are localized on a graph Γ. The positively and negatively charged components of the graph are disjoint Γ = Γ + ∪ Γ − . We denote by n + , n − the number of components and by + , − the number of loops of Γ + , Γ − respectively. The total charge on Γ ± is again constrained to vanish. Excising tubular neighbourhoods of Γ ± we obtain M 3 with associated boundaries Σ ± . By where r counts the number of negative loops which are independent in homology when embedded in M 3 \ Γ + . The chiral index is then computed to be It solely depends on the charge configuration Γ and is independent of the number r. A chiral spectrum is therefore easily generated. Multiple charged graphs will give rise to the same spectrum. Another point to note here is that a non-trivial chiral index will only arise if for some sign of the charge, the number of loops and components is different, i.e. the charge distribution is not localized solely on a disjoint union of circles. This will later on give hints as to how to deform the Higgs bundles for TCS G 2 -manifolds.

BPS-Configurations, Super-QM and Morse-Bott Theory
In the discussion above we were not interested in any particular features of the harmonic function f on M 3 and the computation of the spectrum is in fact valid for any such f . In this section we first specialize to the case where the Higgs field φ = df has isolated, nondegenerate zeros, which is the same as requiring f to be a Morse function -this is the case already studied in [29]. We show that massless chiral matter is localized at the zeros of the Higgs field. We then generalize this to the case where f can have critical loci of dimension one, in which case it is Morse-Bott. The latter will be essential for the TCS geometries.
The main tool here is reformulating the problem of finding the kernel of the Laplacian ∆ f in terms of a supersymmetric quantum mechanics and Morse theory. This approach is useful as it lends itself to the generalized Morse-Bott setup that we are interested in.

Matter, Morse and (Super-Quantum-) Mechanics
Let us consider again the abelian case where φ = df with f harmonic in the decomposition (4.1), which counts the fermionic zero-modes transforming in the representation R q , that are in the kernel of The twisted Laplacian ∆ f can be interpreted as the Hamiltonian of a supersymmetric quantum mechanics (SQM) with the target space M 3 where the supercharges are given by the operators D andD [39]. In section 2.3 we have shown that (due to the partial topological twist) the state space is identified with the space of differential forms on M 3 . However, since M 3 is now a manifold with boundary, we have to restrict the state space to forms satisfying the boundary conditions given in (4.13), which we denote by H = ⊕ i Ω i b (M 3 ). The subscript b indicates that the forms satisfy the boundary conditions. The function f now plays the role of a superpotential and the kernel of ∆ f characterizing the true zero modes in the reduction to 4d is now enumerating the supersymmetric ground states of the SQM [39]. In summary: is discrete and all points p ∈ N are non-degenerate. A critical point p ∈ N is called nondegenerate if its Hessian H f (p) is non-degenerate as a bilinear map. In this case p ∈ N is assigned a number µ(p) called the Morse index given by the number of negative eigenvalues In the case of manifolds with boundary, we further assume that there are no critical points of f on ∂M 3 . Note that this is true in our case, since the normal derivative of f at the boundary is non-zero (see section 4.1). For more details on Morse theory we refer the reader to [47,48].
We can choose a normal coordinate system in which f and the metric g on M 3 take the form where we assumed that p = 0 and c i are the eigenvalues of the Hessian, which due to the harmonicity of f sum to zero. This means that only points with Morse index 1 and 2 can occur. Expanded in these coordinates ∆ tf reduces to the Hamiltonian of a supersymmetric harmonic oscillator with Solving for the ground states of the harmonic oscillator locally, near a critical point p of Morse index µ(p), we find a unique solution given by a differential form of degree µ(p). The zero modes of ψ, which are identified with 1-forms in (2.44), localize at critical points of Morse index 1. For c i with signature (−, +, +), the solution to leading order is Likewise the fermions in R −q are obviously counted by replacing f with −f .

Exact Spectrum from SQM
The perturbative calculation in the previous section does not necessarily give the exact spectrum of the full theory. On the SQM side this is due to the fact that quantum mechanical instanton corrections can cause perturbative ground states to acquire a mass and be lifted in the full theory [39,47]. We now complete the dictionary between the 4d effective theory of 7d SYM and SQM by showing that masses of perturbative zero modes in the 4d theory arise precisely from instanton corrections on the SQM side.
We start our analysis with the off-shell action in (B.9) and split the complex 1-form ϕ = ϕ 0 + δϕ into its background ϕ 0 = tdf and fluctuations δϕ. The 7d fields are expanded in terms of a basis of perturbative ground states of the twisted Laplacian as where (x, y) ∈ R 1,3 × M 3 . Here the sum runs over the charged representations, R q and R −q , and all critical points p a of Morse index 1 with respect to the relevant Morse function, f and −f respectively. The fermionic field ψ (a,q) (x) carries the anti-commuting, gauge and 4d spinor structure while ψ (a,q) (y) is a 1-form on M 3 annihilated by the twisted Laplacian in perturbation theory. In leading order in t these are (5.6) or the CPT conjugate of (5.7). The decompositions for δϕ are of analogous structure.
A mass term in 4d descends from the 7d interaction which for an abelian Higgs background yields the mass matrix This precisely computes the instanton corrections between the perturbative ground states in SQM theory and is simply the matrix element Let us briefly summarize the classic results on these instanton corrections, see [39,47] for a detailed treatment. The (Euclidean) action of the SQM with target space M 3 is given by a standard sigma-model action where g ij is the metric on M 3 , D the covariant derivative and R ijkl the curvature tensor.
Canonically quantizing this action, one gets the SQM we have described in the previous section [39]. The matrix element (5.11) now has the the following path integral expression which is the gradient flow equation. With this the mass matrix is evaluated in [47] to leading order in 1/t as Here the sum runs over all ascending gradient flow lines γ starting at p b and ending at p a . The contribution from a flow line γ is weighted by a sign n γ = ±1, which arises from a choice of orientation on the moduli space of gradient trajectories. The precise derivation from the SQM context is intricate and is given in [49,Appendix F]. The main takeaway is that perturbative ground states form a complex, where the coboundary operator is given by It is possible that the boundary operator of the Morse-Witten complex is trivial. This is equivalent to a vanishing of the mass matrix M ab = 0, i.e. all perturbative ground states are true ground states. In this case the Morse function f is called perfect. This is precisely the We can consider these mass terms also in the M-theory picture. In section 2. M2-instanton wrapped on such a three-sphere S 3 γ reduces to SQM [51,52]. The stationary points of the M2-brane action, which correspond to associative three-cycles, hence become a fibration of the vanishing cycle of the ALE-fiber over the gradient trajectories γ(f ) determined by the Morse function f . These associatives then give a non-perturbative correction to the superpotential [51,52] which is of the form ∆W = n γ exp i We thus find a rather simple criterion to determine whether the true ground state spectrum of the theory is chiral or not: Two examples are shown in figure 5. This result is of course recovered from the more general charge distributions discussed in section 4.3 upon setting the number of loops l + and l − to zero. In particular for generic placements of the n + + n − charges one has Each critical point thus constitutes a true ground state and we recover (4.29). This is made explicit in figure 5. If flow lines between critical points exist, they always do so in pairs with n γ = ±1. Hence the corresponding ground states are not lifted.

Generalized Critical Loci and Morse-Bott Theory
The setup studied in [29] and in the last section assumes that the critical loci of the function f are isolated points. Although this is the generic situation, it will be important to relax this assumption and consider the generalized setup in which the critical locus of f can be one-dimensional, which happens for the recent TCS constructions of G 2 -manifolds. Functions f with critical loci of dimension greater than zero whose Hessian at its critical closed submanifold is non-degenerate in the normal direction are called Morse-Bott functions. An example is given in figure 6. For further background on this see [47,53]. The case where N is a point has been analyzed in section 5.1.
If N = S 1 we can proceed analogously. As N has index 1, f is locally of the form in a suitable normal coordinate chart centered at a point p ∈ N . In this coordinate system x 3 is the coordinate tangential to N and the Hessian H f is diagonalized with the eigenvalues c and −c. In these coordinates the twisted Laplacian (5.1) now takes the form The analysis of perturbative ground states thus splits into normal and tangential parts relative to N . In the normal direction we get a single 1-form solution ψ ⊥ given by Here we have split ψ ⊥ into a 4d Weyl spinor ψ (n,q) carrying the anti-commuting, gauge and spinor structure and its internal profile normal to N . In principle ψ ⊥ is defined only locally on N . However, observe that ψ ⊥ is a volume form on the fiber of ν − N . Hence, assuming that the negative eigenbundle ν − N is orientable, the local solutions can be patched together to a global form on N . Since f is constant on N the tangential equation reduces to a Laplace equation on S 1 . Let the coordinate on the circle be θ. Then we obtain two solutions we find the localized perturbative ground states on every critical circle contributing to the perturbative spectrum to assemble to two chiral multiplets transforming in R q and R −q .

Generalized Critical Loci and SQM
We now turn to the computation of the exact spectrum from the perturbative solutions in the Morse-Bott case, where the critical loci of f consist of points and circles. While it is possible to compute the SQM instanton correction in much greater generality [47,53], the applications for TCS local models allow us to consider only the set-up with this restriction.
The instanton calculation in this case effectively reduces to the one considered in section 5.2.
To find the exact spectrum, we again want to compute the matrix element (5.15) between perturbative zero modes localized at critical submanifolds we use the analogous SQM computation. Let N m denote the disjoint union of critical submanifolds of Morse-Bott index m (recall that this is the dimension of the negative eigenspace of the Hessian matrix). In our case, m can take the values 1 or 2. For m = 2, all of the components of N 2 must be points, whereas N 1 can contain points as well as circles.
Recall that among the ground states localized at critical circles there are chiral multiplets transforming in the representation R q and R −q . As already discussed in section 5.4, this is because perturbative ground states are of the form with deg(ψ ⊥ ) = 1 and α a harmonic form on N 1 . When N 1 is a circle, α can be a function or a one-form. Consider again the matrix element Here we again use the indices a and b to enumerate all the perturbative ground states of total degree 2 and 1 respectively. However, note that for Morse-Bott functions the index is no longer in one-to-one correspondence with critical loci since there are two perturbative ground states localized at each critical S 1 ⊂ N 1 . For the following we will require the assumption that there are no ascending gradient flow lines between connected components in N 1 . 6 To compute M ab we need to consider three cases. First, bothψ (a,q) and ψ (b,q) may be localized at points in which case the discussion of section 5.2 applies verbatim. We now turn to the second possibility, where the ground states are both localized at the same circle critical circle S 1 ⊂ N 1 . The matrix element is then given by the integral where we have used the explicit expression of for such ground states given in (5.30). Using the expression for ψ ⊥ in (5.29) one can see that dθ ∧ * (df ∧ ψ ⊥ ) = 0 and also dθ ∧ ψ ⊥ ∧ * dψ ⊥ = 0.
This implies that the matrix element M ab is zero, ifψ (a,q) and ψ (b,q) are both localized at the same circle.
The third possibility is thatψ (a,q) is localized at a point p a in N 2 and ψ (b,q) is localized at a circle S 1 b ∈ N 1 . To keep track of all of the gradient curves between critical loci of f , we 6 In this case f is said to be weakly self-indexing. This assumption can be avoided at a cost of making the exposition much more technical [53].
introduce the moduli space of gradient trajectories between N m and N n Here, the quotient is taken with respect to the remaining reparametrization invariance of the gradient flow: γ(t) → γ (t) = γ(t + δt). The moduli space M(N m , N n ) is a smooth manifold, and it follows from simple dimensional analysis that its dimension is m−n−1. An illustrative example is given by S 2 with the Morse-Bott function f (x, y, z) = z 2 , see figure 6.
For our purposes, the only relevant case is m = 1 and n = 2 in which case the moduli space is a finite set of points. This means that there are finitely many gradient trajectories connecting N 1 and N 2 and there are finitely many ascending gradient flow lines connecting S 1 b and p a . We can now continue with the computation. In terms of the SQM path integral we have the expression . This is nearly the same expression as in (5.13), with the only difference being that we integrate over all curves with γ(−∞) ∈ S 1 b . However, the same localization argument as before applies. As we have seen above, the number of gradient trajectories is still finite and the result of the path integral computation has exactly the same form as for points, i.e. (5.15). The expression for the operator D also remains unchanged The exact spectrum is given as the cohomology of D, which acts on the following complex This complex is a convenient way to arrange all the perturbative ground states of degree p in C p . It is a specific instance of a Morse-Bott complex for f , which can be defined for f with critical loci of arbitrary dimension [53]. If f is a solution to the electrostatics problem in section 4.1, the Morse-Bott cohomology again recovers the relative cohomology of a pair (M 3 , Σ − ).

Chiral Index from Spectral Covers
We close this section by introducing yet another picture for counting the perturbative zero modes, namely using the spectral cover introduced in section 3. For certain configurations it is possible to read off the exact spectrum using the spectral cover, this was already observed for the U (1) case in [29].
For simplicity let us begin by recalling the statement for the rank 1 Higgsing in (5.18) where G = SU (n + 1). There we turned on a single abelian Higgs background parametrised by the Morse function f via φ = df . The spectral cover C in this case is simply the graph of To compute the perturbative spectrum we thus just need to count the intersections of the different sheets with their signs as in the rank 1 case above. Let C i ⊂ C denote the sheet of a spectral cover C with Graph(df Q i ) = C i then where the notation is as in (5.38). Similarly we compute the chiral index to Perturbative zero modes transforming a representation R Q i which is not associated by λ i = df Q i to a sheet of this spectral cover are enumerated by the intersection of the different sheets The chiral index again given by the difference This is pictorially most clear in the case of A n singularities. In this case the ALE-fiber is given by a circle fibration over R 3 and the eigenvalues λ i , which are characterised by the sheets of the spectral cover, correspond to the points at which the circle collapses. A vanishing sphere is stretched between any pair of these points and collapses whenever they come together, i.e.
when the sheets intersect. This enhances the spectrum and constitutes an additional ground state.

Yukawa Couplings and Higher-Point Interactions
In this section we discuss the interactions of bulk and localized matter. It will be useful to consider the case of a fully factored spectral cover, in which case we can compute the zeromodes. In M-theory interactions between localized matter fields come from M2-instantons wrapped on calibrated 3-spheres of the local ALE-fibration. This is simply a generalisation of the results of section 5.2, where we interpreted non-perturbative mass terms as arising from M2-instantons wrapping three-cycles which connect two critical points over a gradient flow line. For higher point interactions these three-cycles project to gradient flow trees on M 3 and studying the moduli space of these constrains the corresponding interactions in 4d.
Corrections to these couplings are obtained from integrating out states with masses induced by M2-instantons as discussed in section 5.2.
Consider again the background as in section 4.2 globally on M 3 : The matter content is summarized in table 2.

Bulk-Localized-Matter Interactions
First consider the bulk-localized-localized interactions. These interactions are present for rank 1 abelian Higgs background. The contribution to the N = 1 4d superpotential is canonically derived by expanding the partially twisted 7d action  expansions in ordinary harmonics of the bulk fields, denoted by Φ k ,Φ k . The interaction term is then We now turn to interactions of the localized matter fields only.

Yukawa Couplings
For Yukawa couplings we need a rank n = 2 Higgs bundle (or higher). There are two Morse functions f 1 and f 2 and the combination f Q = q 1 f 1 + q 2 f 2 . From the effective field theory we obtain this coupling by expanding the action (6.2) in perturbative zero modes where (a, p) refers to the internal profile of the perturbative zero mode localized at the critical point p a transforming in R Qp . The Yukawa couplings arise from M2-instantons wrapping associative three-cycles. To characterize the three-cycles consider the Morse functions which describe an SU (3) ALE-fibration over the base M 3 . Each of the functions f i controls the volume of a corresponding two-sphere α i in the ALE fiber, which satisfy in the homology of every fiber. Recall that α i shrinks to zero volume precisely over the points p i where df i = 0. To every gradient trajectory γ(f i ) starting at a point p i we can associate a 3-chain, which is given by tracing out the corresponding α i in the ALE-fibration.
Given three sufficiently generic Morse functions f i , there will be finitely many gradient flow trees connecting the three critical points p i (see figure 7). Adding the associated three-chains produces a three-cycle, the boundary of which is given by i α i in the ALE fiber. We may produce a closed three-cycle with the topology of a three-sphere by adding a three-cycle β such that ∂β = α 1 + α 2 + α 3 . Moreover, this S 3 is in fact associative, since it projects to the tree of gradient trajectories and hence minimizes the volume among all the three-cycles which project down to trees connecting p 1 , p 2 and p 3 . Wrapping an M2-brane on such a cycle gives rise to Yukawa couplings between modes localized at the critical points of f 1 , f 2 and The strength of these interactions is governed by the choice of functions f i . The threesphere giving rise to the Yukawa coupling is a supersymmetric rigid homology sphere within the G 2 -manifold and its contribution to the superpotential is again given by (5.17). The sign n γ = ±1 arises in the same manner and is given by an orientation on the moduli space of gradient flow trees. As the Higgs field φ i and the gauge field W i are identified with the periods of the supergravity 3-form C and associative 3-form Φ the integral is evaluated as Here, we have used that we can gauge the background for the gauge field W i to zero. Evaluating the final integrals and using that the homological relation between the α implies 3 (6.9)

Associatives and Gradient Flow Trees
The generation of Yukawa couplings and mass terms from associative three-cycles which project to flow trees on M 3 has a natural generalization [54], which in the effective theory realizes higher point couplings.
We consider a setup in which G ⊥ = S[U (1) k ], so that the Higgs background is described  11) and there are thus finitely many gradient flow trees connecting k points of Morse index 1. An example of a gradient flow tree for the case of k = 5 is shown in figure 9.
As before, we can construct a three-cycle by fibering the two-sphere α associated with the Morse function f over each segment γ(f ). This both guarantees that we end up with an associative, and also that α collapses at the end-points of the flow tree. Furthermore, the fact that we have a tree in M 3 implies that the resulting associative three-cycle has the topology of a three-sphere, so that it contributes to the effective superpotential. Using the same manipulations as in (6.8), we can compute the volume of such a 3-sphere γ as Vol so that the resulting contribution to the superpotential is The scale M φ is set by the vev of φ (see our discussion of this in section 8.1). Note that there can in general exist several flow trees connecting matter localized at the same loci p i , which can cancel out. induced by associates over flow lines between two points, we may integrate out the corresponding massive fields, thereby generating higher-dimensional operators in the effective field theory. As M Inst M φ these corrections are dominant compared to the couplings induced between the same fields by associatives. An example is shown in figure 10.
For non-generic choices of the charge distribution the moduli space of gradient flow lines may increase and (6.11) is no longer valid. In this case the moduli space of gradient flow lines is not discrete but of dimension 1 and isomorphic to a circle. In the ALE geometry this corresponds to a continuous family of associative submanifolds. In [52] it is shown that the contribution of such a family C of associatives is proportional to χ(C).
In the more generic set-up of unfactored spectral cover of rank n the Higgs background can only be diagonalised locally as in (4.20). The source terms ρ in the BPS equations are now oriented arbitrarily along G ⊥ breaking the gauge symmetry to its commutant G in G.
Assuming we can diagonalising the Higgs field in a tubular neighbourhood T of the singularity α is a 0, 1-form and localize at these critical points of Morse index 1, 2 respectively. We are further assuming that the states α ∧ ψ ⊥ and η i are of unit norm. After this perturbation, the previous analysis applies. The true ground state corresponding to α ∧ ψ ⊥ is Finally, let us apply these observations to the TCS constructions. In [26] a chain of string dualities was used to argue for the existence of infinitely many associatives on a class of TCS G 2 -manifolds, and this result was recovered in an orbifold limit in [28]. These associatives are furthermore in one-to-one correspondence with elements of the lattice E 8 ⊕ E 8 . The local limits of these models must be such thatG = G ⊥ = E 8 × E 8 , and the associatives argued for in [26,28] must correspond to flow trees in these local models. In fact, the description of the associatives in terms of string junctions in [26,28] is already deceptively close to our description in terms of flow trees. It would certainly be interesting to flesh out this correspondence in detail.

Higgs Bundles and Twisted Connected Sum G 2 -manifolds
In this section we consider local models associated with twisted connected sum (TCS) G 2manifolds, which form the largest known class of examples of compact G 2 -manifolds [16,17].
The TCS construction has by now been covered extensively in the literature, so we will only briefly recapitulate the main points and refer the reader to [16,17,24] for further details. In a nutshell, the power of the TCS construction is that it shows how compact G 2 -manifolds can be glued from simpler building blocks, which can in turn be constructed using algebraic geometry.
Although this makes finding examples relatively straight-forward, TCS G 2 -manifolds appear to be a rather special class within the set of all G 2 -manifolds [24]. Our analysis of local models for G 2 -manifolds allows us to move away (at least in local models) from the TCS description and explore how to connect TCS G 2 -manifolds to G 2 -manifolds giving rise to chiral spectra.

TCS G 2 -Manifolds
The basic ingredient for the twisted connected sum construction is a pair of algebraic threefolds Z ± , which each admit a K3 fibration with generic K3 fiber S ± . The manifolds Z ± have to satisfy i.e. the first Chern class of Z ± must be equal to the class of a generic K3-fiber. With some further assumptions on the topology (see [17,Definition 3.5]) Z ± are then called the building blocks. Excising a generic fiber S 0 ± from Z ± one obtains a pair of non-compact threefolds X ± = Z ± \ S 0 ± , fibered over a punctured Riemann sphere, which are asymptotically cylindrical (aCyl) Calabi-Yau threefolds. Away from a compact submanifold, the X ± have the topology of the cylinder R + × S 1 b,± × S 0 ± and the Ricci-flat metrics on X ± asymptote to the Ricci-flat product metric on this cylinder. The situation is sketched in figure 11.
To form a compact G 2 -manifold, the aCyl Calabi-Yaus X ± are then multiplied by an extra circle S 1 e,± and glued together along a their cylindrical regions. The diffeomorphism used for the gluing exchanges the 'internal' circles S 1 b,± with the 'external' circles S 1 e,∓ and identifies the K3 surfaces S 0 ± by a diffeomorphism which induces a hyper-Kähler rotation, or Donaldson matching, Re(Ω 2,0 ± ) = ω ∓ Im(Ω 2,0 + ) = − Im(Ω 2,0 + ) . Here, Ω 2,0 ± and ω ± are the holomorphic (2, 0) forms and Kähler forms on S 0 ± which are induced by the complex structures on X ± . The compact topological manifold J resulting from this gluing then admits a metric with holonomy G 2 , which is close to the Ricci-flat metrics on X ± × S 1 e,± . More precisely, there exists a limit, which we will call 'Kovalev limit', in which the cylindrical region becomes arbitrarily long and in this limit the Ricci flat G 2 -holonomy metric approaches the Calabi-Yau metrics on X ± . In compactifications of M-theory, modes localized only on X + ×S 1 e,+ or X − ×S 1 e,− give rise to subsectors with enhanced N = 2 supersymmetry in four dimensions. These subsectors are coupled such that they mutually only preserve N = 1 supersymmetry, and we may think of the length of the cylindrical region as the inverse of their coupling [22,25].
As both X ± × S 1 e,± are fibered by K3 surfaces and the gluing acts separately on the fiber and base, J is (topologically) fibered by K3 surfaces as well. The base of this fibration is a three-sphere S 3 glued together from two solid tori. We can hence think of the local models associated with TCS G 2 -manifolds as describing an ALE space which is cut out from the K3 fiber over a base space M 3 which is S 3 . To engineer non-abelian gauge groups, every ALE fiber of the local geometry and hence every K3 fiber of the associated compact G 2 -manifold must be singular. It is straightforward to construct acyl Calabi-Yau threefolds in which every K3 fiber has a singularity of ADE type and the work of [24,25] suggests that gluing such singular three-folds indeed results in a singular G 2 -manifold.
Let us consider this in more detail. Denote the image of by N ± . The Donaldson matching implies an idenfication of H 2 (S 0 + , Z) H 2 (S 0 − , Z). Using this map, every element of gives rise to an associated harmonic two-form on J: the Poincaré dual cycle to such a form is algebraic for both S + and S − , so that its fibration over the whole base S 3 of J is trivial and it sweeps out a five-cycle, which is Poincaré dual to a two-form on J. The number of independent such two-forms on J is simply given by the rank of g [17]. In compactifications of M-theory on J, there are hence |g| massless U (1) vectors from the Kaluza-Klein reduction of the three-form C 3 8 .
The hyper-Kähler structure on S 0 ± is forced by the Donaldson matching to be such that the integral of both Ω 2,0 ± and ω ± vanishes for every cycle contained in g. This means that whenever there is a root, i.e. a lattice vector of length −2, contained in g, the K3 fibers S 0 ± are singular. As g sits inside of the polarizing lattices 9 of the algebraic families X ± , this implies that every single K3 fiber has a singularity. The type of singularity can be read of by finding the sublattice g root ⊂ g generated by the roots of g. This sublattice must be a (sum of) ADE root lattice(s) and its type determines the corresponding singularity and the resulting simply-laced 10 non-abelian gauge group upon compactification of M-theory. 8 There are in general further massless U (1) vectors associated with classes in the kernel of ρ± [17], which associated with the irreducible components of reducible fibers of the K3 fibrations on Z±. 9 The polarizing lattice of a family of K3 surfaces is the sublattice of H 2 (K3, Z) which is orthogonal to Ω 2,0 for all members of the family. 10 While this data is sufficient to find the singularities associated with simply-laced gauge groups, it is slightly more tricky to find non-simply laced gauge groups. Their emergence in TCS G2-manifolds parallels their emergence in F-theory [55] in that the exceptional divisors of resolutions of ADE singularities of S± may globally become a single divisor in X± [24]. In terms of lattice data, this can be expressed by saying that a cycle of self-intersection n < −2 contained in g can force an ADE singularity in every fiber if it is a linear combination of −2 curves in S+ or S− which are all in the polarizing lattices of the families S+ and S−. The difference between the polarizing lattices and N± determines the 'folding' of the ADE Dynkin diagram.
The matter loci in these models arise as the degeneration loci of the singular K3-fibration i.e. where the singularity worsens. This happens over points in P 1 ± , each of which gets multiplied with a circle in the TCS construction. This implies that in M-theory compactification on a TCS manifold J, matter is localized along circles. This is true at least in the Kovalev (stretched neck) limit in which the metric on the J is well approximated by the metrics on each of the building blocks, which can be thought of as being contained inside J (more precisely, the products X ± × S 1 ± are in J).

Higgs Bundles of TCS G 2 -manifolds
We start by considering the local models of the two building blocks individually. As the discussion is the same for both sides, we will drop the ± subscripts. The first step is to replace the K3 fibration with a local ALE model. The precise details of this local limit depends on the ADE group corresponding to the type of ADE singularity, and are well known in the literature [56,57]. Besides an ADE singularity, every ALE fiber contains a number of compact cycles, the volumes of which vary over the base. Such cycles may collapse over points in the base C.
At these loci the singularity present in the generic fiber is enhanced and matter is localized.
be such cycle which vanishes at (some) of the corresponding points in the base C. Let us denote the hyper-Kähler triple by Θ = (ω I , ω J , ω K ). In terms of (7.4) the hyper-Kähler structure is simply ω = ω I Ω 2,0 = ω J + iω K . (7.8) After taking the local limit and integrating over σ we get a meromorphic function φ on C: with zeros precisely where σ shrinks to zero volume. The poles and zeros of φ are located away from ∞. Moreover, we can identify φ with the meromorphic (1,0)-form as φ dz. Since the base C is contractible φ = df , where f is now a Morse function with critical points of index 1 and singular loci corresponding to the poles of φ. After taking a product with the circle we trivially get a Morse-Bott function.
If (unit) charges are placed at points a i ∈ C, the function φ will be of the form Therefore φ can have at most n − 1 critical loci, which is generically the case. If we impose charge conservation on each side, there can be at most n − 2 critical loci of φ.
With this information let us now consider how the Higgs bundle for TCS manifolds. After gluing C + × S 1 e,+ with C − × S 1 e,− , the base manifold is M 3 = S 3 . In the Kovalev limit, the critical locus of the harmonic Morse-Bott function f consists of a disjoint union of m circles of Morse index 1. As before, we may engineer such an f by an appropriate configuration of charges Γ on S 3 . On C ± × S 1 ± , these charges will simply be given by a collection of points on C ± times the circle S 1 ± . From the above discussion we only need the simple observation that matter loci in TCS G 2 -manifold, at least in the Kovalev limit, are circles. Suppose that there are m matter circles and no points. Using the results of section 5.5 we see that the Morse-Bott complex is 11) and the cohomology gives just We find that every perturbative ground state constitutes a true ground state, the Morse-Bott function f is thus perfect. As each circle gives rise to a pair of chiral and conjugate-chiral zero modes upon Kaluza-Klein reduction, the spectrum associated to this Higgs field configuration φ = df is non-chiral We can use this result to derive constraints on the function f . By the above results the relative Euler characteristic χ(M 3 , Σ − ) = 0 vanishes and by Lefschetz duality we find that this implies χ(M 3 , Σ + ) = 0. We obtain the topological constraint There has been a recent attempt to modify the TCS construction to yield singular G 2manifolds with codimension 6 singularities by Chen [40]. Instead of smooth bulding blocks Chen takes the Z + building block to be a threefold with isolated nodal singularities, which means that the non-compact aCyl G 2 -manifold X + × S 1 e+ has singularities in codimension 6. However, the standard TCS gluing argument does not work in this case; rather it is conjectured [40] that if circles of nodal singularities are replaced by pairs of isolated conical singularities it is possible to glue to a G 2 -manifold with conical singularities using a modified version of the connected sum construction. In terms of the local model, the collapse of circles into points corresponds to deforming the Morse-Bott function to a generic Morse function, where the same collapse of critical circles to critical points occurs (recall that critical points correspond precisely to isolated singularities of the total space of the G 2 -manifold). However, even if this conjecture is true, such G 2 -manifold will still give rise to a non-chiral spectrum by the arguments above.
Finally, let us discuss the spectral covers for a TCS G 2 -manifold which is given to us in terms of building blocks X ± and a gluing map where S ± have ADE singularities over U ± ⊂ C ± . This gluing of the K3 fibers in the TCS geometry also implies a consistent gluing map for the ALE-fibrations associated with the local model. In general, to be able to glue two given ALE-fibrations together, the two ALE-spaces need to be of the same type,G + =G − , and furthermore the periods of the ALE-fibers must satisfy a matching condition. By Torelli theorem for ALE-spaces [58], the structure of an ALE-space is completely determined by the periods of the hyper-Kähler structure forms over the 2-cycles in the root lattice of the algebra of G. Explicitly, the matching condition is where σ j are the 2-cycles generating the root lattice. Note that this implies that the nonabelian part of the group G, i.e. the type if ADE singularity, must be the same on both sides.
Each X ± furthermore has a local model, which is a Higgs bundle φ (±) over C ± , and a corresponding spectral cover C (±) . The asymptotic values the Higgs fields φ (±),0 are similarly related by which induces a gluing of the spectral covers.
Let us explain the origin of U (1) gauge symmetries in glued spectral covers. Each cover C (±) can have a factorization structure, which defines two-forms (five-cycles) and locally U (1) symmetries. This can be detected by the restriction of the map (7.5) to the ALE-fibrations over C ± . Factorization of the spectral cover C over S 3 after gluing C (±) can likewise be detected by (7.6), and only those two-forms in the image that lie in the intersection will globally give rise to a two-form and thereby a U (1) symmetry. An example is shown in figure C + C - Figure 12. Each of the building blocks X ± defines a spectral cover over C ± , and these are then glued to a spectral cover over S 3 . In the example shown here, the cover on C + factors into two components and the cover on C − factors into three components. These covers are glued such that the resulting spectral cover over S 3 has two components with two sheets each. Hence the resulting model has and there is a single unbroken U (1).
12, where C (±) → C ± each is a four-sheeted cover. However C (±) is factored into two (three), and thus locally gives one (two) U (1) symmetries. The gluing is such that the spectral cover C → S 3 has only two factors, and thus only gives rise to a single U (1) symmetry. The scale at which the other U (1) is broken is set by the size of the neck region of the TCS-construction.
Besides an analysis via Higgs bundles, the matter spectrum of M-theory on TCS G 2manifolds can also be found using a purely geometric reasoning. The geometry in the vicinity of each matter locus is that of a Calabi-Yau threefold times a circle. The local Calabi-Yau geometry is that of a fibration of an ADE singularity over C with a further degeneration at a point. Using the usual dictionary between singularities and gauge theory for M-theory or type IIA on Calabi-Yau threefolds, the Cartan generators and weight vectors can be identified with exceptional divisors and curves in the resolved Calabi-Yau geometry [5,43,59,60]. Our analysis of Higgs bundles now implies that the multiplicities must be such that each matter locus gives rise to a single vector-like pair of representations. Furthermore, we may determine the U (1) charges by simply integrating the two-forms in g which give rise to the U (1)s over the exceptional curves of the resolution associated with the matter.

Deformation of TCS Higgs Bundles
Given the local model for TCS G 2 -manifolds we now consider the behavior of the physics under deformations of the Morse-Bott function. We have seen above that circular critical loci arise in the non-generic S 1 -invariant distributions of charges in S 3 which are present in the Kovalev limit. The natural question is what happens if this invariance in is broken by a slight deformation. The strategy we will use to describe deformations is to exploit the A configuration of charges which produces the Morse-Bott function associated with a TCS G 2 -manifold in the Kovalev limit must of course be finely tuned, as a generic configuration of charges will always result in critical loci of dimension zero. Let us discuss this in a simple example -see figure 13: consider a charge distribution of two equally charged coplanar and cocentric circles in R 3 . This setup has rotational symmetry and correspondingly the critical locus is another coplanar and concentric circle. A generic deformation will destroy the rotational symmetry and lead to critical points instead of a circle. Consider e.g. deforming the charges to ellipses while preserving coplanarity. This collapses the critical locus to two points of Morse index 1 near the vertices of the ellipses and two points of index 2 near the co-vertices.
More generally, the function f will become Morse with isolated critical points for a generic deformation. However, since the topology of Σ ± does not change, we still have Physically this means that any deformation of f will give rise to chiral spectrum if under the deformation the topology of Σ − remains unchanged. Denoting the number of points with Morse index i by m i , the Morse inequalities for manifolds with boundary imply Equally, every deformation of the local model of a TCS G 2 -manifold that has an associated charge distribution which consists of a number of circles will satisfy n ± = l ± , so that the resulting spectrum is seen to be non-chiral.

Chirality and Singular Transitions
It is not at all surprising that TCS G 2 -manifolds do not give rise to chiral spectra and that small deformations do not change the chiral index. However the result we have found already has fairly interesting geometrical implications: for a generic small deformation of a TCS G 2manifold, the loci at which matter is localized are no longer one-dimensional but become point-like. This of course implies that the product structure of X ± × S 1 e,± must be broken and the periods of the hyper-Kähler triplet on the K3 fiber must have a non-trivial dependence along S 1 e,± . Although such small deformations will not yield G 2 -manifolds giving rise to chiral spectra, the crucial ingredient, which are point-like singularities, is already present for small deformations of TCS G 2 -manifolds.
Engineering the ALE-fibration from a Morse function which in turn is determined by a configuration of charges allows us to make the key observation for how to deform TCS G 2manifolds to situation with chiral spectra: we need to make a transition after which n ± = l ± no longer holds. The simplest way to do so is to bring two circles of equal charge together and then deform them to an object with l = 2. For TCS G 2 -manifolds, there are essentially two different ways to achieve this.
The first option is to take e.g. a positive charge on C + and another positive charge on C − , bring them together, and fuse them as shown in figure 14. As now l + − n + = 1 while l − − n − = 0, the resulting spectrum must be chiral. In a generic situation in which f is Morse, i.e. f only has isolated critical points, the critical locus of f hence consists of an odd number of points now. As we started from a non-chiral configuration with an even number In fact, the type of transition we have just sketched can also be anticipated from the heterotic duals of TCS G 2 -manifolds, which are given by compactifications on the Schoen Calabi-Yau threefold with different vector bundles [24]. Such models always have non-chiral spectra and a singular transition connecting the Schoen Calabi-Yau threefold to a different Calabi-Yau threefold (together with appropriate vector bundles) is needed to find a chiral spectrum. The Schoen Calabi-Yau threefold can be described as a fiber product of two dP 9 s, and it allows singular transitions in which a singular fiber of one dP 9 is collided with a singular fiber of the other dP 9 . As discussed in [24,26], the duality to a TCS G 2 -manifold implies that the singular fibers of these two dP 9 s are separated into disjoint regions of the common P 1 base. A collison between singular fibers from both ends translates to a collision of singular K3 fibers coming from the two separate ends X + and X − of the dual TCS G 2 -manifold.
The second option is to change the charge configuration corresponding to a TCS G 2manifold by colliding two circles of equal charge which are both located in the same building block. The picture of such a deformation will be similar to the one in figure 14, however initially the charged circles will be unlinked. Again, it is clear that this signals a departure from a TCS G 2 -manifold (and must result in a singular transition on the heterotic side as well): after the transition e.g. X + × S 1 e,+ must become a non-compact G 2 -manifold without the structure of a product.

Higgs Bundles for G 2 s: A User's Manual
We will now give a user-friendly summary of how to build local Higgs bundle models for G 2 manifolds, stripping off most of the mathematical baggage and condensing it to the essentials, which might be useful for the practitioners in the field.

Scales
Let us briefly discuss the mass scales in the problem, and specify what scale separation gives rise to the decoupling of gravity. For this purpose consider the gauge theory on the associative three-cycle M 3 on which the gauge degrees of freedom are localized. The compactification geometry determines the size of the cycle M 3 and the volume of the G 2 -manifold J and the characteristic size R ⊥ of the directions transverse to M 3 in J is hence In terms of these length scales, the Plank masses M 4 in 4D and M 7 in 7D are given by The scale M φ at whichG is broken to G is set by the an appropriate average of the Higgs background. As the volumes of the compact cycles in the ALE fiber over M 3 are set by φ , the local limit corresponds to the limit in which we can decouple gravity from the gauge theory degrees of freedom: Approximating the effective physics by the gauge degrees of freedom is only valid if the two scales M φ and M 7 can be decoupled and permit the limit M 7 → ∞ while keeping M φ constant.
Keeping φ fixed, this limit is equivalent to shrinking M 3 inside of J. Finally, the coupling of the 4D gauge theory G is given in terms of

Matter Content and Interactions
Here we summarise the construction of some simple backgrounds and the resulting matter content and interactions.
1. Choose the rank n of the Higgs bundle, which is equal to the number of Cartan generators t i along which a Higgs field background df i has been turned on. Each of these n abelian directions is sourced by a charge distribution ρ i of different type which determines the background function f i : M 3 → R completely, and must integrate to zero on M 3 . The Higgs bundle is therefore given by 2. This background breaks the gauge symmetry G → G × U (1) n and determines the count of representations in ⊕ Q R Q , where Q = (q 1 , . . . , q n ). To count the zero modes in R Q we need the effective charge distribution and its corresponding potential At every point in M 3 where df Q = 0, there is a localized chiral multiplet transforming in R Q . For every flow line governed by f Q between a pair of such points, there is a mass term for the associated chiral multiplets. Given a charge configuration ρ Q , the resulting massless spectrum can be described in terms of the numbers n Q ± of positively and negatively charged component, and the total number Q ± of loops.
The massless spectrum is counted by Here r Q denotes the number of negatively charged loops which are independent in homology when embedded into M 3 \ ρ + Q .
3. Interactions between three 4d chiral fields localized at points p s transforming in R Qs can only arise if Q 1 + Q 2 + Q 3 = 0 and if there exists a trivalent gradient flow tree between them. In general, there can be several such flow trees and cancellations between them can occur. Furthermore, there are mass terms for some of the localized zero modes and one needs to integrate out those massive fields to find the Yukawa couplings between the massless fields.

Retro-Model-Building 1: Top Yukawa
We close this section with two retro-inspired model building applications. First we will consider the top Yukawa coupling in an N = 1 SU (5) toy model. We take M 3 = S 3 and G = E 6 , i.e. two abelian directions parametrised by f a and f b are turned on. The corresponding decompositions read The effective Morse functions are 10 = 4f a .
(8.10) and they are determined by two independent charge distributions ρ a and ρ b . Using the linear combinations ρ 1 = −ρ a − 3ρ b and ρ 2 = 4ρ a we can write the charge vectors Q s as In our model, we distribute n  (10,10,5)  The top Yukawa was already engineered in the local G 2 spectral cover in section 8.3. To break the GUT group the standard model gauge group can e.g. be achieved by turning on discrete Wilson lines. We will not discuss this here, but it would be interesting to incorporate this into the Higgs bundle framework, see also [44] for a discussion.
There are numerous ways in generating the chiral matter content of an SU (5) GUT model in the Higgs bundle setup. We will choose one that is minimal and requires only point charges. Furthermore to incorporate the superpotential couplings we will consider a Higgs field with G = E 8 and φ taking values in The Higgsing that we consider is a special case of  in which The way to parametrize the charges is in terms of the embedding of a U (1) 4 into the Cartan subalgebra (CSA) of SU (5), where we have five generators t i , which satisfy Any U (1) α can be parametrized now as a linear combination The charges under the CSA generators (t 1 , · · · , t 5 ) of the matter fields of the SU (5) GUT are Matter Field U (1)-Charges under CSA Spectral Cover We also include the loci where these matter fields are localized in M 3 in terms of the λ i , To realise this spectrum via electrostatic charge distributions we translate the above into the language of section 5 and its higher rank generalisations. Note from (8.15) that 4 factors of U (1) have been broken off, this gives us 4 types of charge with which to build the model. Due to this special abelian background the decomposition in (8.14) decomposes further and the fundamental weights of SU (5) ⊥ are now associated with fundamental charge vectors Q F i such that (8.19) now becomes Upon a redefinition of U (1) generators we may take these charge vectors to be

Conclusions and Outlook
The main result of this paper is a study of the gauge sector of M-theory compactifications on G 2 -holonomy manifolds to 4d N = 1 supersymmetric gauge theories. There are various concrete directions to build on the present work: 1. Using the analogy with electrostatics not only allows to (implicitely) construct abelian Higgs field backgrounds, but can furthermore be used to find the zero mode spectrum.
In the case of genuinly non-Abelian Higgs field, i.e. in situations in which the spectral cover does not factor completely, this method can not be straightforwardly applied.
For model building applications one hence needs both a concise way to specify such solutions and a efficient way to determine the resulting spectrum of zero modes. bundles on the SYZ fiber as a generalization of [63]. For TCS G 2 -manifolds this has been done in [24].
We realise the gamma matrices as Here σ 0 = Id 2 and the 4d gamma matrices are where 4d signature is R 3,1 . The chirality, B-matrices, charge conjugation matrices and Lorentz-generators will be denoted by: respectively for 4d, 7d, 10d gamma matrices. and As one of the defining properties of the B-matrices is B * B = 1 there strictly speaking does not exist a matrix B 7 . However the definition of the B-matrices as product of all imaginary gamma matrices can be extended to odd dimensions. We thus define B 7 and its corresponding charge conjugation matrix C 7 by does not exist  The three B-matrices fit together as We collect the relations satisfied by the above matrices in where indices run as µ = 0, . . . , 3 and k, l = 1, 2, 3 andî, = 1, 2, 3.
The 10d Majorana-condition B 10 λ = λ * leads to a symplectic Majorana-constraint on the 7d spinors and a Majorana-constraints on the 4d spinors. We trace through the decomposition of the spinor representation as detailed in (2.1), (2.4) and (2.6) and make these constraints explicit.
We next turn to the spinors λ ααα in 4d of the untwisted symmetry group. By (A.9) these are required to satisfy the Majorana-constraint where the indices run as α, α,α = 1, 2. There are no further constraints on the spinors.
After performing the twist we find the Dirac spinors (λ 0 , λ i ), which carry twisted indices, to be constrained as with i = 1, 2, 3. Decomposing these Dirac spinors into Weyl spinors λ 0 = (iχ α , iξα) and λ i = (ψ iα ,ζα i ) the conditions are rewritten explicitly with the charge conjugation matrix as 16) which are, due to the introduction of a factor of i, nothing but two canonical Majorana- Using this we can rewrite the 4d dimensionally reduced action of 10d SYM interms of 1+3 unconstrained Weyl spinors χ, ψ i in 4d. The resulting action is given in (2.13).

B.1 Off-shell action for 7d SYM
Following [30] we rewrite the twisted 7d SYM action into a mostly off-shell form. We introduce complex auxiliary fields H i packaged in 3 chiral multiplets and a real auxiliary scalar field D contained in a vector multiplet in WZ gauge The N = 1 supersymmetry algebra can not be implemented off-shell for a vector multiplet in WZ gauge. Either manifest gauge invariance or manifest supersymmetry need to be demitted. It is possible to maintain the reduced algebra and manifest gauge invariance at the same time. This is achieved if we choose the bosonic variations to be of the form and fermionic variations (2.15) to be of the form where δ α ( · ) = [Q α , · } andδα(·) = [Qα, · }. The variations of the auxilary fields H i , D are then given by We next derive the action invariant under the above transformations which after integrating out the auxiliary fields reduces to the on-shell action previously obtained via dimensional reduction and twisting (2.13). To obtain an action annihilated by δ α ,δα it suffices to construct a real action which is annihilated by either of the two operators. The other operator will annihilate the action as the operators are conjugate to another. A real action must thus be built from exact and closed terms with respect toδ 2 orδα which are annihilated byδα. give rise toQ-exact and closed terms Tr − 1 4 F µν F µν − D µφk D µ ϕ k + iD µ χσ µχ − iD µ ψ k σ µψk (B.12) Note that the term giving the Pontryagin density can be dropped as it is topological, therefore supersymmetric on its own and thus not necessary for the supersymmetric invariance of the action. The auxiliary fields are eliminated via their equations of motion This returns the on-shell action (2.13) and the corresponding supersymmetric variations (2.14),(2.15).

B.2 Effective 4d Action
In section 2.3 we discussed the zero modes upon reduction to 4d and enumerated these in section 4. In section 6 we studied the 4d N = 1 superpotential terms governing the interaction of localized modes of the 4d N = 1. In this appendix we give more details on the reduction to 4d. We begin by reproducing the assumptions on the Higgs background of section 6.
We take the Higgs background to be turned on along n abelian given by the Cartan which are uncharged under all factors of U (1) and will be referred to as uncharged or bulk fields.
Then there are b 1 (M 3 ) Ad G and Ad U (1) n chiral multiplets respectively denoted as (B.18) We further have chiral multiplets valued in every representation R Qr which we will referred to as charged or localized fields. These we denote as Φ (a,r) = ϕ (a,r) + √ 2θψ (a,r) + θθH (a,r) + · · · , (B. 19) where the index denotes the representation via the corresponding charge vector Q r the multiplet transforms in and a denotes the critical point of Morse index 1 of the function f Qr to which the multiplet is associated. The index a omits critical points of Morse index 2. We think of all manipulations in the gauge algebra as embedded in Ad G.
The 4d Lagrangian organizes itself into three parts. The first part are the standard kinetic terms in 4d N = 1 superspace, the second are the superpotential terms to which bulk fields contribute and the third are the superpotential terms to which only the localized fields contribute. The latter two were focus of section 6.
Expanding the action (B.12) in zero modes the first part of the 4d Lagrangian is found to be The second part of the 4d action can be inferred from the terms and its Lagrangian is found to take the form We conclude with a comment on scales. The internal zero modes along M 3 are dimensionless as can be read off from the 7d action derived in section B.1. The 7d gauge coupling g 7 thus contributes the mass scale and we find M ab Qp , which denotes the mass terms for the conjugate representations R Qp , R −Qp in analogy to (5.32), to be associated to M φ (exponentially suppressed by the volume of three-cycles responsible for the mass term). Of course Y abc pqr is dimensionless. The masses M ab Qp thus contribute a new mass scale M Inst which is linked to M φ via the average volume Volγ of the 3-cycles in the ALE geometry as To obtain the true effective theory we need to integrate out all massive modes. As M Inst is in general much lighter than the lightest state of the KK tower we can indeed constrain considerations to modes whose masses are induced by M2-branes as these will yield the largest corrections. Thus integrating these states out will alter the Lagrangian (B.32) to S 4d = d 4 xTr L kin + L Zuk + L Yuk + · · · + 1 M k Inst L k Int + . . . .

(B.34)
where we have made the mass scale M Inst suppressing the higher point interactions explicit.

C Boundary Conditions
In this appendix we provide some details underpinning the computation of the D-cohomology groups in section 4.1. The argument was originally introduced by Witten [39] for closed manifolds. The extension to the manifolds with boundary includes the additional subtlety of choosing the appropriate boundary conditions, which allow us to apply the Hodge theory arguments on manifolds with boundary.
Following [45] we consider the case where we have a single boundary component. We first and ∆ 2 , with the domain We use the indices to formally distinguish between the two operators, based on their different domains. With these domains, the two Laplace operators are self-adjoint and we have the Hodge decompositions [45,65] (C. 5) Here, N and R denote the nullspace and range respectively. Superscripts denote the degree of forms on which the operators act. Observe that this implies The quotient on the right side is precisely the p-th cohomology group of the operator d i . As already noted above, the cohomology of operator d 1 is the de Rham cohomology H p (M 3 ) by definition. To identify the cohomology groups of the operator d 2 , observe that α t | Γ = 0 is the same as saying α| Σ = 0. This implies that the cohomology of d 2 is isomorphic to the relative cohomology H p (M 3 , Σ) [46].
Let us now consider how the deformation of the complex by a smooth function f relates to the above. Recall that D is in fact the deformed de Rham differential D = e −qf de qf .
Denote D i = e −qf d i e qf and D † i = e qf d † i e −qf i.e. we impose the same boundary conditions on D i (resp D † i ) as we did on d i (resp. d † i ). Let ∆ f,i = D i D † i + D † i D i denote the twisted Laplacian from (4.6). We set D(∆ f,i ) = D(∆ i ). We have the same Hodge decomposition for the operators D i , D † i , ∆ f,i as in (C.5) as well as the identification condition we obtain from the electrostatics problem. However, the above discussion clearly outlines the general structure, which persists in the case of mixed boundary conditions, which we consider in 4.1.