arxiv:gr-qc/ v1 2 May 2005

Transcription

1 Spatially Compact Solutions and Stabilization in Einstein-Yang-Mills-Higgs Theories arxiv:gr-qc/ v1 2 May 2005 Péter Forgács and Sébastien Reuillon Laboratoire de Mathématiques et Physique Théorique CNRS-UMR 6083 Université de Tours, Parc de Grandmont Tours, France November 5, 2018 Abstract New solutions to the static, spherically symmetric Einstein-Yang-Mills-Higgs equations with the Higgs field in the triplet resp. doublet representation are presented. They form continuous families parametrized by α = M W /M Pl (M W resp. M Pl denoting the W-boson resp. the Planck mass). The corresponding spacetimes are regular and have spatially compact sections. A particularly interesting class with the Yang-Mills amplitude being nodeless is exhibited and is shown to be linearly stable with respect to spherically symmetric perturbations. For some solutions with nodes of the Yang-Mills amplitude a new stabilization phenomenon is found, according to which their unstable modes disappear as α increases (for the triplet) or decreases (for the doublet). By now a lot is known about classical non-abelian particle-like and black hole solutions in Einstein-Yang-Mills (EYM) and in Einstein-Yang-Mills-Higgs (EYMH) theories, we refer to [1] for a review and references. The existence of such objects had not been anticipated and their study has revealed some important non-perturbative features of the interaction between the gravitational and gauge fields. In addition to particle-like solutions where the corresponding spacetime is asymptotically flat, a rather different class has been found in the standard SU(2) EYMH theory (i.e. without a cosmological constant) with the Higgs field being in the triplet representation [2]. They form continuous families parametrized by the gravitational coupling strength, α = M W /M Pl, andthedifferent families arecharacterized bytwointegers, thenumber of nodes of the Yang-Mills amplitude, m 1, resp. of the Higgs field, n 0. The spacetimes corresponding to these solutions have spatially compact sections, i.e. they are topologically R S 3, representing static non-abelian Einstein-universes. The solutions of Ref. [2] generalize those found in EYM theory with a positive cosmological constant Λ (EYM-Λ) [3]. Intuitively the self-interaction potential of the Higgs field acts as a dynamical cosmological constant providing the necessary energy density to support a spatially compact spacetime. The very existence of most of these classical solutions is quite surprising and their study is important to understand fundamental phenomena in gauge theories coupled to gravitation. The important problem of stability of the EYM and EYMH asymptotically flat solutions has been studied mostly in linear perturbation theory [4, 5], with the disappointing conclusion that, with the exception of the gravitating analogues of the t Hooft-Polyakov magnetic monopoles, all of them are unstable. The stability analysis of the compact solutions in the EYM-Λ theory also shows the existence of unstable modes [6, 7]. Interestingly it has been recently found that, allowing for a negative cosmological constant, the EYM-Λ equations admit nodeless solutions with asymptotically anti-de Sitter (AdS) spacetimes which have been shown to be stable [8, 9]. 1

2 In this Letter we present globally regular, static and spatially compact solutions in the SU(2) EYMH theory with the Higgs field being in the doublet representation, i.e. the nonabelian bosonic subsector of the standard model of electroweak interactions. Similarly as for the triplet Higgs case [2], there exist distinct continuous families parametrized by α and indexed by the number of nodes of the YM resp. Higgs field, (m,n). We exhibit a new class of compact solutions for both the doublet and triplet case where the Yang-Mills amplitude has no nodes. These new classes have no analogues in EYM-Λ theories and their existence is related to an interplay of the scalar potential and the gravitational field. They could also be interpreted as gravitationally bound (or rather confined) monopole-antimonopole (triplet) or sphaleronantisphaleron (doublet) pairs. We have also investigated the stability properties of the compact solutions. The spectrum of small fluctuations can be divided with respect to parity into even ( gravitational ) and odd sectors. A detailed numerical investigation strongly indicates linear stability of these nodeless solutions, at least against spherically symmetric perturbations for both the triplet and the doublet case in the gravitational sector of perturbations which makes their further study worthwhile. In the odd parity sector of perturbations an analytic argument shows that for the triplet case the number of unstable modes is simply equal to the number of nodes, m, of the Yang-Mills field, establishing that the nodeless solutions have no unstable modes at all. In the doublet case, however, a suitable adaptation of the result of Ref. [10] leads to the result that there is always an unstable mode independently of m. This unstable mode does not couple directly to the metric fluctuations and it might just indicate that the static doublet solution becomes time-dependent, possibly of oscillatory type since 3-space is compact and there are no radiative modes as in an asymptotically flat space-time. In addition, we have found yet another interesting phenomenon stabilization, which occurs in the even parity sector. Our stability analysis has shown that for the triplet case, the (m = 1, n = 0) family of compact solutions discovered in [2] has two unstable modes in the gravitational sector. Quite surprisingly as α increases the two unstable modes actually merge, and for larger values of α no unstable modes are found. This phenomenon can beinterpreted so that the solutions undergo gravitational stabilization since they loose their unstable modes for increasing values of the gravitational coupling strength. To the best of our knowledge this is the first example that a system gets stabilized when the gravitational coupling gets stronger. In fact exploiting the (somewhat unphysical) limit α we have strong numerical evidence that no instabilities occur for arbitrarily strong values of the gravitational coupling. Remarkably there is also a family of solutions (m = 2, n = 0) in the doublet case which looses its two unstable modes in the gravitational sector of perturbations. This occurs, however, for decreasing values of α. The classical solutions we find exist in the regime of strong gravitational coupling at energies close the Planck scale, and therefore they would be mostly relevant for semiclassical computations and to gain more understanding of the classical phase space of EYMH theories. They seem to be quite distinguished by their stability properties and also by the fact that they are present in the standard EYMH theories (i.e. without the introduction of new interactions, or negative vacuum energies) making them of quite some interest. The most general spherically symmetric line element is ds 2 = e 2ν(R,t) dt 2 e 2λ(R,t) dr 2 r 2 (R,t)dΩ 2, (1) and the most general spherically symmetric Ansätze for the gauge resp. for the Higgs fields can be written as A=(A 0 dt+a 1 dr+cosθdφ)t 3 +(W T 1 +KT 2 )dθ +(W T 2 KT 1 ) sinθdφ, Φ=H.1+i F T 3 (doublet), Φ = H T 3 (triplet) (2a) (2b) 2

3 where T a (a = 1,2,3) are the generators of SU(2) and A i,w,k,h,f are functions of R and t. To look for compact solutions we use the minimal Ansatz compatible with the field equations, corresponding to A 0 = A 1 = K = F = 0. Following the notations of [2], we define N := ṙ, κ := N +r ν, where f := df/dσ := e λ df/dr, and the static field equations can be written as Ṅ = κ N N r r drḣ2, (3a) κ=[1 κ 2 +2Ẇ2 β2 r 2 2 (H2 α 2 ) 2 2C]/r, (3b) Ẅ =W W2 1 r C 2 W Ẇ, (3c) r Ḧ = 1 ( β 2 d 2 H(H2 α 2 )+ 1 ) C r 2 κ+n Ḣ, (3d) H r where β = M H /M W, M H being the mass of the Higgs field, and d = 2, C = (W 1) 2 H 2 /2 resp. d = 1, C = W 2 H 2 for the doublet resp. for the triplet case. It should be noted that for the doublet case α is M W / 2M Pl. Additionally there is a constraint equation relating {r,n,κ,h,w} to the first derivatives of the matter fields. Eqs. (3) admit a family of local solutions regular at the origin where r = σ +O(σ 3 ) and N = κ = 1 + O(σ 2 ). For the triplet case, the behaviour of the YM resp. Higgs field at the origin is W = ±(1 bσ 2 ) + O(σ 4 ) resp. H = aσ + O(σ 3 ) with a, b being free parameters. One can impose without restriction W(0) = +1. For the doublet case, however, because of the asymmetry of the coupling term, C, between the Higgs and YM fields with respect to W W, there are two distinct families of local solutions corresponding to the sign of W(0): either W = (1 bσ 2 ) + O(σ 4 ) and H = aσ + O(σ 3 ), or W = +(1 bσ 2 ) + O(σ 4 ) and H = a+o(σ 2 ). For H(σ) 0 the EYMH eqs. reduce to those of an EYM-Λ theory with the cosmological constant given as Λ = α 4 β 2 /4. As foundin Ref. [3] for a special set of values of the cosmological constant {Λ = Λ m,m = 1...}theEYM-Λ theory admits regularcompact solutions. Them = 1 solution is known analytically [11]: r = 2sin(σ/ 2), N = κ = W = cos(σ/ 2), (4) with Λ 1 = 3/4. Due to the compactness of the manifold (S 3 in this case) there is an extremum of the radial function, ṙ(σ e ) = 0, to which we refer as an equator. Following Ref. [2] a natural strategy to search for solutions of Eqs. (3) is to look for solutions bifurcating with those of the EYM-Λ theory. At the bifurcation point the bifurcating solutions are normalizable zero modes of the linearized Higgs field Eq. (3d) in the corresponding EYM-Λ m background solution. In the doublet case the equation for the linearized Higgs field amplitude, h, takes the form: ( (Wb 1) 2 ḧ = 2r 2 b Λ m α 2 ) h κ b +N b r b ḣ, (5) where r b,n b... etc. stand for the background solution. The bifurcation parameter, α, is to be varied to findglobally regular solutions of Eq. (5). In the analytically known m = 1 background (4) the bounded zero modes of the Higgs-field, h n, can be explicitly constructed. They are characterized by an additional node number, n 0, and by the values of the bifurcation parameter α 2 = 6/((2n+1)(2n+5) 2), n = 0,1,..., (yielding 1,n α2 = 2, 1,0 α2 1,1 = 6/19,...). Due to the asymmetric form of C for the doublet, h n is neither symmetric nor antisymmetric around the equator, and the solutions of the full system are asymmetric for odd values of the Yang-Mills node number. For odd values of m the gauge field amplitude interpolates between W = +1 (at the north pole) and W = 1 (at the south pole) of the compact manifold. 3

4 For even values of m, there are two classes of solutions, both symmetric around the equator distinguished by the values of W(0) = ±1 at the origin (referred to as type + resp. ). We give a detailed description of these families in [12], here we exhibit a particularly interesting class of type W(0) = 1 with m = 2 and n = 0, denoted as (m = 2, n = 0). This solution bifurcates with the m = 2 compact solution of the EYM-Λ theory at α 2,0 = with Λ 2 = It exists for a finite range of values of α [α min,α 2,0 ] with α min 0.54(9). For decreasing values of α, the two zeros of W move towards the equator and they merge there for α = For α < α, these solutions become nodeless. One can easily show that this cannot happen for the (m = 2,n = 0) + solutions. For α α min the spatial sections of the limiting solution will no longer be compact (σ ). For the triplet case, the corresponding analysis has been performed in [2] where various oneparameterfamilies ofsolutions ofeqs.(3), indexedbythenodenumbersm 1, n 0, have been found. These solutions are symmetric or antisymmetric around the equator depending on the parity of m and n. Here we exhibit a new type of family of solutions in the triplet case with m = 0, n = 0. This family does not bifurcate with any known EYM-Λ compact solution and it seems to exist for α [α 0, ) where α (7). In the strong coupling limit α, Eqs. (3) reduce to those of an EYMH system with β = 0 and with a cosmological constant Λ = 1/4. In this limit the equator, and hence the zero of W, tend to infinity, and the spatial sections of the solutions will no longer be compact (σ and r 2). The limiting metric has the asymptotic form ds 2 = σ 2 dt 2 dσ 2 2dΩ 2, corresponding to the product of a 2-d Minkowski space (in Rindler coordinates) and a 2-sphere. In the limit α α 0 the solution interpolates between a regular origin and a non-abelian fixed point of Eqs. (3) (see Ref. [13]) with non-compact spatial sections. The nodeless solutions for both the triplet and doublet cases are displayed on Figs. 1, r,κ N W,H α = β = σ Figure 1: A nodeless triplet solution for α close to α 0. Stability Analysis. We shall restrict our stability analysis to spherically symmetric fluctuations around the solutions. In many cases instabilities show up already in the spherically symmetric sector. The complete stability analysis of the asymptotically AdS EYM-Λ solutions performed in Ref. [9] has demonstrated that solutions stable under spherically symmetric perturbations remain stable under non-spherically symmetric perturbations. The general stability analysis is considerably more involved, and it is further complicated for spatially compact solutions due the existence of an equator. The most general spherically symmetric fluctuations split into two sectors according to spatial parity reflections. Making a mode decomposition of the 4

5 W,H r,κ N α = 0.7 β = σ Figure 2: A nodeless doublet solution. type f(r,t) = f(r)+δf(r)e iωt, where f(r) denotes a static solution, the linearized EYMH Eqs. fall into two sets. In the even parity ( gravitational ) sector, the matter perturbations (δw, δh) couple to the gravitational ones (δν, δλ, δr). The odd parity ( sphaleronic ) sector involves only the matter perturbations δa i, δk and δf (for the doublet). In the even parity sector, there exists a particular gauge the Schwarzschild gauge defined by δr = 0, in which the metric perturbations can be eliminated and one obtains a closed system of Schrödinger equations for the perturbation amplitudes Φ := (δw, δh) of the form: e ν (e ν Φ) = (V +ω 2 )Φ, (6) where the 2 2 potential matrix, V, can be expressed explicitly in terms of the functions {r,n,e ν,h,w} and their derivatives. For the limiting solutions (α ) the potential matrix is bounded and one can use a generalized version of the Jacobi criterion [14] to detect bound states with negative eigenvalues (ω 2 < 0), corresponding to unstable modes. According to this criterion, the vanishing of a suitable determinant of the solutions of Eqs. (6) in (0, + ) implies the existence of at least one bound state with a negative eigenvalue. We have numerically integrated (6) in the background of the limiting (m = 0, n = 0) solution and found that it has no unstable mode. For solutions with compact spatial sections (i.e. topologically S 3 ), the potential matrix, V, in Eqs. (6) is singular due to the presence of an equator, see Refs. [6, 7]. One can avoid such problems by working in a globally regular gauge, e.g. δλ = 0. In this gauge, however, one obtains a somewhat unusual eigenvalue problem, which can be nevertheless dealt with numerically. Starting at the bifurcation points where the numerical values of the negative eigenvalues are known [7] it is possible to determine their dependence on α. We remark that the zero mode of the Higgs field at the bifurcation points, becomes an unstable mode for the (m = 1, n = 0) triplet solution, comp. Fig. 3, while it goes into the positive spectrum for the (m = 2, n = 0) doublet solution, comp. Fig. 4. The negative eigenvalues behave quite similarly to what has already been observed for certain families in the EYM-Λ system [7]. It is quite remarkable that in the present cases the solutions loose all of their unstable modes. For the nodeless triplet solutions, we have undertaken a careful numerical search for unstable modes (scanning in the parameter space) and the lack of any evidence for such modes, together with the result of the Jacobi criterion in the α limit strongly suggests that this solution has no unstable modes for any α [α 0, ). 5

6 0-0.1 ω ω α 1, α Figure 3: Stabilization for the (m = 1, n = 0) triplet solution. The bifurcation point is α 1,0 = In the odd parity sector, the stability analysis is simplified by the fact that the metric and matter perturbations decouple. By a suitable adaptation of the main ideas of Refs. [4, 6] we have proved the following: a) in the triplet case, solutions with m nodes for the YM amplitude have exactly m unstable modes; b)in the doublet case, all solutions have at least one unstable mode. Therefore we conclude that the nodeless triplet solutions are linearly stable against spherically symmetric perturbations in both sectors, while the (m = 1, n = 0) triplet solutions stabilize in the even parity sector but still possess one unstable mode in the odd parity sector. The (m = 2, n = 0) doublet solutions are stable for α in the even parity sector and have at least one unstable mode in the odd parity sector. It is not unlikely that the gravitational stability persists at the non linear level, and no singularity forms. Then, the instability of the static solution indicates merely that it becomes time-dependent (of oscillatory type). References [1] M. S. Volkov and D. Gal tsov, Phys. Rept. 319, (1999). [2] P. Breitenlohner, P. Forgács and D. Maison, Phys. Lett. B489, 397 (2000). [3] M. S. Volkov, N. Straumann, G. Lavrelashvili, M. Heusler and O. Brodbeck, Phys. Rev. D54, 7243 (1996). [4] M. S. Volkov, O. Brodbeck, G. Lavrelashvili and N. Straumann, Phys. Lett. B349, 438 (1995); [5] N. Straumann and Z. H. Zhou, Phys. Lett. B237, 353 (1990);ibid. 243, 33 (1990); Nucl. Phys. B360, 180 (1991); G. Lavrelashvili and D. Maison Phys. Lett. B343, 214 (1995); N. E. Mavromatos and E. Winstanley Phys. Rev. D53, 3190 (1996). [6] M. S. Volkov, N. Straumann, G. Lavrelashvili, M. Heusler and O.Brodbeck, Phys. Rev. D54, 7338 (1996). [7] P. Forgács and S. Reuillon, Phys. Lett. B568, 291 (2003); 6

7 1 0 ω ω ω α α 2,0 Figure 4: Stabilization for the (m = 2, n = 0) doublet solution. The bifurcation point is α 2,0 = [8] J. Bjoraker and Y. Hosotani, Phys. Rev. Lett. 84, 1853 (2000); E. Winstanley, Class. Quant. Grav. 16, 1963 (1999); [9] O. Sarbach and E. Winstanley, Class. Quant. Grav. 18, 2125 (2001); ibid. 19, 689 (2002); [10] P. Boschung, O. Brodbeck, F. Moser, N. Straumann, M. S. Volkov Phys. Rev. D50, 3842 (1994) [11] J. Cervero and L. Jacobs, Phys. Lett. B78, 427 (1978); M. Henneaux, J. Math. Phys. 23, 830 (1982); Y. Hosotani, Phys. Lett. B147, 44 (1984). [12] P. Forgács and S. Reuillon, in preparation. [13] P. Breitenlohner, P. Forgács and D. Maison, Nucl. Phys. B 442, 126 (1995). [14] H. Amann and P. Quittner, J. Math. Phys. 36, 4553 (1995). 7