Analytic Baby Skyrmions at Finite Density

We introduce a consistent ansatz for the baby Skyrme model in (2+1)-dimensions which is able to reduce the complete set of field equations to just one equation for the profile function in situations in which the baby baryon charge can be arbitrary. Many analytic solutions both with and without the inclusion of the effects of the minimal coupling with the Maxwell field are constructed. Linear stability and other physical properties are discussed. These analytic gauged baby Skyrmions generate a persistent $U(1)$ current which cannot be turned off continuously as it is tied to the topological charge of the baby Skyrmions themselves. In the simplest non-trivial case of a gauged baby Skyrmion, a very important role is played by the Mathieu equation with an effective coupling constant which can be computed explicitly. These configurations are a very suitable arena to test resurgence in a non-integrable context.


Introduction
The (gauged) baby Skyrme model in (2+1)-dimensions is one of the most interesting models admitting solitons without the huge analytic difficulties present in the (3+1)-dimensional case (see [1], [2], [3][4][5][6]; for two detailed reviews see [7] and [8]). Due to the fact that the (gauged) baby Skyrme model shares many relevant physical features with its (3+1)-dimensional progenitor (see [9] and references therein), it has been deeply analyzed as it can shed light on relevant non-perturbative features of a field theory with non-trivial topological sectors. On the other hand, despite the fact that it is simpler than the (gauged) Skyrme model in (3+1)-dimensions, it is manifestly a non-integrable system for which tools like the Lax pairs or the bi-Hamiltonian formalism cannot be applied. For these reasons, with few exceptions, most of the multi-solitonic solutions in the baby Skyrme model in (2+1)-dimensions have been constructed numerically (see [9] and references therein). In the present manuscript we want to construct analytic multi-solitonic configurations of the (gauged) baby Skyrme model in (2+1)dimensions. These enable one to disclose novel non-perturbative features of the theory.
We are especially interested in the following two important issues. Firstly, we would like to understand how finite-size effects affect multi-solitonic configurations. Secondly, we would like to understand the electromagnetic properties of these solitons when finite-size effects cannot be neglected.
The great physical interest of finite density effects arises from the fact that a theoretical description of cold and dense nuclear matter as a function of baryon number density is still lacking. In this sense, the (gauged) baby Skyrme model in (2+1)-dimensions is a very interesting model which is not integrable (but still easier than the (3+1)-dimensions (gauged) Skyrme model) in which one can analyze finite density effects with non-perturbative methods. In contrast to what happens in condensed matter physics, exact analytic results on the phase diagram of the low energy limit of QCD at finite density and low temperatures are extremely rare. The non-perturbative nature of low energy QCD explains why the very complex and interesting structure of its phase diagram (see [10], [11], [12][13][14][15] and references therein) is mainly analyzed with numerical and lattice approaches. A remarkable phenomena in the QCD phase diagram (which appears at finite baryon density) is the appearance of ordered structures like crystals of solitons (as it happens, for instance, in condensed matter theory with the Larkin-Ovchinnikov-Fulde-Ferrell phase [16]). Many results by now support the existence of these ordered structures at finite density (see, for instance, [17][18][19][20], and references therein). The very few analytic results have been found either in (1+1)-dimensions [18] or when some extra symmetries (such as SUSY) are included (see [9] and references therein). That is why it is important to shed more light on these intriguing phenomena as, often, even the numerical approaches are in trouble when the topological charge is very high: the analytic tools developed here below disclose novel features of multi-solitonic solutions in the (gauged) baby Skyrme model in (2+1)-dimensions.
The analytic approach which will be developed here will also shed light on a novel relation between (gauged) multi-solitons in the (gauged) baby Skyrme model and resurgence. Resurgence (see [21][22][23][24][25] and references therein) provides us with a concrete hope to make sense of the divergent character of perturbation series in QFT. This approach (using non-perturbative informations related with the non-trivial saddle points of the path integral) works, roughly, as follows: one is instructed to use Borel summation in the complex g-plane (g being a relevant coupling constant of the theory) to give a meaning to the factorially divergent series of perturbation theory. The initial divergent series becomes a finite expression. This expression possesses ambiguities which manifest themselves along suitable lines in the complex g-plane. At a first glance, one has just traded one problem for another (maybe even worse). In fact, when the theories one is interested in possess non-perturbative sectors labelled by suitable topological charges, the perturbative expansions in these topologically non-trivial sectors (which, usually, are also factorially divergent) allow a remarkable simplification (as shown for the first time in the physical literature in [25], [26]). The "Borel ambiguities" arising in the topologically non-trivial sectors exactly cancel those of the perturbative sector: in this way, by including all the nonperturbative sectors on the same footing, one gets a well-defined answer. Starting from the nice results in [27], in recent years this topic has attracted a lot of attention. The theories which have been mostly analyzed from the resurgent point of view have been quantum mechanics, topological strings, and integrable quantum field theory in low dimensions (see [21], [28], [29], [30], [31], [32], [33] and references therein). One of the main building blocks of the resurgent scheme is a proper analytic understanding of the non-trivial saddle points of the path integral. Hence, it is of fundamental importance to deepen our knowledge of the classical solutions of field theory with non-trivial topological charges. This task is very urgent for non-integrable field theories in which the theoretical tools coming from SUSY and the saturation of the BPS bound are not effective. As explained in the next sections, the framework developed in the present manuscript is useful from a resurgent perspective. Finite density effects are also interesting from this perspective. In particular, how the properties of multi-solitonic configurations depend on the typical size of a compact spatial domain in which these multi-solitonic configurations live, is especially interesting from a resurgent perspective (see [21][22][23][24][25][26][27][28][29][30][31][32][33] and references therein).
In the present paper, the methods of [34], [35], [36][37][38][39][40], [41] will be generalized to construct analytic solutions for baby Skyrmions at finite density. We will organize the paper as follows. In Section 2 we begin by discussing the nature of our ansatz. In Section 3 we analyse the results of a specific ansatz on ungauged baby skyrme models and in Section 4 we show how this can be extended to the gauged case, providing also some links to the topic of resurgence. We will conclude with a discussion of results in Section 5.

Beyond the spherical hedgehog ansatz
The action of the gauged baby Skyrme system is where λ is a Lagrange multiplier implementing the geometric constraint Φ · Φ = 1 in isospin space, a j , e are coupling constants and F µν F µν is the usual Maxwell Lagrangian. It is convenient to parametrize the fields as follows where and we orient the minimum of the potential in the position n = (0, 0, 1) .
In this parametrization the replacement of ∇ µ Φ with D µ Φ of the usual partial derivative with the covariant U (1) derivative defined above corresponds to just replace ∇ µ G with (∇ µ G + A µ ) in the action while ∇ µ F remains unaltered.
The parametrization in Eqs. (3),(4) is the most general parametrization of a gauged baby Skyrmion. The system is described by two scalar degrees of freedom, as Φ satisfies one scalar constraint, plus the degrees of freedom of the Maxwell field. Moreover, we can always rotate the Isospin vector n along the third internal axis. Hence, one can write the field equations either in terms of Φ (which satisfies Φ · Φ = 1) and A µ or in terms of F , G and A µ : the two ways to write the field equations are completely equivalent. In order to build a good ansatz, it is better for us to work with F , G and A µ .
With the above parametrization, the gauged baby Skyrme action reads where Obviously, the Lagrangian in Eqs. (5)-(6) is equivalent to the one in Eqs. (1)- (2). In other words, the Lagrangian in Eqs. (5)-(6) correctly describes the two scalar degrees of freedom of Φ as well as the degree of freedom associated to the gauge potential A µ . The gauged baby Skyrme field equations can be obtained by taking the variation of the action with respect to F , G and A µ . The easiest way to take into account finite-density effects is to introduce the cylinder flat metric where φ is an angular coordinate with range and L (which has dimension of length) is the compactification radius. The range of the coordinate r will be chosen as follows: where R represents the height of the cylinder. Physical constraints on R will appear when discussing the boundary conditions. Given the metric (8) and requiring the positivity of the energy, we choose for the coefficients a 1 , a 2 and a 0 the conditions The most difficult task is to find a good ansatz which has 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.
The topological charge of the configuration for A µ = 0 is where Σ is a two-dimensional space-like surface (for example t = const surfaces) and, 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 two independent functions) can be interpreted as saying that such two functions "fill a two-dimensional spatial volume" at least locally. 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 one to fix some of the parameters of the ansatz.
Following the strategy of [36][37][38][39][40] and [41], one can consider the following ansatz Note that the baby Skyrme field depends on the function G only through sin G and cos G so that this solution is periodic in time. Note also that, when ∂ t G = 0, the corresponding configurations do possess a non-trivial winding also along the time-like surface {φ = const}. Such configurations will have a time-period which "cannot be undone" as they wrap around a time-cycle in a topologically non-trivial way. The periodic time dependence in Eq. (16) is a natural way to avoid Derrick's no-go theorem on the existence of solitons in non-linear scalar field theories and corresponds to a time-periodic ansatz such that the energy density of the configuration is still static. The present ansatz defined in Eqs. (3), (16) has exactly this property. Indeed, the energy density does not depend on time as one can check with a direct computation. Moreover, the present ansatz also possesses a non-trivial topological charge corresponding to a winding along a t = const surfaces. On the other hand, it is worth to emphasize that the static solutions are the easiest to interpret. Thus, we will consider first in Sec. 3 static solutions in which ω = 0. The importance of introducing a time-dependence will be apparent in Sec. 4 on gauged baby Skyrmions where ω will have a specific relation with p.
When the gauge potential is non-vanishing A µ = 0 an extra term must be included in the definition of the topological charge in order to ensure both the topological and the gauge invariance. In analogy with the gauged Skyrme model in (3+1)-dimensions (see [2] and [44]) the topological charge of the gauged baby Skyrme model in (2+1)-dimensions is given by: Indeed, it is a direct computation to check that the topological charge in Eqs. (17) coincides with the usual definition of topological charge of the gauged baby Skyrme model in the literature [9]. The definition (17) can be decomposed into the sum of the ungauged topological density plus a total derivative From this form, it is evident that the topological charge is a boundary term.

Ungauged baby Skyrme model
In this section we will consider first the baby Skyrme model without the U (1) gauge field (A µ = 0 or formally e → 0). The variation of the baby Skyrme action with respect to F leads to the equation of motion The variation of the baby Skyrme action with respect to G leads to the equation of motion In the above equations c 3 = a 2 a 1 . The above choice of the ansatz (16) the field equations simplifies dramatically (as shown below) due to the fact that, with the function G in Eq. (16), One can check with a direct computation that the ansatz defined in Eqs. (75) and (16) does indeed reduce the full system of field equations to just one integrable equation for F . With such ansatz, the topological charge density (13) reads The above ansatz satisfies the first condition to be topologically non-trivial in Eq. (14) as F and G in Eq. (16) are independent. The boundary condition on F and the corresponding baby Baryon charge are The integers n and m appearing in Eqs. (24), (23) do not appear in the topological charge but, together with p, they are relevant to classify the map. The meaning of these integers needs some comments. If the integer p corresponds to the number of windings from the angular variable θ of the cylinder into the azimuthal angle G of the sphere, the integers n or m are connected with the "polar" windings of the sphere. A way to visualize the map (16) with the boundary conditions (24) or (23) can be thought following two steps. Firstly, we map the circumference of the cylinder at r = 0 p-times into the G angular coordinate at F = 0. Then, we extend this map along the r-direction on the cylinder wrapping the sphere along the polar F -coordinate. If we reach the value of F = π at r = R (the case n = 0 in (24)), the sphere is completely wrapped and the topological charge is B = p. However, if F goes over π, the map starts to unwind the sphere along the F -direction till F = 2π at r = R. In this case (m = 1 in (23)) the sphere is wrapped and unwrapped along the F -direction so that the topological charge is B = 0. Increasing again the value of F the map alternates the windings and unwindings of the target space along the "polar" direction. For this reason, if F at r = R is fixed to be an even number of π, the topological charge will always be zero; otherwise in the odd case the total sum of positive "polar" windings will be always one and the total topological charge B = p, regardless of n. A more rigorous mathematical proof of (24) and (23) can be performed after some comments on the expression (13). For our aim, it is necessary that the topological density ρ B is well defined for every interval of coordinates with particular attention on F when it exceeds the value of π. As a check, we expect a change of sign in the measure of the sphere (13) when F passes from the interval [0, π] to [π, 2π] (if we keep the same orientation along the G-direction). This change of sign must reflect the two different directions in the path from the north pole to the south pole of the sphere. We expect the right sign of the density (13) is provided by the presence of sin F for every interval of F and G. The explicit calculation of the topological charge is therefore as anticipated in (24) and (23). Nevertheless, even if n and m do not contribute to the total topological charge, they play an important role. As we discussed before, the integers n or m classify the number of windings and unwindings along the "polar" coordinate of the sphere. Therefore, even without knowing the exact form of the functions F and G, we can predict that these solutions will oscillate from positive to negative topological charge density. The number of these oscillations will depend on n and m (it will be 2n + 1 for the case (24) and 2m for (23)). These integers thus indicate the presence of solitonantisoliton (baryon-antibaryon) inside a solution of total charge B = p. We call it a topological soliton with charge p "dressed up" with n kinks. In a sense, the integers n and m can be thought of a second "charge" of the configuration.

Kink-dressed topological solitons
When one plugs the ansatz in Eq. (16) (remember that in this section the ansatz is static) into the baby Skyrme field equations in Eqs. (19)- (20), they reduce to only one integrable equation for F (r): When the exponent k in the potential (7) is an integer, some of the formulas below simplify (see the next discussion). It is a trivial computation to see that the field equation for G in Eq. (20) is identically satisfied with the ansatz in Eq. (16). One can write Eq. (26) as follows: and then So Eq. where and E 0 is a constant to be fixed by the boundary conditions (24), (23) integrating equation (30) as Using the ansatz in Eq. (16) all the baby Skyrme field equations in Eqs. (19) and (20) are satisfied if Eq. (28) is satisfied. The integral equation (32) represents a sort of state equation for the system that, once given the length R and the number of "bumps" (with F (R)), fixes uniquely the value of E 0 . In other words, the profile function F (r) can be found not only solving the second order differential equation (19) with the two boundary conditions in (24) or in (23) but even solving the first order equation (28) with a single boundary condition and the constant E 0 . In Figures 1, 2, 3, we show the plots of the profile function F (r) and the topological charge density ρ B (r) for different baby Skyrmion solutions for different boundary conditions, respectively n = 0, m = 1, n = 1 as defined in (24) and (23). The set of parameters used for all the plots are a 0 = 1 Using equation (28) to solve the profile function F (r) with the fixed parameters (33), the constants E 0 for the solutions n = 0, m = 1, n = 1 are respectively The plots obtained with equation (19) or (28) are totally equivalent. As anticipated in this section and as shown in the plots, the solutions with n > 0 or m > 0 have positive and negative oscillations in the topological charge density. Moreover, the number of "bumps" in the charge density is respectively 2n + 1 and 2m for configurations of the type (24) and (23). Every bump has charge p or −p so that these kinds of solutions form a chain of soliton-antisoliton or, precisely, baryon-antibaryon. In this sense, they represent a 2-dimensional example of crystal-like structures in baryonic matter at finite density.
where we define The above equation for E 0 always has a positive real solution. 2 Moreover, one can see that ∂ r F > 0 and that, when n is large, both η (F, E 0 ) and E 0 are of order n. The above integral I (n, p, E 0 , K, k; c j ) belongs to the family of (generalized) elliptic integrals at least for integer values of k in the potential (7). Let us consider first the above integral in Eq. (19) (which plays a very important role, as explained in a moment) when k = 1 and k = 2. In both cases, I (n, p, E 0 , K, k; c j ) in Eq. (38) has the form where the real coefficients N i and can be related with K, c j and E 0 of the theory. Then, both polynomial of fourth order in sin (F/2) can be factorized in their four roots (δ j for the denominator and ν j for the numerator, some of them can be complex conjugated): These types of integral are called "generalized elliptic integrals" in the literature (see [59,60] and references therein). They satisfy recursion relations and their asymptotic expansions are known explicitly. For simplicity we do not include these solutions explicitly here, the interested reader can simply solve equation (28) directly using the DSolve function in Mathematica. On the other hand, when k is too large, it is not possible in general to find explicitly the roots of the polynomial E 0 −W (F ) appearing in the denominator of Eqs. (38) and (39). Thus, in these cases, the theory of generalized elliptic integrals will not help to reduce the integral I (n, p, E 0 , K, k; c j ) in Eq. (38) to elementary functions. Equations (37) and (38) (which implements the boundary condition in order to have topological charge p) can be seen as a relation between the height of the cylinder R, the integration constant E 0 , the topological charge p, the kink number n and the coupling constants of the theory. In other words, Eq. (37) represents a sort of equation of state for the system which must be satisfied in order for the solutions to exist.
The conclusion is that the ansatz in Eq. (16) in the flat cylindrical metric in Eqs. (8) and (9) in which the profile F is given in closed form in Eqs. (29), (30) and (37) gives rise to exact and topologically non-trivial solutions (with baby baryonic charge p) of the field equations for any integers n and p.

Energy density
From the energy-momentum tensor the energy density T 00 with the ansatz in Eq. (16) reads where We show in Figure 4 the plots of the energy density for different solutions with the same parameters used in (33). It is worth to note that since, in this section, we are considering static solutions, the field equations for F in Eq. (26) can also be obtained considering the energy-density as a one-dimensional effective Lagrangian for the static profile F (r). Moreover, (replacing (∂ r F ) 2 with η (F, E 0 ) 2 using Eq. (30)), one can get an expression of the energy-density which does not depend on F : where the integration constant is fixed by equation (37). The plots of the energy density as a function of r from 0 to R reveals 2n + 1 or 2m bumps. Thus, this family of solutions with charge p is dressed with 2n + 1 or 2m kinks. Note also that in this case the energy density of a given field F depending only on one coordinates (r in our cases) can be written as (where H j (F ) are functions of the field F ) and therefore BPS-like tricks are available. Indeed, so that one can "complete the squares" in BPS style, the last term in Eq. (46. ) being a total derivative. The integral of this total derivative defines a topological charge Q(n) which depends on n (or m), i.e. the integer which determines the number of bumps defined above). It is obvious that the above expression defines a boundary term: to see this, it is enough to define a function of F (let us denote it by Ω (F )) with the following property such a function always exists. Therefore, one gets: hence, such new "topological charge" is sensitive to the integer n (or m) appearing in the boundary conditions. From the physical point of view, this integer n (or m) measures the number of bumps (local maxima) in the energy density in the r direction. The corresponding BPS bound is It is worth to emphasize that, for the reasons explained above, this "topological charge" is different from the usual winding number of baby Skyrmions. One can see that there exist configurations which saturate the bound: such configurations satisfy Note that this equation is the same of (28) in the special case of E 0 = 0. Remarkably enough, also in this case, the saturation of the above BPS bound does imply the general field equations (compare Eq. (49) with Eqs. (28) and (30)). Thus, the present analytic solutions have two topological labels: the baby baryonic charge p and Q(n) (or simply n or m). Obviously, the configurations which saturate the above bound are stable. This happens when the integration constant E 0 in Eqs. (28) and (30) vanishes (which is the case when there is only one bumps in the energy density along the r direction). However, as long as the integration constant E 0 is such that F does not change sign (see below), also the solutions with higher n are stable (at least under r-dependent perturbations of the profile F ). As a last remark, the total energy for this solutions can also be written in closed form: Due to the fact that the field equation for F is integrable, one can compute explicitly the above integral over r. From Eqs. (28) and (30) so that taking into account Eqs. (43), (43) and (44) (as well as the boundary conditions on F ) we get Also the above explicit integral (which describes how the total energy depends on both topological charges, on the size of the cylinder as well as on the other parameter of the model) belongs to the family of generalized elliptic integrals described above (see [59], [60] and references therein). However, if the k in the potential term is too big, the roots of the polynomial cannot be computed explicitly. On the other hand, the above formula is very useful since it allows to compute explicitly the derivatives of the total energy with respect to relevant parameters.

A remark on the stability
In this section we will look more carefully at the stability of our solutions. A complete treatment of the stability requires the study of the most generic fluctuations around a background static solution. Using F and G as variables of a vector Ψ, we expand the fields as where F 0 and G 0 are solutions of the equations of motion and δF and δG are generic perturbations.
Expanding the static energy in term of the perturbation δΨ, we arrive at where E 0 (Ψ 0 ) represents the energy of the baby Skyrmion and the ellipses indicate terms with higher powers of Ψ that can be neglected. The operator O is Hermitian and consists of a 2 × 2 matrix whose components are with This operator refers to the case of k = 1 for the potential (7). To simplify the expression (51), the fluctuation δΨ can be properly decomposed as where the set of functions {δΨ n } is a complete and orthonormal basis representing the set of the eigenfunctions of the operator O O δΨ n = λ n δΨ n , with {λ n } eigenvalues. Using the property of the basis {δΨ n }, the energy (51) reduces to for a whatever perturbation δΨ, identified uniquely by the set of coefficients {c n }. The stationary point Ψ 0 in the space of the field configurations is a local minimum of the energy, i.e. a stable field configuration, if and only if λ n ≥ 0 ∀n. Therefore, the problem of the stability, at least under small perturbations, reduces to a problem of eigenvalues for the operator O.
In order to resolve the equation for the eigenvalues, we need to specify the boundary condition for the perturbation δF and δG. These conditions must be chosen in such a way to not change the total topological charge (9) and, moreover, to not change the number of "bumps", i.e. the number of baryons and anti-baryons, on the cylinder. In particular, for background solutions of the type (23) we request periodic boundary conditions for the fluctuations δF (r, 0) = δF (r, 2π) , δG(r, 0) = δG(r, 2π) , The periodicity along the θ-coordinate is due to the cylindrical geometry. Otherwise, the condition along r depends on the periodicity of the background configuration due to (23). In the case of background with condition (24), we need instead δF (r, 0) = δF (r, 2π) , δG(r, 0) = δG(r, 2π) , Although numerical calculations are needed for the calculation of the eigenvalues {λ n }, analytically we are able to calculate the zero modes (λ n = 0) of the operator (52). One of these zero modes arises from the SO(2) exact symmetry of the Lagrangian and can be written explicitly as where α is a constant angle. This solution respects both the boundary conditions independently from the background field. Otherwise, the periodicity of the background solution F 0 built with (23) gives rise to a translational zero mode F 0 (r) → F 0 (r + ) , and then to the eigenvector δG(r, φ) = 0 , respecting the conditions (60). This last eigenvector solution belongs to that kind of fluctuations of the form i.e. fluctuations involving only the field F with the propriety of keeping the structure of the ansatz (16). Reducing the space of perturbations only to the perturbations of the type (65), the eigenvalues system of equations (58) reduces to a 1-dimensional equation. Using the oscillation/nodal theorem, we identify the eigenfunction (64) with the groundstate of the system being such a solution with no node. Since in this case the groundstate has λ 0 = 0, then λ n ≥ 0. As a remark, we can affirm that the baby Skyrmion solutions with boundary (23) are stable under any fluctuation of the type (65). As anticipated before, beyond the analytic calculations, the set of eigenvalues for a generic fluctuation that respects the boundary conditions (60) or (61) has been calculated numerically. When possible, we used the Mathematica built-in function NDEigensystem and when this failed we proceeded by discretizing the coordinate space of the operator manually. This discretization is achieved by using a uniform 2-dimensional lattice of size R × 2πL in the θ and r directions. We vary the number of points in order to check the convergence of our procedure. The derivatives are discretized using second-order finite differences, adapted to backwards and forward differences in the r boundaries, and periodic conditions on the θ ones. This generates a large matrix of the discretised operator, which is then passed through the Eigenvalues function in Mathematica, to extract the required eigenvalues.
In this paper, we examine the stability of the static solutions characterized by n = 0, m = 1, n = 1 in (24), (23) with parameters (33). The study of the small perturbations around the solution m = 1 has been realized with the use of NDEigensystem, calculating the smallest 400 eigenvalues in order of absolute value for the operator (52). Our analysis reveals no negative eigenvalues among this set of values. In the case of solutions n = 0 and n = 1, the more complicated boundary conditions (61) need the implementation of the eigenvalue problem on a lattice. For these solutions, the numerical procedure reveals no negative eigenvalues up to numerical accuracy O(10 −3 ). We therefore conclude that these solutions are stable.
As a remark, since the topological charge does not depend on m or n, all the solutions with m > 0 or n > 0 have to be considered metastable. In other words, in every topological sector the configurations with m > 0 or n > 0 represent local minima of the energy. The stability (or meta-stability) of a sequence of baby Skyrmions and anti-baby skyrmions represent a rather surprising result of this work. From one side, it constitutes an effect of the finite density approach since the system is forced to be confined in a finite box (a finite cylinder in this case). Moreover, the baby skyrmions built on a cylinder with the ansatz (16) seem to have different features compared with the solutions built on the R 2 plane. For example, on the plane a baby and an anti-baby skyrmion can be distinguished from the different orientation of the map G in (3) with respect to the polar angle θ. Otherwise, in the ansatz (16) the map G does not change and the sign of the topological density depends on the orientation of the map F on the sphere. In general, the different geometry between the plane and the cylinder leads to different characteristics for the solitons. As a consequence, we expect a different kind of interaction among them. The analysis of these new interactions can be an interesting point for a future work.

The case of potential
In this section, we consider the baby Skyrmion theory defined on a cylinder with potential (7). Using again the ansatz (16), the equations of motion with the new potential are the same as Eq. (19) and (20) replacing V with the expression (66). All the considerations discussed before about the topology of these solutions hold and again the baby Skyrmions are classified by the topological charge B = p and by the integers n and m. In the following, we show the plots of some solutions calculated with this set of parameters a 0 = 1 The plots obtained with equation (19) or (28)  The study of the stability follows the same analysis of section 3.3 in which the operator O ij has the same expression of (52) replacing O 11 with The procedure and the conclusions of our analysis about the eigenvalues of this new operator are equivalent to the case of section 3.3. About the solution m = 1, using NDEigensystem, we do not find any negative eigenvalues among the smallest 400 eigenvalues in order of absolute value for the operator (52) with (71). In the case of solutions n = 0 and n = 1, the numerical computation on the lattice once again reveals the presence of no negative values up to numerical error. In the case of the two-vacuum potential is possible to build the baby skyrmion on an infinite cylinder keeping the energy finite. To this end, in the simple n = 0 configuration, we must impose the boundary conditions for F as This solution has been discussed in [62] with the name of "Skyrmion wall". A further analysis of these kind of solutions and their stability can be extended to the case of n > 0 or m > 0 in a future work.

Gauged baby Skyrmions
We are finally ready to consider the system coupled to the gauge fields. Despite its high interest, very few analytic solutions have been found in this system. In this subsection we will consider the gauged baby Skyrme action in Eq. (1) and we set the coupling e = 1. This action is invariant under the SO(2) U (1) gauge transformation The Maxwell equations corresponding to the above action read The approach of [36][37][38][39][40][41] leads to the ansatz (3) with and for the gauge field The components of the tensor F µν define the electric-like and magnetic-like fields as Note that we chose a particular chiral relation between ω and p. This will be crucial in order to achieve the simplifications of the equations. With the ansatz (3) and (75) the following terms vanish: Quite remarkably, with the above ansatz for the gauge potential A µ and for the baby Skyrmion, it is possible to reduce the complete set of coupled gauged baby Skyrme model field equation in Eqs. (19), (20) together with the corresponding Maxwell equations in Eq. (74) to a single integrable equation for the profile F as well as to only one linear Schrodinger-like equation in which the effective periodic potential can be computed explicitly in terms of the profile F itself.
The reason behind this surprising simplification is the following. All the terms which can couple A µ to − → Φ in the gauged baby Skyrme field equations involve (at least) one of the following scalars Υ j (j = 1, 2, 3 and 4): Hence, the choice in Eqs. (75) and (76) kills all the terms (see Eq. (78)) which, generically, would couple A µ with − → Φ . Obviously, the system is still coupled since the U (1) current corresponding to the gauged baby Skyrmion (which appears on the right hand side of the Maxwell equations) depends explicitly on − → Φ . However, the present approach greatly simplify the coupled field equations of the Maxwell gauged baby Skyrme system since (thanks to Eqs. (75)-(78)) one can first solve the field equations for − → Φ (in which A µ disappeared for the above reasons and then, once − → Φ is known, one is just left with the Maxwell equations which, in fact, reduce to only one linear equation for η(r) in which − → Φ plays the role of an effective potential. Note that here no approximation will be made: namely we are dealing with the complete set of five coupled non-linear field equations for the gauged baby Skyrme model in (2+1)-dimensions in which both the back-reaction of the solitons on Maxwell field and viceversa are explicitly taken into account in a self-consistent way.
More concretely, the problem to construct analytically gauged baby Skyrmions problem can be divided into two steps.
Therefore, the above ansatz is able to reduce in a self-consistent way (without any approximation) the gauged baby Skyrme equations with a non-trivial gauge potential to just one equation for F which is integrable in the sense that it can always be reduced to a quadrature. When k = 1 or k = 2, Eq. (79) can be integrated explicitly in terms of inverse elliptic functions (Jacobi Amplitude functions [61], as the reader can quickly verify solving the differential equation below) observing that it is equivalent to the following first order equation The integration constant E 0 can determined requiring to have the correct boundary conditions (24), (23) in order to obtain a non-vanishing integer topological charge. Second step: the three Maxwell equations with the current in Eq. (74) corresponding to the gauged soliton reduce consistently to only one equation for η(r) (defined in Eq. (76)): where Note that, since the equation for F is solvable (see Eq. (80)), both sin 2 (F ) and (F ) 2 above are known explicitly. It is useful to rewrite Eq. (81) as Here it is also useful to write explicitly the U (1) conserved current J µ in the general form that due to the ansatz in Eqs. (75), (76) reduces to: The components of J µ represent the source of the electric and magnetic field where ρ e is the electric charge density and J r , J φ the currents. The new energy-momentum tensor with the gauge field A µ reads and the energy density ρ with the ansatz (75), (76) reduces to The expression for the topological charge B differs from the ungauged case due to the presence of A µ as anticipated in (17). The explicit calculation for B with the ansatz (75) gives As expected, the topological charge depends on the boundary conditions of the fields as (recall we always take F (0) = 0 as in (24), (23)) In order to keep the topological charge unchanged, the last result imposes one boundary condition for the field η η(R) = 0 .
This condition is a natural choice since we do not expect that the Maxwell field A µ introduces a baryon charge in the system. However, mathematically this condition can be relaxed if we simply ask the charge B to be an integer independently of which field "carries" the baryon charge. In this case, the most general condition for η becomes The discussion of this point is beyond the aim of this paper and only the condition (90) will be used. Now, let us consider the simplest non-trivial example of analytic gauged baby Skyrmion which can be constructed with vanishing potential a 0 = 0. In this case, using the ansatz defined in Eqs. (75), (76), one can easily see that the profile F must be linear: (as in this way Eq. (79) is satisfied) where k 1 must be chosen in such a way to satisfy the boundary conditions in Eq. (23) With the above expression for the profile F (r), the remaining Maxwell equation (82) reduces to In order to solve this equation we need to specify the boundary conditions for Ψ. One condition is fixed only for the case of F with boundary (23) by Eq. (90), i.e.
Otherwise, for the condition at r = 0 there is not a unique choice. The different boundary conditions that we impose to Ψ fix the electric charge and the current inside the system, leading to different physical configurations. In this paper, we show and consider only three cases but other choices are allowed. In the following, we use for simplicity the same conditions (97), (96) also for the case of F with boundary (  the energy density ρ(r). In all the cases, the electric field E r = −B as shown in (77), the electric charge density ρ e = J φ due to (84) and E φ = 0, J r = 0. As it is clear form the plots, the periodic boundary conditions for Ψ give rise to a positive angular current J φ which generate a magnetic field with opposite sign on the borders of the cylinder. In the same way, a positive electric charge exists and it generate an electric field E r that spreads in opposite directions. The energy density ρ is modified by the presence of the electric/magnetic field and it is mainly concentrated at the boundaries of the cylinder.
In Figure 12, we show the case of anti-periodic boundary conditions for Ψ, i.e. Ψ(0) = −Ψ(R) = −ω, only for the solution n = 0 using again the parameters (98) (the solutions m = 1 and n = 1 have a similar profile as it happens for the previous case). This configuration is characterized by a negative current on one half of the cylinder and a positive current, flowing in the opposite φ-direction, on the other half. The magnetic field is therefore negative everywhere. Consistently, the presence of a negative and positive electric charge leads to an electric field E r which goes in the same direction everywhere along the cylinder. The energy density has the identical profile of the periodic case of figure 9. The periodic and anti-periodic Ψ-configurations are therefore degenerate in energy.
In Figure 13 we impose a vanishing boundary condition for Ψ at one side of the cylinder, i.e. Ψ(0) = 0 and Ψ(R) = ω, for the n = 0 solution with the set of parameters (98). In this case a positive current flows on the cylinder along the φ−direction, similarly to Figure 9. However, here the magnetic field is unexpectedly only negative. Looking at the electric side, we have a positive electric charge whose electric field goes only to one direction instead of spreading its flux in the two opposite sides.
We can also investigate what happens in the more complicated case of a 0 = 0. In particular, we set a 0 = 1 for equation (80) and the other parameters are the same of (98). We also choose k = 1 for the potential. The solutions are calculated imposing periodic boundary conditions for Ψ, i.e. Ψ(0) = Ψ(R) = ω, respectively for n = 0 and m = 1 in figures 14 and 15. These configurations follow the same considerations that we discussed for figures 9 and 10 with the difference, for the case n = 0, that the potential breaks the symmetry with respect to R/2.
The conclusion of this section is that the coupled system of 5 coupled non-linear field equations made by the 2 field equations of the gauged non-linear baby Skyrme model and the 3 Maxwell equations

Gauged baby Skyrmions and resurgence
We now discuss an interesting interplay between the present framework and the resurgence scheme. A very difficult problem in the application of resurgence in non-integrable theories with multi-solitonic configurations is the correct identification of the coupling constant which is responsible for the resurgent behavior itself. When finite density effects are present, the combination of parameters of the theory which plays the role of effective coupling constant can be far from obvious. Besides the "obvious" coupling constant (which plays the role of in the action), other important parameters could be the "size" of the region in which the multi-solitonic configurations live, the soliton density, and so on. The correct resurgent parameter could be a combination of those. This situation should be compared with the many applications of resurgence in the literature in which the identification of the coupling constant is clean. An example which will be relevant for us is the Mathieu equation whose resurgent behavior has been deeply investigated (see [32] and references therein): In the above equation the only possible resurgent parameter is while E plays the role of an eigenvalue.
Here we show an explicit example of a solitonic system in which we are able to clearly identify the effective coupling constant that can be used for the resurgence study. This identification is possible comparing the baby Skyrme equations with the Mathieu equation (99). The example of system that we analyse is the gauged baby Skyrmion discussed in the above section with a 0 = 0. Using a trigonometric equality, the equation (94) can be rewritten as Now, let us define the new variable x and the negative parameter a 1 as follows: The remaining Maxwell equation (100) reduces to Hence, the above computations (together with a comparison with Eq. (99)) suggest that in this simplest non-trivial case the role of 2 2 is played by: Thus, in particular, the resurgent parameter is not only related to the "obvious" coupling constants a 1 and/or a 2 but actually is related with the suitable combination of both of them defined here above (which also involves the size of the system through R). It is worth to emphasize that only when analytic gauged solitonic solutions are available it is possible to discover the "proper coupling constant" (which, in the present case, is defined in Eq. (102)): this shows that the present analytic tools can be extremely useful from a resurgent perspective. Eq. (101) is a particular case of Eq. (99) analyzed, for instance, in [32] in which the eigenvalue has been fixed to -1. On the other hand, the spectrum of the electromagnetic perturbations of these gauged solitons defined in Eqs. (75), (76), (92), (93) with a 0 = 0 in the gauged baby Skyrme model minimally coupled to Maxwell theory is closely related with the spectrum of the operator To prove this, it is sufficient to expand the field Ψ around a background solution Ψ 0 (keeping F and G unchanged) as Ψ = Ψ 0 + δΨ δΨ = c n δΨ n (103) where {δΨ n } is a complete basis. In this case, the energy of the system with the electromagnetic perturbations reads where λ n is the set of eigenvalues of the operator − 8π 2 g|a 1 |R 2 ∂ 2 ∂x 2 + cos (x) + 1 δΨ n = λ n δΨ n .
This operator is simply calculated linearizing equation (101). Defining for simplicity λ n = λ n − 1, we obtain exactly the Mathieu equation (99) with the effective resurgence parameter. Thus, quite interestingly, the full resurgent analysis of the Mathieu equation in [32] is relevant to determine the spectrum of the operator O which, in its turn, plays an important role to determine the (spectrum of the) electromagnetic perturbations of these analytic gauged solitons. We hope to come back on this interesting issue on the resurgent analysis of the perturbations of these gauged solitons in a future publication.
One may also ask what happens in the more general cases in which the potential term in Eq. (79) does not vanish a 0 = 0. In this case, the profile can be integrated explicitly in terms of the quadrature in Eq. (80). When k = 1 or k = 2, the quadrature in Eq. (80) is explicitly solved in terms of inverse elliptic functions (Jacobi Amplitude functions [61], call them A(r)). Thus, the remaining Maxwell equation in Eq. (82) looks similar to a deformed Mathieu equation with the differences that, firstly, the argument "r" of the trigonometric function has been replaced with A(r) and, secondly, the coupling constant a 1 g has been replaced by a new function involving the derivative of A(r). Schematically, one gets: ∂ 2 ∂r 2 + a 1 g sin 2 π R r Ψ = 0 ⇒ ∂ 2 ∂r 2 + a 1 1 − A 2 sin 2 (A(r)) Ψ = 0 .
Thus, besides the explicit appearance of the coupling constants as in the previous simpler case, now the potential depends in a more implicit way on the boundary condition (through the integration constant appearing in the function J(r)). We hope to come back in a future publication on the resurgent analysis of these Schrodinger-like problems arising in the description of the multi-solitonic solutions of the gauged baby Skyrme model in (2+1)-dimensions.

Persistent residual currents
Another interesting characteristic of the gauged baby Skyrmions, which is valid also in the present case, is that they generate a current with many features in common with superconductive currents. Let us consider, once again, the simplest non-trivial case in the previous subsection of eq. (92), (93) in which a 0 = 0 (no potential term). One can notice the following characteristic of the current in Eq. (84).
1) The current does not vanish even when the electromagnetic potential vanishes (A µ = 0 ).

2)
This left-over current J (0)µ cannot be turned off continuously. One can try to eliminate J (0)µ either deforming F to integer multiples of π (but this is impossible as such a deformation would kill the topological charge as well: see Eq. (13)) or deforming G to a constant (but also this deformation cannot be achieved for the same reason). Note also that, as happens in the construction of superconducting currents in [42], G is defined modulo 2π (as − → Φ depends on cos G and sin G rather than on G itself). This implies that J (0)µ defined in Eq. (106) is a persistent current supported by the present gauged baby Skyrmions. Needless to say, the persistent character of a current is the cleanest hallmark of superconductivity. Moreover, just recently in [33], the first glimpse of resurgent structures appeared in integrable many body models of superconductivity. The present results open the intriguing possibility to analyze superconductive currents with the tools of resurgence in a non-integrable model such as the gauged baby Skyrme model in (2+1) dimensions.
A final remark on the stability of these solutions has to be made. The eigenvalue analysis presented above for the ungauged case is applicable in this case aswell, with the extensive complexification of having to deal with all types of gauge perturbations possible. This type of analysis is computationally hard to achieve as the size of the pertubation operator becomes large. If we reduce the space of perturbations to only the ones deforming the field profiles as we did above, then we expect a similar stability analysis resulting in no negative eigenvalues. The full perturbation space is however much larger and we leave this full stability analysis as an open problem.

Conclusions
Analytic topologically non-trivial solutions of the complete set of baby Skyrme field equations of arbitrary charge describing topological solitons dressed by n kinks (with ordered structure) living at finite density have been constructed explicitly. ... These configurations are characterized by two integers: n (which is the number of kinks dressing the baby Skyrmion) and p (which is the charge of the Skyrmion). All these configurations pass a non-trivial stability test. Moreover, the first analytic examples of gauged baby Skyrmions with nonvanishing topological charge in (2+1)-dimensions in the gauged baby Skyrme model minimally coupled to Maxwell theory have been constructed. The complete set of 5 coupled non-linear field equations can be reduced (using a judiciously chosen ansatz) in a selfconsistent way to one linear Schrodinger-like equation with an effective periodic potential keeping alive, at the same time, the baby baryonic charge. These analytic gauged baby Skyrmions generate a persisten U (1) current which cannot be turned off continuously as it is tied to the topological charge of the baby Skyrmions themselves. It appears that the present gauged multi-Solitonic configurations are a very suitable arena to test resurgence in a non-integrable context as these solitons allow a clean identification of the proper coupling constant which is responsible for the resurgent behavior. We hope to come back on the relations of these gauged baby Skyrmions and resurgence in a future publication.