Ordered arrays of Baryonic tubes in the Skyrme model in (3+1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$3+1$$\end{document}) dimensions at finite density

A consistent ansatz for the Skyrme model in (3+1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$3+1$$\end{document})-dimensions which is able to reduce the complete set of Skyrme field equations to just one equation for the profile in situations in which the Baryon charge can be arbitrary large is introduced: moreover, the field equation for the profile can be solved explicitly. Such configurations describe ordered arrays of Baryonic tubes living in flat space-times at finite density. The plots of the energy density (as well as of the Baryon density) clearly show that the regions of maximal energy density have the shape of a tube: the energy density and the Baryon density depend periodically on two spatial directions while they are constant in the third spatial direction. Thus, these topologically non-trivial crystal-like solutions can be intepreted as configurations in which most of the energy density and the baryon density are concentrated within tube-shaped regions. The positions of the energy-density peaks can be computed explicitly and they manifest a clear crystalline order. A non-trivial stability test is discussed.


Introduction
A theoretical description of cold and dense nuclear matter as a function of baryon number density [1,2] is still lacking. Analytic solutions in the low-energy limit of QCD able to describe crystals-like structures in Baryonic matter at finite density would be extremely helpful in this respect. Until quite recently, the only theoretical arguments supporting the appearance of Baryonic crystals at finite density could be found in two-dimensional effective toy models (see [3,4] and references therein) as well as in the case of the so-called magnetic Skyrmions [5]. Clever approximations related to either the so-called product ansatz or to Yang-Mills instantons may lead to analytic results in (3+1)-dimensions [6][7][8][9][10][11][12][13][14][15] which help to clarify many properties of these multi-solitonic a e-mail: canfora@cecs.cl configurations. 1 In fact, it is still an open problem to build crystals of Baryons explicitly in the low-energy limit of QCD. For this reason, different schemes to describe this phase have been introduced.
A common approach (based on the AdS/CFT correspondence [16]) is in the references [17][18][19][20][21][22]. The main difficulty in these cases is to calculate the required solitons: only single baryon solutions have been discussed in details [23,24]. Numerical solutions at finite density in full Sakai-Sugimoto model are extremely difficult (see [25,26]).
Despite that, it has been argued that baryons in such a regime are necessarily in a solid crystalline phase [20,21,[27][28][29]. More complex ordered structures are also expected to appear at finite densities. It is widely accepted that the socalled nuclear pasta phase do exist (see [30][31][32][33] and references therein). In this phase, the energy-density plots in the references just mentioned show, for instance, that ordered configurations of "Baryonic tubes" (nuclear spaghetti) in which most of the Baryonic charge is gathered in tube-shaped regions are observed (nuclear lasagna and nuclear gnocchi phases have also been observed). It is usually assumed that the available theoretical tools are unsuitable to describe such complex structures in the low energy limit of (only numerical results are available in the nuclear pasta phase).
An explicit analytic approach to describe the predictions above in the context of the low energy limit of (3 + 1)dimensional QCD is the main goal of the present paper.
The appearance of crystal-like structures in the Skyrme model (which was predicted on the basis of the product ansatz introduced in [8]) is well established numerically (see [38,48] and references therein). Within the rational map approach [48] one can construct numerically configurations in which the number of "bumps" in the energy density is related with the corresponding Baryon charge. Based on the fact that, in the Skyrme model, the solution representing one spherical Skyrmion (introduced by Skyrme himself) must be found numerically, it has been always assumed that it is impossible to find analytic solutions representing ordered pattern of solitons without approximations. Consequently, interesting modifications of the original Skyrme model have been considered which do allow to find analytic solutions representing Skyrmions at finite Baryon density (see, in particular, [49][50][51][52] and references therein). However, only the original Skyrme model (which is directly related with the low energy limit of QCD) will be considered here.
The generalized hedgehog ansatz introduced in [53-58] [59][60][61] allowed the construction of the first analytic multi-Skyrmions at finite density (however, these Skyrmions at finite density do not have crystal-like structure).
In the present paper, the methods in [59][60][61] will be generalized to construct ordered arrays of Baryonic tubes with crystalline structure living at finite density. One then observes transitions in which for fixed total topological charge, configurations with taller tubes but a smaller number of tubes are energetically favored over configurations made of a lower tubes but with a higher number of tubes. This paper is organized as follows: in the second section the Skyrme model is introduced and the Skyrme field equations are written in two equivalent ways: this is very important in order to check that the present crystals are really solutions of the full Skyrme field equations. In the third section, the method to go beyond the spherical hedgehog ansatz is introduced and the analytic crystal-like Skyrmions are derived. In the fourth section the appearance of fractional Baryonic charge is discussed. In the final section, some conclusions will be drawn.

Beyond the spherical hedgehog ansatz
The action of the SU (2) Skyrme system is where K and λ are the coupling constants, 2 1 2 is the 2 × 2 identity matrix and the t j are the basis of the SU (2) generators (where the Latin index j corresponds to the group index). The three coupled Skyrme field equations can be obtained by taking the variation of the action with respect to the SU (2)-valued field U : The energy density (the 0 − 0 component of the energymomentum tensor) reads The following parametrization of the SU (2)-valued scalar U (x μ ) will be adopted where, generically, α, F and G are functions of the four space-time coordinates.
It is worth to emphasize here a relevant technical point. The number of Skyrme field equations is equal to the dimension of the corresponding Lie Algebra (in the present case the SU (2) case will be considered so that the task is to solve three coupled non-linear field equations). One can write down the complete set of Skyrme field equations in many different equivalent ways. In this paper two natural representations of the field equations will be considered. The first one is described in Eq. (2). In this matrix representation the three field equations (namely, E j = 0, j = 1, 2, 3) are obtained as variation of the action with respect to the SU (2)-valued scalar field U . An equivalent way (which will be described in the next subsection) corresponds to take the most general representation of SU (2)-valued scalar field U in the terms of three scalar degrees of freedom (α, F and G defined in Eqs. (4), (5)) and (6)) and replace such general representation of U into the action (which will become an action for the three interacting scalar fields α, F and G). In this second representation, the three field equations (which will be derived explicitly in the next subsection) will be simply Obviously, the two representations are equivalent (in particular, it will be shown explicitly in both representations that the analytic crystals constructed in the present paper are solutions of the complete set of Skyrme field equations). The only warning is that, in general, the field equations in one representation are linear combinations of the field equations in the second representation. The most difficult task is to find a good ansatz which keeps alive the non-trivial topological charge and, at the same time, allows for a crystal-like structure in the energy density without making the field equations impossible to solve analytically. It is worth to note that with the usual hedgehog ansatz in [34][35][36] (as well as in its finite density generalization [59,60]) the energy density in Eq. (3) (due to the trace on the internal indices) only depends on the Skyrmion profile. Until now, this fact prevented one from describing analytically multi-solitonic configurations with crystal-like order. It is also interesting to note that the main building blocks of these ordered multi-solitonic structures are not necessarily single nucleons. Indeed, it is widely accepted that the so-called nuclear pasta phase may appear (see [30][31][32][33] and references therein). In this phase one can have, for instance, ordered arrays of "Baryonic tubes" (nuclear spaghetti) in which the shapes of the regions of maximal energy density and Baryon density are tubes. To describe such configurations one necessarily needs an ansatz in which both the Baryon density and the energy density depend not only on the profile α but also on an extra spatial coordinate.
The Baryon charge of the configuration reads In terms of α, F and G, the topological density ρ B reads so that a necessary condition in order to have non-trivial topological charge is From the geometrical point of view the above condition (which simply states that α, F and G must be three independent functions) can be interpreted as saying that such three functions "fill a three-dimensional spatial volume" at least locally. In other words, dα, d F and dG can be used as 3D volume form. Hence, the condition in Eq. (10) ensures that the configuration one is interested in describes a genuine three-dimensional object. On the other hand, such a condition is not sufficient in general. One has also to require that the spatial integral of ρ B must be a non-vanishing integer: Usually, this second requirement allows to fix some of the parameters of the ansatz. From now on, as it is customary in the literature, the terms Skyrmionic crystals (and also crystals of solitons) will refer to smooth regular solutions of Eq. (2) with the properties that both the topological charge (defined in Eqs. (7), (8), (9) and (10)) is non-vanishing and larger than 1 and, moreover, that the local maxima of energy density and Baryon density manifest a crystal-like ordered structure.

Explicit parametrization
It is useful to write down the Skyrme field equations explicitly in terms of the three scalar degrees of freedom α, F and G. In this way, it is possible to check directly that the novel ansatz proposed in the present paper is consistent. It is convenient to introduce the following functions Y 0 (x μ ) and In this way, the most general element U of SU (2) can be written as These functions are useful to write down the Skyrme action explicitly in terms of α, F and G. Introducing the tensor the Skyrme action is then defined as Now, we are in the position to write down the general Skyrme field equations in terms of α, F and G (which obviously are equivalent to the Skyrme field equations in matrix form in Eq. (2)). The variation of the Skyrme action with respect to α leads to the equation of motion The variation of the Skyrme action with respect to F leads to the equation of motion The variation of the Skyrme action with respect to G leads to the equation of motion

Example: the original Skyrme ansatz
Here we give a simple and well known example showing that the above equations (16), (17) and (18) are suitable to devise a strategy to find good ansatz which reduce the full Skyrme field equations to only one consistent equation for the profile in a non-trivial topological sector. The original Skyrme ansatz is defined by the choice in the flat metric in spherical coordinates: One can check that, due to the identities and to the fact that both F and G are linear functions in the chosen coordinates system, the three Skyrme field equations (16), (17) and (18) (16), (17) and (18) one can see directly that Eqs. (17) and (18) are identically satisfied and that Eq. (16) reduces to the usual equation for the Skyrme profile α s (R) which can be found in all the textbooks (see, for instance, [38]). In the following, a different trick to get a consistent equation for α in a topologically non-trivial sector will be introduced.

Comparison with the rational map ansatz
The idea of the rational map approach [48] is to obtain an approximate description of multi-Skyrmionic configurations replacing the isospin vector n i for the spherical hedgehog defined in Eqs. (5), (19) and (20) by a more general map between two spheres. Even if the rational map ansatz is only an approximation, it provides one with a very clear intuitive picture of multi-Skyrmionic configurations. The usual starting point is where z is a complex coordinate which, using the stereographic projection, can be identified with the coordinates on the 2-sphere: The rational map (z) (which gives the name to the whole approach) reads where p and q are polynomials in z with no common factor. The degree N of the rational map is defined as It is well-known that where n p and n q are the degrees of p (z) and q (z) respectively. Correspondingly, the Baryon number B of the configuration in Eqs. (22), (23) and (24) reads where The winding number is the product of the contribution coming from the profile function α times the degree N of the rational map . One can think that n is the number of "bumps" associated to α while N are related to the bumps in the directions orthogonal to α. When (and only when) the rational map ansatz in Eqs. (22), (23) and (24) reduces to the original spherical hedgehog ansatz discussed in the previous subsection and it gives rise to a (numerical) solution of the complete set of Skyrme field equations. On the other hand, in all the higher degrees cases N > 1 the rational map solutions are only approximated solutions of the complete set of Skyrme field equations. The strategy 3 of this approach is to minimize the total energy with respect to both α and the rational map . This framework is tied to situations in which there is an approximated spherical symmetry: in particular, the crystals-like configurations typical of nuclear pasta (see [30][31][32][33] and references therein) cannot be analyzed in this approach. The main advantage of this rational-map framework is that it disentangles the radial coordinate R from the angular coordinates. The main disadvantage is that, as already emphasized, such a disentangling is only an approximation.
In terms of the explicit parametrization introduced in Eqs. (12) and (13) (in which the complete set of Skyrme field equations are (16), (17) and (18)) one can interpret the rational map approach as follows. The first step of the procedure (in which one minimizes the energy functional with respect to the rational map) gives rise to an effective equation for the radial profile which is very similar to the original one of Skyrme with the difference that some parameters are rescaled (see, in particular, the discussion in [58] and references therein). Thus, in a sense, one can think that in the rational-map approach Eqs. (17) and (18) are solved "on average" (freezing α) and then the corresponding Eq. (16) for α is solved (more often than not) numerically.
Instead, in the present paper it will be introduced an ansatz with three main advantages.
Firstly, the complete set of Skyrme field equations will be solved exactly (and not only approximately). In other words, the ansatz described in the next sections will solve Eqs. (16), (17) and (18) without any approximation.
Secondly, the solutions will be analytical and not numerical.
Thirdly, the ansatz is not tied to situations with spherical symmetry.

Skyrme crystals
The main physical motivation of the present work is to study finite density effects. The easiest way to take into account finite-density effects is to introduce the following flat metric where 4π 3 L A is the volume of the box in which the Skyrmions configurations to be constructed below will live. The adimensional spatial coordinates r , θ and φ have the range Following the strategy of [59,60], the boundary conditions in the θ direction here will be chosen to be Dirichlet while in the r and φ directions they can be both periodic and antiperiodic. These are not the only allowed boundary conditions: however this is the simplest choice as it will be discussed in a moment.
The ansatz for the Skyrme crystal is The periodic time dependence 4 in Eq. (31) is natural for three very good reasons. Firstly, the most popular way to avoid Derrick's famous no-go theorem on the existence of solitons in non-linear scalar field theories corresponds to search for a time-periodic ansatz such that the energy density of the configuration is still static, as it happens for boson stars [62] (in the simpler case of U (1)-charged scalar field: see [63] and references therein). The present ansatz defined in Eqs. (4), (5) and (31) has exactly this property. Indeed, the energy density does not depend on time (see Eqs. (44), (45) and (46) below) so that these configurations possess a static energy density.
Secondly, unlike what happens for the usual Bosons star ansatz for U (1)-charged scalar fields, the present ansatz for SU (2)-valued scalar field also possesses a non-trivial topological charge.
Thirdly, with the above choice of the ansatz the complete set of Skyrme field equations simplifies dramatically (as it will be shown below) due to the fact that, with the function G in Eq. (31), ∇ μ G∇ μ G = 0.
With the ansatz defined in Eqs. (4), (5) and (31) the topological density in Eq. (8) reads The above ansatz satisfies the first condition to be topologically non-trivial in Eq. (10) as the three function α, F and G in Eq. (31) are independent and "fill a three-dimensional space" (see the discussion about Eq. (10)). As far as the boundary conditions are concerned, one can check directly that if one chooses, for instance, periodic boundary conditions in θ (taking q to be an even number in Eq. (31)) the topological charge vanishes (even if the topological density does not) so that the second condition in Eq. (11) would be violated. Hence, the simplest choice to get topologically non-trivial configurations is to take q to be an odd integer. I hope to come back on the analysis of other possible boundary conditions in a future publication. The boundary condition on α and the corresponding Baryon charge are Thus, the present ansatz cannot be deformed continuously to the trivial vacuum U 0 = 1 2x2 since the third homotopy class (which represents the Baryon charge) is non-vanishing.

Tube-shaped regions
In order to understand the shape of the regions of maximal Baryon density (which will coincide with the shape of the regions of maximal energy density: see Eqs. (44), (45) and (46) below) one can proceed as follows. The baryon density ρ B in Eqs. (8) and (32) is periodic in r and θ while it is constant in the φ-direction. 5 One observes that the positions of the maxima of the Baryon density are defined by the conditions which define a periodic 2-dimensional lattice in r and θ . The contour-plots in Fig. 1 for the energy-density have a similar pattern. Since both the Baryon density and the energy density are constant in the third spatial direction φ, a threedimensional contour plot 6 would show that the shapes of the regions of maximal energy density (and similarly for the Baryon density) are tubes of length 2π L.
In other words, in order to understand the three-dimensional shape of the regions of maximal energy and Baryon density one has simply to take the clear-yellow spots in Fig. 1 and move them along the φ-direction (which is orthogonal to the r − θ plane of Fig. 1). Hence, such regions are (ordered arrays of) parallel tubes. The "φ-constant" sections of these tubes correspond to the contour plots in Fig. 1.

The explicit solution
When one plugs the ansatz in Eq. (31) into the complete set of three coupled Skyrme field equations in Eq. (2), they reduce to only one integrable equation for α(r ): 5 The fact that the Baryon density is periodic both in r and in θ can be seen taking into account the explicit solution for α in Eqs. (37), (38) and (39). The same is true for the energy-density: see Eqs. (44), (45) and (46) below. 6 A three-dimensional contour plot can be easily done since all the present solutions are known in explicit analytic form. However, I have been unable to reduce such 3D contour plots to a reasonable size (less than 10 MB). On the other hand, the two-dimensional contour plots should be enough to understand the shapes of the regions of maximal energy density.
where E 0 is an integration constant to be fixed requiring the boundary condition 7 in Eq. (33). In other words, with the ansatz in Eq. (31) all the three Skyrme field equations in Eq.
(2) E j = 0 have the common factor P [α] defined in Eq. (36) so that if P [α] = 0 then all the Skyrme field equations are satisfied. This is the reason why the ansatz in Eq. (31) has been chosen. An alternative way to see that the above ansatz reduce the complete set of Skyrme field equations to just one equation for the profile α keeping alive the topological charge is to use the explicit parametrization of the Skyrme action in term of α, F and G in Eqs. (12 ), (14) and (15). In this way, one can read directly the field equations for the three functions α, F and G. The relevant properties of the ansatz in Eq. (31) are If one takes into account Eqs. (40), (41) and (42) together with the fact that both F and G depend linearly on the coordinates defined in Eqs. (29) and (30), one can see easily that Eqs. (17) and (18) 8 in Eq. (41). All in all, the three coupled Skyrme field equations Eq. (2) with the ansatz in Eq. (31) reduce to a simple quadrature which can be integrated using elliptic functions. The boundary condition in Eq. (33) reduces to: The above equation for E 0 always has a positive real solution 9 . Moreover, one can see that ∂ r α > 0 and that, when n is large, both η (α, E 0 ) and E 0 are of order n. The conclusion is that the ansatz in Eqs. (4), (5), (6) and (31) in the flat metric in Eqs. (29) and (30) in which the profile α is given in closed form in Eqs. (37), (38), (39) and 7 The positive sign of the square root of (∂ r α) 2 has been chosen. 8 It is worth to note that in the usual case of the spherical hedgehog introduced by Skyrme himself (see Eqs. (19) and (20)) ∇ μ G∇ μ G = 0. 9 The left hand side of Eq. (43) as function of E 0 increases from very small values (when E 0 is very large and positive) to very large values (when E 0 is close to zero but positive). Thus, there is always a value of E 0 which satisfies Eq. (43).
(43) gives rise to exact and topologically non-trivial solutions (with Baryonic charge np) of the complete set of Skyrme field equations 10 for any integers n and p and for any odd integer q.
Firstly, in that references, the factor sin 2H (where H is the profile defined in Eqs. (9) and (10) of [59]) appears linearly in the Baryon density in Eq. (16) of [59]. Thus, in that references, it is not possible to increase the Baryon charge by increasing the number of "bumps" in the profile H .
Secondly, the energy density (defined in Eq. (15) of [59]) only depends on H (and, consequently, it only depends on one spatial coordinate). This prevents one from describing explicitly crystal-like structures in which the number of bumps in the energy-density is related with the Baryon number (as, in that references, only one bump in H is allowed).
In the present case both problems are solved since in the Baryon density in Eq. (32) sin α appears quadratically (so that by increasing the number of "bumps" in α one also increases the Baryon charge) and the energy density in Eqs. (44), (45) and (46) (even after taking the trace over the SU (2) indices) depends non-trivially both on the profile α and on the spatial coordinate θ (so that a clear crystal-like structure in its peaks emerges).
The two contour plots of the energy-density 11 in the r −θplane below in Fig. 1 (with p = 1, q = 1 and n = 25 on the left and p = 1, q = 3 and n = 25 on the right) show the crystal-like pattern of the bumps (whose positions can be determined explicitly maximizing the energy-density): p = 1, n = 25 and q = 1 p = 1, n = 25 and q = 3 (b) (a) Fig. 1 The left panel (a) shows the contour plot of the energy density with p = 1, n = 25 and q = 1 (where the highest energy density is in clear yellow while the lowest is blue). The right panel (b) shows the contour plot of the energy with p = 1, n = 25 and q = 3 (in both cases it has been assumed that K = 2 and λ = 1; note that the horizontal axis corresponds to the r direction while the vertical axis to the θ direction).
The patterns of the peaks clearly show a crystal-like structure. It is clear to see that, when passing from q = 1 to higher values of q (keeping fixed n and p), the bumps in the energy density (which, when q = 1, carry one unit of Baryon charge) are divided into q smaller bumps each of which carries 1/q of Baryon charge. a p = 1, n = 25 and q = 1. b p = 1, n = 25 and q = 3 These configurations are ordered arrays of Baryonic tubes (namely, tubes in which most of the energy density and Baryonic charge is concentrated within tube-shaped regions). It is easy to recognize that this is the case since both the energy density and the Baryon density do not depend on the third spatial direction φ (while they are periodic both in r and in θ ). A three-dimensional contour plot would reveal that the clear-yellow regions in Fig. 1 of maximal energy density would become, in 3D, clear-yellow tubes. The positions of these Baryonic tubes in the r − θ plane (defined by the clearyellow spots in Fig. 1) manifest a clear crystalline order. The similarity with the spaghetti-like configurations found (numerically) in the nuclear pasta phase (see the plots in [30][31][32][33] and references therein) is quite amazing.
Consequently, the interpretation of the above results together with the plots in Fig. 1 tell that, when q = 1, the integer n corresponds to the number of Baryonic tubes in the configurations while p is the Baryonic charge of each tube. The role of the integer q will be analyzed in the next section.
The total energy E tot of the configuration can be obtained in a closed form (using Eqs. (38) and (39)). One is left with an integral of an explicitly known function of α which contains all the relevant informations: where B = np is the Baryon number, n is the number of bumps associated with the profile α in the r direction. The above equation represents the explicit expression for the total energy as a function of the parameters of the system. As it was already emphasized in the discussion below Eq. (34), the total energy does not depend on p and L separately but only on the ratio p/L. It is convenient from the energetic point of view neither to have very small p/L nor to have very large p/L. This can be seen as follows: let us denote the variable p/L as ξ . Then, the total energy per-Baryon (namely, E tot /B from Eqs. (47) and (48)) as function of ξ (assuming that n is large so that where c 1 and c 0 (which can be read from the explicit expressions in Eqs. (47) and (48)) depend on all the other parameters of the theory. Thus, the total energy per-Baryon as function of ξ has a non-trivial local minimum: This fact means that if one increases the value of p one also has to increase the length of the tube L in such a way to keep p/L equal to its optimal value in the above equation. Therefore one can interpret p as the length of the Baryonic tubes in unit of f m. A further natural question is the following: given a total Baryon number B, is it energetically more convenient to have higher p and lower n or the other way around?
The answer to this question in the low energy limit of QCD can be found just by analyzing the function (α; p, q, n) defined in Eq. (48) which contains all the relevant informations. For instance, can see explicitly that if A decreases it is better to have larger values for p.
In other words, if one keeps B = np fixed and A decreases, then it is better to have taller tubes (and a smaller number of tubes to keep B fixed).
On the other hand, if one increases A it is better to have a bigger number of smaller tubes (to keep B fixed). Although a plot of the total energy (or total energy per-Baryon as function of all the parameters is not very useful), if one keeps some of the parameters fixed and only considers pairs of related parameters (such as p and L or A and n) the pattern is always similar to the one in Eqs. (49) and (50). I hope to come back on the properties of the function (α; p, q, n) defined in Eq. (48) in a future publication.

A remark on the stability
A remark on the stability of the above crystals is in order. In many situations, when the hedgehog property holds (so that the field equations reduce to a single equation for the profile) the most dangerous perturbations are perturbations of the profile which keep the structure of the ansatz (see [64,65] and references therein). In the present case these are It is a direct computation to show that the linearized version of Eq. (36) around a background solution α 0 (r ) of charge B = np always has the following zero-mode: u (r ) = ∂ r α 0 (r ). Due to Eqs. (38), (39) and (43) u(r ) has no node so that it must be the perturbation with lowest energy. Thus, the present solutions are stable under the above perturbations. It is also worth to remark that isospin modes U A = A −1 U 0 A (where U 0 is the solution of interest and A is a generic SU (2) matrix which only depends on time) have positive energies (if the energy of U 0 is positive, as in the present case) [42]. The effective action for these modes is the one of a spinning top and the energy of the corresponding perturbations is positive definite. Isospin modes are "transverse" to the perturbations of the profile in Eq. (52) as they do not touch α(r ): the reason is that both U A = A −1 U 0 A and U 0 have the same profile. This stability argument holds for any n, p and q in Eqs. (4), (5), (31) and (33). It is worth to note that in the original spherical ansatz of Skyrme with radial profile α s discussed in the previous sections, the topological density is also proportional to sin 2 α s so that one can increase the winding increasing the "spherical bumps" in the energy density (which in the case of the original spherical ansatz of Skyrme only depends on the radius). As it is well known, one can construct (numerically) these "higher charges spherical Skyrmions" but all of them are unstable. The instability arises from zero modes of the form ∂ R α s with nodes (so that there are perturbations with negative energies and, indeed, the only stable solution of this family is the spherical Skyrmion of charge 1).

Bumps with fractional Baryonic charge
The role of the odd-integer q (which does not enter directly in the Baryon charge in Eq. (33)) has not been discussed. A clear hint is the comparison of the energy-density contour plots above for two such crystal-like structures with the same Baryon charges, p and n but different q (for instance, p = 1, q = 1 and q = 3). In these plots, there are n bumps in the first structure with q = 1 which correspond exactly to the Baryon charge of the layers: thus, each bump carries one unit Baryon charge as expected. In the second structure with q = 3 there are 3n bumps in the r − θ plane but the Baryon charge has not changed. Thus, each bump carries 1/3 (in general 1/q) of the unit topological charge. Thus, the integer q is related to the fractional Baryonic charge carried by the bumps of the crystal. In the generic case in which p = 1 the interpretation is the following. The number of Baryonic tubes N tubes in the solution is N tubes = qn. On the other hand, each tubes carries a Baryonic charge equal to p q so that the total Baryonic charge is still B = np.

An example of topologically trivial solution
Here it is useful to discuss the following natural question: Since both the energy-density and the Baryon density only depend on two spatial coordinates, are the solutions constructed above really three-dimensional?
In fact, it is more useful to give an example of a configurations which looks similar to the ones analyzed in the previous section but, in fact, has vanishing topological density.
It is easy to see that the topological density vanishes (since G = 0). On the other hand, the above ansatz still gives rise to nontrivial Skyrme field equations for α and F. In the explicit parametrization in Eqs. (16), (17) and (18) one can see that, with the choice in Eq. (53), Eqs. (17) and (18) are identically satisfied while Eq. (16) reduces (once again) to Eq. (36). Despite the fact that both in this trivial case and in the topologically non-trivial configurations analyzed previously the equations for the profile α are the same, the corresponding configurations are totally different.
One obvious difference is that there is no topological argument which prevents the configuration defined in Eq. (53) from decaying into the trivial vacuum as the corresponding Baryon charge vanishes. Another difference is that the total energy of the trivial configuration defined in Eq. (53) which reads (assuming the same boundary conditions for α and F so that the definition of η (α, E 0 ) is the same as in the previous sections): trivial (α; q, n) = sin 2 α η (α, E 0 ) The differences between this expression and the ones in Eqs. (47) and (48) can be ascribed to the genuine threedimensional nature of the Baryonic tubes discussed in the previous sections and to the effective two-dimensional nature of the topologically trivial configurations defined in Eq. (53). In particular, comparing the expressions in Eqs. (47) and (48) with Eqs. (54) and (55), one can observe that to have Bary- 12 As shown both in the matrix form in Eq. (2) and in the explicit parametrization in Eqs. (16), (17) and (18). onic charge is quite expensive from the energetic point of view.
The example in this subsection discloses neatly the differences between configurations with and without Baryonic charge.

Conclusions
Analytic topologically non-trivial solutions of the complete set of Skyrme field equations of charge np describing ordered arrays of Baryonic tubes with crystalline structure living at finite density have been constructed explicitly. These configurations are characterized by three integers: n, p and q (q being odd). When q = 1, the Baryonic charge of each tube is p while n is the number of tubes so that the total Baryonic charge is B = np. On the other hand, when q is greater than 1, the Baryon charge of each tube is p/q while the total number of tubes is nq (so that the total Baryonic charge is still B = np). The positions of the peaks in the energy-density can be computed explicitly by maximizing the energy density. These configurations pass a non-trivial stability test 66.