Analytic non-Abelian gravitating solitons in the Einstein-Yang-Mills-Higgs theory and transitions between them

Two analytic examples of globally regular non-Abelian gravitating solitons in the Einstein-Yang-Mills-Higgs theory in (3+1)-dimensions are presented. In both cases, the space-time geometries are of the Nariai type and the Yang-Mills field is completely regular and of meron type (namely, proportional to a pure gauge). However, while in the first family (type I) $X_{0} = 1/2$ (as in all the known examples of merons available so far) and the Higgs field is trivial, in the second family (type II) $X_{0}$ is not 1/2 and the Higgs field is non-trivial. We compare the entropies of type I and type II families determining when type II solitons are favored over type I solitons: the VEV of the Higgs field plays a crucial role in determining the phases of the system. The Klein-Gordon equation for test scalar fields coupled to the non-Abelian fields of the gravitating solitons can be written as the sum of a two-dimensional D'Alembert operator plus a Hamiltonian which has been proposed in the literature to describe the four-dimensional Quantum Hall Effect (QHE): the difference between type I and type II solutions manifest itself in a difference between the degeneracies of the corresponding energy levels.


Introduction
Non-Abelian solitons and instantons are fundamental pillars in our understanding of gauge theories beyond perturbation theories (see [1] [2] [3] [4] [5] and references therein). In many situations (such as in the cases of topological defects in the early universe), the effects of gravity cannot be neglected (see [6] and references therein). From a genuine General Relativistic viewpoint, non-Abelian gravitating solitons represent severe tests for well-known conjectures such as the no-hair conjecture. Consequently, it is a mandatory task to shed further light on these types of gravitating solitons. The gravitating solitons which have been analyzed in more details in the eighties and nineties are asymptotically flat and (almost always) numerical. However, it is worth noting that the requirement of asymptotic flatness is somewhat "foreign" to the requirement of regularity in the sense that while gravitating solitons must be regular by definition, genuine non-Abelian and globally regular gravitating solitons need not to be asymptotically flat 1 . The first example was discovered by Bartnik and McKinnon (BK) [7]. Soon after the BK gravitating solitons, genuine non-Abelian black holes were also constructed numerically in [8] [9] [10]. These results closely related to the no-hair conjecture and to the black holes uniqueness theorem attracted a lot of attention since then (see, for instance, [11] [12] [13] [14] [15] [16] [17] and references therein). In the present paper, we will be mainly interested in the construction of analytic non-Abelian gravitating solitons. To the best of our knowledge, the only known analytic example in (3+1) dimensions of a globally regular gravitating soliton with "bona fide" non-Abelian gauge field has been found in [18] (many nice numerical examples are described in [11] [12] [13] [14] [15] [16] [17] and references therein). Such remarkable analytic solution has been found in the Einstein-Yang-Mills-dilaton system: so far it has not been possible to extend this solution to Einstein-Yang-Mills theory without dilatonic coupling. The main goal of the present manuscript is to construct two types of globally regular analytic non-Abelian gravitating solitons in the Einstein-Yang-Mills-Higgs system in (3+1)-dimensions without any dilatonic coupling and to analyze the transitions between these families as well as their remarkable physical properties.
An obvious question is: why should one insist so much on finding analytic solutions if these equations can be solved numerically? Since the pioneering works of Bartnik and McKinnon mentioned above, very powerful numerical techniques have been proposed in the literature to construct non-Abelian gravitating solitons (see, for instance, [11], [19] [20] [21] and references therein). There are really good reasons to strive for analytic solutions nevertheless.
Besides the obvious fact that a systematic tool to construct analytic gravitating solitons can greatly enlarge our understanding of these configurations, the main reason is that the availability of analytic gravitating non-Abelian solitons allows to disclose phenomena that would be very difficult to see otherwise. In particular, the present formalism discloses the possibility to have transitions between these solitons and a very surprising analogy with the four-dimensional Quantum Hall Effect (QHE) which will be discussed in a moment.
As we are seeking for genuine non-Abelian effects, meron types configurations (introduced in [22]) are really good candidates as such configurations can only appear in non-Abelian gauge theories (see [2] and references therein). A gauge potential of meron type can be defined as proportional to a pure gauge: such a gauge potential is non-trivial due to the presence of the commutator in the non-Abelian field-strength. On the other hand, in Abelian gauge theories, a gauge potential which is proportional to a pure gauge is itself a pure gauge 2 and therefore is trivial. Consequently, merons are bona fide non-Abelian configurations. Since the pioneering works [23] [24], the important role of these types of configurations in understanding the non-perturbative phase of QCD has been widely recognized (see [25] and references therein).
All the known examples of merons so far have X 0 = 1/2; thus, a first question that we will answer (affirmatively) is: is it possible to have non-trivial merons configurations with X 0 = 1/2 ?
The second and most difficult issue related to merons is the following. Merons on flat spaces are necessarily singular: they play an important role (see [23]  The paper is organized as follows: In Section II, we present the model and give our ansatz and the corresponding field equations. In Sections III and IV, we construct the two families of the analytic regular solutions of the Einstein-Yang-Mills-Higgs theory in (3+1)-dimensions. In Section V, we introduce a standard coordinate system of the Nariai class of spacetime to work out the surface gravity and the temperature of the system. In Section VI, we present the entropy functions of the two types of the solutions and provide some useful plots of the entropy functions to see which configuration is favored for given sets of the physical parameters. In Section VII, we analyze the Klein-Gordon equation to disclose the difference between two types of solutions by making use of the four-dimensional quantum hall effects. In the last Section, our conclusions are drawn.

Action and ansatz
The starting point is the action of the (3+1)-dimensional Einstein-Yang-Mills-Higgs field: where R, Λ, and V (Φ * Φ) are the Ricci scalar of the space-time, the cosmological constant, and the self-interacting Higgs potential, respectively. The dimensionless constant e is the gauge coupling and κ = 8πG for the gravitational constant G. The field strength of the Yang-Mills field A µ is F µν = ∂ µ A ν − ∂ ν A µ + [A µ , A ν ] and the gauge-covariant derivative is D µ = ∇ µ + [A µ , · ]. In this notation, the Einstein's equations are written as where the energy-momentum tensor is given by The equations for the Higgs and Yang-Mills fields are The Higgs potential is given by where λ is the self-interacting coupling constant of the Higgs fields. We consider this Einstein-Yang-Mills-Higgs (briefly, EYMH) system in the space-time with the metric given by where L and F 0 are constant (without loss of generality F 0 can be assumed to be positive). As has been already emphasized, the Yang-Mills field is assumed to have the meron form where X 0 is a constant such that X 0 = 0, 1. The SU (2)-valued scalar field U is parametrized as where t A = iσ A for the Pauli matrix σ A The Higgs field is given in an adjoint representation by where the constant W 0 has to be determined solving the Higgs field equations which (with the ansatz defined above) reduce to the single algebraic equation The Yang-Mills equations (4) also reduce to the following single algebraic equation; From Eqs. (6) and (7) it is clear that there are two families of solutions.
The first family (which will be called type I) corresponds to the usual meron solution together with the condition that the Higgs field profile W 0 is in the "VEV": The second family (which will be called type II) corresponds to the conditions The explicit form of the type II solution will be discussed in the next sections.
3 Type I family: standard Meronic configuration It follows from Eq. (6) that the value of W 0 is precisely the vacuum expectation value v 0 if and only if the configuration of the gauge field is a meron: In this case, the Eqs. (6) and (7) are automatically satisfied and the Eq. (10) fixes the value of F 0 L 2 , which is the size of the S 2 , in terms of the cosmological constant Λ and of the other parameters: This equation admits one or two positive roots for F 0 L 2 and the number of roots depends on the range of physical parameters as follows:

Analytic Solutions of type II
The Higgs equation (6) has the solution only when the constant W 0 lies in the range of When W 0 = v 0 the configuration of the Yang-Mills field becomes that of a standard meron X 0 = 1/2 (and the Higgs field becomes trivial as it completely disappears from the energy-momentum tensor) so the solutions become of type I. In this section, we will discuss W 0 such that W 2 0 < v 2 0 (the type I solutions will be discussed in the next section). The Yang-Mills equation (7) gives Using (13) and (15) in Eq. (10), we obtain a quadratic equation for This equation admits one or two positive roots for W 2 0 within 0, v 2 0 . Let us examine the corresponding configurations of the physical system in order.

Type II configurations
There are three cases that must be considered separately.
Option 1) Eq. (16) can have two positive roots for W 2 0 (in this case type II configurations can be divided into two sub-families, one for each positive root of Eq. (16)).
Option 2) Eq. (16) can have one positive roots for W 2 0 (in this case, there is just one family of type II configurations).
Option 3) Eq. (16) has no positive roots for W 2 0 (in this case, there is no family of type II configurations). Whether (for instance) option 1 is realized instead of options 2 or 3 depends on the values of the parameters of the models. Especially relevant are the Higgs coupling constant λ and the VEV v 2 0 . Option 1 is the most interesting one from the thermodynamical viewpoint since, in this case, there is a competition between three types of solutions: the standard meron X 0 = 1/2 and W 2 0 = v 2 0 (type I), the type II solution corresponding to the larger positive root of Eq. (16) and the type II solution corresponding to the smaller positive root of Eq. (16). As it will be discussed in the next sections, this opens the very intriguing possibility of multiple transitions between these three types of solutions. Option 2 is the second most interesting case since there is a competition between two types of solutions: the standard meron X 0 = 1/2 and W 2 0 = v 2 0 (type I), and the only viable type II solution corresponding to the unique positive root of Eq. (16). On the other hand, when option 3 is realized, no transition is possible since the only viable non-Abelian gravitating solitons belong to type I. In the discussion below, it will be convenient to introduce the following boundary values of Λ: The Eq. (16) has two different positive roots where .
when one of the following sets of conditions is satisfied: 1. When 0 < λ < 2/e 2 : 2. When 2/e 2 < λ: The equation (16) has one positive root W 2 0 ∈ 0, v 2 0 and one negative root W 2 0 ∈ − ∞, 0 when the conditions given below are satisfied: In this case, the physical solution is given in Eq. (17).
The Eq. (16) has a positive double root when the following condition is satisfied: In this case, the cosmological constant Λ and the other coupling constants are related by The Eq. (16) has no positive root within 0, v 2 0 when one of the following sets of the conditions is satisfied: 1. When 0 < λ < 2/e 2 : 2. When λ > 2/e 2 :

Equation for y(u, v) and Ricci Scalar (both type I and II)
In all of the above cases, the only non-trivial differential equation is the θθ-component of the Einstein's equations, which can be written as 1 where The solution to this equation is found to be Correspondingly, the Ricci scalar R s is given by The gauge field for both type I and type II are regular and free of singularity, since the following scalars invariant under coordinate transformation are found to be One can also find that the space-time of this solution belongs to the Petrov type D.

Thermodynamics in Nariai Coordinates
Without loss of generality, one can rescale the coordinates u and v to choose the constants C 1 = 0, C 2 = 1, and Then, the metric becomes we can express the metric (36) becomes Again, introducing new coordinates t and r given by the metric (37) can be transformed to the following metric representing a product space dS 2 × S 2 : This metric clearly shows that our space-time belongs to the Nariai family of space-times. With this coordinate system, it is straightforward to compute the surface gravity κ κ = 1 r 0 , so that the temperature T is found to be It should be noticed that (as it will be shown in the next section) the entropies of the two different families are inversely proportional to r 2 0 , so that the ratio of the entropies (which determines which family prevails) is independent of r 2 0 . In the next section, we will compute the entropy functions of type I and type II and find the range of the physical parameters where one of these type dominates the other.

Thermodynamics and Transitions between families
One of the most interesting outcomes of the present analysis is that there are different types of gravitating non-Abelian solitons in the Einstein-Yang-Mills-Higgs in the sector described above. The space-time geometry is of Nariai type, while the Yang-Mills field can be either a standard (X 0 = 1/2) meron with a trivial Higgs field or a non-standard meron with a non-trivial Higgs field. Hence, the natural question is: which one of the gravitating solitons described above will prevail? The analysis shows that the answer to this question depends in a crucial way both on the Higgs coupling and on the "VEV" of the Higgs field itself. There are two related tools to answer this question: the first is the computation of the "entropy" of the solutions, the second is the computation of the Euclidean action of the different configurations (as it will be shown, these tools give consistent answers).

Entropy function
The entropy function(al) of space-times of the form AdS 2 ×S D−2 (associated to the near-horizon geometries of extremal black holes) for an arbitrary dimension D was studied by A. Sen [40]. The entropy function of this class of space-times is the product of 2π with the Legendre transform of the Lagrangian density integrated over S D−2 . Since the electric fields associated with the gauge fields play the role of configuration variables, this entropy is a Routhian density over AdS 2 rather than a Hamiltonian density. A direct application of the Sen's method to Nariai class shows that the entropy of a Nariai space-time is "minus" Routhian density over dS 2 [41]. The first step corresponds to the on-shell evaluation of the Lagrangian of the matter field (while the second step corresponds to integrate the on-shell Lagrangian for the matter field over S 2 ). The on-shell Lagrangian reads and the Ricci scalar was given by (35). In the Nariai coordinates, we have sin θ .
Since the system has no electric field, the Routhian density H becomes where L is the Lagrangian of the system. The computations given in [40,41] tells us that the entropy function S will be given by where H max is the maximum value with respect to F 0 y 0 and F 0 L 2 (note that H max is negative so that the entropy is positive definite as it should). The entropy function of the system is found to be

Euclidean action
In order to compute the Euclidean action for both families, it is convenient to define the "sizes" of the dS 2 and S 2 as respectively. Then, the Euclidean action can be written as The on-shell conditions that extremize the Euclidean action are found to be One can easily check that these equations are equivalent to the equations of motion (6), (7), (10), and (33). Thus, it follows from (41) that the on-shell Euclidean action is precisely equal to the entropy of the system. The overall factor 2π should be understood as the circumference of the imaginary time.

Useful Plots
The previous analysis showed that a critical parameter to determine which of the solutions is thermodynamically favored is the VEV of the Higgs field v 2 0 . Here below, we will include the plots which clarify the comparisons between the two families of gravitating solitons in the three options defined in the previous sections (depending on the number of roots in the equation for W 2 0 ).

Option 1 Plots
In the "Option 1 case" defined in the previous sections, Eq. (16) has two different positive real roots (let us call them Hence, using the results from the previous subsection, we get where Here we plot together S I , S + II and S − II as function of x in three different cases (case 1 : e = 100, λ = 1; case 2 : e = 1, λ = 1; case 3 : e = 1, λ = 100) in order to show the differences between the situations in which the Higgs and Yang-Mills coupling are equal, and one small and one large. In these three plots, both Λ and κ will be fixed to 1. As one can see from these figures, the most preferred configuration is determined in a sensitive way depending on the physical parameters. For the parameters used in these three figures, the most favorable configuration is the standard meron. The last figure shows that S − II wins S + II for some certain values of the parameters.

Option 2 Plots
In the "Option 2 case" defined in the previous sections, Eq. (16) has one positive real roots (let us call it just W 2 0 ). In this case, transitions may appear as the thermodynamics is determined by the comparison of two solutions: the type I and the type II with root W 2 0 . These two solutions will be characterized by their own entropy (Euclidean actions): let us call S I and S II the entropy of the type I solution and of the type II solution respectively. As in the option 1 case, S I and S II depends on all the parameters of the model, but we will emphasize the dependence on v 2 0 . Using the results from the previous subsection, we get the entropy function of the type II with a double root S (D) The entropy function for the case with one positive and one negative roots equals to S + II . Also in this case, we plot together S I and S II as function of x in two different cases (case 1 : e = 1, λ = 1; case 2 : e = 1, λ = 100) in order to show the differences between the situations in which the Higgs and Yang-Mills coupling are equal, and λ is larger than e. In these three plots, r 0 , Λ, and κ will be fixed to 1. It should be emphasized that for a given set of the parameters {e, λ, κ, Λ, r 0 }, the most preferred configuration can change as the square of the VEV varies. Interestingly, such a configuration changes continuously in Fig. (3a) whereas it changes in a discontinuous way, as can be seen in Fig. (3b) and (3c). and λ = 100, the interval 0.29 < x < 1.71 is forbidden. In Fig. 3a, the type I is favored for x < 1.359 and the type II is favored for x > 1.359. In Fig. 3b, the type II always prevails, and in Fig. 3c, the type I is always preferred.

Spin from Isospin and 4DQHE
The analysis of the Klein-Gordon equation reveals some crucial differences between the present non-Abelian versions of Nariai space-times and the well known Nariai space-times, which are solutions of the Einstein-Maxwell field equations [43,44].
Let us begin this section with a very short review of the spin-from-isospin effect [45] [46] [47]. Roughly speaking, in the case of topologically non-trivial configurations that are not spherically symmetric but whose energy-momentum tensor is spherically symmetric the lack of spherical symmetry under spatial rotation is compensated by an internal transformation. This leads to a modification of the definition of the angular momentum, which in the usual example of non-Abelian SU (2) monopoles reads: where − → L is the orbital angular momentum and − → τ are the SU (2) generators entering in the ansatz of the gauge field.
The new term in the definition of − → J is exactly related to the infinitesimal internal rotation needed to compensate the lack of spherical symmetry under spatial rotation. As it has been discussed in the original references [ where 2D is the two-dimensional D'Alembertian in the t and r coordinates and c 0 is a constant which depends on the representation of the test scalar field. Hence, Eq. (53) explains in a very simple way the effect of the "need to compensate" the lack of spherical symmetry with an internal rotation: the centrifugal barrier (which is the term that goes as 1/r 2 for large r) is modified. Obviously, it is precisely from the centrifugal-barrier like term that one usually reads the spin of the "test fields". A quite trivial (but useful as we will now show) observation is the following: the 1/r 2 factor which multiplies − → J · − → J arises because of the r 2 term in front of the two-sphere in the metric in Eq. (52).

Now, we are ready for an important question:
What happens if in the spherically symmetric space-time of our interest (sourced by a non-Abelian soliton) in front of the two-sphere we have just a constant instead of the coordinate-dependent factor r 2 ? What happens to the spin-from-isospin effect?
Eq. (53) suggests an intriguing answer: the (modified) centrifugal barrier becomes a term that does not depend on r at all: such a term possesses discrete degenerate energy levels (as it is proportional to − → J · − → J ) with an energy gap proportional to the (homogeneous) magnetic flux. In other words, when in front of the two-sphere we have a constant instead of r 2 the "spin-from-isospin" term is replaced by the typical Hamiltonian which is used to describe the QHE in higher dimensions 4 (see in particular, [30] [31]). We will now show that this is indeed the case, and that the degeneracy of the energy levels of the non-Abelian Klein-Gordon equation changes dramatically when passing from the type I to the type II solutions: such an effect is a genuine non-Abelian fingerprint of the present families of gravitating solitons.
Let us consider scalar test fields Ψ a (a = 1, 2), which transforms like a two-component vector, in our background space-time. The field equation is given by where m is the mass of the scalar fields. Note that the Yang-Mills field associated with our solution satisfies The Klein-Gordon equation can be written as Explicitly, it can be written as whereL =L i τ i for the standard angular momentum operatorsL i , and ∆ (u,v) is the 2 dimensional D'Alembertian operator given by The eingenfunction of this operator has the form of with the eigenvalue −2ab/F 0 y 0 , for arbitrary constants a and b, and a constant doublet Ψ 0 .

Differences between type I and type II non-Abelian solitons
In the case of the type I non-Abelian gravitating solitons one can compute the sum of the orbital and isospin angular momenta through a standard procedure. In particular, the eigenvalues of the spin-orbit coupling can be obtained by which gives the eigenvalues of the coupling operator L · τ Thus, the part of the Hamiltonian which results in a spin-isospin effect has the same degeneracy of the physical system as the usual spin-orbit couplings when X 0 = 1/2. It gives a very similar behavior 5 of the Hamiltonian of the four dimensional quantum hall effect proposed in [30], Here, F 0 L 2 , the denominator of (56), plays the same role as the radius of the orbit in the system with QHE.
On the other hand, in the gravitating solitons of type II (In the case with X 0 = 1/2 and a non-trivial Higgs field) the (would be) total angular momentum − → J becomes with a real coefficient X 0 different from 1/2 (as, generically, X 0 is not even a rational number: see Eqs. (13) and (16)). The most dramatic effect manifests itself in the degeneracy of the energy level associated to the operator (56). This can be understood easily looking at the standard manipulations in Eq. (55) in the case in which X 0 is a generic real number: Although one may still hope to find a rigorous definition of the "total angular momentum − → J " with X 0 a real number, it is clear that one should not expect that the spectrum of the above "spin-from-isospin operator" H spin-isospin in Eq. (56) is still related to the combination j (j + 1) , The reason is that Eq. (57) suggests that the eigenvalues of H spin-isospin depend on l (related to the eigenvalues of L) and s (related to the eigenvalues of τ ) separately (and therefore the degeneracy is expected to be lower than in the case with X 0 = 1/2 when the eigenvalues only depend on j). Consequently, the non-Abelian Klein-Gordon equation associated to the type II gravitating merons should have energy levels with a different degeneracy than in the case of the type I gravitating solitons.
This has the following potential consequence. A multi-Fermionic system (charged under the gauge group) living in the type I gravitating solitons (in the approximation in which these Fermions can be considered as test fields) would be subject to a Hamiltonian with many of the features of the 4D QHE (as it has been explained here above). The same multi-Fermionic system would perceive a Hamiltonian with very different degeneracies in the type II gravitating solitons. Therefore, if there is a semiclassical transition from one family to the other 6 the multi-Fermionic system would suddenly be subject to a different QHE-like Hamiltonian with completely different degeneracies. We hope to come back on the fascinating physical properties of these scenarios in a future publication.

Conclusions
In this article, we constructed the first two analytic families of globally regular non-Abelian gravitating solitons in the Einstein-Yang-Mills-Higgs theory in (3+1)-dimensions with the Higgs field in the adjoint representation (however, the case in which the Higgs field is in the fundamental is very similar). The space-time metric is of Nariai type in both cases. The Yang-Mills fields are of meron type (namely, proportional to a pure gauge: A µ = X 0 U −1 ∂ µ U for some parameter X 0 ) for both families. On the other hand, while in the first family (called type I in the main text) of non-Abelian gravitating soliton X 0 = 1/2 (as in all the known examples of merons available so far) and the Higgs field is trivial, in the second family (type II) X 0 = 1/2 and the Higgs field is non-trivial (to the best of the authors knowledge, this is the first example of meron solutions with X 0 = 1/2). We have compared these two families of globally regular gravitating solitons by computing the Euclidean action of both types. This allows to determine when type II solitons (with a non-trivial Higgs and X 0 = 1/2) are favored over type I solitons and viceversa. It turns out that the most favored configuration is determined in a sensitive way depending on the parameters of the model. Even for a given set of the parameters other than x = v 2 0 , the most preferred configuration changes continuously or discontinuously as x varies. In order to disclose the differences between type I and type II gravitating solitons we analyzed the non-Abelian Klein-Gordon equation for a test scalar field minimally coupled to the non-Abelian fields sourcing the gravitating solitons themselves. The Klein-Gordon equation is able to detect very neatly the difference between type I and type II solitons (despite the fact that the space-time metric is similar in both cases). The Klein-Gordon equation can be written as the sum of a two-dimensional D'Alambert operator plus one of the Hamiltonians which has been proposed in the literature to describe the four-dimensional Quantum Hall Effect (QHE): the difference between type I and type II solutions manifests itself in a difference between the degeneracies of the corresponding energy levels.
This opens the very intriguing perspective to analyze the "many-body" wave functions of multi-Fermionic systems minimally coupled to the regular meronic fields of type I and type II solutions (in the test field limit in which these Fermions do not modify substantially the space-time metric). The idea would be to distinguish type I and type II solutions by looking at the degeneracy of the corresponding Landau Levels. We will return to this very interesting issue in a future publication.