On the global existence of spherically symmetric hairy black holes and solitons in anti-de Sitter Einstein-Yang-Mills theories with compact semisimple gauge groups

We investigate the existence of black hole and soliton solutions to four dimensional, anti-de Sitter (adS), Einstein-Yang-Mills theories with general semisimple connected and simply connected gauge groups, concentrating on the so-called"regular"case. We here generalise results for the asymptotically flat case, and compare our system with similar results from the well-researched adS $\mathfrak{su}(N)$ system. We find the analysis differs from the asymptotically flat case in some important ways: the biggest difference is that for $\Lambda<0$, solutions are much less constrained as $r\rightarrow\infty$, making it possible to prove the existence of global solutions to the field equations in some neighbourhood of existing trivial solutions, and in the limit of $|\Lambda|\rightarrow\infty$. In particular, we can identify non-trivial solutions where the gauge field functions have no zeroes, which in the $\mathfrak{su}(N)$ case proved important to stability.


Introduction
Therefore, it seems logical to perform the same experiment on the asymptotically adS, spherically symmetric EYM system for a general compact semisimple gauge group, and to see how many features are present in both the general case and the specific su(N ) case. Also strongly motivating this work is the possibility of exploring a very wide class of matter theories, both for the sake of CMP, and for further refinement of the "no-hair" theorem (see Sect. 9) which is relevant to gravitational physics. For the regular case at least, which is the main case considered in the literature so far, we see that it is not even necessary to know the YM one-form connection explicitly in order to obtain the field equations-all the information one needs is essentially in the Cartan matrix of the Lie algebra of the structure group G which represents the gauge field, making it easy to apply to a wide spectrum of EYM theories.
The outline of this paper is as follows. First, in Sects. 2 and 3 we will describe how we use our ansätze to carve down the general field equations for four dimensional adS EYM theory with a general compact gauge group in the case of the 'regular action', which we will describe later; and we show that in doing so, it coincides with the principal action-this allows us to simplify the field equations considerably. In fact, they become very similar in form to the field equations for su(N ) [9]. In Sect. 4, we consider the boundary conditions needed for our solutions to be regular at r = r h (or r = 0) and as r → ∞. In Sect. 5, we examine the asymptotic limit of the field equations r → ∞ in a 'dynamical systems' sense, which turns out to be much simpler than it was for asymptotically flat space. Then in Sect. 6 we identify some trivial embedded solutions, which are important to our final results.
In Sect. 7, we prove the existence of solutions locally at the boundaries, which are unique and analytic in their boundary parameters. Finally, in Sect. 8, after proving that solutions may be regularly integrated out from the initial boundary into the asymptotic regime, we finish by establishing our main results: that global nodeless black hole and soliton solutions may be found in a neighbourhood of some trivial solutions found in Sect. 6, which are everywhere regular and uniquely and consistently specified by their boundary conditions; and that nodeless black hole and soliton solutions can be found in the limit | | → ∞ (Sect. 8.2), anticipating a later investigation into the stability of these solutions. In Sect. 9 we present our conclusions.

Spherically symmetric, purely magnetic Yang-Mills connections for asymptotically adS spacetime
For asymptotically flat space, it is found [34] that we can reduce our attention from considering all possible conjugacy classes of bundle automorphisms by restricting focus to those for which the YM fields decay sufficiently fast at either boundary (r → ∞, and/or r = 0 if the solution is a soliton). These are called 'regular models' in [36] and correspond to the 'zero magnetic charge' case in [37]. A conjugacy class of SU (2) bundle automorphisms is characterised by a generator W 0 which is an element of the Cartan subalgebra h-for regular models, W 0 must be an A 1 -vector, i.e. the defining vector of a sl(2)-subalgebra of g. There is a remarkably wide variety of such actions for the case of su(N ), as noted by Bartnik [36]; and such A 1 -vectors are finite and have been tabulated [38,39].
The presence of a non-zero does not directly affect the automorphism classes on the bundle structure, and therefore some similar results to [34] will here be derived, as we describe how to express the field equations for these regular models. But does make a difference asymptotically, and so we find a big difference in the regularity requirements for solutions in the limit r → ∞ (as may be expected from previous treatments of su(N ) [9]); as such, we note that the definition of 'regular models' as given above must be amended a little for asymptotically adS space.
Let G from here on be a compact semisimple connected and simply connected gauge group with Lie algebra g. To consider spherically symmetric EYM connections is to consider principal SU(2) automorphisms on principal G-bundles E with base manifold M (our spacetime), such that the automorphisms project onto isometry actions in M whose orbits are diffeomorphic to 2-spheres. Since there is no natural action of SU (2) on E, we must consider all conjugacy classes of such automorphisms. These conjugacy classes are in one-to-one correspondence to integral elements W 0 of a closed fundamental Weyl chamber W ( ) belonging to a base of the roots of g with respect to a chosen Cartan subalgebra h [35,36,40].
Let g 0 be the (real) Lie algebra of the structure group G of the bundle E, so that g = (g 0 ) C , its complexification. Also, let {τ i }, i ∈ {1, 2, 3} be the standard basis of su(2) defined using the Pauli matrices, with commutator relations [τ i , τ j ] = i jk τ k , for i jk the Levi-Civita antisymmetric symbol. Then W 0 may be chosen such that where λ is the homomorphism from the isotropy group I x 0 of the SU(2)-action on M at the point x 0 ∈ M, determined by k ·π 0 = π 0 ·λ(k), ∀k ∈ I x 0 if π 0 ∈ π −1 (x 0 ), where π −1 (x 0 ) is the fibre above x 0 and the central dot notation denotes the adjoint action.
The subject of possible classes of connections over principal bundles has been covered in the literature by Wang et al. [41][42][43]. For instance, it is known that we may write the metric in common spherical Schwarzschild-type co-ordinates (t, r, θ, φ) as ds 2 = −μS 2 dt 2 + μ −1 dr 2 + r 2 dθ 2 + sin 2 θ dφ 2 . ( Note that we here consider only static solutions, meaning all field variables are functions of r alone. In addition, Brodbeck and Straumann [35] show that in this case a gauge may always be chosen such that the Yang-Mills one-form potential is locally given as A ≡ A μ dx μ =Ã + W 1 dθ + (W 2 sin θ + W 3 cos θ ) dφ.
In the above,Ã is a one-form defined on the quotient space of the manifold which is entirely parametrised by the (t, r ) co-ordinates, representing the 'electric' part of the connection. Here we consider the purely magnetic case, and hence we setÃ ≡ 0. We note that for = 0 this sector is not available in regular models [34]; it is available for < 0 but we find in the su(N ) case that the conditionÃ = 0 still yields a rich space of solutions [9]. Also, we have W 3 = − i 2 W 0 as the constant isotropy generator, and (2) gives us constraints on W 1 , W 2 (both also functions of r ), which we refer to as the Wang equations [42]. However, we still have a countably infinite number of possible actions of SU(2) on E: one for each element in W ( ) ∩ I , the intersection of the closed fundamental Weyl chamber and the integral lattice defined by I ≡ ker(exp | h ). Now for regular models, we require the YM fields to be non-singular at the centre r = 0 (for solitons) and asymptotically as r → ∞.
In the case of = 0, this implied that and/or where we define for i ∈ {1, 2}, k ∈ {0, ∞}. That is to say, for asymptotically flat space, in at least one of these limits (if they exist) there has to exist a Lie algebra homomorphism from su(2) into g 0 ; and if both limits exist, there also must exist a homomorphism between ∞ i and 0 i . The reason for the constraints (6) and (7) is that in asymptotically flat space, the values of the gauge field functions ω j at r = 0 and as r → ∞ (taken in a particular basis that we will describe) must be equal to a particular set of constants {λ j } that depend on the Cartan matrix of the reduced subalgebra in question. This implies that the soliton solutions have no magnetic charge, according to [37]. The constraints on the boundary values of the gauge fields are necessary so that the tangential pressure p θ and energy density e (see Sect. 3) remain regular at infinity.
However, for < 0 we have a different scenario. As we shall see, the values of the gauge field functions at the centre r = 0 are still highly constrained, reflecting the singular nature of that boundary, and thus (6) still holds; but asymptotically, the "fall-off" conditions required to force the gauge field to be regular are much laxer than for = 0, and thus the gauge field functions and their derivatives will in general approach arbitrary asymptotic values. Again this is due to the nature of the asymptotic system considered in a dynamical systems sense.
Our investigation in Sect. 5 will show that this lack of asymptotic constraints on the YM field is to do with the nature of the variable change that we perform to render the asymptotic field equations autonomous, which in the case of asymptotically flat space 123 necessitates the trajectory of every regular solution to end at a critical point (which we'll call * i , i = 1, 2) in the phase plane of the system. The critical points of the field equations are thus ω * 2 j = λ j for j = 1, . . . , L, where L = rank(g); the important point here being that for = 0, one is forced to have ∞ i ≡ * i (i = 1, 2), whereas for < 0, ∞ i = * i (i = 1, 2) in general. Hence, (7) does not have to hold for our solutions, and as we will see, this is manifested in the fact that for adS space, no constraints are placed on the gauge field functions or their derivatives as r → ∞, and we are allowed solutions with a global magnetic charge fixed essentially by the Cartan matrix of the reduced subalgebra, for which the tangential pressure and the energy density remain regular asymptotically. (Of course, (6) and (7) will both be trivially satisfied by embedded Schwarzschild anti-de Sitter solutions (see Sect. 6), and so for this solution at least, there must also exist a Lie algebra homomorphism from 0 i into ∞ i [34].) It must be noted though, it is still obviously true from the field equations that for regularity we must have Thus, for asymptotically adS space, the system itself still will possess the constraints (9) at the critical point * i , but solutions will not reach the critical point of the system in general, freeing the asymptotic solution parameters from the constraints that are seen in the = 0 case. This is what is responsible for the much larger space of black hole solutions in the su(N ) case, which we see need obey neither (6) nor (7); though we also emphasise that at the origin, regular solutions must still obey (6). Thus, as in the case of su(N ) for adS, we may expect the local existence proofs to be straightforward for r = r h and r → ∞ and much more involved at the origin r = 0. Now since W 3 is constant, (6) and (9) represent constraints also on W 3 , and hence on W 0 which must be the generating vector of an A 1 -(i.e. sl(2)-) subalgebra of g. However the set of such so-called A 1 -vectors is finite, and have been tabulated by Dynkin [38] and Mal'cev [39] using what they call "characteristics", which are in one-to-one correspondence with finite ordered sets of integers chosen from the set {0, 1, 2}. These strings of integers then represent the value of the simple roots on W 0 , the defining vector of the A 1 -subalgebra, chosen so that it lies in W ( ); and the tables of Mal'cev and Dynkin therefore give us a classification of all possible spherically symmetric, purely magnetic EYM models which obey the correct regularity conditions asymptotically and at the centre, for any compact semisimple simply connected gauge group.

Field equations in the case of the 'regular' action
To proceed, we can note that out of all the possible actions classified by Dynkin and Mal'cev [38,39], these exists a privileged class of actions which corresponds to a principal A 1 -vector in Dynkin's terminology, which Oliynyk and Künzle [34] called principal actions. There exists a slightly larger class of actions called 'regular' in [34] (and 'generic' in [37]), for which the defining vector lies in the interior of a fundamental Weyl chamber. (The other irregular case involves the defining vector being on the boundary of a Weyl chamber.) In this section we will show that for < 0, as it was for = 0, all models with a regular action can be reduced to those with the principal action, for any semisimple gauge group. In terms of the field variables, this means that the YM potential can be chosen to be composed of real functions due to a gauge freedom, and that there are L of such functions where L = rank(g). We also have two metric functions governed by the Einstein equations: m (the mass function) and S (the lapse function). Then the field equations are determined by L + 2 real functions of the radial co-ordinate r alone (for static, spherically symmetric solutions), and possess singularities at the centre r = 0, the event horizon r = r h and as r → ∞.
A more convenient basis to use here for the Wang equations (5) in place of the generators W 1 and W 2 is in which case equations (5) become Then W ± (r ) are g-valued functions, W 0 is a constant vector in a fundamental Weyl chamber of h, and {W 0 , W ± } is a standard su(2) triple in the limit r = 0 and at the critical points of the system. Also, h is the Cartan subalgebra of the complexified form of the Lie algebra, i.e. h = h 0 + ih 0 , for h 0 the real Cartan subalgebra of g 0 , which in turn is the real compactified form of g. Naturally, we introduce a complex conjugation operator c : g → g with convention This implies that W − = −c(W + ).
Therefore the solutions will only depend on the functions m(r ), S(r ) and the complex components of W + (r ). The field equations in the case = 0 are well-known [34,35]. It is not difficult to use the general adS Einstein and YM field equations to derive the analogous forms for < 0. These general field equations are also well-known: where g μν is the metric tensor defined using (3), G μν is the Einstein tensor, F λ μ is the mixed anti-symmetric field strength tensor defined with

123
A μ represents the YM one-form connection (4), and the energy-momentum tensor T μν is given by We note that Tr is the Lie algebra trace, we have used the Einstein summation convention where summation occurs over repeated indices, and we have rescaled all units so that (for the gauge coupling constant q).
Using (3), (4) and (14), we may show that the field equations for < 0 become with ≡ d/dr, and , the adS radius of curvature, given by only valid for < 0. In (19), ( , ) is an invariant inner product [relating to the Lie algebra trace in (16)] on g determined up to a factor on each simple component of a semisimple g, which induces a norm | | on (the Euclidean) h and therefore also on its dual. These factors are chosen so that ( , ) is a positive multiple of the Killing form on each simple component. We may calculate the energy density e, the radial pressure p r and the tangential pressure p θ . As we mentioned in Sect. 2, these are important quantities which help us assess the physicality of our solutions. First we note that since c(F) =F, and X | Y ≡ −c(X ), Y ) is a Hermitian inner product on g, then G ≥ 0 and P ≥ 0. Then, we have [in our units (17)] e = r −2 (μG + r −2 P), p r = r −2 (μG − r −2 P), p θ = r −4 P.
Now we describe how to reduce the field equations down to the case of a regular action as described above. We select a Chevally-Weyl basis for g. Let R be the set of roots on h * and = {α 1 , . . . , α L } be a basis for R (where L is the rank of g). We also define Then {h i ≡ h α i , e α , e −α | i = 1, . . . , L; α ∈ R} is a basis for g, and induces the decomposition for R + , the set of positive roots expressed in the basis . For this decomposition, we adopt the conventions From the commutator relations defining an sl(2)-subalgebra span{e 0 , e ± } of g, i.e.
Hence, writing W 0 in the basis then equations (11) imply that where we have defined λ , a set of roots depending on the homomorphism λ (or equivalently the constants λ i ), as In a similar way we find that 123 for functions α (r ), but given that the complex conjugation operator c maps h i → −h i , e α → −e −α , we easily see that Therefore, the system is determined by two real functions m(r ), S(r ) and L complex functions ω α (r ), ∀α ∈ λ . It is noted in [34] that we may naïvely proceed by substituting the expansion (29) into the field equations and calculate the various Lie brackets using (24), but this may produce many more equations that unknowns, and in addition there is still some gauge freedom left in the connection A. However we may simplify the system a great deal by considering only the so-called regular case, where W 0 is a vector in the open fundamental Weyl chamber W (S) [37]. We begin with a theorem due to Brodbeck and Straumann: Theorem 1 [35] If W 0 is in the open Weyl chamber W ( ) then the set λ is asystem, i.e. satisfies: and is therefore the base of a root system R λ which generates a Lie subalgebra g λ of g spanned by If W 0 is an A 1 -vector then W ⊥ 0 = 0 (though h ⊥ λ need not be trivial).
This allows us to rewrite the field equations in a much simpler form -in fact, in a form that renders them very similar-looking to the well-studied su(N ) case. First we can consider W + to be a g λ -valued function, and write where we now take {α 1 , . . . ,α L λ } as the basis for λ and defineẽ j ≡ eα j . This means that using (24), (18d) becomes implying that the phase of ω j (r ) is constant and can be set to zero using a gauge transformation. Hence we can conclude that the ω j (r ) may we taken as real-valued functions. We note that in = 0, this is only possible for the regular case [34]. Also using this basis, we may define the Cartan matrix of the reduced subalgebra g λ as noting that by definition this is a symmetric and positive operator.
The results in Sect. 3 of [34] depend only on the root structure of the reduced subalgebra, and therefore we may also apply the same logic when reducing the field equations (18) to the regular case. Finally then, dropping tildes from α j and losing the λ index from g et cetera for clarity, we can show that the field equations become with The final step is to determine the values of the constants λ j , which involves determining the subalgebra g λ for a given A 1 -vector W 0 in the open fundamental Weyl chamber. For a semisimple group, for which the Cartan subalgebra splits into an orthogonal sum h = k h k , the orthogonal decomposition given in Theorem 1 splits into analogous decompositions of each of h k . Hence we only need consider the regular actions of simple Lie groups. However, we note that the A 1 -vector in the Cartan subalgebra h of a Lie algebra g is uniquely determined by the integers 123 which integers are chosen from the set {0, 1, 2}. In [38], this is referred to as the characteristic. From (30), it is obvious that for the principal action, for h λ . A 1 -vectors satisfying this define principal su(2)-subalgebras, and hence principal actions of SU(2) on the bundle. As in [34], we may rely the following theorem: Theorem 2 [34] (i) The possible regular su (2) The essence of this theorem is that the regular action here coincides with the principal action. This finally allows us to determine an expression for the constants λ j , derived by using (38b), (38e), (40), and (41):

Boundary conditions
In order to get a sense of the possible term dependencies in the power series expansions of the field variables near the boundary points, and thus decide what methods we will need to prove local existence, it is very enlightening to calculate the lower order terms in the power series expansions of the field variables nearby the boundaries r = 0, r = r h and r → ∞. We do this below, in anticipation of the later proofs of local existence at these points in Sect. 7.
In the black hole case, i.e. for the boundaries r = r h and r → ∞, we find that the situation is relatively uncomplicated. For r = r h , the lower order terms show that the solutions can be characterised entirely by the values of Asymptotically, we find that the solution is parametrised entirely by the values of the limits of m(r ), ω j (r ) and r 2 ω j (r ) ( j = 1, . . . , L) as r → ∞. We find no constraints on the boundary values of the field variables asymptotically, and near r = r h , we merely find a couple of constraints on the metric function μ(r ) that must be satisfied, which are physically necessary to ensure a regular and non-extremal event horizon.
In the soliton case however, i.e. at r = 0, the situation is much more complicated, as it was in the su(N ) case [9,31]. There, we had to solve a tridiagonal matrix equation by using expansions in the eigenvectors of the matrix in question; for this we used Hahn polynomials, an orthogonal class of polynomials defined using hypergeometric functions [44]. In that case, as in this, appears at O(r 2 ) and above in the field equations (18a)-(18c), and therefore near r = 0 we do not expect the appearance of the cosmological constant to make any appreciable difference.
In light of all of this, we now review the boundary conditions we expect in each case.

Origin
Near r = 0 we may simply use the independent variable r , and hence we expand all field variables and quantities as for a general function f (r ). Thus we obtain the following recurrence relations for m k+1 , S k and ω j,k+1 : Here, A ≡ A i j is the matrix defined by δ i j is the Kronecker symbol; and the left-hand side of (43c), the vector , is a complicated vector expression involving the coefficients of the field variable expansions.
We can see that these equations are identical to the su(N ) case [41], and so again, we may solve (43a) and (43b) and obtain a solution with L free parameters on condition that the recurrence relations (43c) can be solved. This in turn is conditional upon the vectors b k lying in the left kernel of the matrix A. As we noted, b k is a complicated expression and so this is difficult to prove in general. In Sect. 7.1, we generalise proofs in [34] which depend directly on the root structure of the Lie algebra g treated as an sl(2, C) submodule.
We note here that G k = P k = 0 for k < 2. For the lower order terms, we find: 123 Table 1 This table shows For the classical Lie algebras the table shows k j for j = 1, . . . , L, L = rank(g). Note that k = 1 belongs to all Lie algebras, thus 1 ∈ E always The Eq. (43) are identical to those we found in the su(N ) case, therefore we expect a similar situation to occur here, in that the higher order terms of the power series expansions near the origin will in general display a complicated interdependence. This reflects the fact that r = 0 is a singular point of the field equations. At this boundary, the higher order coefficients which remain arbitrary occur at the orders r k for which k(k + 1) is an eigenvalue of the matrix A. But in fact, the eigenvalues of A can happily be shown to be k(k + 1) for a series of integer values of k, which series depends on the Lie algebra in question. (For su(N ), this series of integers is simply the natural numbers from 1 to N − 1 inclusive.) For all the simple Lie algebras, we may calculate the spectrum of eigenvalues from the Cartan matrix by using the definition (44)-see Table 1 for this information. The proof for the classical Lie algebras then follows from the properties of the root structure and the results at the end of Sect. 7.1.1.
We will see in Sect. 7.1.2 that in some neighbourhood of r = 0, the relevant field variables have the following behaviour: Here, Q i j is a non-singular matrix, k j are integers andû j are some functions of r -all of these we will define later. Also, m 3 is fixed by (43a), S 0 is fixed by the requirement that S → 1 as r → ∞, and ω 2 j,0 = λ j . Therefore altogether we have L free solution parameters here in total, namelyû j (0) for each j.

Event horizon
For a regular non-extremal event horizon, we require μ h to vanish and μ h to be finite and positive. This severely restricts the solution parameters here and hence reduces the degrees of freedom of any solution, which makes boundary conditions easy to find.
Using the notation f h ≡ f (r h ) and transforming to a new variable ρ = r − r h , we find that where The constraint μ h = 0 implies that with The condition μ h > 0 places a bound on m h : with Therefore, it is clear that fixing r h and , and regarding S h as fixed by the requirement that the solution is asymptotically adS, the solution parameters are given by the set {ω j,h }. Thus, as at the origin, we have L solution degrees of freedom for solutions existing locally at the event horizon.

Infinity
We assume power series for all field variables which are good in the asymptotic limit, i.e. of the form f (r ) = f ∞ + f 1 r −1 +· · · . It is easy to see that this implies G = O(r −4 ), 123 meaning that examining (18b), S must be of the form S(r ) = S ∞ + O(r −4 ). We also use the basis W + (r ) = L j=1 ω j (r )e α j . Therefore, we find that the expansions near infinity must be The power series expansions here are a lot less complicated than for the asymptotically flat case. No constraints appear on ω j,∞ or c j . Similarly, no constraints are placed on S ∞ or m ∞ , so we rescale to S ∞ = 1 and let m ∞ = M [the constant Arnowitt-Deser-Misner (ADM) mass] so that the solution asymptotically is the SadS solution (or pure adS space if M = 0). We find that each new term we calculate in the expansions is entirely determined by previously calculated terms, and this trend continues for higher order terms. For instance, the lower order terms are Therefore we anticipate that proving the existence of unique solutions to the boundary value problem will be a lot less involved than in the case of = 0. In summary, our solution parameters here are {M, ω j,h , c j } and thus we have 2L + 1 degrees of freedom in total.

Asymptotic behaviour of the field equations
As we saw, the asymptotic boundary conditions (53) imply that any regular solutions in this limit will have gauge functions which are characterised entirely by the arbitrary values ω j,∞ and c j , with all higher order terms in the expansions determined by these parameters. This is in opposition to the = 0 case, where the asymptotic values of the gauge field have to approach particular values, and the higher order terms display complicated interdependence related to the intercoupling of the gauge functions caused by Eq. (43c).
Therefore what we wish to do now is take the asymptotic limit of the field equations, transform the independent variable r so that the system becomes 'autonomous' in the dynamical systems sense, and examine the nature of the phase plane of the system. As we will see, it is not so much the asymptotic field equations themselves which give us the difference in behaviour between the = 0 and < 0 cases-it is the form of the parameter we must transform to which dictates the asymptotic behaviour of the field variables, and which gives us an infinitely more plentiful space of regular solutions.
First, we note that as r → ∞, μ ≈ 1 + r 2 2 . Noting also (53), the YM field equations (18c) become asymptotically Using the parameter τ = r −1 , we find that (55) becomes In the more explicit basis (29) using the regular action, defined in Sect. 3 where the field equations become (37), this is equivalent to It is easy to see that the critical points ω * j of this autonomous system satisfy Noting that C i j is of full rank, this gives us two sets of critical points: either ω * j = 0, or ω * j = ±λ Eigenvalue analysis shows these (for each j) to be a centre and a pair of saddles, respectively. We noted that the analysis of the asymptotic boundary conditions (53) implied no such constraints on the asymptotic value of ω j (r ), though the autonomous asymptotic equations (56) are identical to those for = 0.
We may resolve this apparent discrepancy by noting that for < 0, the trajectory of a solution in the phase plane ω j , dω j dτ will not in general reach its critical point. This is due to the nature of the parameter we used to render the equations autonomous. In the case of = 0 the parameter used was τ ∝ log r , so that the range r ∈ [r 0 , ∞) (r 0 = r h for black holes, or r 0 = 0 for solitons) corresponds to τ ∈ (−∞, ∞), and hence any trajectory for a regular solution in the limit r → ∞ will be destined to end at a critical point.
For < 0 however, we use τ ∝ 1/r , meaning that the range r ∈ [r 0 , ∞) corresponds to the range τ ∈ [0, r −1 0 ). Therefore, as we take the asymptotic limit r → ∞, the corresponding trajectories in terms of τ will shrink and only traverse a short distance in the phase plane. Hence the trajectories, and therefore the values of the gauge field functions and their derivatives, will in general approach arbitrary values asymptotically. We note that this is precisely the same as in the su(N ) case [9].
In summary then, our investigation has shown that we need not be concerned with the behaviour of the field equations for r arbitrarily large-as long as we can integrate into the asymptotic region, the solution will remain regular until reaching the (arbitrary) boundary conditions at r → ∞. We will return to this point in Sect. 8.
It may finally be noted that since we are not concerned with the nature of the critical points, we could have stopped at Eq. (56); so this argument therefore applies also to the irregular case, i.e. if the defining A 1 -vector W 0 lies on the boundary of a Weyl chamber.

Embedded solutions
Our argument in Sect. 8 will rely on the existence of embedded (or 'trivial') solutions, as we will prove the existence of global solutions to the field equations (37a) to (37c) in some neighbourhood of these. Therefore, we here review some easily obtainable embedded solutions to our field equations.

Reissner-Nördstrom anti-de Sitter (RNadS)
Here we let ω j (r ) ≡ 0. In that case, we find that G = F = 0 and therefore S becomes a constant, which we scale to 1. The metric function μ(r ) becomes where M is the ADM mass of the solution, and the magnetic charge Q is defined with Therefore we have obtained the embedded Reissner-Nördstrom anti-de Sitter solution, which only exists with this value of Q 2 , and coincides with the su(N ) case [9], using (41) and the su(N ) Cartan matrix.
To summarise, the RNadS solution is given by

Schwarzschild anti-de Sitter (SadS)
Here we let ω 2 j (r ) ≡ λ j , ∀r, ∀ j = 1, . . . , L. Then from (38) we find that P = G = F = 0, implying the following. From (37a), we get m (r ) = 0, so that m(r ) is a constant which we again set to the ADM mass M. From (37b) we have S (r ) = 0, so that S is a constant which we scale to 1 for the asymptotic limit. Finally, the YM equations (37c) are automatically satisfied. Since P = 0, this solution carries no global charge, and can be identified as the embedded Schwarzschild anti-de Sitter solution.
Summarising this solution:

Embedded su(2) solutions
Noting that we can embed SU(2) isomorphically into any semisimple gauge group G, then there must always exist trivial embedded su(2) solutions to the field equations (18a) to (18c). We may show this by a simple rescaling.

Proposition 3
Any solution to the field equations (18a)-(18c) can be rescaled and embedded as a solution which satisfies the field equations for su(2) adS EYM theory.
Proof Consider the gauge group G, fixing the symmetry action such that W 0 is regular. Select any basis such that the set {W 0 , + , − } spans su (2), with c( + ) = − − . We rescale the field variables as follows: with Q 2 given in (60). Then the field equations (18a)-(18c) become with These equations are identical to those for the su(2) adS case, for which the existence of (nodeless) solutions has been proven [22].
It is interesting to note that the scaling involves the magnetic charge itself, which can possibly be put down to the fact that the RNadS solution for su (2), embedded in the su(2) equations, only exists where the magnetic charge Q 2 = 1.
Finally, it should also be noted that using the definition of the Cartan matrix for su(N ), i.e.
and normalising so that the length of the long roots |α k | 2 = 1 ∀k, the field equations (18a) to (18c) yield exactly the su(N ) adS EYM equations [9].

Local existence proofs at the boundaries
Now we have much information about the behaviour of the solutions to the field equations nearby the boundaries of our spacetime, enough to prove local existence at those boundaries. To do this, we rely on a well-known theorem of differential equations [10], generalised to the appropriate case by [34].
Theorem 4 [34] The system of differential equations where μ i , ν j ∈ Z >1 , f i , g j are analytic functions in a neighbourhood of (0, c 0 , 0) ∈ R 1+m+n , and the functions h j : for |t| >r for somer > 0 if |c − c 0 | is small enough. Moreover, the solution depends analytically on the parameters c i .
Essentially, the proof of this theorem proceeds from the requirement that formal power series may be found for the field variables at the boundaries in question. We now consider those boundaries one by one.

Existence at the origin: r = 0
As we hinted in Sect. 4, we do not expect much of a difference between the asymptotically flat and asymptotically adS cases nearby the origin, because as r → 0, the terms in the field equations involving the cosmological constant become negligible. Hence we may proceed along very similar lines to those in [34]. Therefore, we now collect all necessary results from [34] needed to prove local existence of solutions near r = 0. The general idea is to consider the root structure of sl(2, C) taken as a Lie algebra submodule of g. Note that the results in this section are only necessary for this boundary, and hence only for solitons.

Necessary results for local existence at r = 0
First we introduce our conventions. We begin by defining a non-degenerate Hermitian inner product | : g × g → C, such that Then | is a real positive definite inner product on g 0 , since c : g → g is the conjugation operator determined on the compact real form g 0 . It is elementary to show that | satisfies for all X, Y, Z ∈ g. Now we introduce a positive definite, real inner product | : Let be the norm induced by (71), i.e. X 2 = X | X ∀X ∈ g. Then we can easily verify the following properties of | : for all X, Y, Z ∈ g. Let + , − ∈ g be two vectors such that Then span C {W 0 , + , − } ∼ = sl(2, C). We again use a central dot notation · to represent the adjoint action, i.e.
But since W 0 is a semisimple element, ad(W 0 ) is diagonalisable, and so from sl (2) representation theory we know that the eigenvalues are integers. Therefore we define V n as the eigenspaces of ad(W 0 ), i.e. with It also follows from sl(2, C) representation theory that if X ∈ g is a highest weight vector of the adjoint representation of span C {W 0 , + , − } with weight n, and we define X −1 = 0, X 0 = X and Now we are ready to state a series of results proven in [34] which will help us to prove existence locally at r = 0. Essentially, these are necessary because we find that the term F in the YM equation (18c) is the only term which resists our rearrangement of the field equations in a form appropriate to Theorem 4, and it is necessary to argue that certain lower order term of F (in a power series sense) are zero. Hence we proceed.

Proposition 6
The R-linear operator A : g → g defined by is symmetric with respect to the inner product This shows that the operator A restricts to V 2 : we therefore denote this operator by Now we label the set of integers k j from Proposition 5 as follows: where we define the series of integers J 1 = 1, J k + m k = J k+1 for k = 1, . . . , I and J I +1 = − 1. To ease notation we define κ j ≡ k J j , for j = 1, . . . , I.
As noted in Proposition 5, the set {ξ Then due to Proposition 6, A is symmetric, and so also is A 2 , and hence A 2 must be diagonalizable. Then the following Lemma is true.
We now define the spaces and Then E 0 = ker(A 2 ) and E l + is the eigenspace of A 2 corresponding to the eigenvalue κ j (κ j + 1). Also, from Proposition 5 (iv) we see that Proposition 14 Let W 0 be regular. Then if + ∈ α∈ λ Re α , E + = α∈ λ Re α .

Proof of local existence at the origin (r = 0)
Now we use Theorem 4 and the results of Sect. 7.1.1 to prove the existence of solutions, unique and analytic with respect to their boundary parameters, in some neighbourhood of the origin. We begin by introducing some necessary notation, which will be used throughout this section. First, we define the set for κ j given in (81); and a set of projection operators between the spaces defined in (84) where for our purposes, > 0 is small. Using Proposition 14 and Eq. (35), we know that the solution W + (r ) of Eq. (18c) is completely characterised by the condition We noted previously that Eq. (18b) decouples from the others, so that once we have solved Eqs. (18a) and (18c) for μ and W + , we may easily solve (18b) to give S. However, for completeness, we shall include S in our analysis. We now have everything we need to state our Proposition: Above, Q i j is a non-singular matrix for which the jth column is the eigenvector of the matrix A (44) with eigenvalue k j (k j + 1), andû j (r ) are some functions of r . Each solution is entirely and uniquely determined by the initial valuesû j (0) ≡ β j , for arbitrary values of β j . Once these are determined, the metric functions m(r ) and S(r ) are entirely determined.
Proof Since W + (r ) ∈ E + , we introduce new functions u k (r ) with with + = W + (0) and u s+1 (r ) ∈ Es + ∀r, ∀s ∈ E. This transformation is clearly invertible since E + = I q=1 E q + . Define Then we may write (92) as W + (r ) = + + ∞ k=0 χ k u k (r )r k . Substituting this into the YM equations (18c), we find: for some N 1 ∈ Z, and The need for the results of Sect. 7.1.1 becomes apparent if we examine those results alongside the forms of (94) and (95). Now since A 2 (u k+1 ) = k(k + 1)u k+1 , (94) becomes We proceed by defining new variables v s+1 ≡ u s+1 , ∀s ∈ E. The YM equations (18c) become 123 Now we apply projection operators pk + (88) to equations (97) for each k ∈ E, giving for all k ∈ E. The main obstacle to writing this equation in the correct form for Theorem 4 is the final term, as was the case for su(N ) [9,41]. As written it contains terms of much lower order than we want, i.e. terms of order r −s where s > 0. Happily we may rewrite the final term using the following equality: We make the derivation of this plain by using the results from Sect.
because if k ∈ E, then k = κk and so if s < k = κk, thens <k, proving (99). Using (99) in (98) and rearranging gives Using the properties of | and the fact that A 2 (u 2 ) = 2u 2 , we can show that there exist analytic functionŝ with P = r 4 u 2 2 + r 5P (u, r ), and where u = s∈E u s+1 , v = s∈E v s+1 , and X 2 = X |X .

Now we rewrite the Einstein equations (18a, 18b). We introduce a new mass variable
(We know that u 2 is always defined since κ 1 = 1 always and hence 1 ∈ E.) Then (18a, 18b) become We make one last variable change: pk + ( f k+2 ). (106) We proceed by fixing a vector X ∈ E + and definev = s∈Ev s+1 . Then from (101, 104, 106), we can show there exists a neighbourhood N X of X ∈ E + , some > 0, and a sequence of analytic maps Also, with (105, 106) and using v s+1 = u s+1 , there exist analytic maps such that where Y s = ps + (Y ). It is helpful to note that from the definition of M (104), we can show that m(r ) = O(r 3 ), and so in this regime, 1 (98)]. Also, it is easy to see from (102, 106, 111) that From the results of Sect. 7.1.1, there must exist an orthonormal basis {w j | j = 1, . . . , } for E + consisting of the eigenvectors of A 2 , i.e. A 2 (w j ) = k j (k j + 1)w j . So we introduce new variables in this basis: From Proposition 5, we know that = | λ |, so we can write λ = {α j | j = 1, . . . , }; and from Proposition 14, we find that {e α j | j = 1, . . . , } is also a basis for E + . Therefore we can write With this definition of the matrix Q i j , it is clear that the columns of Q i j are the eigenvectors of A 2 . Now we expand + and W + (r ) in the same basis: Then Eqs. (92, 113, 114, 115) imply that with ω 2 i,0 = λ i . Finally, from (111) and (113) we obtain with β j (Y ) ≡ w j |Y . Therefore, we obtain the expansions (91).

Proof of local existence at the event horizon r = r h
Here, the situation is again quite similar to the asymptotically flat case [34]. Therefore, as was the case in [34], we have no need of the results in Sect. 7.1.1. In particular, the space E + that we will use does not have to be of the form defined in (85)-we may replace E + everywhere in the following with α∈ λ Re α , and it is not necessary to know that E + = α∈ λ Re α (which is the essence of Proposition 14). Thus, we use the notation E + purely for convenience. We begin by introducing the variable so that for r → r h we are considering the limit ρ → 0. Keeping in mind the boundary conditions in Sect. 4.2, we prove the following Proposition:

Proposition 16
In a neighbourhood of the event horizon r = r h = 0 (i.e. ρ = 0), there exist regular black hole solutions to the field equations (18a)-(18c), analytic and unique with respect to their initial values, of the form where μ h > 0.
Proof Along with (118), we introduce some new variables: forλ, V + functions of ρ, and ν some constant yet to be determined. Immediately we have and it is clear that there exist analytic mapsF : E + → E + ,P : E + → R, witĥ Define an analytic mapĜ : 123 Then we can see that G =Ĝ(V + ,λ). Using these we can rewrite the EYM equations (18a) to (18c) as In order to cast the equations in the form necessary for Theorem 4, we introduce some final new variables:λ We continue by defining an analytic map γ : it is obvious that γ (Y, 0) = 0. Therefore, define an open neighbourhood D of (Z , 0) ∈ Then from (121, 124, 125) we can show there must exist some > 0 and analytic maps such that It can be seen that equations (130) are in the form applicable to Theorem 4. Hence there is a unique solution To gain a more explicit solution, we expand Z , W + in the basis {e α j | j = 1, . . . , }, as follows: Noting (131a), this yields Finally, it is easy to show from (120a, 125a, 131c) that and hence Therefore, we have obtained the expansions (119).

Proof of local existence as r → ∞
The behaviour of solutions in the asymptotic limit is the biggest difference between the asymptotically flat and adS cases. Because of the constraints on the asymptotic values of the gauge functions for = 0, the proof followed a similar route to the local existence at the origin. However for < 0, our situation is much more similar 123 to the local existence at the event horizon, so we follow a similar method to that used in Proposition 16 from Sect. 7.2. Hence, the same comments apply as at the beginning of Sect. 7.2: we do not need any of the results of Sect. 7.1.1 here, and thus we use the notation E + out of utility.
To deal sensibly with the limit r → ∞ we transform to the variable whence we are now dealing with the limit z → 0. We state our Proposition:

Proposition 17
There exist regular solutions of the field equations in some neighbourhood of z = 0, analytic and unique with respect to their initial values, of the form for arbitrary constants ω j,∞ , c j ; where in order to agree with the asymptotic limit of adS space, we have let m ∞ = M, the ADM mass of the solution, and S ∞ = 1.
Proof As well as (136), we introduce also the following new variables: We immediately find that and it is clear that there exist analytic mapsF : E + → E + andP : E + → R witĥ Also we find that which means that For λ and v + , it can be shown that It is useful to note that in the asymptotic limit, μ ∼ 1 + 1 z 2 2 , from which we may see that 1 μz 2 2 − 1 = O(z 2 ), and Examining the number of degrees of freedom we expect at this boundary, we fix two vectors X, C ∈ E + . Then from results (138b)-(144), it is clear that there exists an > 0 and analytic maps with (noting that G ∞ is just the map v + → −v + ). Now we are at the stage where we may apply Theorem 4; and hence it is clear that these equations possess a unique solution However, noting (136) and (138b), we may integrate (147d), choosing the constant (vector) of integration to agree with (147c). This combines (147c) and (147d), yielding To gain an explicit solution in terms of the components of X , C and W + , we expand them all in the same basis: 123 No constraints are placed on the constants ω α,∞ or c α . Then it is clear that near z = 0, the gauge field functions have the form Finally, noting that we expect our solution to approach adS space in the asymptotic limit, we set λ ∞ ≡ 2M, S ∞ ≡ 1, and thus recover the expansions (137).

Global existence arguments
Now we turn our attention to proving the existence of global solutions to our field equations. Here we have a choice of approaches. We considered using the more novel approach of Nolan and Winstanley [29] who let the initial conditions and embedded solutions reside in appropriate Banach spaces, and then recast the field equations so that they could apply the Implicit Function Theorem, hence proving that non-trivial solutions exist in some neighbourhood of embedded solutions. However, it appears to be necessary to their argument that m(r ) is constant for the embedded solution, something we have not been able to get around yet, meaning that we could only identify solutions in a neighbourhood of the embedded SadS solution.
Alternatively, the traditional argument that has been used in this case is the 'shooting argument' (used in e.g. [22,26]), which basically involves proving the existence of solutions locally at the boundaries, and then proving that solutions which begin at the initial boundary r = r h (r = 0) near to existing embedded solutions can be integrated out arbitrarily far, remaining regular right into the asymptotic regime, where they will 'meet up' with solutions existing locally at r → ∞; and that these neighbouring solutions will remain close to the embedded solution. While this seems somehow less elegant, there are no restrictions on the embedded solution we may use, and hence the proof we are able to create is more general and hence more powerful. Therefore, we resign ourselves to using the more traditional techniques.
We begin by noting that we have already considered the behaviour of the field equations in the asymptotic limit and shown that solutions will in general remain regular in this regime (Sect. 5), so we must now make sure that any solution which begins regularly at the initial boundary r = r h (r = 0) can be integrated out arbitrarily far while the field variables remain regular. We also note that as in Sect. 5, we here do not require W 0 to be regular: we use the original field equations (18), and so this proof applies to both the regular and irregular actions.

Proposition 18
If μ(r ) > 0 ∀r ∈ [r h , ∞) for black holes, or ∀r ∈ [0, ∞) for solitons, then all field variables may be integrated out from the boundary conditions at the event horizon (or the origin) into the asymptotic regime, and will remain regular.
Proof Define Q ≡ [r 0 , r 1 ) andQ ≡ [r 0 , r 1 ], where r 0 = r h for black holes and r 0 = 0 for solitons, and r 0 < r 1 < ∞. Our strategy is to assume that all field variables are regular on Q, i.e. in a neighbourhood of r = r 0 , and then show using the field equations that as long as the metric function μ(r ) > 0 ∀r ∈ [r 0 , ∞), then they will remain regular onQ also, i.e. at r = r 1 ; and thus we can integrate the field equations out arbitrarily far and the field variables will remain regular.
First notice that G, P > 0 by the definitions (19). This means using (18a) that m (r ) > 0 ∀r and thus m(r ) is monotonic increasing, as expected for the physical mass. This means that (if it exists), The same applies to (ln |S(r )|) [see (18b)], showing that ln |S(r )| and hence S(r ) is monotonic increasing too, so that (again, if we can prove that S is finite onQ) The condition μ(r ) > 0 ∀r ∈ [r 0 , ∞) gives us our starting point, since this implies that giving us an absolute upper bound to work with. This in turn implies that m(r ) is bounded onQ [and so (151) holds], and thus also that μ(r ) is bounded onQ. Thus we may define μ min ≡ inf{μ(r ) | r ∈Q}. Now we examine (18a). It is clear that and integrating, we can show that which implies from (18b) that ln |S| and hence S is bounded onQ. Equation (155) also implies that G is bounded onQ, and since then again by integrating and using the Cauchy-Schwartz inequality, and hence r 1 The left hand side is bounded, and the right hand side is a sum of positive terms and hence bounded below by 0. Thus W + and hence W + is bounded onQ. Since W 0 is constant and W − = −c(W + ), this also means thatF and hence F and P are similarly bounded onQ (see (19)). Finally, we may rewrite the YM equations (18c) as Integrating and rearranging gives and since all functions on the right hand side are bounded onQ (see (19)), as are μ and S, then we can finally conclude that W + is bounded onQ.

Global existence of solutions in a neighbourhood of embedded solutions
Finally, we may prove the major conclusions of our research, which hinge on the following Theorem. The gist of it is that global solutions to the field equations (37a)-(37c), which we have proven are uniquely characterised by the appropriate boundary values and analytic in those values, exist in open sets of the initial parameter space; and hence that solutions which begin sufficiently close to existing solutions to the field equations will remain close to them as they are integrated out arbitrarily far into the asymptotic regime, remaining regular throughout the range. It can be noted that this argument is quite similar to those we have used for the su(N ) case [9,31]. Proof Assume we possess an existing solution to the field equations (37a) to (37c), where each gauge function ω j (r ) has n j nodes and initial conditions ω j,0 = 0 in general. Proposition 18 and the analysis in Sect. 5 show that as long as μ(r ) > 0 we may integrate this solution out arbitrarily far into the asymptotic regime to obtain a solution which will satisfy the boundary conditions as r → ∞. For the rest of the argument, we assume that is fixed and so is r h for black holes and that each gauge function ω j has n j nodes. From the local existence results (Propositions 15, 16 and 17), we know that for any set of initial values, solutions exist locally near the event horizon for a black hole, or the origin for a soliton, and that they are analytic in their choice of initial conditions. Again we use the notation r 0 = r h for black holes and r 0 = 0 for solitons. For an existing solution, it must be true that μ(r ) > 0 for all r ∈ [r 0 , ∞). So, by analyticity, all sufficiently nearby solutions will also have μ(r ) > 0 for all r ∈ [r 0 , r 1 ] for some r = r 1 with r 0 < r 1 < ∞. By Proposition 18, this nearby solution will also be regular on [r 0 , r 1 ]. Now, let r 1 >> r 0 , so that for the existing solution, m(r 1 )/r 1 << 1. Let {ω j,0 } be a different set of initial conditions at r = r 0 for gauge fieldsω j , such that {ω j,0 } are in some small neighbourhood of {ω j,0 }; and letm(r ) be the mass function andμ be the metric function of that solution. By analyticity (as above),μ(r ) > 0 on this interval, so this new solution will also be regular on [r 0 , r 1 ]; and since the two solutions must remain close together, the gauge functionsω j will also each have n j nodes.
Also it is then the case thatm(r 1 )/r 1 << 1, and since r 1 >> r 0 we consider this the asymptotic regime. Provided r 1 is large enough (and hence τ 1 is very small), the solution will not move very far along its phase plane trajectory as r 1 → ∞ (see Sect. 5). Thereforem(r )/r remains small, the asymptotic regime remains valid, and the solution will remain regular for r arbitrarily large.

Corollary 20
Non-trivial solutions to the field equations which are nodeless, i.e. for which ω j (r ) = 0 ∀r , exist in some neighbourhood of both existing trivial SadS solutions (described in 6.2), and embedded su(2) solutions (proven in Proposition 3).

Existence of solutions in the large | | limit ( → 0)
So far we have proven the existence of global black hole and soliton solutions in some neighbourhood of existing solutions, for fixed r h and . But there is a further consideration, revealed by investigations into su(N ). On the one hand, we discovered numerically that as N increases, regions of the parameter space in which we may find nodeless solutions shrink in size [9,45]; on the other, for | | large enough, all solutions we found were nodeless. In addition, when we investigated the linear stability of these solutions [25], we were only able to prove stability in the limit | | → ∞, due to terms arising in the gravitational sector.
In view of the similarities between the case under consideration and the su(N ) case, it is sensible to investigate this limit in the case of a general compact gauge group. Our strategy is to transform the field variables such that we may sensibly find a unique solution to the equations at = 0. Then, noting that it is only in the asymptotic limit that the influence of is felt, we modify Proposition 17 using our new variables, and show that the arguments used in Sect. 8 may be easily adapted to serve in a neighbourhood of = 0.
We must emphasise that we cannot prove the existence of global non-trivial solutions at = 0, since in that case the asymptotic variable we used in Sect. 5 becomes singular and therefore that part of the proof breaks down. Proof We'll take the black hole case to begin with, noting that we fix r h for the rest of the argument. Let us change to the variables The field equations (18a)-(18c) then become Taking the (now allowed) limit → 0: The first of these is easily integrated to givem constant, which we therefore set tō m(r ) =m h . We also notice that sincē then we must havem(r ) = r 3 h 2 at = 0. The second integrates to S constant, which we set to 1 in agreement with the asymptotic limit. The third is readily integrated to give for X a constant of integration. However this is singular at both r = r h and as r → ∞ unless we take X = 0, giving X + (r ) ≡ 0. Examining (161b) and noting that we will want to vary this solution away from = 0 to small non-zero values of , we see that W + (r ) is also a constant, for which we are forced to take W + (r ) ≡ W + (r h ). Hence using an appropriate basis for W + (r ) (29), the unique solution obtained is We note that this is identical to the su(N ) case. Now we take Proposition 17 and re-purpose it to the case at hand. Defining new variablesλ the field equations (143) become and the equation for S is unchanged. But the structure of the field equations is unaltered, and so the proof given in Sect. 7.3 is unchanged. Then, for arbitrarily small , we may find solutions that exist locally in the asymptotic limit.
The argument that proves that non-trivial global solutions exist for small is very similar to Proposition 19. We fix r h , take the existing solution (166), and consider varying {ω j,h }, and varying away from 0. Note that for the embedded solution (166), all gauge fields will be nodeless. We then choose some r 1 r h so that we can consider r 1 in the asymptotic regime. Proposition 16 confirms that for sufficiently small we can find solutions near the existing unique solution which will begin regularly near r = r h and remain regular also at r = r 1 , and that those solutions will have nodeless gauge field functions due to analyticity. Finally, since we are now in the asymptotic regime, we can use the logic in Sect. 5 and Proposition 18 to ensure that solutions will remain regular as r → ∞ and that all ω j will be nodeless.
The corresponding proof for solitons is similar to that for black holes, though we must be more careful about how we take the limit → 0. The parameter τ ∝ r −1 that we use in the asymptotic regime is fine for black holes since min{r } = r h so τ is bounded and thus r −1 remains regular throughout the range [r h , ∞); but this is clearly no longer the case for solitons as min{r } = 0 so that τ becomes singular.
We follow the clues in the su(N ) case [9] and rescale all dimensionful quantities: In addition, we find it best to work with the gauge functionsû j (r ) which we defined in the proof of local existence at the origin, Proposition 15, using and working with the field equations in the form (37a)-(37c). Substituting (169, 170) into the field equations, again we find thatm(x) and S(x) must be constant, which due to boundary conditions we are forced to set equal to 0 and 1 respectively. We also see that if = 0, all gauge functions ω i (x) ≡ ω i,0 , and the solution reduces to the SadS case where ω j ≡ ±λ 1/2 j , which are manifestly nodeless.

123
However it is important to examine the behaviour of the equations for small but non-zero. When = 0, the YM equations (37c) decouple to produce the following: where we have used results (94, 96, 99). Fortunately, though not necessarily unexpectedly, this is also very similar to the su(N ) case [9] (set k j ≡ k in the above) in that the term containing F vanishes in both cases when = 0. Therefore our more general case has a very similar unique solution in this limit: for j = 1, . . . , L, and where the integers k j for the group G in question are given in Table 1. The constant of proportionality above is simply β j from Proposition 15. It can be seen that this is regular at x = 0, and due to the properties of hypergeometric functions, that it satisfies the required boundary conditions (53). We proceed in a very similar fashion to the black hole case. Proposition 17 adapts in a very obvious way, similar to the above (161a, 161b). So we take the existing solution (172) with arbitrary β j , and consider varying {β j } and varying away from 0. Note again that for the embedded solution (166), all gauge fields will be nodeless. We then choose some r 1 >> 0 so that we can consider r 1 in the asymptotic regime. Propositions 15 guarantees that for fixed sufficiently small we can find solutions near the existing unique solution which will begin regularly near r = 0 and remain regular in the range (0, r 1 ], and that those solutions will have nodeless gauge field functions due to analyticity. Finally, once we are in the asymptotic regime, we can again use Proposition 18 and the logic in Sect. 5 to ensure that solutions will remain regular as r → ∞, and that furthermore all these nearby ω j will be nodeless. tions 15 and 16) carried over similarly from the asymptotically flat case [34], with some minor alterations. The biggest difference in the analyses was in the asymptotic behaviour of solutions (Proposition 17). There, we found that the gauge functions and their derivatives were entirely specified by the arbitrary values they approach at infinitythis differs greatly from the = 0 case, in which the gauge field was specified by higher order parameters in the power series, and these parameters were intercoupled in a complicated way. This difference is explained in Sect. 5, where it is noted that due to the parameter we use to render the equations autonomous, the solutions to this system (in terms of dynamical systems) need not reach their critical points, which was what forced the asymptotically flat system to be so tightly constrained as r → ∞.
Due to this difference, it became possible in Sect. 8 to prove the existence of global solutions to the field equations in some neighbourhood of embedded solutions, of which we found three separate cases (Sect. 6). We proved that as long as μ(r ) > 0 throughout the solution range, then if we begin at the initial boundary (r = r h for black holes or r = 0 for solitons) and integrate the field equations out arbitrarily far, the field variables will all remain regular (Proposition 18). We recall that we already established in Sect. 5 that general solutions will remain regular in the asymptotic regime. Therefore, we were able to argue the existence of black hole and soliton solutions which begin regularly at their initial conditions and can be regularly integrated out arbitrarily far, where they will remain regular as r → ∞ (Theorem 19). We finally considered the limit of | | → ∞, which we explained was necessary in the su(N ) case to guarantee nodeless and hence stable solutions, and proved that nodeless non-trivial solutions exist in this regime too, which are similarly globally regular and analytic in their boundary parameters (Theorem 21).
Our main results are the proof of global non-trivial solutions to the field equations (18a)-(18c), both nearby trivial embedded solutions, and in the limit of | | large. It is remarkable to see how many of the general features of this model carry across to the specific case of su(N ) [9]. These include the forms of the field equations themselves, the embedded solutions we find, the qualitative behaviour of the solutions at the various boundaries, and the existence of solutions both near embedded solutions and in the limit | | → ∞. This is very pleasing, since it may be noticed that the field equations (18a)-(18c) may easily be adapted to any gauge group without precise knowledge of the gauge potential itself, the construction of which for a given gauge group is a non-trivial task. This quite general system, even restricted to solely the regular case, could thus prove to be a powerful analytical model which may give insight into a range of different matter field theories.
There are many future directions that this work could take. Considering the work in [46], a logical next step might be to consider the 'irregular' case, where W 0 lies on the boundary of a fundamental Weyl chamber, and the situation is more intricate. For instance, for = 0 it is known that this means the gauge functions ω j will in general be complex. An analysis of that case, in combination with the results here presented, would cover an existence analysis for black holes and solitons in all possible static, spherically symmetric, purely magnetic EYM adS models with a compact semisimple gauge group.
Another obvious thing to do is to consider the question of the stability of the solutions that we have found. In [18], Brodbeck and Straumann give a proof of instability for a general compact gauge group in asymptotically flat space, for the case of the regular action; but here we find that we are able to establish solutions which fulfil the same conditions which guaranteed stability in the case of su(N ). This would be very enlightening to investigate. In addition, there is the issue of extending this work to higher dimensions, though due to the fact that we would now be dealing with essentially SU(3) principal bundle automorphisms for the isometry group of S 3 , and the higher order Cherns-Simons terms in the action needed to obtain finite-mass solutions [47,48], this is likely to be highly technical.
The main impact of this research is on some outstanding questions in gravitational physics. For instance, we consider Bizon's modified "no-hair" theorem in light of this work, which states: Within a given matter theory, a stable black hole is characterised by a finite number of global charges. [49] Since this work concerns a general gauge group, it opens up the interesting possibility of verifying the no-hair theorem for a large class of gauge structure groups, given some further work. In addition, Hawking very recently raised the interesting possibility that hairy black holes may be used to resolve the 'black hole information paradox' [50]. The possibilities that this research opens up for our field are as yet unknown but potentially significant, and it would be of great interest to know if our recent work may be able shed any light on this long-standing problem.
Finally, there is the important question of whether this research will open up new insights into the adS/CFT correspondence. It is known that for black hole models there are observables in the dual CFT which are sensitive to the presence of hair (see [51] for a discussion of non-Abelian solutions in the context of adS/CFT), and correspondences to CMP problems have been found relating to both superconductors [52,53] and superfluids [54]. Therefore, it is possible that within the class of models considered in this paper, there exist many more applications to QFT phenomena, and this could be a rich and worthwhile vein of study.