Gauged Aloof Baby Skyrme Model
Abstract
We present a study of gauged modification of the 2+1 dimensional planar Skyrme model with a particular choice of the symmetry breaking potential term which combines a shortrange repulsion and a longrange attraction. In the absence of the gauge interaction the multisolitons of the model are aloof, they consist of the individual constituents which are well separated. Peculiar feature of the model is that there are usually several different stable static multisoliton solutions of rather similar energy in a topological sector of given degree. We investigated the pattern of the solutions and find new previously unknown local minima. It is shown that coupling of the aloof planar multiSkyrmions to the magnetic field strongly affects the pattern of interaction between the constituents. We analyse the dependency of the structure of the solutions, their energies and magnetic fluxes on the strength of the gauge coupling. It is found that, generically, in the strong coupling limit the coupling to the gauge field results in effective recovering of the rotational invariance of the configuration.
1 Introduction
The study of topological solitons in field theory can be traced back to the seminal paper by Skyrme [1] where valued nonlinear model for atomic nuclei was suggested. The Skyrme model can be derived from the expansion of the QCD low energy effective Lagrangian in the large limit [2], then the topological charge of the multisoliton configuration is set into correspondence to the physical baryon number. Further, under certain assumption the semiclassical quantization of rotations and isorotations of the Skyrmions allows us to get a good approximation to the isospinning nuclei and describe the corresponding excitations which are associated with pions [3, 4]. Thought, the Skyrme model has a limited success, there are several problems with description of the nuclear masses since the interaction energy of the Skyrmions is much higher than the corresponding experimental data for nuclei. Recent experimental observation of the heavy pentaquarks [5], which could be considered as a sort of baryonmeson bound state [6] also cannot be explained in the conventional Skyrme model with the usual pion mass term.
Several modifications of the Skyrme model in 3+1 dimensions were suggested recently [7, 8, 9, 10], mostly related with modification of the potential of the model. Furthermore, the contribution of the Coulomb electromagnetic energy is necessary to get a good agreement between the binding energies of heavy nuclei and the predictions of the reduced BPS Skyrme model [7]. Therefore it is physically natural to extend the model by gauging it to describe various electromagnetic processes of nucleons.
The gauged Skyrme model was originally proposed in [11], later the axiallysymmetric gauged Skyrmions were considered in [13, 14]. It was noticed that the gauging of a subgroup may stabilize the solitons even if the Skyrme term is dropped [12], furthermore, in the gauged Skyrme model the topological energy bound becomes saturated.
The planar reduction of the nonlinear sigma model is known as baby Skyrme model [15, 16]. This (2+1)dimensional simplified model resembles the basic properties of the genuine Skyrme model in many aspects. Furthermore, the baby Skyrme model has a number of applications on its own, e.g., in condensed matter physics where Skyrmion configurations were observed experimentally [17], in the description of the topological quantum Hall effect [18, 19], or in brane cosmology where the solitons of the model induce warped compactification of the 2dimensional extra space [20]. Also it was found that the restricted baby Skyrme model in 2+1 dimensions has BPS soliton solutions saturating the topological bound [22].
A peculiar feature of the planar Skyrme model is related with the particular choice of the potential term which is necessary to stabilize the solitons^{1}^{1}1 Note that a special choice of the parameters of the baby Skyrme model allows to evade the potential term and construct static soliton solutions in the reduced model [25], however dynamical properties of these special configurations are still unknown.. In lowdimensional systems the effect of this term becomes more significant than in the original Skyrme model, for example there are different choices of the potential related with various ways of symmetry breaking [26, 27, 28, 29]. In particular, a suitable choice for the potential term allows us to separate the individual constituents of the planar Skyrmions, each of them being associated with a fractional part of the topological charge of the configuration [28]. Another possibility is to combine a shortrange repulsion and a longrange attraction between the solitons [29, 30]. Then the multisoliton configuration consists of aloof constituents.
Clearly, coupling the model to the electromagnetic field yields another possible channel of interaction between the solitons. Analysis of the gauged baby Skyrmions [23] reveals very interesting features of the corresponding solutions which carry a nonquantized nontopological magnetic flux. Further, if the ChernSimons term is additionally included in the Lagrangian, the planar Skyrmions become electrically charged [21]. Recently, the properties of the soliton configurations in the gauged BPS baby Skyrme model were investigated [24]. An interesting observation is that in the strong coupling limit the magnetic flux becomes quantized though there is no topological reasons for that [23, 31].
The aim of this paper is to discuss the structure of the soliton solutions of the gauged aloof baby SkyrmeMaxwell system with a rotationally invariant potential term analogous to that used in [29, 30]. Our calculations are performed for multisoliton solutions up to charge 10. We study numerically the dependence of masses of these configurations and the corresponding magnetic fluxes on the gauge coupling constant, both in perturbative limit and in the strong coupling limit without any restrictions of symmetry.
2 The model
We consider a gauged version of the model with the Skyrme term in dimensions [16] with a Lagrangian
(1) 
Here denotes a triplet of scalar fields, which is constrained to the surface of a sphere of unit radius: . We introduced the usual Maxwell term with the field strength tensor defined as . Clearly, the scaling dimensions of this term and the Skyrme term, which is quartic in derivatives, are identical.
The coupling of the scalar field to the gauge field is given by the covariant derivative [12, 23, 32]
(2) 
where is the gauge coupling constant. The localized field configuration has finite energy if , and as .
Topological restriction on the field is that it approaches its vacuum value at spacial boundary, i.e. . This allows a onepoint compactification of the domain space to and the field of the finite energy solutions of the model is a map which belongs to an equivalence class characterized by the topological charge . Explicitly,
(3) 
A special feature of the solutions of the planar Skyrme model is that their structure strongly depends on the particular form of the potential term . The most common choice is so called "old potential" [16]
(4) 
which is an analogue of the standard pion mass term in (3+1)dimensional Skyrme model. The symmetry is broken via the potential to and there is a unique vacuum . The corresponding solitons of degree are axially symmetric [16] however the rotational symmetry of the configurations of higher degree becomes broken [34].
In the model with double vacuum potential (or "easyaxis" potential) [15, 33]
(5) 
the multisoliton solutions are rotationally invariant over entire range of values of the mass parameter . The most general case of the oneparametric potential
(6) 
with was considered in [27]. Since the "old" potential (4) corresponds to the attractive force acting between the solitons, while the "holomorphic" potential is repulsive [35, 36], the parameter in the potential (6) is responsible for the balance of the repulsive and attractive interaction between the Skyrmions. Further, one can consider the linear combination of the "old" and "holomorphic" potentials [29, 30]
(7) 
which corresponds to a shortrange repulsion and a longrange attraction between the solitons. Following [29, 30], we restrict our consideration to the case .
The resulting multisoliton configuration is no longer rotationally invariant, each constituent of unit charge is clearly separated from other, they form a cluster structure. Further, in the absence of the gauge interaction there are usually several different static multisoliton solutions of rather similar energy in a topological sector of given degree [29]. Note that this feature is in common with Hopfion solutions of the FaddeevSkyrme model [39]. Indeed, the structure of both models looks similar, the corresponding Lagrangian (1) includes the usual sigma model term, the Skyrme term, which is quartic in derivatives of the field, and the potential term. Further, the field of the FaddeevSkyrme model in 3 spacial dimensions is also a threecomponent unit vector restricted to the unit sphere . However, the domain space of the latter model is the compactified three dimensional sphere and the Hopfion solutions are classified by the linking number associated with the homotopy group . Thus, in some sense the baby Skyrme model can be considered as a planar reduction of the FaddeevSkyrme model [40].
In 2+1 dimensions the electric field of configuration is vanishing everywhere [23] and we can consider purely magnetic field generated by the Maxwell potential
(8) 
where the gauge fixing condition is used to exclude the component of the vectorpotential. Thus the magnetic field is orthogonal to the plane: . Note we do not use here the rotationallyinvariant parametrization of the fields where the gauge fixing condition is imposed [16, 15, 33].
The complete set of the field equations, which follows from the variation of the action of the baby SkyrmeMaxwell model (1), can be solved when we impose the boundary conditions. As usually, they follow from the regularity on the boundaries and symmetry requirements as well as the condition of finiteness of the energy and the topology. In particular we have to take into account that the magnetic field is vanishing on the spacial asymptotic. Explicitly, on the spacial boundary we impose
(9) 
Here is the usual radial variable.
As initial guess for further computation we used various combinations of single solitons with preassigned phases and positions. These field configurations can be constructed via usual parametrization of the unit scalar triplet in terms of the complex field via the stereographic projection
(10) 
Let , , be the field of a single Skyrmion of degree one with an arbitrary position and phase, then the function , allows us to obtain the field of multisoliton configuration of degree with various relative phases and separations. For example, the function
(11) 
yields the input configuration in the sector of degree , which represents the chain of baby Skyrmions with opposite relative orientations and separation between the constituents .
3 Numerical results
In this section we outline our method for calculating the multisoliton solutions of the gauged baby Skyrme model. The numerical calculations are mainly performed on a equidistant square grid, typically containing lattice points and with a lattice spacing To check our results for consistency we also considered the lattice spacings .
To construct multisoliton solutions of the model (1) we minimize the corresponding rescaled energy functional of the static configuration
(12) 
with the aloof baby Skyrmions potential (7). The gauge coupling is a parameter, each of our simulations began at at fixed value of and , then we proceed by making small increments in . For comparative consistency with [29] in most of our calculations we choose and .
The numerical algorithm employed was similar to that used in [37]. Wellchosen initial configurations of given degree were evolved using the Metropolis method to minimize the energy functional (12). In our numerical calculations we introduce an additional Lagrange multiplier to constrain the field to the surface of unit sphere. Simulations were considered to have converged to local minima if the quantity , where is the time of computation in minutes, was less than . We also verify that the evaluated topological charge of the configuration is in agreement with the input integer value.
Charge  g  

1  0  0.4  1  2 
2  0  0.15  0.3  1 
3  0  0.3  0.6  2 
4  0  0.4  0.6  1 
5  0  0.4  0.6  1 
6  0 (global)  0.4  0.6  2 
6  0 (local)  0.4  0.6  2 
Charge  g  

7  0 (global)  0.4  0.6  1 
7  0 (local)  0.4  0.6  1 
8  0  0.4  0.6  1 
8  0  0.4  0.6  1 
9  0  0.4  0.6  1 
10  0  0.4  0.6  1 
Another check of the correctness of our results was performed by finding corresponding solutions of the EulerLagrange equations which follow from the Lagrangian (1) subject to a set of boundary conditions (9). The relative errors of the solutions we found implementing the NewtonRaphson fourth order finite difference method, are of order of or smaller.
First, we recover the usual solution. Its energy for is this is within of the previously known result [29]. This similarity is another nice validation of our numerical algorithm.
Since there are variety of configurations of planar aloof Skyrmions, in order to label different these multisolitons we are using the following notation: where is the topological charge of the planar configuration or its building block, and is the dihedral group of corresponding symmetry. For example, in the sector of degree 4 we have found the (which is the global minimum at ), and solutions.
Note, that among solutions we constructed using the rational map input (10), there are new configurations, which were not considered in [29] since they represent local minima. In Fig. 1 we exhibited the energy density plots of some of these new solutions, , and , respectively.
Constructing these configurations we vary relative orientation of the constituents in the internal space. Indeed, the total potential of interaction between two neighboring baby Skyrmions is proportional to [16], it is repulsive if the solitons are in phase () and the most attractive channel corresponds to opposite orientations of the solitons, .
However, the interaction energy is still negative for the solitons which are not completely out of phase. If the number of constituents is even, they always form pairs with opposite relative orientation, however for a system of solitons, where is an odd number, another possibility exist. The solitons may form a circular necklace where each of baby Skyrmions is rotated by angle with respect to its neighbor. Hence, in this case we considered composite solutions with soliton species, not only binary species like in [29]. Thus, for odd the total angle of internal rotation along the ring is , where is integer. Particular examples of this type, and solutions are displayed in Fig. 1. Clearly, configurations of that type represent local minima.
Another interesting example of configuration, which also was not discussed in [29], is the "tristar" configurations of degree . Here the relative phase of the 3 solitons with respect to the central constituent is , however, with respect to each other, they in phase, . Thus, the central component provides strong attraction to the outer solitons binding the configuration together.
As the gauge coupling increases from zero, the energy of the gauged aloof Skyrmions decreases since the magnetic flux is formed. Now each individual Skyrmion is coupled to a magnetic flux and the electromagnetic interaction modifies the usual pattern of scalar interaction of constituents. Effectively, the asymptotic field of the gauged soliton represents a triplet of dipoles, two scalar dipole moments are associated with asymptotic form of the scalar fields in the plane and the third magnetic dipole is orthogonal to this plane [23].
Note that effectively, using the Maxwell equation , one can set the magnetic field into correspondence with an effective circular electric current . On the other hand the term in the total Hamiltonian of the system effectively contributes to the mass of the scalar field, thus the Yukawa interaction between the gauged baby Skyrmions becomes stronger.
In Fig. 2 we exhibited the contour energy density plots of the gauged baby Skyrmions in the model (1) for at some set of values of the gauge coupling . First, we observe that for relatively small values of the coupling constant , two neighboring solitons with opposite orientations form pairs in accordance with the a binary species model [29].
Let us consider a few particular configurations. In Fig. 3, left plot we have plotted the dependency of energy of the static gauged static baby Skyrmion of unit charge and the chains as functions of . On the right plot we also displayed the soliton’s binding energy as function of the gauge coupling. Here is the energy of the onesoliton configuration at the corresponding value of the gauge coupling and we used the normalized units of energy per unit charge. From these plots it is clear that the energy of the configurations decreases as the gauge coupling becomes stronger while the binding energy is increasing. Thus, the balance between the repulsive and attractive forces in the ungauged model with the potential (7) is shifted towards the attraction and, as coupling becomes strong enough, the rotational invariance of the multisoliton configuration is restored. Further increasing of the gauge coupling makes the solitons width increasingly localized.
Note that as the coupling remains smaller than one, the electromagnetic energy is increasing, however in the strong coupling limit its contribution begin to decrease as continue to grow, as seen in Fig. 4, left plot.
We can understand this effect if we note that the conventional rescaling of the potential leads to . Thus, the very large gauge coupling effectively removes the Maxwell term leaving the limiting configuration of gauged planar Skyrmions coupled to a circular magnetic vortex of constant flux. Apparently, in such a limit the strong coupling with a vortex yields an effective potential term which supplements the potential (7).
Let us remind that the particular choice of the parameters of the aloof potential (7) is related with the condition of balance between a shortrange repulsion and a longrange attraction, decreasing of the parameter increases the repulsive component. Coupling to magnetic field modifies the structure of interactions, as the coupling increases, the attractive part of interaction becomes stronger. As a result, we can approach the limit providing the gauge coupling is nonvanishing. In such a limit, variation of the coupling allows us to manipulate the attractive part of interaction between the Skyrmions. Further, at some critical value of the gauge coupling there is a transition from the aloof structure of the multisolitons to the rotationally invariant configuration, for example in the sector of degree it happens as and .
Thus, the configurations carry total magnetic flux , which is in general, nonquantized. The flux of the gauged baby Skyrmions is associated with the position of the solitons, it is orthogonal to the plane [23]. In the usual model with rotationally invariant potential (5), or in the strong coupling limit of the gauge model (1), there is a single magnetic flux through the center of the soliton, in the aloof system, where rotational invariance becomes violated, each unit charge constituent of the multisoliton configuration is coupled to a flux.
An interesting observation is that as the gauge coupling becomes stronger, the magnetic flux of the degree baby Skyrmions grows from 0 to , i.e. in the strong coupling regime the magnetic flux is quantized thought there is no topological reasons for it [23]. Indeed, in Fig. 4, right plot, we display the results of our numerical calculations of the integrated magnetic field of the gauged planar Skyrmion through the plane. In the limit the magnetic flux is vanishing, in the weak coupling regime the fluxes are attached to the individual partons of the multisoliton configuration.
As gauge coupling increases further, the radius of the each vortex is getting smaller and the magnitude of the magnetic field increases significantly. For the flux tends to be quantized in units of and the configuration becomes rotationally invariant. The energy density distribution in the strong coupling limit is localized near the center of the configuration approaching a singular stringlike distribution in the limit [23].
4 Conclusion
The main purpose of this work was to present new type of gauged solitons in the planar SkyrmeMaxwell theory. In the model with aloof potential (7) the individual solitons are composed out of the constituents with unit topological charge. Peculiar feature of the model is that, similar to the solutions of the FaddeevSkyrme model, a number of local energy minima of various types exists in each topological sector. Coupling of the solitons to the gauge sector yields additional attractive interaction which modifies the structure of configurations of given degree. It has been shown that the rotational invariance of the multisoliton configuration is restored in the strong coupling regime.
Similar to the corresponding solutions in the gauged Skyrme model and FaddeevSkyrme model [38], the planar configurations are topologically stable and in the weak coupling regime they carry a nonquantized magnetic flux which is orthogonal to the plane and penetrates the Skyrmion. In the strong coupling limit the total magnetic flux, associated with the partons, becomes quantized in units of topological charge.
We confirm that the mass of the static configuration decreases when the electromagnetic coupling constant is increased, thus a baby Skyrmion can lower its mass by interacting with the electromagnetic field. Also when g is increasing, the binding energy of the solitons increases.
Finally, note that the planar Skyrmions appear as quasiparticles in various systems, in particular they are natural objects at the description of the integer quantum Hall effect [19]. The vortices coupled to the planar Skyrmion constituents may appear in multicomponent superconductors [41], thus as avenues for further research, it would be interesting to extend the solutions in this work to the effective condensed matter systems. Another interesting possibility is to consider the gauged baby Skyrme model with aloof potential and ChernSimons term [21] without restrictions of symmetry, this would allow us to include the electric field into consideration.
Acknowledgments
We thank Igor Bogolubsky, Derek Harland, Martin Speight and Andrzej Wereszczyński for useful discussions and valuable comments. This work is supported in part by the A. von Humboldt Foundation in the framework of the Institutes linkage program and by the JINR HeisenbergLandau program (Y.S.).
References
 [1] T. H. R. Skyrme, Proc. Roy. Soc. Lond. A 260 (1961) 127.
 [2] E. Witten, Nucl. Phys. B 223 (1983) 422; ibid 223 (1983) 433.
 [3] G.S. Adkins, C.R. Nappi and E. Witten, Nucl. Phys. B 228 (1983) 552.
 [4] G. E. Brown and M. Rho (Eds.) ‘‘The Multifaceted Skyrmion’’, World Scientific Publishing, 2010.
 [5] R. Aaij et al. [LHCb Collaboration], Phys. Rev. Lett. 115 (2015) 072001
 [6] M. Karliner and J. L. Rosner, Phys. Rev. Lett. 115 (2015) 12, 122001
 [7] C. Adam, C. Naya, J. SanchezGuillen and A. Wereszczynski, Phys. Rev. Lett. 111 (2013) 23, 232501
 [8] P. Sutcliffe, JHEP 1104 (2011) 045
 [9] V. B. Kopeliovich, B. Piette and W. J. Zakrzewski, Phys. Rev. D 73 (2006) 014006
 [10] S. B. Gudnason and M. Nitta, Phys. Rev. D 91 (2015) 8, 085040
 [11] C. J. Callan Jr. and E. Witten, Nucl. Phys. B 239 (1984.) 161.
 [12] B.J. Schroers, Phys. Lett. B 356 (1995) 291.
 [13] B. M. A. G. Piette and D.H. Tchrakian, Phys. Rev. D 62 (2000) 025020.
 [14] E. Radu and D.H. Tchrakian, Phys. Lett. B 632(2006) 109.

[15]
A.A. Bogolubskaya and I.L. Bogolubsky,
Phys. Lett. A 136 (1989) 485
A.A. Bogolubskaya and I.L. Bogolubsky, Lett. Math. Phys. 19 (1990) 171. 
[16]
B. M. A. G. Piette, W.J. Zakrzewski, H.J.W. MuellerKirsten and D.H. Tchrakian,
Phys. Lett. B 320 (1994) 294
B. M. A. G. Piette, B.J. Schroers and W.J. Zakrzewski, Z. Phys. C 65 (1995) 165.  [17] X.Z. Yu et al, Nature 465 (2010) 901.
 [18] A. Neubauer et al, Phys. Rev. Lett. 102 (2009) 186602.
 [19] S.M. Girvin, "The quantum Hall effect: novel excitations and broken symmetries" In "Aspects topologiques de la physique en basse dimension. Topological aspects of low dimensional systems" (pp. 53175). Springer Berlin Heidelberg (1999).
 [20] Y. Kodama, K. Kokubu and N. Sawado, Phys. Rev. D 79 (2009) 065024.
 [21] A.Yu. Loginov, JETP , 118 (2009) 217.
 [22] C. Adam, T. Romanczukiewicz, J. SanchezGuillen and A. Wereszczynski, Phys. Rev. D 81 (2010) 085007.
 [23] J. Gladikowski, B. M. A. G. Piette and B. J. Schroers, Phys. Rev. D 53 (1996) 844.
 [24] C. Adam, T. Romanczukiewicz, J. SanchezGuillen and A. Wereszczynski, JHEP 1411 (2014) 095
 [25] J. Ashcroft, M. Haberichter and S. Krusch, Phys. Rev. D 91 (2015) 10, 105032
 [26] R.S. Ward, Nonlinearity, 17 (2004) 1033.
 [27] I. Hen and M. Karliner, Nonlinearity 21 (2008) 399.
 [28] J. Jäykkä, J.M. Speight and P. Sutcliffe, Proc. Roy. Soc. Lond. A 468 (2012) 1085.
 [29] P. Salmi and P. Sutcliffe, J. Phys. A 48 (2015) 3, 035401
 [30] P. Salmi and P. Sutcliffe, arXiv:1511.03482 [hepth].
 [31] Y. M. Shnir, Phys. Part. Nucl. Lett. 12 (2015) 4, 469.
 [32] C. Adam, C. Naya, J. SanchezGuillen and A. Wereszczynski, Phys. Rev. D 86 (2012) 045010.
 [33] T. Weidig, Nonlinearity, 12 (1999) 1489.
 [34] B. M. A. G. Piette, B.J. Schroers and W.J. Zakrzewski, Nucl. Phys. B 439 (1995) 205.
 [35] R.A. Leese, M. Peyrard and W.J. Zakrzewski, Nonlinearity 3 (1990) 773.
 [36] P. Sutcliffe, Nonlinearity 4 (1991) 1109.
 [37] M. Hale, O. Schwindt and T. Weidig, Phys. Rev. E 62 (2000) 4333
 [38] Ya. Shnir and G. Zhilin, Phys. Rev. D 89 (2014) 105010.

[39]
L.D. Faddeev, Quantization of solitons, Princeton preprint IAS75QS70 (1975)
L.D Faddeev and A. Niemi, Nature 387, 58 (1997); Phys. Rev. Lett. 82, 1624 (1999).  [40] M. Kobayashi and M. Nitta, Nucl. Phys. B 876 (2013) 605
 [41] E. Babaev, J. Jaykka and J.M. Speight, Phys. Rev. Lett. 103 (2009) 237002.