Yang-Mills solutions and Spin(7)-instantons on cylinders over coset spaces with G2-structure

We study g\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathfrak{g} $$\end{document} -valued Yang-Mills fields on cylinders ZG/H=ℝ×G/H\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ Z\left(G/H\right)=\mathbb{R}\times G/H $$\end{document}, where G/H is a compact seven-dimensional coset space with G2-structure, g\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathfrak{g} $$\end{document} is the Lie algebra of G, and Z(G/H) inherits a Spin(7)-structure. After imposing a general G-invariance condition, Yang-Mills theory with torsion on Z(G/H) reduces to Newtonian mechanics of a point particle moving in ℝn\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {\mathbb{R}}^n $$\end{document} under the influence of some quartic potential and possibly additional constraints. The kinematics and dynamics depends on the chosen coset space. We consider in detail three coset spaces with nearly parallel G2-structure and four coset spaces with SU(3)-structure. For each case, we analyze the critical points of the potential and present a range of finite-energy solutions. We also study a higher-dimensional analog of the instanton equation. Its solutions yield G-invariant Spin(7)-instanton configurations on Z(G/H), which are special cases of Yang-Mills configurations with torsion.


Introduction
Higher-dimensional Yang-Mills theory naturally appears in superstring theories [1,2], for example as part of the low-energy limit in the presence of D-branes and as a sub-sector of heterotic supergravity. Moreover, the massless spectrum of closed strings contains an JHEP03(2016)038 JHEP03(2016)038 the R-direction. In some cases, the flow equations are supplemented by constraints, which can be algebraic or first-order differential equations. For eq. (1.2), the flow is secondorder and determined by the gradient of a quartic potential, whereas for eq. (1.1), the flow is first-order and determined by the gradient of a cubic superpotential. For this setup, both numerical and analytical solutions have been found in various dimensions and for different kinds of G-structures [9][10][11][12][13][14][15][16]. In this paper, we consider the instanton equation and the Yang-Mills equation with torsion on cylinders R × G/H, where G/H is a compact seven-dimensional coset space with G 2 -structure and the cylinder has an induced Spin(7)-structure.
The paper is organized as follows. In the next section, we set up our conventions on cylinders of the general form R×G/H. The specialization to G/H being seven-dimensional with G 2 -structure is the subject of section 3. After reviewing the key properties of G 2structures in seven dimensions and Spin(7)-structures in eight dimensions, we compute the instanton and Yang-Mills equations. We discuss some simple solutions common to all such coset spaces. However, exploring the full space of solutions requires explicit knowledge of the coset space in question. We thus study the instanton and Yang-Mills equations on a case-by-case basis for three explicit G 2 -structure coset spaces, namely the Berger space SO (5) 1)) and Q pqr = (SU(2) × SU(2) × SU(2))/(U(1) × U(1)). We conclude in section 5 with a summary of the main results and a discussion of possible future directions. Throughout the main text, we frequently need the structure constants of the coset spaces. For reference and to make our presentation self-contained, we thus include in appendix A a list of the structure constants of all the coset spaces appearing in this paper.

Instantons and Yang-Mills equation with torsion on R × G/H
We shall now briefly review the specialization to the case M = R × G/H, following [9]. Throughout this paper, G/H is assumed to be a compact coset space. That is, G denotes a compact connected Lie group and H a closed subgroup, such that G/H is a compact reductive homogeneous space.
The metric g on the cylinder Z(G/H) := R × G/H with a coordinate τ on R can be written as g = δ µν e µ ⊗ e ν = dτ ⊗ dτ + g G/H = dτ ⊗ dτ + δ ab e a ⊗ e b , (2.1) where {e µ } = {e 0 = dτ, e a } is a local orthonormal basis of T * (R × G/H). Note that Greek and lower case Latin indices label the coordinates of R×G/H and G/H, respectively. That is, µ = (0, a) = 0, 1, . . . , (d − 1). The connection A and curvature F can also be expanded
Since G/H is assumed to be a coset space, the Lie algebra g of G can be decomposed as g = h ⊕ m. Here, m is the orthogonal complement of the Lie algebra h of H in g. The generators {I A } A=1,...,dim g spanning g can be divided into two sets {I i } i=dim m+1,...,dim m+dim h and {I a } a=1,...,dim m . The latter are generators of h and m, respectively. They satisfy the following commutation relations, 3) The generators are normalized such that the usual Killing-Cartan metric on g can be written as Throughout this paper, we use the adjoint representation, (I A ) C B = f C AB , for the generators of g. The Killing-Cartan metric can also be succinctly expressed as (g g ) AB = − tr(I A I B ) = −f C AD f D BC = δ AB by virtue of (2.4). Making use of the above coset space construction, we choose a G-invariant ansatz for the gauge connection A [11], A = e i I i + e a X a , (2.5) with X a = X b a I b and {e i } being left-invariant 1-forms on G/H dual to {I i }. The e i can be expressed as linear combinations of the 1-forms {e a } via e i = e i a e a , with real functions e i a . For the precise definitions of e a and e i we refer to [9][10][11], and for a more pedagogical introduction to [17]. For the calculations carried out in this paper, it suffices to keep in mind that the 1-forms {e A } = {e a , e i } are always well-defined on a contractible open subset U of G/H, and moreover that they obey the Maurer-Cartan equations, (2. 6) In order to ensure the G-invariance of the ansatz (2.5), the matrices X a need to satisfy (2.7) In component notation, this condition reads X b a f c ib = f b ia X c b . The curvature F of the connection (2.5) becomes where the dot denotes the derivative with respect to τ .

JHEP03(2016)038
3 Seven-dimensional coset spaces with G 2 -structure In this section, we study the instanton equation and the Yang-Mills equation with torsion on cylinders over compact seven-dimensional coset spaces G/H admitting a G 2 -structure. Before considering Yang-Mills theory in this set-up, we begin by reviewing some important facts about G 2 -structures in seven dimensions and Spin (7)-structures in eight dimensions.

G 2 -structures in seven dimensions
The G 2 -structure on a compact seven-dimensional coset space G/H is defined by a 3-form P . The 4-form dual to P is denoted Q := * 7 P . In a suitably chosen basis, the canonical forms (P, Q) are given by the standard expressions [18,19], Note that in our conventions, the orientations of the seven-dimensional coset space G/H and of the cylinder R × G/H are chosen such that ε (7) 1234567 = +1 and ε (8) 01234567 = +1, respectively. Consequently, the Hodge star operator * on the cylinder is in our conventions related to the Hodge star operator * 7 on the coset space via * ω (7) p = −( * 7 ω (7) p ) ∧ dτ , * (dτ ∧ ω (7) p ) = * 7 ω (7) p . (3.6) Here, ω p is a p-form with legs only in the coset space directions. G 2 -structures can be distinguished and classified by means of the intrinsic torsion. This is analogous to the case of six-dimensional SU(3)-structures, reviewed for example in [3]. In the present case, it holds that where T ab c are the components of the torsion tensor, Λ 1 is the space of 1-forms and so(7) = g 2 ⊕ g ⊥ 2 . The g 2 piece drops when acting on G 2 -invariant forms. The corresponding torsion is called intrinsic torsion, denoted T 0 . It can be decomposed into irreducible G 2 representations [19][20][21][22] according to The tensors τ p are p-forms (p ∈ {0, 1, 2, 3}), and they correspond to the four torsion classes appearing in the exterior derivatives of the structure forms, dP = τ 0 Q + 3 τ 1 ∧ P + * 7 τ 3 , dQ = 4 τ 1 ∧ Q + τ 2 ∧ P .
On the other hand, the Spin(7)-structure on Z(G/H) is locally conformal, if and only if τ 2 = 0 and τ 3 = 0, which corresponds to a locally conformal G 2 -structure on G/H. This includes the special case of a nearly parallel G 2 -structure, which plays an important role in the present work. In this case (and after setting τ 0 = 4, without loss of generality), the Spin(7)-structure equation becomes dΨ = 4 dτ ∧ Q, which can be re-written as implying that Ψ then indeed defines a locally conformal Spin(7)-structure [27]. Before closing this subsection, we mention that the cone over a nearly parallel G 2structure manifold and the cylinder over a G 2 manifold (that is, a manifold with parallel G 2 -structure) are examples of Spin(7) manifolds (that is, manifolds with parallel Spin(7)-structure).

Spin(7)-instantons
For the special case of a Spin(7)-structure in eight dimensions, the instanton equation (1.1) with Q M = Ψ becomes This is equivalent to the statement that F ∈ spin (7). Further specializing to the case where the eight-dimensional manifold is a cylinder over a G 2 -structure manifold and using eq. (3.12), leads to In our case, the second equation is identically satisfied by means of the identity (3.4). The remaining equation, F 0a = 1 2 P abc F bc , can be rewritten on Z(G/H) aṡ after inserting (2.8).

JHEP03(2016)038
Expressed in components, X a = X b a I b , eq. (3.21) reads as follows, The decomposition into two equations is due to the fact that F has a part valued in m, a part valued in h and both need to satisfy (3.21) separately. The first equation in (3.22) comes from the m-part of F and is a first-order ordinary differential equation. The second equation in (3.22), coming from the h-part of F , yields a set of purely algebraic conditions known as quiver relations [28]. The origin of the latter can be traced back to the reduced ansatz for X a which was taken to be m-valued only, X a = X b a I b . For the full g-valued X a = X b a I b +X i a I i , the algebraic conditions instead become first-order differential equations in the extra components X i a . In general, to proceed further, we require explicit knowledge of the groups G and H. This can be seen for example from eq. (2.7), which involves the structure constants. We will thus study the instanton equation in section 3.5 on a case-by-case basis on cylinders over explicit seven-dimensional coset spaces G/H that admit a nearly parallel G 2 -structure.
On the other hand, there exists a specialized ansatz, which has already been considered in [9] and which can be solved explicitly without the need to specify G/H, provided some additional assumptions are imposed. The idea is to trivially solve eq. (2.7) by reducing the number of degrees of freedom to just a single real scalar field φ(τ ). The G-invariance condition (2.7) is obviously satisfied upon setting (3.23) Consequently, the gauge connection A and the corresponding curvature F are given by Now, we assume 2 that the structure constants with all indices lowered, f ABC := δ CD f D AB , are totally anti-symmetric, f [ABC] = f ABC , and that they can be used to construct a welldefined G 2 -structure on G/H. That is, we set P abc = σf abc for some σ ∈ R \ {0}. The second equation in (3.25) is then identically satisfied by virtue of (2.4).
In addition, we assume that the structure constants satisfy the following two equivalent relations

JHEP03(2016)038
implies α = 6/σ 2 . Note that the Killing-Cartan metric is assumed to be δ AB in this subsection. Otherwise it must be reinstated in all expressions where indices of structure constants are manipulated. The instanton equation (3.22) reduces in the above set-up to a first-order ordinary differential equation for the single function φ(τ ), It has two static solutions, φ = 0, 1, and the well-known interpolating kink solution [9] φ(τ ) = Here, τ 0 ∈ R is an arbitrary integration constant fixing the position of the instanton in the τ direction.

Yang-Mills equation with torsion
We now turn to the Yang-Mills equation with torsion, (1.2), on a cylinder over a compact seven-dimensional coset space G/H admitting a G 2 -structure. It is clear that the instanton solutions discussed in the previous section also solve (1.2). Our aim in this section is to analyze whether there are additional analytical solutions.
To start, it is useful to write out (1.2) in components (for a more detailed derivation see [16,17]), Here, ω µ νρ are the components of the affine spin connection ω µ ρ = ω µ νρ e ν on Z(G/H) with torsion T µ = 1 2 T µ νρ e νρ = de µ + ω µ ν ∧ e ν , and H µνρ are the components of the 3-form H = 1 3! H µνρ e µνρ defined in section 1. For the cylinder (product) metric (2.1), one finds ω 0 0b = ω a 0b = ω 0 cb = 0 , and Given the coset space construction, it is natural to relate the components of the torsion tensor of the affine spin connection on G/H to the structure constants [9], parameterized by κ ∈ R. The affine spin connection on G/H becomes

JHEP03(2016)038
Note that κ drops out of the combination 1 2 T ν ρσ − ω ν ρσ and contractions thereof. Instead, it enters in the Yang-Mills equation (3.29) only via the choice (3.33). In sections 3.5.3 and 4, we will encounter examples where P abc ∝ f abc . In such a situation we still use (3.31)-(3.32) but H a bc ∝ T a bc and hence κ does not enter in the Yang-Mills equation. Throughout this paper, ω, T and H never acquire components along the τ -direction. Under this assumption, we find that the ν = 0 component of (3.29) reduces to the Gausslaw constraint on the matrices X a [11], a [X a ,Ẋ a ] = 0 , (3.34) after inserting (2.5). This is a consequence of the gauge fixing A 0 = 0 and becomes nontrivial when X b a is non-diagonal. Carrying over the assumptions on the structure constants from section 3.3 and using (3.30)-(3.32), we find that the remaining ν = a component of (3.29) reduces to a second-order equation given bÿ where repeated indices are to be summed over. If P abc ∝ f abc holds, we may assume, in addition, (3.33) and obtain [11] X a = 1 2 Note that in obtaining (3.35) and (3.36) from (3.29), we used the fact that the terms involving the functions e i c vanish by virtue of the Jacobi identity (for details on this calculation, we refer to [9,17]). We also remark that the choice κ = 1 corresponds to the canonical connection introduced in [6]. With (3.26), κ = 6 and α = 1/5 (i.e. σ 2 = 30), one finds that (3.36) is the τ -derivative of the instanton equation (3.21).
Beyond this point, we again need explicit knowledge of the groups G and H, as in the case of the instanton equation. We will study the Yang-Mills equation with torsion in more detail later. For now, we will follow the same idea as in section 3.3 and consider the single-field ansatz (3.24). This special case has already been studied in [9] and can be solved explicitly. We find a point particle equation of motion of the form after inserting the single-field ansatz (3.24) into (3.36) and assuming that (3.26) holds. Different choices of the parameters (α, κ) correspond to different types of solutions of (3.37). This has been analyzed in detail in [9]. For example, the choice α = 0 leads tö which can be integrated toφ = ± 1 2 (1 − φ 2 ). The solutions, φ = ± tanh τ −τ 0 2 , are known as the φ 4 kink (+) and anti-kink (−), respectively. Here, τ 0 ∈ R is an arbitrary integration JHEP03(2016)038 constant fixing the position of the kink/anti-kink in the τ direction. On the other hand, for the choice (α, κ) = (3/5, 1) we find This can be integrated toφ = ± 1 √ 5 φ(φ − 1). Hence, in this case we recover the tanh-kinktype instanton solutions presented in section 3.3.

Cylinders over the Berger space SO(5)/SO(3) max
For the coset space 4 SO(5)/SO(3) max , we use the conventions of [29], but with the generators and structure constants rescaled as I A → 1 This normalization is necessary in order for (2.4) to hold. For reference, the rescaled structure constants used in this subsection can be found in full detail in appendix A.1.
We observe that the structure constants with all indices lowered, f ABC := δ CD f D AB , are totally anti-symmetric, f [ABC] = f ABC , and that f abc defines a nearly parallel G 2 -structure on SO(5)/SO(3) max . We set P abc = √ 30f abc to ensure the correct normalization such that eqs. (3.2)-(3.5) are satisfied. In the terminology of section 3.3, that means σ = √ 30. Moreover, we note that this coset space satisfies (3.26) with α = 1/5.
Instanton and Yang-Mills equations. In order to study the instanton and Yang-Mills equations, we first need to solve the G-invariance condition (2.7), where now G = SO(5). We find X a = φ(τ )I a and hence this case has already been covered by the general analyses in sections 3.3 and 3.4. In particular, the instanton equation yields, besides constant solutions, the standard tanh-kink-type solutions (3.28) with (α, σ) = (1/5, √ 30). Before moving on to the next coset space, we remark that it is not possible to consistently deform the G 2 -structure on SO(5)/SO(3) max away from being nearly parallel. Indeed, after solving eq. (3.10) for this example, the only remaining admissible deformation parameter is an overall volume rescaling, which does not affect the nature of the G 2 -structure. For the coset space Sp(2) × Sp(1)/Sp(1) 2 , we adopt the conventions of [33]. We rescale the generators and structure constants according to
Also for this coset space, the structure constants with all indices lowered, f ABC := δ CD f D AB , are totally anti-symmetric, f [ABC] = f ABC , and f abc defines a nearly parallel G 2structure. We fix P abc = √ 30f abc to ensure the correct normalization such that eqs. (3.2)-(3.5) are satisfied. In the terminology of section 3.3, that means σ = √ 30. The rescaled structure constants satisfy (3.26) with α = 1/5.

Instanton equation.
We now insert f A BC , P abc and X b a into the instanton equation (3.22). The quiver relations reduce to the quadratic constraint With φ 1 = ±φ 2 ≡ ±φ, the first equation in (3.22) again yields the standard tanh-kink-type equation (3.27) and solution (3.28) with (α, σ) = (1/5, √ 30). It is possible to deform the nearly parallel G 2 -structure on Sp(2) × Sp(1)/Sp(1) 2 in an Sp(2) × Sp(1)-invariant fashion. Solving eq. (3.10), we find the deformed quantities g G/H andP , with two real deformation parameters R 1 , R 2 . The deformed G 2 -structure satisfies (3.2)-(3.4), after replacing P →P , Q →Q, δ ab → (g G/H ) ab . The expressions for the exterior derivatives of the deformed G 2 -structure forms are given by The symbol * 7 denotes the seven-dimensional Hodge star operator with respect to the deformed metricg G/H . In the general classification of G 2 -structures, this is called a cocalibrated (or semi-parallel) G 2 -structure (see section 3.1). For the special case R 1 = R 2 ≡ R, this reduces to the original nearly parallel G 2 -structure, up to an irrelevant volume rescaling. This bears a striking resemblance to the analogous case in one dimension lower, namely the deformed nearly Kähler structure on the sixdimensional coset space Sp(2)/SU(2) × U(1) considered, for example, in [34,35].
In the deformed case, the quiver relations imply On the other hand, the differential equation in (3.22) turns intȯ (3.49) By differentiating eq. (3.48) and inserting eq. (3.49), we arrive at After inserting eq. (3.48) into eq. (3.50), we immediately learn that non-constant solutions only exist for which corresponds to the original (that is, undeformed) nearly parallel G 2 case.

Yang-Mills equation.
To compute the second-order equations of motion for φ 1 and φ 2 , we start with (3.29) and insert (3.42). We follow the philosophy of section 3.4 and choose the same torsion (3.31), affine spin connection (3.32) and 3-form H (3.33). Restricting to the undeformed nearly parallel G 2 case, we then find two differential equations, and one algebraic equation, As a non-trivial consistency check, we remark that the same set of equations can also be obtained by inserting (3.41) into eq. (2.27) of ref. [11].

Cylinders over Aloff-Wallach spaces SU(3)/U(1) k,l
It is possible to define a U(1) subgroup of SU(3) given by matrices of the form where 0 ≤ ϕ ≤ 2π. The integers k, l parameterize the embedding of U(1) into SU(3) and as a mnemonic we write U(1) k,l . For relatively prime integers k and l, the coset spaces SU(3)/U(1) k,l are simply connected manifolds known as Aloff-Wallach spaces. One can show that it is always possible to choose an orthonormal coframe {e a } on U ⊂ SU(3)/U(1) k,l together with a connection one-form e 8 defined as the dual of the generator I 8 (a rescaled version of (3.59)) of the group U(1) k,l , such that the coset space metric g G/H = δ ab e a ⊗ e b is Einstein for a connection with a torsion 3-form P given by P = e 135 − e 245 − e 146 − e 236 + e 127 + e 347 + e 567 .
(3.60) Furthermore, the holonomy group of this connection is contained in G 2 and the 3-form P defines a G 2 -structure on SU(3)/U(1) k,l . Here, we adopt the conventions 5 of [12]. Note 5 When comparing expressions in this subsection and in appendix A.3 with those in [12], it should be kept in mind that we adjusted some notation in order to fit seamlessly into this paper. In particular, our [12], respectively.

JHEP03(2016)038
that in this subsection, the orientation is flipped compared to the rest of the paper. Hence, ε With the structure constants as defined in appendix A.3, we find that P abc ∝ f abc , even though the structure constants with all indices lowered, f ABC := (g g ) CD f D AB , are totally anti-symmetric, f [ABC] = f ABC . We also note that generically (3.26) does not hold in this case. The Killing-Cartan metric is given by and it is used to raise, lower and contract Lie algebra indices. Here, ∆ 2 := 2(k 2 + l 2 ) and ζ 1 , ζ 2 , ζ 3 , µ are four real parameters corresponding to the metric deformations allowed by (3.10). They can be incorporated as rescalings of the generators {I A } and 1-forms {e A } as explained in [12]. (Note that for k = ±l, there are four additional metric deformations, which will, however, not be considered in this paper.) For the exterior derivatives of the structure forms, we find (3.63) In general, P thus defines a cocalibrated (or semi-parallel) G 2 -structure on SU(3)/U(1) k,l . However, for all (k, l), it is possible to find an appropriate set of parameters ζ 1 , ζ 2 , ζ 3 , µ ∈ R such that τ 3 = 0 (see [36]). This justifies why we include SU(3)/U(1) k,l in the list of examples of coset spaces with nearly parallel G 2 -structure. We emphasize however that most of the analysis carried out in this subsection not only applies to the special locus of parameter space where P defines a nearly parallel G 2 -structure, but also to the general case where τ 3 = 0. We will solve the G-invariance condition (2.7), where now G = SU(3), in a slightly more general fashion as compared to the rest of this paper. Namely, we allow X a to also contain a piece valued in the Lie algebra h = u(1). That is, Usually, we set X i a = 0, but in this subsection we consider X i a = 0, thereby ensuring full compatibility with [12]. We then find from the G-invariance condition (2.7) eight τdependent scalar degrees of freedom, arranged into three complex fields φ α (α = 1, 2, 3) and two real fields χ i (i = 7, 8). The only non-vanishing components of X B a are where Re(z) := (z +z)/2 and Im(z) := (z −z)/(2i) for any z ∈ C. We use the notation(·) and (·) * interchangeably to denote complex conjugation. Due to the complex nature of the fields φ α , it is beneficial to distinguish between upper and lower field indices in this subsection. This distinction is not necessary for the other coset spaces considered in this paper. The structure of (3.64) suggests to consider a complexified version of the Lie algebra g = su(3), namely su(3) ⊗ C. As a basis for this Lie algebra, we take can be expanded in the complex basis as (3.66) After solving the G-invariance condition (2.7), we find which is the complex version of (3.64). Here, we introduced the notationφᾱ := (φ α ) * to denote the complex conjugated field.
Instanton equation. This case has already been analyzed in [12]. For completeness, we briefly review the key results relevant to the present paper. From the instanton equation (3.22), one obtains a set of coupled first-order ordinary differential equations, The set of equations (3.68) can be succinctly expressed by means of a superpotential W , schematicallyφ

JHEP03(2016)038
where K AB is the inverse of K AB := − tr(J A J B ), which is the Killing-Cartan metric in the complex basis. The superpotential W is a cubic function of the scalar fields {φ α , χ i }, Eq. (3.68) is a complicated system of coupled, non-linear first-order ordinary differential equations and finding the general solution is a considerable task. Nonetheless, a particularly simple yet important special solution is readily found: Here, c 1 , c 2 ∈ R are constants of integration. For C 7 = 0, c 1 = 0, this solution is stationary and corresponds to the abelian (rescaled, if C 8 = 0) canonical connection on a line bundle over SU(3)/U(1) k,l . The same is also true for C 7 = C 8 = 0 with the rescaled canonical connection corresponding to the case c 2 = 0. This is arguably the simplest example of a G 2 instanton on SU(3)/U(1) k,l . Another special case, which allows a more thorough analysis, is the choice k = l = 1. In addition, we fix ζ 2 1 = ζ 2 2 = 2ζ 2 3 = 2α 2 , µ = 1/α 2 , where α ∈ R is a residual free parameter. 6 The instanton equation then becomes When searching for the critical points of W , one may use its symmetries to argue that it suffices to consider the case φ α ∈ R, φ 1 ≥ 0 and φ 2 ≥ 0, without loss of generality.
As an example, we consider the case α = ±2. In this case, SU(3)/U(1) 1,1 admits a cocalibrated (or semi-parallel) G 2 -structure. In other words, P satisfies (3.62) with τ 0 = 0 and τ 3 = 0. The two signs correspond to the choice of orientation. For α = +2, the critical points of W are listed in the following table: On the other hand, for α = −2, the table of critical points of W is given by: where c ± = 11 ± √ 105/6. There is a kink solution for α = −2 flowing from W = −16 at τ = −∞ to W = 0 as τ → ∞. The numerical solution for this case is shown in figure 1. The shape of these curves resembles that of a hyperbolic tangent type kink. The maximal deviation of a fitted hyperbolic tangent from the numerical solution is of the order of 2%.
Yang-Mills equation. The second-order equations of motion for {φ α , χ i } can be worked out from (3.29), taking over the ansätze (3.31)-(3.32) for T and ω from section 3.4. For H we choose with Q M = Ψ = P ∧ dτ − Q and λ ∈ R being some new free parameter. This ansatz for H differs from (3.33) and is a consequence of P abc ∝ f abc . This is a precursor of the situation encountered in section 4, where the appropriate choice of H will be discussed and motivated further. We now insert (3.31)-(3.32), (3.74) together with (2.5), (2.8) and (3.64) into (3.29). First, we obtain the Gauss-law constraint (3.34), which now becomes,

JHEP03(2016)038
Taking suitable linear combinations yields In addition, we obtain three complex and two real second-order equations expressible as a gradient system of the form, for some potential V defined below. These are the Euler-Lagrange equations derived from the action (1.3), with the replacement Q M → λQ M . Indeed, the action S in this case turns into The energy density E = T + V comprises the kinetic energy T = K αβφ αφβ − 1 4 K ijχ iχj and the potential V given by A full analysis of this large system of equations, for generic k, l, ζ 1 , ζ 2 , ζ 3 , µ and λ, is beyond the scope of the present paper. However, we shall now discuss some special cases, where simplifications occur and thus analytical methods apply. First, we study the case λ = 1. For this value of λ, it is possible to write V as the square of the superpotential W defined in (3.70), where W α = ∂W/∂φ α , Wβ = ∂W/∂φβ and W i = ∂W/∂χ i . The gradient of V thus becomes (3.82) From this and (3.81) we learn that critical points of the superpotential are both zeros and critical points of the potential. We are interested in solutions with finite total energy, Such solutions necessarily interpolate between zero-potential critical points. For λ = 1, they are critical points of the superpotential, and the second-order equations of motion (3.78) can be integrated to the first-order instanton equations (3.69). Hence, we have reduced this case to the construction of instanton solutions discussed in [12] and reviewed above. 7 Second, we consider the single-field reduction φ 1 = φ 2 = φ 3 = χ 7 = 0. The only remaining degree of freedom is χ 8 and the Gauss-law constraint (3.77) is trivially satisfied. The equations of motion (3.78) collapse to a single linear equation, The solution reads Here, c 1 , c 2 ∈ R are constants of integration. For λ = 1, the coefficients A and B factorize, A = 4(C 7 ) 2 , B = 4C 7 C 8 , and (3.84) can be integrated tȯ This corresponds to the instanton case, whose solution was given in (3.71). We remark that the only finite-energy solution is χ 8 = const., which corresponds to the abelian (rescaled, if C 8 = 0) canonical connection on a line bundle over SU(3)/U(1) k,l . This requires either Our last example is the choice k = l = 1, together with ζ 2 1 = ζ 2 2 = 2ζ 2 3 = 2α 2 , µ = 1/α 2 , α = −2, λ = 3/2, φ 2 = φ 3 = χ 8 = 0, χ 7 = 11/4 and φ 1 real. This solves the Gauss-law constraint (3.77), and the equations of motion (3.78) reduce tö   In this section, we consider compact seven-dimensional coset spaces G/H with SU(3)structure. Following the useful references [38][39][40], we will review here some key properties of this slightly unusual set-up, in order to equip us with the necessary knowledge to study the instanton equation and the Yang-Mills equation with torsion on (cylinders over) these spaces.

JHEP03(2016)038
Here, we introduced the symbol to denote the interior product of differential forms on G/H. It is defined as ω η := * 7 ( * 7 ω ∧ η). It is possible to show that a number of additional useful relations hold, namely The intrinsic torsion T 0 ab c decomposes as follows under SU (3), From this we can read off the torsion classes. The five singlets are conveniently arranged into one real and two complex 0-forms, denoted R (real) and W 1 , E (both complex), respectively. The five 1-forms are denoted V 1,2 , W 0 and W 4,5 . The four 2-forms are arranged into two real ones, denoted T 1,2 , and a complex one, denoted W 2 . Finally, there are two 3-forms, denoted W 3 and S.
When arranged in such a way, the exterior derivatives of the forms (V, J, Ω) can be succinctly expressed in terms of these torsion classes [38][39][40], The numerical coefficients in these expressions are fixed by demanding compatibility with the seven-dimensional SU(3)-structure relations (4.1). We remark that the torsion classes have been arranged such that the subset {W 1 , . . . , W 5 } precisely equals the set of torsion classes of a six-dimensional SU(3)-structure (see, for example, [3]). This can be understood as the SU(3)-structure that fixes the geometry of the six-dimensional base-manifold X 6 corresponding to the foliation defined by V .
To make contact with the previous section, we note that SU(3) ⊂ G 2 and thus an SU(3)-structure automatically implies the existence of a G 2 -structure on the sevendimensional manifold. Actually, an SU(3)-structure on a seven-dimensional manifold can be used to define two independent G 2 -structures via (4.7) The intersection of the two G 2 -structures P ± is precisely the SU(3)-structure. Henceforth, we will mostly be working with the G 2 -forms P ± , instead of the SU(3)-forms (V, J, Ω). It should be kept in mind though that the two formulations are equivalent.

Explicit coset space constructions
To proceed further, we will now consider cylinders over explicit seven-dimensional coset spaces G/H that admit an SU(3)-structure. Specifically, we shall examine four cases, namely cylinders over SO (5) As a starting point, we use the conventions of [40] in order to define the coset space SO(5)/SO(3) A+B . Although the structure constants are not written down explicitly, they can be extracted straightforwardly from eq. (D.12) in [40]. We work with rescaled generators and structure constants, to ensure the correct normalization of (2.4). Our rescaled structure constants are listed in appendix A.4 for completeness.
The structure constants with all indices lowered, f ABC := δ CD f D AB , are totally antisymmetric, f [ABC] = f ABC , as before. However, we will now encounter a crucial difference in comparison to the coset spaces studied in sections 3.5.1 and 3.5.2, namely that P abc ∝ f abc does not lead to a well-defined G 2 -structure on SO(5)/SO(3) A+B . This assertion can be verified inter alia by computing the so-called associated metric (see for example [41][42][43][44]), (g P ) ab = B ab det(B) −1/9 , with B ab = − 1 144 P ac 1 c 2 P bc 3 c 4 P c 5 c 6 c 7 ε (7)c 1 ...c 7 . (4.8) For P abc ∝ f abc this metric is singular, that is det g P = 0.
To resolve this problem, recall that there is an SU(3)-structure on this coset space. Hence, we may instead use (4.7) to define a G 2 -structure. On SO(5)/SO(3) A+B and with the correct normalization such that eqs. (3.2)-(3.5) are satisfied, we then find  The corresponding 4-form Q ± := * 7 P ± is explicitly given by
The G-invariance condition (2.7), where now G = SO (5), is solved by (4.12) Hence, we need to deal with a priori five real τ -dependent scalar degrees of freedom. The full expressions for the gauge connection A and the corresponding curvature F can be obtained straightforwardly from (2.5), (2.8) and are rather lengthy. We thus omit them here.

JHEP03(2016)038
Yang-Mills equation. We will now turn to the second-order equations of motion for {φ 1 , . . . , φ 5 }. Our starting point is the Yang-Mills equation in components (3.29). First, we need to choose ansätze for the torsion T , affine spin connection ω and 3-form H. For T and ω we take over the expressions from section 3.4, that is (3.30) and depending on the real parameter κ ∈ R.
The ansatz for H is more delicate, however. In most of section 3, we took H a bc = −T a bc . Since H a bc ∝ T a bc ∝ f a bc ∝ P a bc , this choice ensured that we could recover the instanton solutions, for which it is required that H ∝ P (see footnote 3), as special cases of the second-order equations of motion. Here, P abc ∝ f abc and hence this logic breaks down. Nevertheless, we aim at keeping the property that the instanton solutions can be recovered as special cases. We are thus led to setting where λ ∈ R is some new free parameter. The instanton case is then contained as a special case, namely λ = 1. After inserting Q M = Ψ ± = P ± ∧ dτ − Q ± into (4.17) and using eqs. (4.9)-(4.11), we obtain 8 Note that, although allowed by (4.17), H ± has no components with legs in the τ -direction. This is a consequence of dQ ± = 0. We now insert the above ansätze (4.16), (4.18) together with (2.5), (2.8) and (4.12) into (3.29). First, we obtain eqs. (4.13)-(4.14) again, provided λ = 0 and irrespective of the choice P + or P − . Hence, (φ 3 , φ 5 ) = (±φ 1 , ∓φ 4 ), as in the case of the instanton equation. Second, we obtain a system of differential equations. With the reduction of the number of degrees of freedom taken into account, the system of differential equations becomes The sign ambiguity in eqs. (4.19)-(4.21) is not related to the choice P + or P − , but rather a consequence of the identification (φ 3 , φ 5 ) = (±φ 1 , ∓φ 4 ). Note that one may switch between 8 An alternative way to compute H from (4.17) is to use * d * Q M,λ = (1/3!)∇ A (Q M,λ )ABCDe BCD and (4.16). This yields H ± = −λ κ+1 4 f e ab (Q ± ) a cde e bcd , which shows that our choice of H is very similar to those considered in related works (for example, see [16] eq. (3.8)). As a side result, we learn that κ = 1 in order to achieve compatibility with (4.18).

JHEP03(2016)038
the two sign choices by sending φ 2 → −φ 2 . Since this does not lead to qualitatively new solutions, we will henceforth ignore the sign choice and restrict to the upper sign, without loss of generality.
Eq. (4.22) is a consequence of the Gauss-law constraint (3.34). It is solved by φ 4 = c φ 1 for some c ∈ R (or φ 1 = 0, which is however equivalent to the case c = 0, due to the symmetry of eqs. (4.19)-(4.22) under interchanging φ 1 ↔ φ 4 ). The remaining eqs. (4.19)-(4.21) form an a priori over-determined system of non-linear second-order ordinary differential equations. We refrain from solving the most general case here. Instead, we shall now examine the two obvious special choices c = 0 and c = 1.
Upon setting c = 0, that is φ 4 = 0, we are left with the following reduced system of differential equations in the two remaining fields {φ 1 , φ 2 }, with energy density 9 E given by (4.26) plus a total derivative, which is omitted. Here, we used the fact that tr(I A I B ) = f C AD f D BC = −δ AB , as explained below (2.4), to resolve the trace appearing in the action (1.3). We also replaced Q M by Q M,λ in the action to ensure compatibility with (4.17).
Eqs. (4.23)-(4.24) can be re-written as a gradient system of the form, where ∂ α := ∂/∂φ α , α = 1, 2, and the potential V is determined in (4.26). In order to analyze this gradient system further, we study the critical points of V , that is field values for which ∂ α V = 0 holds. We find seven distinct critical points and list them, together with some important properties, in table 1. Below we are interested in solutions with finite total energy, Such solutions necessarily interpolate between zero-potential critical points. As an example, consider the case φ 2 = 1 and λ = −1/2. Then, (4.24) becomes trivial, whereas (4.23) reduces to (3.38), implying that φ 1 is a φ 4 kink or anti-kink. The solution in this case is a JHEP03(2016)038 Table 1.
Solutions connecting other pairs of zero-potential critical points generally have both φ 1 and φ 2 non-constant, and thus lie beyond the scope of analytical methods. We will not present the respective numerical solutions here, albeit they may be constructed straightforwardly (see [11] for analogous constructions in one dimension lower). Another option is to further reduce the number of fields by setting φ 2 = ± √ 3φ 1 and λ = −1. The single remaining field φ 1 is governed by the following equation of motion, The corresponding potential V derived from (4.26) cannot be written as a square V (φ 1 ) = f (φ 1 ) 2 . Hence, (4.30) cannot be integrated straightforwardly to a first-order equatioṅ φ 1 = 2f (φ 1 ). The potential V has three critical points 0, ∓1/ √ 3 and ± √ 3/2, with values V | φ 1 =0 = 1/6, V | φ 1 =∓1/ √ 3 = 1/9 and V | φ 1 =± Finally, we shall examine the case c = 1, that is φ 4 = φ 1 . The action in this case becomes with energy density given by (4.32) Up to a field rescaling of the form φ 1 → φ 1 / √ 2, this expression is identical to (4.26). Hence, the case c = 1 reduces to the case c = 0, which has been discussed above.

Cylinders over N pqr = (SU(3) × U(1))/(U(1) × U(1))
We now turn to studying the coset space N pqr = (SU(3)×U(1))/(U(1)×U(1)) and cylinders thereover. For the description of N pqr , we will closely follow [46], and begin with the observation that the group G = SU(3) × U(1) has a subgroup H = U(1) × U(1) which can be embedded into G in various ways. The different embeddings are parameterized by three integers p, q and r. Without loss of generality, they can be taken to be relatively prime.
To proceed further, we adopt the conventions of [46]. We now encounter an important difference compared to the coset spaces studied above. Previously, we used, where possible, an orthonormal basis of g (orthonormal in the sense of (2.4)), which is then neatly compatible with the choice of metric (2.1). In the present case, the Lie algebra g = su(3) ⊕ u(1) u(3) is not semisimple and hence, there is no orthonormal basis.
However, there is a basis in which at least the coset space components of the Killing-Cartan metric are equal to the Kronecker delta, (g g ) ab = δ ab . This can be achieved by taking the conventions of [46] with the following rescalings of the generators I A and structure JHEP03(2016)038 where ζ := 3p 2 + q 2 + 2r 2 and η := 3p 2 + q 2 . We will henceforth work in this basis. The full set of structure constants is listed in appendix A.5, in order for this paper to be self-contained. We emphasize that the cylinder metric is still taken to be (2.1). However, it is necessary to use the full Killing-Cartan metric, (4.34) in order to raise, lower and contract Lie algebra indices. We will highlight places where this is relevant in due course. For example, the structure constants with all indices lowered, Note that it is closed, dQ ± = 0, whereas the exterior derivative of P ± can be expressed as, with τ ± 0 = −(±1) √ 3/2 and (2 √ 3ζ 2 ) τ ± 3 = −(2qη ± ζ 2 )e 127 − (3pη + qη ∓ ζ 2 )e 347 − (3pη − qη ± ζ 2 )e 567 . Hence, the SU(3)-structure defines two cocalibrated (or semi-parallel) G 2structures on N pqr .
The G-invariance condition (2.7), where now G = SU(3) × U(1), is solved by Hence, the dynamical degrees of freedom are in this case a priori seven real τ -dependent scalar fields. By means of (2.5), (2.8), it is straightforward to compute the explicit expressions for the gauge connection A and the corresponding curvature F .

JHEP03(2016)038
Instanton equation. In order to compute the system of equations for {φ 1 , . . . , φ 7 } imposed by the instanton equation, we insert f A BC , P ± and X b a into (3.22). Independent of the choice P + or P − , the quiver relations yield This pair of equations has two branches of solutions, r = 0 and r = 0. For r = 0, it becomes (2φ 2 , thereby removing one degree of freedom. For r = 0, we have two conditions, namely p(2φ 2 This removes two degrees of freedom if pq = 0, and one, otherwise. (Note that p = q = 0 is excluded since the triple (p, q, r) is pairwise coprime, by assumption.) The first-order ordinary differential equations derived from (3.22) are given bẏ The sign ambiguity is due to the choice P + or P − . Finding the general solution of this system of equations and constraints is a formidable task, even for fixed values of (p, q, r). One may however study special sub-sectors of field space, where simplifications occur. For example, let us consider the case φ ≡ φ 1 = φ 2 = φ 3 , φ 4 = 1, φ 5 = φ 6 = φ 7 = 0 and (p, q, r) arbitrary. This solves the quiver relations identically. The differential equations collapse to the standard tanh-kink-type instanton (3.28) for the single remaining degree of freedom φ with ασ = ∓2/ √ 3.
Yang-Mills equation. Our next goal is to compute the second-order equations of motion for {φ 1 , . . . , φ 7 }. We will use the same ansätze for the torsion T , affine spin connection ω and 3-form H as in the case G/H = SO(5)/SO(3) A+B presented in section 4.2.1. In particular, from (4.17) we find the explicit expression (4.48)

JHEP03(2016)038
We have now collected all necessary ingredients to compute the equations of motion following from (3.29). We first obtain two algebraic constraints, as well as two first-order constraints, The first-order equations (4.51)-(4.52) originate from the Gauss-law constraint (3.34). Note that (4.50) and (4.52) become trivial for r = 0. In addition, we obtain seven second-order equations expressible as a gradient system of the form,φ where η αβ := diag(1, 1, 1, 2, 1, 1, 1) αβ and the potential V is defined below. This gradient system corresponds to the Euler-Lagrange equations obtained from the action (1.3) with the reduction to scalar fields, (4.38), inserted. Indeed, the action S in this case turns into where we used tr(I A I B ) = −(g g ) AB and the replacement Q M → Q M,λ . The energy density E = T + V comprises the standard kinetic term T = 1 2 η αβφ αφβ , where η αβ is the matrix inverse of η αβ , and the potential V , given by (4.55) -31 -
Instead of further considering this system of equations in full generality, we shall now briefly discuss a special case. Namely, let us consider the ansatz φ 1 = φ 2 = φ 3 , φ 5 = φ 6 = φ 7 = 0 and λ = 1. The system of equations then enforces φ 1 = 0 or φ 4 = 1. In the former case, we find φ 4 (τ ) = a τ + b and all other φ's vanishing. This solution is unphysical, since the total energy (4.28) is divergent. In the latter case, φ 4 = 1 and we recover the tanh-kink-type instanton solution for φ 1 that has already been found above (see below (4.47)). This case can be studied in close analogy to N pqr , which was presented in the previous subsection. For the sake of brevity, we thus condense the discussion as much as possible and only highlight relevant differences. Of course, the main difference is the group structure, which, for our purpose, manifests itself in a different set of structure constants.

Cylinders over
The group G = SU(3) × SU(2) × U(1) admits various embeddings of the subgroup H = SU(2) × U(1) × U(1). The different embeddings are parameterized by three coprime integers p, q and r. The Lie algebra g = su(3) ⊕ su(2) ⊕ u(1) is not semisimple and hence, does not have an orthonormal basis. Instead, we use a basis in which at least the coset space components of the Killing-Cartan metric are equal to the Kronecker delta, (g g ) ab = δ ab . This can be achieved by taking the conventions of [47] with the following rescalings of the generators I A and structure constants f A BC , Here, we set ζ := 3p 2 + q 2 + 2r 2 , η := 3p 2 + q 2 and ρ := 9p 2 + 2q 2 . The precise expressions for the structure constants in the new basis can be found in appendix A.6. The Killing-Cartan metric, (g g ) AB = −f C AD f D BC , in this basis is given by, (g g ) AB = δ A=a,B=b +3(δ A8 δ B8 +δ A9 δ B9 +δ A,10 δ B,10 )+ 2ρ 2 r 2 ζ 2 η 2 δ A,11 δ B,11 + 3(2p 2 +q 2 ) η 2 δ A,12 δ B,12 It is used to raise, lower and contract Lie algebra indices. The structure constants with all indices lowered, f ABC := (g g ) CD f D AB , are totally anti-symmetric, f [ABC] = f ABC . Next, we need to construct a 3-form P defining the G 2 -structure on M pqr . As before, the 3-form P is, up to discrete sign choices, fully determined by (3.11), (3.2) and the requirement that the associated metric g P (see (4.8)) be positive definite. At this point, we encounter a crucial novelty when compared to N pqr . Namely, a 3-form P on M pqr that solves the aforementioned requirements only exists for p = ±q, r = 0 and without loss of generality, we may assume p = q = 1. This observation is a manifestation of the well-known fact that only M 110 admits an SU(3) × SU(2) × U(1)-invariant G 2 -structure [44,48].
For the rest of this subsection, we thus restrict to M 110 . There are two independent G 2 -structures on M 110 defined by P ± = −e 127 ∓ e 135 ± e 146 ± e 236 ± e 245 − e 347 − e 567 .  It is closed, dQ ± = 0, whereas the exterior derivative of P ± can be expressed as, The induced Spin(7)-structure on the cylinder is fixed by Ψ ± = P ± ∧dτ −Q ± , according to (3.12). Using that τ 0 = 0, τ 3 = 0, dΨ ± = 0 and the results of section 3.2, we immediately conclude that this is a general Spin(7)-structure where both torsion classes W 8 and W 48 are turned on.
The G-invariance condition (2.7), where now G = SU(3) × SU(2) × U(1), is solved by Hence, the dynamical degrees of freedom are in this case a priori five real τ -dependent scalar fields.

JHEP03(2016)038
Instanton equation. We are interested in computing the system of equations for {φ 1 , . . . , φ 5 } imposed by the instanton equation on Z(M 110 ). To achieve this, we insert f A BC , P ± and X b a into (3.22). We find a coupled system of first-order ordinary differential equations,φ supplemented by an algebraic constraint, which follows from the quiver relation. Note that the entire system of equations is independent of the choice P + or P − . Before moving on to the second-order Yang-Mills equations, we mention a special solution which is analytical. This is the static solution, Yang-Mills equation. In this part, we shall compute the equations of motion for {φ 1 , . . . , φ 5 }, using the same ansätze for the torsion T , affine spin connection ω and 3form H as before. In particular, from (4.17) we find the explicit expression We have now collected all necessary ingredients to derive the equations of motion on Z(M 110 ) from (3.29). We obtain one algebraic constraint, as well as two first-order constraints, The latter are a consequence of the Gauss-law constraint (3.34) and are solved by 11 φ 4 = c 1 φ 1 , φ 5 = c 2 φ 2 for some c 1 , c 2 ∈ R. Eq. (3.29) also yields a set of second-order differential equations, which can equivalently be obtained as the Euler-Lagrange equations of the action (1.3) with the reduction to scalar fields, (4.64), inserted. After using that tr(I A I B ) = −(g g ) AB and replacing Q M by Q M,λ , the action becomes and the potential V , given by In order to bring the expressions for T and V in the above form, we performed a field redefinition, . Note that the equivalence of the second-order equation for φ 3 derived from (4.70) and the respective equation obtained from (3.29) holds up to a term proportional to the algebraic constraint, which reads as follows, after the field redefinition. We also remark that all equations starting from (4.68) and below are independent of the choice P + or P − . We shall now consider an interesting scenario, where the constraint can be easily satisfied, namely λ = −1/2. The constraint then implies either φ 2 1 = 2φ 2 2 or φ 3 = 1/(2 √ 2). The first case, φ 2 1 = 2φ 2 2 , can be quickly dealt with, for it leads to φ 1 = φ 2 = 0 and φ 3 = aτ + b with integration constants a, b ∈ R. This solution is unphysical since the value of the total energy is divergent, The second case, namely φ 3 = 1/(2 √ 2), turns out to be more fruitful however. The energy density then reduces to and the equations of motion becomë (4.75) The critical points of the potential V in the field space (φ 1 , φ 2 ) are listed in table 2. In total, there are nine distinct critical points, subdivided into a local maximum at (0, 0), four saddle points (0, ±1/2), (±1/ √ 2, 0), and four local minima at (±1/ √ 2, ±1/2). There are different types of non-constant solutions, depending on whether one allows both fields to have non-trivial τ -dependence or not. In the former case, the equations of motions (4.75) can be integrated tȯ  Table 2. Critical points of the potential V and some of their properties on Z(M 110 ) with λ = −1/2 and φ 3 = 1/(2 √ 2). In the first two columns we list the values of the critical points and of the potential at the corresponding critical point, respectively. In the third column we summarize the signs of the eigenvalues of the Hessian matrix at the respective critical point. In the last row, the signs are mutually independent in the entries of (φ 1 , φ 2 ). Table 3. Different types of non-trivial analytical solutions on Z(M 110 ) with λ = −1/2 and φ 3 = 1/(2 √ 2). Here, "tanh" refers to a rescaled φ 4 kink/anti-kink solution of the form (4.77). The second column holds the value of the total energy E of the respective solution. The last two rows correspond to unphysical solutions with divergent total energy. The solution is a system of two uncoupled rescaled φ 4 kinks/anti-kinks,

Cylinders over
The last seven-dimensional coset space we consider in this paper is Q pqr = (SU(2)×SU(2)× SU(2))/(U(1) × U(1)). The treatment is similar to N pqr and M pqr , which is why we refer to sections 4.2.2 and 4.2.3 for the technical details. Before discussing the instanton and Yang-Mills equations on Z(Q pqr ), we list the coset space's defining data, which is necessary to carry out our computations.
The group G = SU(2) × SU(2) × SU(2) admits various embeddings of the subgroup H = U(1) × U(1). The different embeddings are parameterized by three coprime integers p, q and r. In order to define the coset space, we start with the generators I A and structure constants f A BC from [47], and perform a rescaling, , so as to ensure the correct normalization of (2.4). Our rescaled structure constants are listed in appendix A.7 for completeness.
The structure constants with all indices lowered, f ABC := δ CD f D AB , are totally antisymmetric, f [ABC] = f ABC . However, taking P abc ∝ f abc does not yield a well-defined G 2 -structure on Q pqr (for example, the associated metric g P is singular for this choice). Instead, we construct the 3-form P subject to (3.11), (3.2) and the requirement that g P be positive definite. We find that such a P only exists if p = ±q, q = ±r and without loss of generality, we may restrict to the case p = q = r = 1. This observation is a manifestation of the well-known fact that only Q 111 admits an (SU(2)) 3 -invariant G 2 -structure [44,49].
For the rest of this subsection, we thus restrict to Q 111 . There are two independent G 2 -structures on Q 111 defined by P ± = −e 127 ± e 136 ± e 145 ± e 235 ∓ e 246 − e 347 − e 567 .
The G-invariance condition (2.7), where now G = (SU(2)) 3 , is solved by (4.83) Hence, the dynamical degrees of freedom are in this case a priori seven real τ -dependent scalar fields.

JHEP03(2016)038
Yang-Mills equation. Our next goal is to compute the second-order equations of motion for {φ 1 , . . . , φ 7 }, using the same ansätze for the torsion T , affine spin connection ω and 3-form H as before. In particular, from (4.17) we find the explicit expression With all the necessary ingredients at hand, we are in a position to compute the equations of motion on Z(Q 111 ) following from (3.29). It turns out that they are independent of the choice P + or P − . We obtain two algebraic constraints, as well as three first-order constraints (originating from the Gauss-law constraint (3.34)), (4.88) After taking suitable linear combinations, the latter become and thus, φ 5 = c 1 φ 1 , φ 6 = c 2 φ 2 , φ 7 = c 3 φ 3 for some c 1 , c 2 , c 3 ∈ R (without loss of generality, we may assume that φ 1 = 0, φ 2 = 0, φ 3 = 0, since the equations of motion are symmetric under interchanging φ 1 ↔ φ 5 , φ 2 ↔ φ 6 , φ 3 ↔ φ 7 ). Eq. (3.29) also yields a set of second-order equations, which can equivalently be obtained as the Euler-Lagrange equations of the action (1.3) with the reduction to scalar fields, (4.83), inserted. After using that tr(I A I B ) = −δ AB and replacing Q M by Q M,λ , the action becomes Here, the energy density E = T + V comprises the standard kinetic term and the potential V , given by In order to bring the expressions for T and V in the above form, we performed a field redefinition,
The second case is This is a three-component φ 4 theory invariant under S 3 × (Z 2 ) 3 . The different types of critical points (modulo S 3 × (Z 2 ) 3 -transformations) together with some of their properties are listed in table 4, where we introduced the 3-vector φ = (φ 1 , φ 2 , φ 3 ) . In total, we find 27 distinct critical points, which fall into four classes.

JHEP03(2016)038 5 Conclusions and outlook
In this paper, we studied the instanton equation, (1.1), and the Yang-Mills equation with torsion, (1.2), on cylinders Z(G/H) = R × G/H over seven-dimensional coset spaces G/H with G 2 -structure. The G 2 -structure on G/H lifts naturally to a Spin(7)-structure on Z(G/H).
In order to be able to construct instanton and Yang-Mills solutions, we chose an ansatz for the gauge connection A, which can be written as after gauge fixing. The dynamical degrees of freedom are carried by the τ -dependent 7 × 7matrix X b a (τ ), where τ is the cylinder coordinate. The matrix X b a (τ ) is constrained by the requirement that A be G-invariant, which translates into the condition The simplest solution of this condition is given by a single-field ansatz X b a (τ ) = φ(τ )δ b a . This common solution exists on any G/H and has already been found in [9]. Under certain additional assumptions, the instanton equation reduces tȯ modulo a trivial overall rescaling. Eq. (5.2) has two static solutions, φ = 0, 1, and the well-known interpolating kink solution φ(τ ) = 1 2 (1 − tanh [τ − τ 0 ]), where τ 0 ∈ R is an arbitrary integration constant fixing the position of the instanton in the τ direction.
In the general case, there are multiple τ -dependent scalar fields given by the matrix components X b a (τ ) and restricted by X b a f c ib = f b ia X c b . This restriction involves the structure constants f A BC of G and hence, the solutions of this condition will depend on the choice of coset space G/H. In particular, even the number of independent components of X b a , or in other words, the number of scalar degrees of freedom, is different for different coset spaces. This forbids a unified treatment. Instead, we had to consider the multi-field set-up on a case-by-case basis.

JHEP03(2016)038
In total, we considered seven explicit coset space constructions that fall into two families, namely those with nearly parallel G 2 -structure and those with SU(3)-structure. The first family comprises the Berger space SO(5)/SO(3) max , the squashed 7-sphere Sp(2) × Sp(1)/Sp(1) 2 , and the Aloff-Wallach spaces SU(3)/U(1) k,l , for a co-prime pair of integers (k, l). All three seven-dimensional coset spaces admit a nearly parallel G 2 -structure and an induced Spin(7)-structure on the cylinder thereover.
For G/H equal to SO(5)/SO(3) max , the G-invariance condition forces X b a (τ ) = φ(τ )δ b a . Thus, this case reduces to the single-field set-up already discussed above. In the second example of this group, namely Sp(2) × Sp(1)/Sp(1) 2 , we a priori have two real scalar fields φ 1 (τ ) and φ 2 (τ ), after solving the G-invariance condition. However, the dynamical equations (that is, both, the instanton equation and the Yang-Mills equation with torsion) force the two fields to be proportional to one another or one field to be constant. In effect, also this case reduces to the well-known single-field scenario. The instanton equation on R × SU(3)/U(1) k,l has already been analyzed in [12] and was reviewed in section 3.5.3 for completeness. There we also studied the Yang-Mills equation with torsion on R × SU(3)/U(1) k,l . After solving the G-invariance condition, we showed that it can be reduced to a gradient system in three complex and two real scalar fields. The gradient system is determined by a quartic potential and subject to the Gauss-law constraint. Beyond the instanton constructions, we obtained two further analytical solutions in other special corners of configuration and parameter space. The first one corresponds to a single varying field which however diverges at τ → ±∞ either exponentially or quadratically depending on the values of the parameters. The second one describes a rescaled φ 4 kink/anti-kink.
Then, in each of the four cases SO(5)/SO(3) A+B , N pqr , M 110 and Q 111 , the SU(3)structure induces two cocalibrated (or semi-parallel) G 2 -structures on the coset space and two Spin(7)-structures on the cylinder thereover. After solving the G-invariance condition, we are left with five, seven, five and seven real τ -dependent scalar fields, respectively. The instanton equation (1.1) reduces to a system of coupled non-linear first-order ordinary differential equations in those fields, supplemented by a number of (typically quadratic) algebraic constraints known as quiver relations. The full system of equations can generally not be solved analytically. We did present however several solutions in special corners of fieldand parameter-space. These include the trivial case of constant (that is, τ -independent) solutions, unphysical configurations exhibiting exponential growth, but also solutions that reduce to the well-known tanh-kink-type instanton (5.2).

JHEP03(2016)038
The Yang-Mills equation with torsion for the four cases SO(5)/SO(3) A+B , N pqr , M 110 and Q 111 reduces inter alia to a system of coupled non-linear second-order ordinary differential equations. This system of equations is always expressible (possibly, after performing some suitable field redefinitions) as a gradient system, with a quartic potential V = V ({φ α }) and a constant metric η αβ . In addition, we found a set of quadratic quiver relations and first-order constraints. The latter originate from the Gauss-law constraint (3.34) and are generally of the form f 1 ({φ αφβ }) = 0, f 2 ({φ αφβ }) = 0, etc., where f 1 , f 2 , etc. are some linear functions. In total, we faced an overdetermined system of coupled non-linear algebraic as well as differential equations. Any attempt at solving such a system is a formidable task. It is however possible to consistently restrict one's attention to a reduced set of scalar degrees of freedom such that the algebraic and first-order constraints are trivially satisfied, while retaining the non-trivial multi-dimensional nature of the motion. This is achieved by setting some of the scalar fields equal to one another or to constants and checking that this is compatible with all equations. The remaining equations still form a gradient system as in (5.6), but for a subset of {φ α }.
We applied this procedure to each of the four cases SO(5)/SO(3) A+B , N pqr , M 110 and Q 111 . The resulting gradient systems can be analyzed by standard methods. First, we computed the critical points of the respective potential V . We are interested in nonconstant finite-energy solutions flowing from a critical point at τ → −∞ to a critical point at τ → +∞. Such solutions necessarily interpolate between zero-potential critical points.
In this way, we obtained a multitude of solutions on coset spaces that have not been studied before in this context. Besides constant solutions, linearly diverging solutions and single-field tanh-kink-type solutions, we found a couple of additional solutions which are new and interesting. For example in the case G/H = SO(5)/SO(3) A+B , there exists a setup where the dynamics is described by a single field φ governed by an equation of motion of the form,φ Intriguingly, this equation can in general not be integrated to a first-order equation. However, we have not been able to find a non-constant finite-energy solution. For G/H = M 110 , we showed that in a special corner of field-space the gradient system (5.6) de-couples, thus yielding two independent copies of (5.3) with different values of the parameters. We constructed several different types of analytical solutions for this case and showed their plots in figure 3. Most of them correspond to finite-energy configurations. Finally, for G/H = Q 111 , we set one field to a constant and obtained a three-component φ 4 theory invariant under S 3 ×(Z 2 ) 3 . Non-trivial finite-energy configurations are given by a 3-vector with independent (rescaled) φ 4 kinks/anti-kinks in each entry of the vector. We end this paper by pointing out some open problems and possible future directions. Besides the obvious possibility of extending our methods to even more coset space constructions, it may also be fruitful to analyze in more detail the systems of coupled non-linear JHEP03(2016)038 algebraic and differential equations obtained in this work. In particular, a more systematic and exhaustive exploration of the various parameter spaces may lead to new analytical or numerical solutions.
Another direction in which this work can be continued is to alter the properties of the cylinder direction. For example, one may consider the instanton and Yang-Mills equations on a Lorentzian cylinder iR × G/H, that is sending τ → it, or on S 1 × G/H, that is identifying τ τ + L. The corresponding finite-energy solutions offer an interpretation as dyons (for iR × G/H) and sphalerons (for S 1 × G/H), respectively.
Finally, it would be interesting to examine possible embeddings into string theory or supergravity. A promising candidate for the embedding is provided by heterotic supergravity. One may attempt a construction analogous to [34,35,50,51], where ten-dimensional space-time is taken to be a (possibly) warped product of a four-dimensional external part and a compact six-dimensional SU(3)-structure manifold. With the set-up of this paper, one may instead consider a product of a three-dimensional external part and a compact seven-dimensional G 2 -structure coset space. Some explicit constructions of this type have already been realized, for example, in [6,52]. However, these constructions are based on specific assumptions, which may not exhaust the possible types of embeddings. It is left for future work to determine whether other types of embeddings are also viable. Here, we list the structure constants f A BC of g = so(5) used in section 4.2.1. They were extracted from eq. (D.12) in [40] and rescaled, f C AB → 1 √ 6 f C AB . The only non-vanishing