Some Exact Solutions for Maximally Symmetric Topological Defects in Anti de Sitter Space

We obtain exact analytical solutions for a class of SO($l$) Higgs field theories in a non-dynamic background $n$-dimensional anti de Sitter space. These finite transverse energy solutions are maximally symmetric $p$-dimensional topological defects where $n=(p+1)+l$. The radius of curvature of anti de Sitter space provides an extra length scale that allows us to study the equations of motion in a limit where the masses of the Higgs field and the massive vector bosons are both vanishing. We call this the double BPS limit. In anti de Sitter space, the equations of motion depend on both $p$ and $l$. The exact analytical solutions are expressed in terms of standard special functions. The known exact analytical solutions are for kink-like defects ($p=0,1,2,\dotsc;\, l=1$), vortex-like defects ($p=1,2,3;\, l=2$), and the 'tHooft-Polyakov monopole ($p=0;\, l=3$). A bonus is that the double BPS limit automatically gives a maximally symmetric classical glueball type solution. In certain cases where we did not find an analytic solution, we present numerical solutions to the equations of motion. The asymptotically exponentially increasing volume with distance of anti de Sitter space imposes different constraints than those found in the study of defects in Minkowski space.

1 Introduction In this article we obtain exact analytical solutions for a class of SO(l) Higgs field theories in a non-dynamic background n-dimensional anti de Sitter space AdS n . These field theories admit maximally symmetric p-dimensional topological defects. The world brane of a maximally symmetric p-defect is a q = p + 1 dimensional timelike submanifold Σ q ≈ AdS q that is isometrically embedded in AdS n ; it is the gauge invariant set corresponding to the zero locus of the Higgs field. The value of l is determined by n = q + l.
The search for a maximally symmetric defect solution to the equations of motion requires the Lorentzian submanifold Σ q to admit the largest possible group of isometries. For q-dimensional manifolds, this Lie group has dimension 1 2 q(q + 1). The choices of this q-dimensional manifold are Minkowski space M q , anti de Sitter space AdS q , and de Sitter space dS q . In this article we mostly discuss the anti de Sitter cases. We show that the Minkowski and de Sitter cases do not give a maximally symmetric solution. The anti de Sitter case gives a maximally symmetric solution when the isometric embedding AdS q → AdS n is totally geodesic. In our defect considerations we assume that when we refer any of these maximally symmetric manifolds, we are implicitly considering the simply connected universal covering space.
We need that p ≥ 0 or equivalently that q ≥ 1. The mathematical reason is that the formalism we employ requires an Euclidean signature for the metric of the normal tangent space (T σ Σ) ⊥ for σ ∈ Σ. If we allow q = 0 then AdS 0 is a point and its normal tangent space is the whole tangent space of AdS n at that point, which has Minkowski signature, and our formulas do not apply directly 1 . We are restricted to AdS n with dimensionality n = q + l ≥ 2.
We classify AdS topological defects by their degrees of longitudinal and transverse freedom using the tuple notation (q, l). We develop a universal formalism that allows us to study all values of (q, l). We only need to study in detail three types of defects: kink defects, vortices, and monopoles. Our parlance is that a defect with l = 1 (with one transverse dimension) is kink-like, a defect with l = 2 is vortex-like, and a defect with l = 3 is hedgehog-like or monopole-like. For a review of Minkowski space kinks (l = 1) look in [1,Chapter 6]. The Nielsen-Olesen vortex [2] has l = 2, and p = 1 or equivalently q = 2. The 'tHooft-Polyakov monopole [3,4] has l = 3, and p = 0 or equivalently q = 1. For l ≥ 4, the transverse energy of these spherically symmetric solutions in AdS n diverges. For a comprehensive review of topological solitons we recommend the book by Manton and Sutcliffe [5]. We do not discuss electrically charged defects such as the Julia-Zee dyon [6].
There is a major difference between the study of the equations of motion for maximally symmetric p-defects in Minkowski space M n and in anti de Sitter space AdS n . In Minkowski space, the equations of motion only depend on the transverse dimensionality l, and the study of the solutions is independent of the dimensionality p of the defect. This is not the case for p-defects in anti de Sitter space where the equations of motion depend on both p and l. This forces us to do a case by case analysis as we vary p and l, e.g., see figure 4.
Line vortices in AdS 4 were discussed in reference [7], and line vortices in dS 4 in [8]. In the first reference, approximate analytic solutions for vortices in AdS 4 were found after applying some simplifying approximations to the equations of motion. Vortex holography played a strong part of the discussions in these papers.
The study of magnetic monopoles in AdS 4 , the case with (q, l) = (1, 3), has been around for a while. The earliest work we are familiar with are papers authored by Lugo and Schaposnik [9], and Lugo, Moreno and Schaposnik [10]. We collectively refer to these two articles as LMS. In appendix B we discuss how their work is related to ours. Numerical axially symmetric monopole solutions in AdS 4 are explored in [11]. No exact analytic solutions were found in these references. Approximate analytic and numerical methods are used to discuss multi-monopoles and multi-monopole walls and their importance in the AdS/CFT correspondence in [12] and in [13]. Atiyah [14] had earlier discussed magnetic monopoles in Euclidean hyperbolic 3-space H 3 by exploiting the conformal invariance of the self-dual Yang-Mills equations. His 4-dimensional manifold was S 1 × H 3 with the Euclidean signature product metric. This manifold is conformal to Euclidean space E 4 . He used the observation due to Bogomolny that the self-dual Yang-Mills equations applied to a time independent SU(2) gauge field are equivalent to the Bogomolny equations 2 . The product manifold Atiyah uses is not the Lorentzian manifold AdS 4 , and his monopoles are not AdS 4 monopoles.
In the 'tHooft-Polyakov monopole in Minkowski space M 4 there are two independent length scales determined by the mass of the Higgs scalar m φ and by the mass of the massive vector meson m A . A good way to see this is to observe that after an appropriate rescaling of the fields, the action for the SO(3) Georgi-Glashow model may be schematically written as here φ 0 is the vacuum expectation value of the scalar field. For example, you expect solutions of the classical equations of motion for the rescaled dimensionless scalar field ϕ to depend on m A and on the dimensionless ratio m φ /m A : ϕ = ϕ(m A x, m φ /m A ). Prasad and Sommerfield [15] discovered an exact solution to the equations of motion by considering a non-trivial limit of the equations of motion in which m φ ↓ 0, keeping m A fixed, and enforcing an appropriate topological asymptotic boundary condition of the Higgs field ϕ BPS line  at infinity 3 . The net effect is that the equations of motion only depend on one length scale 1/m A that controls the asymptotic behavior along with the correct boundary conditions imposed manually. The dimensionless parameter m φ /m A ↓ 0 in the Prasad and Sommerfield limit. The Prasad-Sommerfield solution satisfies the first order equations of Bogomolny [16] that guarantee a solution with a saturated lower bound on the energy [16,17]. It is known that the mass of the monopole is given by where the function F satisfies F(0) = 1 and F(∞) ≈ 1.787, see [5, p. 255].
In this paper we also discuss the Prasad and Sommerfield limit m φ ↓ 0 for other values of q and l. We refer to this as the limit of Bogomolny, and Prasad and Sommerfield (BPS) even though there may be no Bogomolny equations. The BPS limit was studied in LMS for the case of monopoles in AdS 4 .
In studying the equations of motion for defects in AdS n , we encounter an additional length scale ρ, the radius of curvature 4 of AdS n . Now there are three independent length scales 1/m φ , 1/m A , and ρ. This leads to a three dimensional parameter space, see figure 1, that can be explored for exact solutions of the equations of motion. The appearance of this additional length scale was already noticed in [9]. The scalar field solution ϕ of the equations of motion depends on ρ and on two dimensionless quantities m φ ρ and m A ρ: ϕ = ϕ(x/ρ, m φ ρ, m A ρ). We attempt to extend the methods of Bogomolny, and Prasad and Sommerfield, and look for regions in the parameter space where we might find exact solutions. We consider a limit for the equations of motion for a maximally symmetric defect where ρ = 0 is fixed, but m φ ρ ↓ 0 and m A ρ ↓ 0. The net effect is that the equations of motion only depend on one length scale ρ that controls the asymptotic behavior along with the correct boundary conditions imposed manually.
N.B. This is a very delicate limit because the action is singular in this limit but the equations of motion are not.
We take the limiting equations of motion as the starting point in our analysis, and we abandon the action. These limiting equations are consistent and do not follow from an action principle. This is analogous to consistent equations of motions such as the selfdual Yang-Mills equations, or the self-dual equations of motion for the 4-form in type IIB supergravity which are not derivable from an action. The limiting equations of motion partially decouple. The one for the gauge field is completely decoupled from the Higgs field and can be solved independently. The gauge field solution can then be inserted into the Higgs field equation of motion which is now linear. In addition, this "double BPS limit" preserves the nonlinear interactions of the gauge fields. The vanishing of the dimensionless parameters, m φ ρ and m A ρ, leads to exact analytic solutions in a variety of cases, see figure 5 for admissible pairs (q, l). For example, there are exact analytic solutions for kink defects (q, l) = (q, 1) given by (5.8), for the Nielsen-Olesen vortex line (q, l) = (2, 2) given by (5.12), and the 'tHooft-Polyakov monopole (q, l) = (1, 3) given by (5.19) and (5.22). We note that the transverse size of these p-defects in the double BPS limit is comparable to the radius of curvature of AdS n .
These double BPS equations of motion are a first step in a perturbative expansion of the full equations of motion where the small parameters are m φ ρ and m A ρ. In this sense, we can make contact with the action again.
We have not studied the stability of these double BPS solutions in anti de Sitter space. Our only attempt at trying to prove stability was to look for a Bogomolny type bound. This bound is used to establish the stability of the Prasad-Sommerfield solution [16,17] in Minkowski space. The curvature of AdS n invalidates some of the Minkowski space arguments as explained in appendix C. We show that in the monopole case there are no first order Bogomolny type equations that imply the equations of motion. Additionally, there is a partial Bogomolny bound relating the energy to the magnetic charge when the magnetic charge density is non-negative. In principle, the underlying symmetries of AdS q and AdS n should greatly aid in the stability analysis of linear perturbations of the equations of motion around these double BPS solutions.
The organization of this article is the following: A formalism is developed to study maximally symmetric p-defects in AdS n in sections 2 and 3, leading to the general maximally symmetric equations of motion (3.15) for all values of (q, l). One useful feature of our formalism is that the normal radial coordinate ν we employ is always the physical distance from the defect world brane no matter what (q, l) pair we are studying. We found it convenient to use an orthonormal frame because the Pythagorean theorem automatically organized the calculation for us. For example, the transverse energy functional consists of five positive semi-definite summands, therefore finiteness of transverse energy follows if each summand is finite. The theorem that requires Σ q to be totally geodesic is relegated to mathematical appendix A. The finite transverse energy constraints are discussed in section 4, and summarized in figure 4. In section 5 we discuss the double BPS limit and the exact solutions we found. In section 6, we briefly explore numerically a portion of the parameter space m A ↓ 0 where the gauge field is explicitly known but the scalar field has to be studied numerically. In appendix B we relate our work to the work of Lugo, Moreno and Schaposnik. Finally in appendix C we restrict to magnetic monopoles in AdS 4 in the BPS limit. We discuss that the BPS equations do not imply the equations of motion. We show the existence of a partial bound on the mass if the magnetic charge density is non-negative.

Defects in constant curvature spaces
We use a SO(l) gauged Higgs type field theory as the model for a topological defect. We first discuss the definition of a maximally symmetric p-dimensional defect in Minkowski space and subsequently generalize the notion to a Lorentzian constant curvature space. A p-dimensional maximally symmetric defect in M n is a topologically stable solution to the equations of motion that is invariant with respect to the action of the subgroup P(q) × SO(l) ⊂ P(n) where q = p + 1 and q + l = n. Here P(n) is the Poincaré group, the isometry group of Minkowski space M n , and M n ≈ P(n)/ SO(1, n − 1). The world brane (time evolution of the defect brane) for the defect is the q-dimensional manifold Σ q . Since the symmetry group of the solution is P(q) × SO(l) we know that the world brane Σ q for the defect is a timelike q-plane. Notice that Σ q is intrinsically flat, and that the invariance group of the solution implies that the defect is static for any choice of time direction in Σ q .
Let M n be a Lorentzian manifold with constant sectional curvature k and with isometry group 5 Isom 0 (M n ) of dimension 1 2 n(n + 1). A p-dimensional maximally symmetric defect in M n is a topologically stable solution to the equations of motion that is invariant with respect to the action of the subgroup Isom 0 (Σ q ) × SO(l) ⊂ Isom 0 (M n ) where Σ q is a maximally symmetric q-dimensional Lorentzian submanifold. A maximally symmetric Σ q is a constant curvature manifold with dim Isom 0 (Σ q ) = 1 2 q(q + 1). Here Σ q is the world brane of the defect core 6 . Finding such a Σ q ⊂ M n requires an appropriate generalization of choosing a plane. The correct "flatness" notion that leads to a maximally symmetric defect is to require Σ q to be a totally geodesic submanifold.
An embedded submanifold Σ q of a general manifold M n is said to be totally geodesic if every geodesic (with respect to the induced metric) on Σ q is also a geodesic on M n . Next, we see how this is related to the differential geometric data. Let D Σ and D M be the Levi-Civita connections on the respective manifolds. A Darboux frame is an orthonormal frame adapted to the orthogonal decomposition T σ M = T σ Σ + (T σ Σ) ⊥ for σ ∈ Σ q . In a submanifold neighborhood of the point σ, we consider an orthonormal framing (ê a ,n i ), where theê a are tangential to Σ and then i are normal to Σ. We use the index convention that latin indices from the beginning of the alphabet a, b, c, d run from 1, 2, . . . , q and latin indices from the middle of the alphabet i, j, k, . . . take l = n − q values from q + 1, . . . , n.
Let u = u aê a and v = v bê b be vector fields tangent to Σ q , then the two connections are related by The symmetric tensor K ab i is called the second fundamental form or the extrinsic curvatures. If u is a tangent vector on Σ then In a totally geodesic submanifold, we would have that D M u u = 0 and D Σ u u = 0 for all geodesics on Σ. This is only possible if the extrinsic curvatures K ab i = 0. Summarizing, totally geodesic submanifolds are those where the extrinsic curvatures vanish. From the viewpoint of standard General Relativity, totally geodesic submanifolds are very desirable because if a test mass in Σ is given an initial velocity tangential to Σ then its motion will remain in Σ.
As shown in detail in appendix A, a totally geodesic q-dimensional submanifold Σ q of AdS n is a constant curvature Lorentzian submanifold with the same sectional curvature k as M n and with a flat normal bundle (T Σ) ⊥ . These are the only submanifolds that admit the possibility of finding a maximally symmetric solution, see remark 1 in the appendix A.3.

The Darboux frame and the spherically symmetric ansatz
To work out the equations of motion for our maximally symmetric defect, it is convenient to use a coordinate system adapted to the geometry of the problem, i.e., an analog of spherical coordinates. The construction is based on the method of Cartan discussed in our previous paper [18]. Let Σ q be a Lorentzian q-dimensional submanifold of AdS n . If σ ∈ Σ q then in a submanifold neighborhood of the point σ consider a Darboux frame (ê a ,n i ). The dual Darboux coframe is denoted by (ϕ a , ϕ i ). Choose a geodesic of AdS n starting at σ with initial normal velocity ν = ν in i , and go a distance ν along the geodesic to a point x ∈ AdS n . The coordinates of the point x are (σ, ν i ). Formally, this is the exponential map exp σ : ν ∈ (T σ Σ) ⊥ → x ∈ AdS n . Note that ϕ i = dν i . Cartan's idea is to extend the orthonormal Darboux coframe by parallel transporting it along the normal geodesics; in this way you construct an orthonormal coframe at (σ, ν) denoted by (ϑ a , ϑ i ). In our previous paper we showed that The velocity of a geodesic is constant, and in this way we denote the orthogonal projector along the velocity vector by (P L ) ij = ν i ν j / ν 2 , and the orthogonal projector perpendicular to the velocity by (P T ) ij = δ ij − ν i ν j / ν 2 . Using these, it is easy to write down the orthogonal decomposition of the extended Darboux frame component In the above (Dν) i = dν i +ω ij ν j , where ω ij is the connection on the normal bundle (T Σ) ⊥ . The formulas above are general. We are interested in the Lorentzian case where AdS q ≈ Σ q → M n is a totally geodesic submanifold. As discussed previously, we know that K abi = 0 and that the normal bundle (T Σ) ⊥ is flat, see appendix A. We can always locally trivialize the normal bundle so we set ω ij = 0, i.e., the normal part of the Darboux framen i is parallel along Σ q . In summary we have The metric on AdS n is given by ds 2 Going to spherical coordinates we see that where ds 2 AdSq = g M ab ϕ a ⊗ ϕ b is the standard constant curvature metric on Σ q ≈ AdS q , ν = ν is the radial distance, and ds 2 S l−1 is the round metric on the unit (l − 1)-sphere. The part of the metric in the square brackets is the pullback of the metric on AdS n to (T σ Σ) ⊥ via the exponential map exp σ . Said differently, this is the induced metric on exp σ (T σ Σ) ⊥ , the image of the normal tangent space under the exponential map. This induced metric is isometric to the standard metric on Euclidean hyperbolic space H l , see the discussion associated with eq. (3.5). We re-emphasize that ν = ν is the physical distance from a point σ ∈ Σ q to the point (σ, ν) ∈ M n . We will use this physical distance to measure the behavior of our fields as you move away from the defect world brane.
A simple model that has a topological p-defect in AdS n is a Higgs model with an SO(l) gauge symmetry, l ≥ 2, where n = (p + 1) + l = q + l. The model has a scalar field Φ I that transforms under the vector representation of SO(l). The uppercase latin indices I, J, K, . . . from the middle of the alphabet will take values from 1 to l. We are looking for a maximally symmetric p-defect that is invariant under the action of Isom 0 (AdS q ) × SO(l), and therefore our fields do not depend on the coordinates σ a of AdS q and only depend on the normal coordinates ν i . A connection that is compatible with the symmetries is The equations of motion are obtained by extremizing the action I = M n L, where L is the Higgs model Lagrangian density. We are looking for p-defect solutions that are maximally symmetric with symmetry group Isom 0 (Σ) × SO(l) and under these conditions the action for the Higgs type model is where ζ Σ is the volume element on Σ. The transverse energy 7 E ⊥ in a local orthonormal 7 E ⊥ is the tension of the associated p-brane. frame for (T σ M ) ⊥ is given by The volume element formula is a special case of a result from our previous paper [18]. We implicitly assumed spherical symmetry for the Lagrangian density to do the angular integrals. Here U is the potential, g is the gauge coupling constant, and V l−1 is the (l − 1)volume of the unit sphere S l−1 . We emphasize to the reader the presence of the hyperbolic cosine factor that would not be there in the case of M n = M n . The origin of this hyperbolic cosine factor is the metric (3.2). You can make a mistake influenced by the familiarity of working in Minkowski space where the transverse energy would be obtained by using the pullback metric to (T σ Σ) ⊥ via the exponential map. In this case for p-defects in AdS n , the metric on the normal tangent space is the metric on Euclidean hyperbolic space H l Said differently, removing the hyperbolic cosine term in (3.4) leads to incorrect equations of motion. The equations of motion arise from varying the action (3.3). Note that the k → 0 limit of (3.4) is the familiar transverse energy for a spherically symmetric p-defect in Minkowski space M n . It is convenient to define the Jacobian factor where the "radius of curvature" of AdS n is ρ = 1/|k| 1/2 .
The SO(l) spherically symmetric ansatz we employ is a generalization and slight variant of the original one used by 'tHooft [3] and by Polyakov [4] in the SO(3) Georgi-Glashow model. The ansatz is where φ and f are functions only of the radius ν = ν . We choose the potential U to be of the general symmetry breaking form such as the one shown in figure 2. The form chosen for Φ is a hedgehog type ansatz with φ(ν) ν→∞ − −−− → ±φ 0 corresponding to the topological winding number ±1 solutions. We want the world brane of the p-defect to be a gauge invariant set thus we require Φ(ν = 0) = 0 which implies φ(0) = 0. The form for A is motivated by the abelian constant curvature spherically symmetric ansatz 8 and by the identification of the SO(l) Lie algebra associated with the gauge group to the one associated with the rotational isometries of (T σ Σ) ⊥ . In this article we only discuss situations where both the functions φ(ν) and f (ν) are non-trivial functions of ν that describe localized defects. We note that expressions (3.4) and (3.7) are actually valid for l = 1 if you interpret V 0 = 2 since S 0 ⊂ E 1 consists of two points {−1, +1}. Ansatz (3.7) for l = 1 is the odd parity kink (domain wall) solution with an automatically vanishing gauge field. A brief computation and implementing (3.1) gives To compute DΦ 2 we exploit that in the orthonormal polar coframe, the longitudinal direction is the radial direction, and that P T and P L are orthogonal projectors to conclude (3.9) Next we compute the curvature by using differential forms F IJ = dA IJ + A IK ∧ A KJ . The computation is a bit more involved but greatly simplifies by using the orthogonal projectors: The expected behavior near the world brane is constant field strength so we expect f (0) = 0. Figure 3: The expected behavior of φ and f as a function of radial distance ν for a topological defect. The region with the gradient shade in yellow represents the core of the defect.
Using the properties of the orthonormal polar coframe we conclude (3.11) The norm on 2-forms is normalized by observing that in an orthonormal coordinate system in E 2 , a constant U(1) field strength is given by F = f dν 1 ∧ dν 2 and F 2 = f 2 . When we write F IJ 2 we mean sum over all I and J. For example in an SO(2) gauge theory, you have F IJ 2 = F 12 2 + F 21 2 = 2 F 12 2 . In general when summing over the spacetime indices i, j you obtain F IJ ij F IJ ij = 2 F IJ 2 . Note that the Yang-Mills lagrangian is quadratic in f in the abelian SO(2) case.
Collating all the terms, we have that the transverse energy E ⊥ of the defect is (3.12) This expression for E ⊥ is valid for l ≥ 1. Equation (3.12) may be simplified by introducing the auxiliary function h defined by The expected defect boundary conditions on h are h(ν) ν→0 − −− → +1 and h(ν) ν→∞ − −−− → 0, see eq. (4.2). In terms of φ and h, the transverse energy E ⊥ looks like: (3.14) With the defect boundary conditions our ansatz does not admit pure gauge solutions. To see this we note that if l ≥ 2 then setting F IJ = 0 in (3.11) leads to f (ν) = 0 or f (ν) = 4/ν 2 , algebraic forms for f (ν) that are incompatible with the defect boundary conditions.
If we restrict to the case of l = 2, then in polar coordinates (ν, ϕ) for the normal bundle we would have that A 12 = f (ν) · 1 2 ν 2 dϕ ν→∞ − −−− → dϕ and thus we conclude that we have nontrivial holonomy because if we integrate along the circle at infinity we have A 12 = 2π. This tells us that the total vortex flux is 2π in our normalization, that may be related to the conventional normalization by A = gA conv .
In the above discussion, the action leads to the standard Laplacian in the equations of motion. The situation is the same for the theory with the conformal Laplacian, which is obtained by adding to I spherically sym a term of the type cR M n Φ 2 where c is some fixed constant. Since the scalar curvature R M n is a constant, that term can be absorbed into the potential energy function U as a correction to the quadratic term in Φ. We still choose the potential U to be of the general symmetry breaking form such as the one shown in figure 2.
Minimizing E ⊥ , we can derive the equations of motion for the defect. The equation of motion for φ is Similarly, the h equation of motion is (3.15b)

Finite transverse energy constraints
We would like the p-defects in AdS n to have finite transverse energy (3.14). To achieve this, we have to study the convergence of the integral for E ⊥ in the two asymptotic limits ν → ∞ and ν → 0 where we require that 0 < ρ < ∞. Of these two limits, the ν → ∞ one is more important. If the transverse energy diverges in the ν → ∞ limit then there is no hope for a defect. A divergence in the integral as ν → 0 may be resolved by considering a different ultraviolet completion of the model. For example, you could add higher derivative terms to the action analogous to what is done in the 4D Skyrme model [19] to stabilize the model. In this section, we consider both limits within the context of expression (3.14) for E ⊥ .

Behavior as ν → ∞
The ν → ∞ behavior is tricky, and has to be carefully analyzed on a case by case basis: l = 1, l = 2, and l ≥ 3. We are interested in solutions that satisfy the boundary conditions The E ⊥ integrand (3.14) consists of five positive semi-definite summands. A finite E ⊥ solution requires that the integral of each summand converge. These conditions impose growth rates on the fields h and φ. We observe that the growth rates that lead to convergent E ⊥ are not necessarily the growth rates given by the equations of motion. We are interested in topological defect solutions with finite transverse energy so the boundary condition that we choose is φ(ν) ν→+∞ − −−−− → +φ 0 for the winding number +1 defect. The form of the potential, see figure 2, tells us that as ν → +∞ we are near a quadratic minimum and we have that We would like for the length scale given by the radius of curvature ρ to dominate the Compton wavelength 1/m φ of the scalar field, so in our later applications we require ρ 1/m φ even though this does not enter into the convergence analysis. Note that the case ρ 1/m φ essentially reduces to the flat space case. For all l ≥ 1, we note that requiring that the scalar field decays exponentially as ν → ∞ for some > 0 guarantees the convergence of the φ (ν) 2 term and the U (φ) term 9 in (3.14). For the future we note that for topological defects in AdS n , the asymptotic behavior of φ and h depend not only on the transverse dimensionality l but also on the dimension p = q − 1 of the defect, see (4.6) and (4.10). This is different than the situation of the familiar defects in M n .
To have a gauge theory, we need l ≥ 2 which implies that n = q + l ≥ 3. In the same spirit, we analyze the asymptotic behavior of summands two and four in (3.14) that are associated with the kinetic energy of the gauge field h. Those integrals will converge if Finiteness of E ⊥ imposes strong constraints if l ≥ 3 because the non-abelian gauge fields' self interactions contribute to the energy via the term In fact the integral will converge with the weaker condition ψ( in eq. (3.14). We know that h(ν) ν→∞ − −−− → 0 and the convergence of the integral may be determined from the asymptotics of the hyperbolic functions ∞ dν e (n−5)ν/ρ . (4.4) This integral will converge if n < 5. Because l ≥ 3 and q ≥ 1, we see that there is exactly one case where the integral converges, namely q = 1 with l = 3. This is the soliton case, a 0-defect, which is the 'tHooft-Polyakov monopole in AdS 4 . The n = 5 case where the integral is linearly divergent may be salvageable via some unknown method (see the flat space discussion below). There you would have an SO(4) soliton (p = 0), and an SO(3) line defect (p = 1).
The situation for the existence SO(l) nonabelian topological defects for n ≥ 6 appears to be quite dire because of the exponentially divergent energy.
The results are different in the flat space case, i.e., ρ = ∞. Here, the asymptotics of (4.3) are given by This integral converges if l = 3 for all values of q. Thus SO(3) topological p-defects are energetically allowed for all values of p ≥ 0, note that n ≥ 4. These are 'tHooft-Polyakov p-defects that have finite energy. The l = 4 case has logarithmically divergent energy and may be salvageable by modifying the model in some way. An analogy is the XY-model 10 where the vortex has a logarithmically divergent energy that can be tamed by adding an abelian gauge field. Another approach to finite energy in the XY model is to just have the scalar field but restrict the field configurations to the topological sector with zero net winding number. The consequence of this is that the field decays faster at infinity and leads to finite energy. In this way you can have finite energy vortex anti-vortex pairs. The vortices have an ultraviolet energy divergence near the world brane. This is usually remedied by imposing a finite radius cutoff in the core, or including higher derivative terms such as in the Skyrme model.

Asymptotic behavior of φ
If we look at (3.15a) we note that as ν → ∞ the term containing h 2 is very small compared to the other terms because h is exponentially small and 1/ sinh 2 (ν/ρ) is also exponentially small 11 . In this case, the asymptotic behavior of the equation of motion (3.15a) is The XY model is the SO(2) nonlinear sigma model in E 2 , i.e., n = 2 flat space. 11 Note that this term is absent in the l = 1 case. and ψ * > 0 is a constant. This decay behavior leads to convergence in the relevant transverse energy integral summands if m 2 φ + [(n − 1)/2ρ] 2 > 0 by (4.1). This condition is always satisfied in our models because n ≥ 2. The contribution to the transverse energy density in the asymptotic region from the purely scalar field part is Thus, we have that The important result here is that the pure φ part contribution to the transverse energy integral converges for all values of n ≥ 2. The topological kink (l = 1) exists for all q ≥ 1.

Asymptotic equation of motion for h
As ν → ∞ with our boundary conditions, the asymptotic equation of motion is To have a gauge field, we need l ≥ 2 which implies that n ≥ 3. The relevant transverse energy terms will converge if [(n − 3)/2ρ] 2 + m 2 A > 0 according to criterion (4.2). To compute the contribution of h to the energy density, we consider the second and fourth terms of (3.14): The asymptotic energy density for the gauge field is 12) and the associated integrand is The asymptotic behavior of h in the case l = 2 The l = 2 case is the abelian vortex and we see that in (3.14) that the non-abelian selfinteraction term vanishes identically. The asymptotic energy density integrand for the gauge field is Notice that the h part contribution to the transverse energy integral converges for all values of q ≥ 2. Therefore, the topological abelian vortex localized at a totally geodesic AdS q → AdS q+2 is energetically possible, i.e., vortex line, vortex sheet, etc.. For q = 1, the vortex soliton requires m A > 0.

The asymptotic behavior of
Expression (4.13) is the contribution from the second and fourth summands of (3.14). We note that since n ≥ 4 this contribution to the transverse energy converges even if m A = 0. The only concern is the fifth summand and we have already addressed it in section 4.1 in the derivation of (4.4). The conclusion is that we do not expect solutions if n ≥ 5 for l ≥ 3. Figure 4 is a summary of values of (q, l) that admit or do not admit solutions.

Behavior as ν → 0
In studying the ν → 0 behavior of φ and h, we will encounter two second order linear ODEs with regular singular points at ν = 0. We look for Frobenius solutions of the form ν α 1 + O(ν 2 ) . There are two real Frobenius indices, α + and α − , with α + > α − . In both the φ and h cases, α + − α − = l ∈ Z. From the theorem of Fuchs, we know that there is a solution of the form ν α + 1 + O(ν 2 ) . The solution involving α − may also have a logarithm. This solution does not satisfy the desired boundary conditions or the finite transverse energy constraint. From now on we only consider the α + solution.
The boundary conditions as ν → 0 are φ(ν) → 0 and h(ν) → +1. First, we examine the ν → 0 behavior of (3.15a). For ν ρ, the equation of motion is approximately First, we discuss the l = 1 case where we immediately see that the solution behaves like φ(ν) ∼ a + bν. We are interested in setting a = 0 because that solution vanishes at ν = 0 and has an odd extension to ν < 0 corresponding to a kink localized at the origin. Next, we look at l ≥ 2, where a Frobenius type solution of the form φ(ν) = Cν α 1 + O(ν 2 ) leads to α + = 1 or α − = 1 − l. Thus we conclude that the small ν behavior of φ is Next we look at the ν → 0 behavior of (3.15b). From (3.14) we see that we have to assume that l ≥ 2 to discuss the h equation of motion. We defineh(ν) = h(ν) − 1, then the boundary condition becomesh(0) = 0. For ν ρ and |h(ν)| 1, the equation of motion is approximately To summarize, the asymptotic behaviors of the fields for this kind of topological defect is given by Near region: Far region: In deriving the asymptotic behavior we explicitly assumed that the radius of curvature satisfied 0 < ρ < ∞. Taking a limit such as (1/ρ) ↓ 0 may be delicate and we have to be very careful. There are no issues with the flat space limit (1/ρ) ↓ 0 if m φ > 0 and m A > 0. On the other hand, you cannot directly apply the asymptotic formulas above when you approach the BPS line in figure 1. This is a delicate limit in which you have (1/ρ) ↓ 0 and m φ ↓ 0. The reason is that the assumptions leading to eq. (4.19) are now invalid. The solution (7.2) of the flat space BPS equations lead to a Coulombic tail asymptotic behavior φ(ν) → φ 0 + O(1/ν) because the flat space φ field is massless. You have to be equally careful when you approach the NS-BPS line. The behavior along the double BPS line is safe because the radius of curvature remains non-zero and governs the asymptotic behavior.
5 Double well potential model with the double BPS limit In this section, we will explore a method of obtaining exact solutions to the equations of motion for some of these defects. We study the equations of motion for a model with potential function It is convenient to rescale to dimensionless variables via ν → ρν, φ → φ 0 φ. In flat space, the mass of the Higgs boson is m 2 φ = λφ 2 0 , and (l − 1) vector bosons acquire mass via the Higgs mechanism with value m 2 A = g 2 φ 2 0 . The transverse energy (3.14) may be written as We note that [ρ] = M −1 = L, [φ 0 ] = M (n−2)/2 , therefore the prefactor of the integral may we written as ρ n−2 φ 2 0 /ρ q , and we see that [E ⊥ ] = M q = M/L p as required. The equations of motion are Note that the equations of motion for a p-defect in AdS n depend on both q and l. The equations of motion (7.1) for a p-defect in M n depend only on l.
In the Minkowski space M n version of the Higgs model, there are two independent length scales 1/m φ and 1/m A that enter the equations of motion. We are studying topological defects in AdS n and there is automatically an independent third length scale ρ, the radius of curvature of AdS n . The spherically symmetric equations of motion (5.3) depend on three mass scales m 2 φ , m 2 A , and 1/ρ 2 if l ≥ 2. The radius of curvature ρ is implicit in (5.3) because the ν coordinate that appears in those equations of motion is actually the dimensionless radial distance ν/ρ. From the form of the equations of motion there are various parameter limits that can be studied, see figure 1. The limit 1/ρ 2 ↓ 0 corresponds to the Minkowski flat space equations of motion. This is the yellow planar region in the figure. At a boundary of this region is the BPS line corresponding to solutions where m 2 φ ↓ 0 but m 2 A = 0. There is also a very delicate non-standard BPS (NS-BPS) like limit in flat space where you consider m 2 A ↓ 0 with m 2 φ = 0; there are no solutions with our ansatz because because the h-field does not satisfy the two point boundary conditions at ν = 0 and ν = ∞. This NS-BPS limit is denoted by the dashed green line. LMS studied the gray planar region given by m 2 φ ↓ 0, and they were interested in the the BPS line boundary but did not consider the double BPS line, see appendix B. We will find explicit analytic solutions with 1/ρ 2 = 0 but m 2 φ ↓ 0 and m 2 A ↓ 0 along what we call the double BPS limit. Finally in the limit m 2 A ↓ 0 we have the pink planar region where we have an analytic form for the h-field, and a numerical solution for the φ-field which are briefly discussed in section 6.
At this point, it is worth the effort to be more explicit in the limit we are taking, as it will appear in much of what is to follow. As previously mentioned, this limit is when m 2 φ ↓ 0 and m 2 A ↓ 0, and can only be taken when there is a third length scale for the physics. In our case, we have one given by the radius of curvature of AdS n . It is also important to understand the way that this limit is taken. Up until this point, we have been using the action of a topological defect embedded in AdS n , this is what leads to the transverse energy integral (5.2) and the equations of motion (5.3). What we do is take the double BPS limit in the equations of motion, allowing the masses to fall towards zero and partially decoupling the ODEs, see (5.4) below. As remarked in the Introduction, the action is singular in this limit but the equations of motion are not. The double BPS equations of motion do not follow from any action. The solutions that we obtain from this method are a good point to start in a perturbative analysis of the full equations of motion that do follow from the action.
In this double BPS limit, the equations of motion become Note that the h equation of motion (5.4b) has decoupled, yet it is nonlinear. Thus we have a standalone second order ODE for h, whose solution can be inserted into eq. (5.4a) to obtain a standalone linear second order ODE for φ. If you are interested in pure Yang-Mills theory then you can just study equation (5.4b) for h and ignore the φ equation of motion.
On the other hand, we see that if a gauge field is present then we get exponential decay in (5.6) if and only if n > 3. The asymptotic analysis for the case of n = 3 with q = 1 and l = 2 is delicate because we lose the exponential decay factor in (5.6) but the expression has a factor of (n − 3) 2 . To assess this case correctly requires more detailed analysis, see section 5.2.1, and we will conclude that there is no acceptable solution with q = 1 and l = 2. This means that in the double BPS limit we lose the vortex solution with q = 1 and l = 2 which is a 0-defect. This discussion is summarized in figure 5.
In this article we consider the first step in a perturbative analysis of eqs.

The double BPS kink-like defects (l = 1)
In the l = 1 case we only have one equation of motion (5.4a) that is easily solved with the required boundary conditions. Since q ≥ 1, the solution to the equation of motion is given by The integrals above may be performed exactly for integer q ≥ 1 and the answer is given in terms of the hypergeometric function:

The double BPS vortex-like defects (l = 2)
The abelian vortex l = 2 in the double BPS limit has various simplifications. First we remark that for a finite energy solution eq. (5.6) requires that the p-defect must have p ≥ 1 or equivalently q ≥ 2. The marginal case with q = 1 is discussed in the next section. The h equation of motion (5.4b) is decoupled and linear since the abelian gauge field has no self-interactions: The solution to this equation with the two point boundary conditions is From this we see that f (0) = q − 1. Inserting this solution for h in the φ equation of motion (5.4a) leads to the linear ordinary differential equation The solution to this equation must satisfy the two sided boundary conditions φ(0) = 0 and φ(∞) = 1. We point out that this ODE has a regular singular point at ν = 0. The behavior near the origin may be studied using the Frobenius power series method which leads to an indicial equation with indices ±1. We consider the regular solution that begins like ν and ignore the singular solution.

Case of q = 1 (no acceptable solution)
The analysis of q = 1 requires the observation that the general solution to (5.4b) is h(ν) = c 1 · 1 + c 2 ln cosh ν. Note that h(ν) = O(ν) as ν → ∞. This function cannot satisfy the two point boundary conditions on the h field. We can require h to have the correct behavior near ν = 0, and we write h(ν) = 1−f 0 ln cosh ν where f 0 is a constant. The asymptotic behavior is given by h(ν) ∼ −f 0 ν as ν → ∞, and this is problematic because the contribution to E ⊥ in (3.14) from second summand diverges cubically while the contribution from the fourth summand diverges linearly. Thus there is no acceptable vortex solution for q = 1 and l = 2 in the double BPS limit.

Case of q = 2
There is an exact analytic solution to (5.11) for the case q = 2 that can be expressed by a hypergeometric function where the normalization factor is given by A plot of the solution is given in figure 7.

Case of q = 3
There is an exact analytic solution to (5.11) for the case q = 3 that can be expressed in terms of Bessel functions where the normalization factor is given by A plot of the solution is given in figure 8.

Case of q ≥ 4
In the case of q = 4, we found an exact solution in terms of confluent Heun functions, see the NIST Digital Library of Mathematical Functions [20, § 31.12]. This (q, l) = (4, 2) solution is not very illuminating because of the unfamiliarity of the confluent Heun functions. For the case q ≥ 5 we were unable to find exact solutions to the linear differential equation (5.11), but we constructed numerical solutions for many q ≥ 2. There are well known subtleties in trying to construct numerical solutions because the natural initial conditions φ(0) = 0 and φ (0) = 1 are numerically unstable due to the regular singular point at ν = 0. To get around this, we compute the terms of the regular power series solution φ ps to (5.11) to O(ν 7 ) with initial conditions φ(0) = 0 and φ (0) = 1: You can verify that the above agrees with the power series expansion of the exact solution (5.12a) for the case q = 2, and (5.14a) for the case q = 3. The idea is to replace the numerically unstable initial conditions at ν = 0 with nearby initial conditions at ν * where Difference between exact and numerical solutions for q=2 0 < ν * ≪ 1. The initial conditions are set by the truncated power series expansion: φ(ν * ) = φ ps (ν * ) and φ (ν * ) = φ ps (ν * ). Typically we chose ν * from 10 −8 to 10 −6 . The numerical solution φ num thus obtained will not have the correct normalization at ν = ∞ but since we have a linear ODE we know that the correct normalized solution will be φ norm (ν) = φ num (ν)/φ num (∞). The difference between the exact solution in q = 2 and the numerical solution is plotted in figure 9. The normalized numerical solutions for q = 2, 3, 4, 8, 24 are shown in figure 10.

Double BPS hedgehog-like defects (l = 3)
The discussion about the finiteness of the energy in section 4.1.4 tells us that we only have to consider the case of a soliton (p = 0 or equivalently q = 1) embedded in AdS 4 . The equations of motion become This ODE admits even and odd solutions, and it has a regular singular point at ν = 0 with Frobenius indices α + = 1 and α − = −2. The α = 1 solution may be taken to be an odd function, and we ignore the α − = −2 solution that may be taken to be an even function. With our initial conditions φ(0) = 0 and φ (0) = 1, we find the power series solution The exact normalized solution is given by the odd function It is easy to verify that the power series expansion (5.21), derived directly from ODE (5.20), agrees with the power series expansion of (5.22) if you remove the normalization factor N (1,3) in order that both functions have derivative 1 at ν = 0. The magnetic charge for the spherically symmetric ansatz is discussed in appendix C. Surprisingly, the exact scalar field profiles for (q, l) = (2, 2) and (q, l) = (1, 3) are very similar, see figure 12.

Static spherically symmetric classical glueballs
Note added after submission to arXiv: On the day our manuscript appeared on the arXiv, there was also a manuscript by Ivanova, Lechtenfeld and Popov [23] where they study finite action and finite energy solutions of Yang-Mills theory in AdS 4 . This work is based on an earlier paper [24]. We did not realize at the time of submission that our double BPS limit   3). In both cases, the field φ satisfies the same two point boundary conditions, also eq. (4.6) tells us that both solutions decay exponentially at the same rate e −3ν/ρ . It is not surprising that both solutions will be similar but the 1% discrepancy is a bit surprising.
monopole solution automatically gives a static spherically symmetric solution of the SO (3) Yang-Mills equations in AdS 4 with finite energy. It should have been obvious. The reason is that in the double BPS limit, eq. (5.18) does not depend on the scalar field φ. This equation is the Yang-Mills equation. All we have to do is ignore the scalar field φ and its equation of motion (5.17). Equivalently, you can look at the transverse energy (5.2) and ignore the scalar field φ. We know the solution is h(ν) = sech ν, see (5.19). The mass of this static spherically symmetric Yang-Mills solution is This solution is a classical glueball in AdS 4 . We do not know if this solution is stable. Also, this solution goes away in the flat space limit ρ → ∞ in accordance with known theorems. We have not explored whether there is a generalization to AdS 4 of the Minkowski space theorems of Coleman [25], and Coleman and Smarr [26] that constrain the existence and the properties of glueballs. It is not immediately apparent how to relate the static Yang-Mills solution in [23, section 6] to the static solution presented here. Another observation is that the discussion of the double BPS limit in the SO(2) case in section 5.2 provides exact maximally symmetric purely electromagnetic p-defects associated with an embedding AdS q → AdS n given by eq. (5.10). The transverse energy density of this solution is given by It is interesting that the transverse energy density is proportional to the dimensionality p of the defect.
6 ρ 2 = 0 and m 2 A ↓ 0 Here we briefly discuss solutions to the spherically symmetric ansatz in the parameter plane defined by m A ↓ 0 with 1/ρ 2 = 0, see the pinkish plane in figure 1. The equations of motion in this limit become Equation (6.1b) already appeared as eq. (5.4b) in section 5 where we analyzed its solutions in detail. There we learned that the solution is given by h(ν) = 1/(cosh ν) q−1 if l = 2 and q ≥ 2, see (5.10). For l = 3, there is only a q = 1 solution h(ν) = 1/ cosh ν, see eq. (5.19).
The strategy is to insert the known h solution into ODE (6.1a), and look for a solution that satisfies the two point boundary conditions φ(0) = 0 and φ(∞) = 1. The analysis is more complicated than in the double BPS limit because the field φ now satisfies a non-linear ODE; the solution can be determined numerically. We use a modification of the power series technique discussed in section 5.2.4 to get around the numerically unstable initial conditions at ν = 0.

Vortex-like defects (l = 2)
The power series solution to the φ equation of motion with initial condition φ(0) = 0 and φ (0) arbitrary is given to O(ν 7 ) by φ ps (ν) = νφ (0) + 1 48 Because the ODE is nonlinear, we see nonlinear behavior on the initial condition φ (0) beginning at O(ν 5 ). We replace the initial conditions at ν = 0 with those given by the power series at a nearby point ν * = 1 × 10 −8 . We use [21] the "shooting method" where we vary φ (0) numerically until we find a solution with approximately the correct asymptotic behavior φ(∞) = 1, see figure 13 and figure 14. A comparison between the numerical solutions and the exact double BPS solution is given in figure 15. in the m A ↓ 0 limit with the Higgs field Compton wavelength 1/m φ = ρ/5. As expected, the decay behavior of the φ field is dominated by the Compton wavelength. As a comparison, the curve in green is the exact double BPS solution where the the Higgs field has infinite Compton wavelength and the length scale is set by the radius of curvature ρ. The shooting method is very sensitive to the initial condition φ (0) due to the potential presence of exponential growth terms; compare the dotted blue curve with the dashed orange curve. Numerical methods become more finicky as m φ ρ increases.
where, once again, = m φ ρ. We see a similar nonlinear behavior to that of the previous case. We use the initial conditions provided by our power series at a sufficiently nearby point. Here we use ν * = 1 × 10 −8 . See figure 16 for a comparison between our numerical solution and the double BPS solution when 1/m φ = ρ/2, and figure 17 for the relative difference between these two solutions.

Flat space equations of motion
For completeness, the flat space equations of motion for spherically symmetric defects are given by taking the ρ ↑ ∞ limit of eqs. (5.3): These equations are well known and have been thoroughly studied over the past 40 years.
The Prasad-Sommerfield monopole solution (l = 3) is given by The scalar field φ has a 1/ν Coulomb tail because m φ = 0.

Conclusions
In this article we studied the equations of motion for maximally symmetric p-defects in AdS n . In the double BPS limit the radius of curvature ρ is the only length scale that appears in the equations of motion, and we were able to find exact analytic solutions in many cases. We also saw that the radial exponential increase in volume in AdS n plays a crucial role, and requires a case by case study for admissible values of (q, l). The method we advocate in this paper is part of a broader strategy to study solutions of the Yang-Mills Higgs system in AdS n . The solutions found in this double BPS limit are the first step in a perturbative expansion in small parameters m φ ρ and m A ρ for the full equations of motion.

A Maximally symmetric submanifolds of maximally symmetric spaces
In this section we derive necessary conditions satisfied by a maximally symmetric submanifold of a maximally symmetric manifold.
Any local orthonormal coframe θ µ on a constant curvature manifold M n with associated Levi-Civita connection ω µν = −ω νµ will satisfy the Cartan structural equations for a manifold of constant sectional curvature k: These constant curvature manifolds are maximally symmetric spaces with dim Isom 0 (M n ) = 1 2 n(n + 1), where Isom 0 (M n ) is the component of the isometry group of M n that is connected to the identity. In fact, eqs. (A.1) are the Maurer-Cartan equations for a Lie group. Next we use the Cartan structural equations and a bit of the theory of exterior differential systems [22], mostly the Frobenius theorem for integrability of a Pfaffian system of equations. Assume Σ q is an isometrically embedded q-submanifold of the constant curvature M n . If we use the index conventions that latin indices from the beginning of the alphabet a, b, c, d run from 1, 2, . . . , q and latin indices from the middle of the alphabet i, j, k, . . . take l = n − q values from q + 1, . . . , n then the structural equations in an orthonormal coframe adapted to the tangent bundle of the submanifold may be written as The submanifold Σ q in this adapted coframe is given by the exterior differential system θ i = 0 and consequently dθ i = 0. Using these conditions in (A.2) we find Applying Cartan's Lemma to (A.3b), we conclude that on Σ q we have where K ab i = K ba i are tensors on Σ q called the extrinsic curvatures or the second fundamental form. When restricted to Σ q , the Cartan structural equation for the intrinsic curvature is dω ab + ω ac ∧ ω cb = 1 2 R Σ abcd θ c ∧ θ d . Thus we obtain the Gauss equation for embedding in a space of constant curvature: where g Σ ab is the induced metric on Σ q due to the isometric embedding.

A.1 Totally geodesic submanifolds
A submanifold Σ q of a general manifold M n is said to be totally geodesic if every geodesic (with respect to the induced metric) on Σ q is also a geodesic on M n . If D Σ and D M are respectively the Levi-Civita connections on the respective manifolds then the definition of the second fundamental form (extrinsic curvatures) tells us that for the tangent vector field X to a curve on Σ we have D M X X = D Σ X X − K i (X, X)n i , where {n i } is a local orthonormal frame for the normal bundle (T Σ) ⊥ . The definition of totally geodesic implies that K i (X, X) = 0 for all X(σ) ∈ T σ Σ and so we conclude that K ab i = 0. It is easy to show that totally geodesics submanifolds exist by constructing an example directly. Fix a point x ∈ M and a q-dimensional vector subspace V x ⊂ T x M . Consider all the geodesics in M n that begin at x with initial velocity in V x . The locus of all these geodesics 12 is a q-dimensional totally geodesic submanifold Σ q ⊂ M n . From the viewpoint of standard General Relativity, totally geodesic submanifolds are very desirable because if a test mass in Σ is given an initial velocity tangential to Σ then its motion will be restricted to Σ.
Next we ask what are the totally geodesic submanifolds of a constant curvature manifold M n . The totally geodesic condition ω ai = 0 means that equations (A.3) restricted to Σ q become Thus we have derived necessary conditions for submanifold Σ q to be a totally geodesic submanifold of the constant curvature manifold M n . The first two equations in (A.6) are the Cartan structural equations for a q-dimensional constant curvature manifold Σ q with the same constant sectional curvature k as M n . These totally geodesic submanifolds are maximally symmetric spaces with dim Isom 0 (Σ q ) = 1 2 q(q + 1). Equation (A.6c) is the statement that the connection on the normal bundle (T Σ) ⊥ is flat.

A.2 Examples of totally geodesic submanifolds
First we consider the example of the sphere of radius a, M n = S n a . The main observation is that the equatorial sphere S n−1 a also has radius a. A geodesic in S n−1 a is a great circle in S n−1 a , and this great circle is also a great circle in S n a . Thus the equatorial S n−1 a is a totally geodesic submanifold of S n a . Repeating the argument we see that the equatorial sphere S n−2 a of S n−1 a is a totally geodesic submanifold of S n−1 a and consequently a totally geodesic submanifold of S n a . We can repeat this argument until get get down to the equatorial S 1 a . Thus we have shown that there exists an "equatorial" S q a that is a totally geodesic submanifold of S n a . To extend the result above to arbitrary k and Euclidean or Minkowski signature for the metric, we repeat the discussion in the previous paragraph in terms of equations. The n-sphere of radius a is the set of points in E n+1 that satisfies the equation (x 1 ) 2 + (x 2 ) 2 + · · · + (x n ) 2 + (x n+1 ) 2 = a 2 . The isometry group of S n is SO(n + 1), the stability group of the North Pole (0, 0, . . . , a) is SO(n), and S n ≈ SO(n + 1)/ SO(n). A polar (n − 1)sphere 13 has radius a and contains the point (0, . . . , 0, a), the North Pole. It is obtained by adjoining to the spherical constraint equation an additional equation x n = 0. To see that this "polar" S n−1 is totally geodesic, we observe that the the vector field ∂/∂x n along this S n−1 is tangent to S n and is the unit normal vector field to the polar S n−1 . This vector is parallel in E n+1 and thus D S n X (∂/∂x n ) = 0 where the vector field X is tangent to S n−1 , and D S n is the connection induced on S n from the Euclidean structure of E n+1 . Therefore, the extrinsic curvature vanishes. You can inductively repeat this argument. A polar (n − 2) sphere containing (0, . . . , 0, a) is obtained by adjoining an additional constraint x n−1 = 0. You can either do an inductive argument or note that the extrinsic curvature vanishes again since the normal bundle is spanned by ∂/∂x n and ∂/∂x n−1 . In this way we can go all the way down to a polar S 1 . This sequence of "polar spheres" is a collection of totally geodesic submanifolds of S n . You can move these totally geodesic spheres to other locations by using a transformation in the isometry group Isom 0 (S n ) = SO(n + 1).
The argument for the vanishing of the extrinsic curvatures is identical in the three cases below and we skip it. The geodesic submanifolds that we construct can be moved by using the appropriate isometry group.
For hyperbolic space H n (Euclidean signature and k < 0), we consider the connected set of points that satisfies (x 1 ) 2 + (x 2 ) 2 + · · · + (x n ) 2 − (x n+1 ) 2 = −a 2 in R n+1 with the signature of the R n+1 metric being (+, +, . . . , +, −) and contains the "North Pole" (0, . . . , 0, a). The isometry group of H n is SO(n, 1), the stability group of the North Pole is SO(n), and thus H n ≈ SO(n, 1)/ SO(n). The "polar" H n−1 is obtained by adjoining the constraint x n = 0. The "polar" H n−2 is obtained by adjoining the additional constraint x n−1 = 0. This argument can be repeated until you get down to a "polar" H 1 ≈ R. The H q constructed this way are all totally geodesic submanifolds of H n .
If we are looking for manifolds of Lorentzian signature then we get de Sitter space dS n (k > 0) or anti de Sitter space AdS n (k < 0).
For de Sitter space dS n , we consider the connected set of points that satisfies −(x 1 ) 2 + (x 2 ) 2 + · · · + (x n ) 2 + (x n+1 ) 2 = a 2 in R n+1 with signature of the R n+1 metric being (−, +, . . . , +, +) and contains the point (0, . . . , 0, a). The isometry group of dS n is SO(n, 1), the stability group of the North Pole is SO(n − 1, 1), and dS n ≈ SO(n, 1)/ SO(n − 1, 1). The "polar" dS n−1 is obtained by adjoining the constraint x n = 0. The "polar" dS n−2 is obtained by adjoining the additional constraint x n−1 = 0. This argument can be repeated until you get down to a "polar" dS 1 . The dS q constructed this way are all totally geodesic submanifolds of dS n .
For anti de Sitter space AdS n , we consider the connected set of points that satisfies −(x 1 ) 2 + (x 2 ) 2 + · · · + (x n ) 2 − (x n+1 ) 2 = −a 2 in R n+1 with signature of the R n+1 metric being (−, +, . . . , +, −) and contains the point (0, . . . , 0, a) 14 . The isometry group of AdS n is SO(n − 1, 2), the stability group of the North Pole is SO(n − 1, 1), and AdS n ≈ SO(n − 1, 2)/ SO(n − 1, 1). The "polar" AdS n−1 is obtained by adjoining the constraint x n = 0. The "polar" AdS n−2 is obtained by adjoining the additional constraint x n−1 = 0. This argument can be repeated until you get down to a "polar" AdS 1 . AdS q constructed this way are all totally geodesic submanifolds of AdS n . We remind the reader that we always have in mind the simply connected universal cover of the quadric hypersurface in question.

A.3 Intrinsically flat submanifolds
Next we show that there are intrinsically flat submanifolds that can be isometrically embedded in negative constant curvature spaces. These manifolds are not totally geodesic submanifolds. We impose the intrinsic flatness condition R Σ abcd = 0. We also require a special form for the extrinsic curvature tensor K ab i = h ab C i , where h ab is the flat metric on Σ. Note that C i is the mean curvature vector. The reason for this form is that we require the extrinsic curvature to be compatible with the isometries of the vector space Σ q . Inserting these conditions into (A.3), we obtain To satisfy the structural equations above, we require from (A.7b) that k = − C 2 ≤ 0. Thus, we have our first result that the space M n must have constant negative curvature 15 since we are restricting our analysis in this paper to k = 0. Next, we observe that since the {θ a } are linearly independent, eq. (A.7c) implies that 0 = DC i = dC i + ω ij C j . Thus the mean curvature vector C = C in i is covariantly constant and we have a preferred direction in the flat normal bundle (T Σ) ⊥ , see condition (A.7d). If M n has Lorentzian signature then the maximal possible symmetry group of a defect is P(q) × SO(l − 1) as opposed to the AdS q case where you can have SO(2, q − 1) × SO(l). There is no spherically symmetric defect in the case of Σ ≈ M q . Said differently, the "transverse part of the Lagrangian" is at most SO(l − 1) invariant. Remark 1. If AdS q is not a totally geodesic submanifold then there is a non-zero mean curvature vector that gives a preferred normal direction. The Gauss equation in this case becomes k Σ q = C 2 + k M n . The symmetry group of the p-defect is maximally Isom 0 (Σ q ) × SO(l − 1). There are no maximally symmetric solutions in this case.
You can verify that the world branes in references [7,9,10] are totally geodesic submanifolds.
It is easy to discuss a basic example that is commonly used in physics. Consider AdS n with embedded submanifold Σ q ≈ M n−1 , i.e., q = n − 1. The metric on AdS n is written in the "upper half space" form If we restrict to the q = n − 1 dimensional submanifold Σ q defined by y = y 0 , y 0 constant, then the induced metric on it, is a rescaled Minkowski metric. The mean curvature vector is parallel to the unit normal to Σ. Any vector subspace of this Σ n−1 is automatically a flat submanifold of AdS n . In figure 18 we compare a totally geodesic constant curvature submanifold and a zero curvature submanifold of AdS n in the "upper half space" representation.

B Relation to the work of Lugo, Moreno and Schaposnik
This is an attempt at relating the work of Lugo, Moreno and Schaposnik (LMS) with ours. They choose a different coordinate system for AdS 4 where the metric is It is important to notice that r is not the radial distance. You can verify that Σ 1 , the "t-axis", is a totally geodesic submanifold. A notational dictionary is in table 2.
A first observation is that we have a different viewpoint on what is BPS. Their philosophy is described in the paragraph following eq. (4.32) in reference [9]. They assume that setting λ LMS = 0 means that there is no potential and the asymptotic value of the Higgs field can float to some non-zero value H(∞) as you solve the equations of motion. This is natural because in their numerical work they are using a relaxation method [21] taking Figure 19: The gray plane above, taken from figure 1, corresponds to the generalized BPS limit where m 2 φ ↓ 0. We have reinterpreted the results of LMS as the dashed straight lines m 2 A = ν • (ν • + 1)/ρ 2 in this plane, where LMS discovered the existence of a power series solution if ν • is a positive half integer. LMS also studied other values of ν • numerically in [10]. The BPS line corresponds to formally setting ν • = ∞, which is the same as sending ρ ↑ ∞. The double BPS line corresponds to formally setting ν • = 0, which is the same as sending m A ↓ 0. The work of LMS seems to indicate that solutions along the dashed lines are special but we have not explored this. Note that those dashed lines interpolate between the two lines where we have exact analytic solutions. LMS found reasonable numerical solutions using a relaxation method with starting point the BPS solution; this is not surprising in light of this figure. On the other hand, this figure suggests that in some sense, the double BPS solutions are far away from the Prasad-Sommerfield solutions.
C The lack of Bogomolny equations and a partial bound

C.1 The magnetic flux
In this appendix we make an unsuccessful attempt to determine the stability of our double BPS monopole solution. We analyze whether ours is a minimum energy static solution. We discuss various topics related to static magnetic monopoles in AdS 4 . We address whether in the λ ↓ 0 Prasad-Sommerfield limit [15] there is an analog of the first order Bogomolny equations [16] that gives solutions of the equations of motion. More precisely, we show that the arguments [16,17] that lead to the Bogomolny bound in the M 4 theory do not generalize to the AdS 4 case. In this section we ignore the potential term because we are always implicitly taking the BPS limit. The value of g is generic.
First, we compute the magnetic flux. The metric on M = AdS 4 may be written in the form ds 2 AdS 4 = − cosh 2 |k| 1/2 ν dτ 2 + dν 2 + sinh |k| 1/2 ν |k| 1/2 ν 2 ν 2 ds 2 S 2 , where the coordinate τ is the proper time along the 1-dimensional defect world line Σ 1 = AdS 1 . The restriction of the AdS 4 metric to the normal tangent space (T σ Σ) ⊥ is given by At σ ∈ Σ 1 , the normal tangent space (T σ Σ) ⊥ is a three three dimensional vector space. Define the magnetic 2-form by B I = 1 2 IJK F JK and observe that d Φ I B I = DΦ I ∧ B I as a consequence of the Bianchi identity DB I = 0. If we fix σ ∈ Σ 1 and we apply Stoke's Theorem we find where S 2 ν → (T σ Σ) ⊥ is the 2-sphere of radius ν. Here we used the definition (TσΣ) ⊥ · · · = lim ν→∞ B 3 ν · · · , where B 3 ν → (T σ Σ) ⊥ is the 3-ball of radius ν and S 2 ν = ∂B 3 ν . As shown by 'tHooft [3] and Polyakov [4], the SO(2) unbroken gauge symmetry magnetic field 2-form is Φ I B I /φ 0 and thus we expect that, for a monopole, the abelian magnetic flux F m , which is a measure of the magnetic charge, is given by: Thus, we conclude that The integrand in the left hand side of the equation may be interpreted as the magnetic charge density. Note that the conventionally normalized magnetic field B conv is related to our magnetic field by B = g B conv , and one has to be careful in analyzing the g ↓ 0 limit. The topological flux quantization condition should hold for all values of g. Next, we want to see how the above computation works within our ansatz for the spherically symmetric hedgehog because this illuminates why there are no Bogomolny equations in this case. Using eq. (3.10) we obtain where ζ S 2 ν is the area 2-form on S 2 ν → (T σ Σ) ⊥ . Note that that S 2 ν has area 4πρ 2 sinh 2 (ν/ρ). Due to its spherical symmetry, the integral is easily computed and we obtain the result found by 'tHooft and Polyakov: where we used the asymptotic boundary conditions. With our normalization, the magnetic flux is F m = 4π, which is the correct topological quantization condition for the case where the SO(3) gauge group breaks down to SO(2) via a vector representation scalar field 16 . Note the cancellation of the hyperbolic sine factors, which was necessary to obtain a finite non-zero magnetic flux. Next, we discuss a couple of attempts at trying to obtain the Bogomolny bound and the associated Bogomolny equations. We implicitly assume that we are always taking the Prasad-Sommerfield limit λ ↓ 0, and maintaining the appropriate boundary conditions.
It is computationally useful to use the Hodge inner product. We give a brief review to establish conventions. Let V be a vector space with inner product, either Euclidean or Lorentzian signature, dim V = n. The normalized volume element determined by an orthonormal basis is ζ. If ω is a k-form, it is convenient to define a norm by ω 2 = ω µ 1 ···µ k ω µ 1 ···µ k /k!. With this normalization each independent term in the k-form only contributes once in the summation. A pointwise bilinear product, •, • , and the Hodge dual operator on k-forms are defined by ω, ω ζ = ω ∧ ω = ω 2 ζ. If s = +1 in Euclidean signature and if s = −1 in Minkowski signature then please note that ω 2 = s ω 2 , and ω = s(−1) k(n−k) ω.

C.2 From the magnetic flux to the energy functional
In this section, we fix σ ∈ Σ and we restrict to the three dimensional normal tangent space (T σ Σ) ⊥ . In particular, the metric on (T σ Σ) ⊥ is given by (C.2), and the Hodge duality operator refers to Hodge duality with respect to this metric. The strategy in this section is to begin with the correct expression for the magnetic flux and to try to get to the energy functional whose variation gives the equations of motion.
In this appendix, we have to use Hodge duality on the normal tangent space (T σ Σ) ⊥ , and on the full manifold AdS 4 . For this reason we introduce a special notation: on the normal tangent space (T σ Σ) ⊥ with metric (C.2), we denote the Hodge duality operation by , and the associated norm by • 2 .
We observe that Here ζ (TσΣ) ⊥ = 1 is the volume element on (T σ Σ) ⊥ . Rearranging terms we find This is the expression [16,17] that is used to show in M 4 that a solution to the Bogomolny equations DΦ I = B I /g is an absolute minimum of the left hand side of (C.8).
Unfortunately, this is not what we need because the functional that has to be minimized to obtain the equations of motion is not the left hand side of (C.8) but the transverse energy functional (3.4), which in Hodge star notation is The hyperbolic cosine factor in the previous equation is necessary to obtain the correct equations of motion. The static equations of motion are obtained by restricting the 4dimensional action, constructed with the AdS 4 metric, to time translationally invariant field configurations. Note that ∂/∂τ is a Killing vector since the metric components in (C.1) are independent of τ . Equation (C.8) is what you would get for an instanton solution in Euclidean hyperbolic 3-space H 3 . Solutions to the BPS equations would be absolute minima of the equations of motion in the Prasad-Sommerfield limit, see also the discussion in Atiyah [14].
The Bogomolny argument works in a situation where the spacetime is a product manifold M 4 = R × V 3 with the product metric ds 2 M = −dτ 2 + ds 2 V . It also works in the Euclidean signature instanton case where the manifold M 3 is Euclidean hyperbolic 3-space, ibid. and footnote 1.

C.3 From the energy functional to the magnetic flux
In this section, we begin with the correct functional that has to be minimized to obtain the static equations of motion and we try to see if we can find a bound for this functional that is related to the magnetic flux . From now on, the metric is the AdS 4 Lorentzian metric (C.1) with associated inner product •, • . The Hodge duality operation is with respect to this metric, the AdS 4 volume element is denoted by ζ = 1, and the normalized timelike 1-form isθ 0 = cosh(ν/ρ) dτ . The action for our static configuration is Next we rewrite the action as To try to get Bogomolny type equations, we observe thatθ 0 ∧ DΦ is a timelike 2-form, and that B is also a timelike 2-form. Note thatθ 0 is orthogonal to DΦ and thus, we can If we set ρ = ∞, i.e., k = 0, then (C.11) leads to the standard Bogomolny type argument for a lower bound on the mass of 4π/g. The bound is saturated by field configurations that satisfy the Bogomolny equations DΦ = B/g. If 0 < ρ < ∞ and • 2 is a positive definite norm, then this equation implies a bound This bound does not appear very useful because the bound is field configuration dependent. You can verify the field dependency by computing the variation of the integral with respect to a variation of the scalar field. If δΦ(ν) is a variation of the scalar field with compact support, then you obtain Note that the right hand side of the equation above vanishes for ρ = ∞, and it exemplifies the topological nature of the magnetic flux (C.5), which is invariant under deformations of the field configuration.

C.4 A bound on the mass for non-negative magnetic charge density
We previously mentioned that the left hand side of magnetic flux equation (C.5) may be interpreted as the magnetic charge density. Here we prove a theorem that states that if the magnetic charge density is non-negative in the Prasad-Sommerfield limit and if 0 < ρ < ∞, then the mass satisfies the strict inequality E ⊥ (ρ) > 4πφ 0 /g. A 3-form α on (T σ Σ) ⊥ is said to be non-negative if there exists a function f : (T σ Σ) ⊥ → R with f ≥ 0 such that α = f ζ (TσΣ) ⊥ . Let C = {(Φ, A)} be the space of admissible field configurations, and let C + ⊂ C be the subset of field configurations where the 3-form DΦ I ∧ B I = d Φ I B I is non-negative. These are the configurations with non-negative magnetic charge density. Since cosh(ν/ρ) > 1 for ν > 0 and 0 < ρ < ∞, the use of (C.12), (C.5) and (C.7) leads to the conclusion that if we restrict to field configurations in C + , then there is a lower bound provided by magnetic flux quantization