Analytic self-gravitating 4-Baryons, traversable NUT-AdS wormholes, flat space-time multi-Skyrmions at finite volume and a novel transition in the SU(3)-Skyrme model

We construct the first analytic self-gravitating Skyrmions with higher Baryon charge in four dimensions for the SU(3)-Skyrme–Einstein-Λ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda $$\end{document} theory by combining the generalized hedgehog ansatz with the approach developed by Balachandran et al. to describe the first (numerical) example of a non-embedded solution. These are genuine SU(3) analytic solutions instead of trivial embeddings of SU(2) into SU(3) and its geometry is that of a Bianchi IX Universe. The Skyrme ansatz is chosen in such a way that the Skyrme field equations are identically satisfied in the sector with Baryon charge 4. The field equations reduce to a dynamical system for the three Bianchi IX scale factors. Particular solutions are explicitly analyzed. Traversable wormholes with NUT-AdS asymptotics supported by a topologically non-trivial SU(3)-sigma soliton are also constructed. The self-gravitating solutions admit also a suitable flat limit giving rise to Skyrmions of charge 4 confined in a box of finite volume maintaining the integrability of the SU(3) Skyrme field equations. This formalism discloses a novel transition at finite Baryon density arising from the competition between embedded and non-embedded solutions in which the non-embedded solutions prevail at high density while are suppressed at low densities.

Many of the important results in both the Skyrme model and the Einstein-Skyrme system have been derived numerically due to the highly non-linear character of the field equations and, until very recently, there were no analytic solutions with Baryon charge in these models. Analytic results are not just of academic interest. For instance, it was known that large isospin chemical potentials lead to Skyrmion instabilities on flat space. However, it is just lately that this critical behavior in the chemical potential has been fully understood thanks to an analytic formula [27][28][29][30][31]. Although the most relevant case corresponds to the SU (3) group, many of the theoretical and numerical works have been performed for the SU (2) case, which is already quite difficult in itself but not as much as the former one that not only exhibits five more generators infinitesimally, but also a non-constant curvature group manifold. However, since genuine features of the SU (3) Skyrme model could have very important physical consequences it is important to pursue any hint contributing to its understanding. In the seminal works [6,32,33] the first numerical example of a non-embedded solution was constructed. This is a genuine feature of the SU (3) Skyrme model since the solution has spherical symmetry and at the same time Baryon charge equal to 2; in contrast to the SU (2) Skyrme model where the only stable solution with spherical symmetry, in the hedgehog sense, has unit Baryon charge. These pioneering ideas have been generalized in [34][35][36][37]. Remarkably enough, in the present paper we provide the first analytic examples of these types. Very recently, following the techniques developed in [38][39][40][41][42][43][44][45][46], the first analytic self-gravitating SU (2) Skyrmions have been constructed in [47]. Here it will be shown that the configurations found in [47] can be generalized to genuine self-gravitating and topologically non-trivial SU (3) configurations. The resulting ansatz is based on the SO(3) subgroup of SU (3) that has been introduced in [6,32,33]. Moreover, we can define a suitable flat limit of these self-gravitating configurations in which they remain at the same time topologically non-trivial and integrable. Such a limit corresponds to confining the 4-Baryon (our analytic configuration with topological charge equals to 4) in a finite volume. This construction discloses a novel finite density transition: the non-embedded solutions prevail at high Baryon density while the trivially-embedded solutions prevail at low Baryon density.
The paper is organized as follows: in the Sect. 2, the Einstein-Skyrme model is introduced. In Sect. 3, the gravitating 4-Baryon is constructed. In Sect. 4, the NUT-AdS wormhole is described. In Sect. 5, 4-Baryons living within a finite flat spatial volume are introduced. In the Sect. 6, the possibility of a phase transition arising from the competition between embedded and non-embedded solutions is discussed. Eventually, some conclusions and perspectives are presented.

The Einstein-Skyrme theory
The Einstein-Skyrme system in presence of a cosmological constant is described by the action where R is the Ricci scalar, , is the cosmological constant, κ the gravitational constant, and the positive couplings K and λ are fixed by experimental data. In our convention c =h = 1 and Greek indices run over the four dimensional space-time with mostly plus signature.
The complete Einstein-Skyrme field equations read (2) where G μν is the Einstein tensor, and the Skyrme energymomentum tensor is defined by The winding number for a given solution is given by When the topological density ρ B is integrated on a space-like surface, B represents the Baryon number of the configuration.

SU(3) self-gravitating 4-Baryon
Before presenting the new results, first we will shortly describe the trivial embedding into SU (3) of the SU (2) selfgravitating solution found in [47].

The SU (2) embedded self-gravitating Skyrmion
The SU (2) generalized hedgehog ansatz reads where n 1 = sin cos , n 2 = sin sin , n 3 = cos and λ j j=1,..., 8 are the Gell-Mann matrices. Notice that {λ 1 , λ 2 , λ 3 } generate the SU (2) subgroup of SU (3). The above scalar functions are chosen in [47] as while the metric is with the range of coordinates uniquely fixed by requiring the regularity of the metric. With this ansatz one can verify that the Skyrme field equations are identically satisfied, while the Einstein equations reduce to [47] Using Eqs. (4) and (8) the Baryon number of this configuration turns out to be B = 1.

The new SU (3) non-embedded self-gravitating 4-Baryon
To construct a self-gravitating solution with higher Baryonic charge we use the remarkable ansatz introduced in [6,32,33] for a diBaryon in flat space-time. This ansatz is constructed with the subalgebra of the Gell-Mann matrices generating T = ·n ,n = (sin cos , sin sin , cos ) , In Refs. [6,32,33] ψ and χ are the radial profiles of the spherically symmetric SU (3) Skyrmion, while and are chosen as the spherical angles. Here, we promote the four functions to new profiles depending on the coordinates of a generally curved space-time, which must be determined by solving the Skyrme field equations.
In the flat case, the Baryon charge can be different from zero even if the profile ψ vanishes [6,32,33] but the field equations require ψ to be non-trivial. On the other hand, in the Einstein-Skyrme case one can take ψ = 0 for the curved space-times we consider as it will be now discussed.
At this point, we need an educated ansatz for the four functions appearing in the above ansatz. A slight modification of the generalized hedgehog ansatz of [47] defined by ψ = 0 and does the job. Indeed, the SU (3) Skyrme configuration (10)- (12) identically satisfies the Skyrme field equations on any metric of the form (7)-(8).
Thus, Einstein equations with the energy-momentum tensor in Eq. (3) corresponding to the Skyrme configurations defined in Eqs. (10)-( 12) reduce tȯ One can write down the most general solution of Eqs. (13) following the analysis of [44]. The above Eq. (13) in the SU (3) case are very similar to the SU (2) one in Eq. (9), but one can notice that the coefficients appearing are different; in fact, the SU (3) energy-momentum contribution is four times that of SU (2). This difference also appears when computing the topological charge density ρ of the SU (3) ansatz, defined in Eqs. (10)- ( 12), and the corresponding ρ SU(2) B of the SU (2) ansatz, defined in Eqs. (5) and (6), Taking into account Eqs. (4) and (8), the topological charge of the SU (3) configuration is B = 4. These are genuine SU (3) configurations: in SU (2) one cannot obtain configurations with topological charge 4 compatible with the metric in Eq. (7). Another interesting quantity is the ratio between the three-dimensional volume of the spatial section of the static gravitating 4-Baryons (corresponding to the static solutions of Eq. (13)) and the static gravitating Skyrmions of [47] (corresponding to the static solutions of Eq. (9)). It reveals the non-trivial strong and gravitational interactions of the system (as, in the case of four noninteracting solitons, one should expect = 4). The above regular solutions of the Einstein--SU (3)-Skyrme system defined in Eqs. (7), (8), (10), (11) and (12), are the first analytic self-gravitating Skyrmions of higher Baryonic charge in (3 + 1)-dimensions.
The fact that the Baryon charge is 4 instead of 2 as in [6,32,33], is related to the compactness of the t = const. hypersurfaces of the metric (7). Here, rather than requiring boundary conditions at spatial infinity as in [6,32,33], one has to require periodic boundary conditions for U B compatible with the compact spatial metric (7)- (8). This charge 4 arises due to the fact that, unlike what happens in [6,32,33], the present gravitating 4-Baryons must wrap around three compact spatial directions instead of two.

Bianchi IX
The above construction can be further extended to a Bianchi IX cosmology. Indeed, the SU (3) Skyrme field equations on the metric ds 2 = −dt 2 + I 2 1 (cos θ dγ + dφ) 2 + I 2 2 (cos φdθ + sin θ sin φdγ ) 2 where I j = I j (t), are still identically satisfied with the same ansatz defined in Eqs. (10), (11) and (12)! The reason behind this quite remarkable fact is that the left-invariant forms that appear in the construction of the most general Bianchi IX metric are proportional to the left-invariant forms that appear when computing the Skyrmion derivatives (11) and (12). Since the a μ characterizing the Skyrmionic configuration play, at the same time, also the role of "drei-bein" of the spatial metric, huge simplifications appear in the field equations. Hence, the complete set of coupled Einstein--SU (3)-Skyrme field equations (2) and (3) for the metric (14) with the ansatz in Eqs. (10 ), (11) and (12), reduce to a consistent dynamical system for the three Bianchi IX scale factors where we have defined I (2) = I 2 1 + I 2 2 + I 2 3 , I (4) = I 4 1 + I 4 2 + I 4 3 , I = I 2 1 I 2 2 + I 2 1 I 2 3 + I 2 2 I 2 3 , I (3) = I 1 I 2 I 3 . As the SU (3) Skyrme field equations are automatically satisfied, the above dynamical system for the scale factors I j (t) describes the dynamical evolution of the four self-gravitating Skyrmions. The analysis of such system can reveal many interesting features about the interplay between the gravitational and the strong interactions of these 4-Baryons. A comparison with the numerical gravitating Skyrmions constructed in [48] is useful. Their numerical solutions with Baryon charge 2 (which are based on the SU (2) subgroup of SU (3)) are asymptotically flat while the present analytic solutions with Baryon charge 4 (which are based on a non-embedded ansatz) either have compact spatial sections (in the Bianchi IX sector) or possess two asymptotically NUT-AdS regions (see the section below).

NUT-AdS Wormhole
In this section we will consider the limit of vanishing Skyrme coupling λ = 0. Following [47], a double-Wick rotation of the ansatz in Eq. (12) leads to ψ = 0 and are identically satisfied. A direct computation shows that the coupled Einstein--SU (3)-non-linear sigma model field equations (2) and (3) with λ = 0 are satisfied for provided < 0. Notice that the square of the NUT parameter Q 2 is four times the result for SU (2) [47]. The above metric represents a traversable Lorentzian wormhole with NUT-AdS asymptotic supported by a reasonable physical source (see [47,49] ); the non-linear sigma models do not violate energy conditions. The well-known no-go results on the existence of wormholes are avoided due to the presence of a NUT parameter. The SU (3) wormhole throat is larger than the one of the SU (2) configuration in [47,49] since the present matter field possesses a higher topological charge.
found by adapting the self-gravitating solution found in [47] to a flat space with finite volume in [27]. A similar construction also works in the present case. The sector explored below describes four low-energy Baryons confined within a finite volume in flat space. Consider the SU (3) Skyrme configuration defined by ψ = 0 and together with Eqs. (10) and (11) in a flat metric of the form where L 0 has dimension of length and represents the size of the box in which the Baryons are confined, and the dimensionless coordinates have the following ranges The topological density is again different from the SU (2) case since Considering the boundary conditions as H (t, 0) = 0 and H (t, 2π) = ± π 2 , the SU (3) topological charge is ±4 (which is four times the charge found in [27][28][29]). Concretely, these Skyrmions confined to a finite volume can only have charges −4, 0, or 4 (while in the SU (2) case defined in [27][28][29], the corresponding Skyrmions can only have charges −1, 0, or 1). This confirms the genuine SU (3) nature of the ansatz defined in Eqs. (10), (11) and (17).
The SU (3) Skyrme field equations on the flat metric (18) for the Baryon-like ansatz (10), (11) and (17) become a single partial differential equation for a scalar a profile Hence, the eight coupled SU (3) Skyrme field equations collapse to a single integrable PDE. This opens the intriguing possibility to analyze many interesting non-trivial properties of these multi-Skyrmions confined to a finite volume using Sine-Gordon theory such as possible phase transitions between the trivial embedding of the SU (2) solutions [27] and the present non-embedded SU (3) solutions living in the same box. At this point it is important to recall that the value 4 of the Baryonic charge cames only from the nature of the SU (3) ansatz and the compactness of the space-time under consideration, both, in the present case and in the previous gravitating configurations. Using the same boundary conditions for the usual SU (2) Skyrme field is not possible to obtain B = 4. In fact, Eqs. (17) and (19) describes a bound state of 4-Baryons whose energy density is localized. In the next section will be show that indeed, the energy of the 4-Baryon state is different from the one of the four Skyrmions on the SU (2) ansatz.

A novel transition at finite Baryon density
In the present section we will show that, within the same box of finite volume, there is a competition between embedded and non-embedded configurations with the same Baryonic charge, the control parameter of the transition being the Baryon density. We show here below that at high Baryon density non-embedded solutions are favoured over embedded solutions (and vice-versa at low Baryon density).
A suitable ansatz to describe SU (3) configurations living in the flat box defined above that are trivial embedding of where p and q are non-vanishing integer numbers (see [50]) and λ 1 , λ 2 , λ 3 the first three Gell-Mann matrices. As we explained in Section V, as well as in [50], both in the 4-Baryon case (Eqs. (10) and (17)) and in the case of trivially embedded solutions (Eqs. (10) and (20)) the consistent boundary conditions for the profile H (t, r ) are The topological densities read We know that the topological charge in the SU (3) case is B SU(3) = 4, while the topological charge in the trivially embedded solutions with the above slightly generalized ansatz is where the integers p and q are the ones appearing in the ansatz in Eqs. ( 10) and (20) (see also [50]). Thus, if we choose pq = 4 , then we will have two inequivalent configurations 1 (the one in Eqs. (10) and (17) as well as the one in Eqs. (10) and (20)) with the same charge living in the same box. The natural question is: Which type of configuration is energetically favoured? Thanks to the remarkable properties of the ansatz described above, one can answer this question explicitly. The first step is to notice that in both cases (the one in Eqs. (10) and (17) as well as the one in Eqs. (10) and (20)) the SU (3) Skyrme field equations in the flat metric (18) reduce to a single partial differential equation of the form with for (17) and (20), respectively.
TABLE 1 2 shows the energy comparison between the SU (2) and the SU (3) configurations with Baryon charge B = 4, while it varies the size of the box (measured by L 0 ). We have fixed K = 2, λ = 1 (which means that we are measuring length in f m). One can see that a phase transition 1 Using the elegant arguments in [6], one can easily see that these two configurations cannot be deformed into each other using global isospin transformations. 2 The notation E SU (2) X,Y means "the energy of the SU (2) configuration in Eqs. (10) and (20) with p = X and q = Y ". In the table XY = 4 as we need to compare configurations with the same topological charge. appears at 3 L * 0 ≈ 1.4. When L 0 < L * 0 then the 4-Baryon are favoured over the trivially embedded solutions (and viceversa when L 0 > L * 0 ). The appearance of this transition is related to the fact that 4-Baryons deal more efficiently with the well known repulsive interaction between Skyrmions. From the plots of the energy densities below for the two types of configurations one can see that, at high Baryon density, it pays off that the 4-Baryon is "less peaked" around its maximum while at low Baryon densities the energy densities of both types of solutions become relatively flat and, in such cases, the trivially embedded solutions prevail. Fig. 1 shows the energy density of the three relevant configurations, for three different values of the size of the box.
At last, we make a small observation about the quantization of the solutions. Since, with the above ansatz, the hedgehog property holds (as the 8 coupled SU (3) field equations reduce to a single PDE for the profile H ). The small fluctuations of the profile H around these solutions are described by the effective action obtained replacing the ansatz itself into the original SU (3) Skyrme action ( [51,52]). Thus, the quantization of the 4-Baryon living in the finite box defined above can be analyzed using known results on the sine-Gordon theory plus the semiclassical quantization of the Isospin degrees of freedom described in [32,33]. This observation also shows that the flat configurations constructed in the manuscript are stable, at least, under the perturbations which keep the hedgehog properties. Namely, the flat solutions considered here are stable under the following type of perturbations (see the analysis in [51] and [52]): We hope to come back on the interesting but rather difficult problem of the full stability analysis (which must be analyzed numerically) in a future publication. We can see that at high Baryon density it is convenient that the 4-Baryon is "less localized" around its maximum, and at low Baryon density the energy density for both types of solutions tend to flatten

Conclusions and perspectives
The first regular analytic self-gravitating Skyrmions of higher Baryon charge have been constructed. The space-time corresponds to a general Bianchi IX cosmology whose three scales factors evolve accordingly, while the Skyrme field equations are identically satisfied in the sector with Baryon charge 4. All these solutions disclose genuine features of the SU (3) Skyrme model. Traversable wormholes with NUT-AdS asymptotic supported by an SU (3) regular solitonic solution of the resulting non-linear sigma model are also constructed. Also, a suitable flat limit in a finite volume of the self-gravitating Skyrme configurations can also be defined. In this case, the eight SU (3) Skyrme field equations become the Sine-Gordon field theory without sacrificing the higher Baryon charge. This flat limit explicitly describes charge 4 Baryons confined in a box. The present formalism discloses the existence of a novel transition between embedded and non-embedded configurations with the same Baryonic charge. The control parameter of this transition is the Baryon density: at high Baryon densities, it is energetically convenient to have non-embedded solutions while, at low Baryon densities, trivially embedded solutions prevail. To the best of authors knowledge, this phenomenon related to the competition between embedded and non-embedded solutions at finite Baryon density is new. The reason is that the present techniques are especially suitable to deal with Skyrmions and 4-Baryons at finite density as these techniques allow to determine how relevant physical properties of these solitons depend on the Baryon density. Without these informations it would have been impossible to disclose such possibility.