Gravitating superconducting solitons in the (3+1)-dimensional Einstein gauged non-linear σ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma $$\end{document}-model

In this paper, we construct the first analytic examples of (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} -dimensional self-gravitating regular cosmic tube solutions which are superconducting, free of curvature singularities and with non-trivial topological charge in the Einstein-SU(2) non-linear σ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma $$\end{document}-model. These gravitating topological solitons at a large distance from the axis look like a (boosted) cosmic string with an angular defect given by the parameters of the theory, and near the axis, the parameters of the solutions can be chosen so that the metric is singularity free and without angular defect. The curvature is concentrated on a tube around the axis. These solutions are similar to the Cohen–Kaplan global string but regular everywhere, and the non-linear σ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma $$\end{document}-model regularizes the gravitating global string in a similar way as a non-Abelian field regularizes the Dirac monopole. Also, these solutions can be promoted to those of the fully coupled Einstein–Maxwell non-linear σ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma $$\end{document}-model in which the non-linear σ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma $$\end{document} -model is minimally coupled both to the U(1) gauge field and to General Relativity. The analysis shows that these solutions behave as superconductors as they carry a persistent current even when the U(1) field vanishes. Such persistent current cannot be continuously deformed to zero as it is tied to the topological charge of the solutions themselves. The peculiar features of the gravitational lensing of these gravitating solitons are shortly discussed.


Introduction
Topological defects are formed in phase transitions when a system goes from a state of higher symmetry to a state of lower symmetry. They can be classified as local or global depending on the fact if a local or global symmetry is broken. Topological defects occur in very different areas of physics like, for example, condensed matter physics, high energy physics, and cosmology. In condensed matter physics perhaps the most famous and intuitive example is the formation of domain walls in ferromagnetic materials, which separate domains with different magnetization. Another famous example is the formation of vortex lines in superfluid helium and line defects (dislocations) in crystals. A topological defect can sometimes be completely regular in which case it is known as a topological soliton. They play a fundamental role in quantum field theory, nuclear physics, and high energy physics [1,2].
The first example of stable topological soliton in three space dimensions was proposed by Skyrme and it is known as Skyrmion [3][4][5]. It has the remarkable property of possessing Fermionic excitations despite the fact that the dynamical field is an SU (2)-valued scalar field. At leading order in the 't Hooft large N expansion [6][7][8][9], the Skyrme model represents a (phenomenologically successful) low energy description of Quantum Chromodynamics (QCD).
The Skyrme model arose from a clever modification of the non-linear σ -model (NLSM henceforth), which is the low energy description of the dynamics of Pions (for nice reviews see e.g. [10,11]). The Skyrme term was added to the NLSM in order to avoid Derrick's no-go scaling argument [12] preventing the existence of static soliton solutions of finite energy. On the other hand, it should be kept in mind that the elegant arguments in [13,14] (see also [15][16][17][18], and references therein) to show that Skyrmions represent Fermions (at least semiclassically) are only based on the existence of stable solitons with non-trivial third homotopy class while they do not use directly the Skyrme term in itself. Thus, it is extremely interesting to search for alternative ways to avoid Derrick's scaling argument in order to achieve "Fermions out of Bosons" with the simplest possible ingredients. The main approaches to avoid the no-go scaling argument in [12] are, first of all, to minimally couple the NLSM to gravity and/or to Maxwell theory. Indeed, using the techniques developed in [19][20][21][22][23][24][25][26][27][28][29][30][31][32][33][34][35], exact self-gravitating NLSM solutions with non-trivial topological charge have been found in [36,37] and [38]. Secondly, it is very helpful to construct time-dependent ansatz for the SU (2)-valued matter field with the property that the energymomentum tensor is time-independent (this idea is the SU (2) generalization of the Bosons star ansatz for a U (1)-charged scalar field: see [39,40] and references therein). Such a generalization has been achieved in [24,25,30,34] and [35].
Further relevant topological solitons which play an important role in high energy as well as condensed matter physics, are the Abrikosov-Nielsen-Olesen vortex line [41,42] and the 't Hooft-Polyakov monopole in the SU (2) Yang-Mills-Higgs system [43,44]. It is worth to mention that the latter at large distances looks like a Dirac monopole, however, the Dirac monopole, which describes a point-like magnetic charge, is singular at the origin whereas the 't Hooft-Polyakov monopole is regular at the origin.
The formation of topological defects are very important in grand unified theories as the actual symmetry group of the standard model is supposed to be a result of a series of spontaneous symmetry breaking of a larger symmetry group. This means that topological defects play a fundamental role from microscopic scale to extremely large scale namely in cosmology due to the fact that the universe in its evolution expanded and cooled down, and therefore went through several phase transitions were topological defects have formed. As in cosmology, the most relevant interaction is gravity, the topological defects must be studied in the context of field theories coupled to gravity.
In cosmology the topological defects which attracted the most attention of the scientific community are the cosmic strings. The simplest exact cosmic string solution is given by an energy-momentum tensor concentrated in a line [45] (for example the z axis) The exact solution of the Einstein field equations associated with this energy-momentum tensor, written in cylindrical coordinates, is locally but not globally flat as the range of the angular coordinate θ is not, as usual, 2π but where G and μ are the gravitational constant and the mass of the cosmic string per unit length, respectively. Therefore, this locally flat space-time has a conical defect, an angular deficit = 8Gμ and has a curvature singularity on the z axis. A possible way to smooth out the singularity is to smear the energy-momentum tensor on a cylinder of finite radius δ and it is possible to find exact solutions [46][47][48].
The principal problem of this procedure is that there is a sharp boundary whose radius is arbitrary and moreover the energy-momentum tensor is not derived from some fundamental action principle.
In order to find a cosmic string solution from a fundamental action principle usually the Einstein-Yang-Mills-Higgs action is used, but unfortunately no exact solutions are known. An explicit gravitating global cosmic string metric (where the fundamental field is a Goldstone boson instead of the Higgs and Yang-Mills fields) was found by Cohen and Kaplan, and has the peculiarity that the matter distribution has no sharp boundary and the angular defect varies with the distance from the symmetry axis [49]. This metric has a curvature singularity at a finite distance from the axis. It was shown that non-singular gravitating global strings can exist if there is an explicit time dependence in the metric [50,51]. For a nice explanation of the properties of the Cohen-Kaplan global string see e.g. [52, pages 196-198].
The existence of cosmic strings has many important cosmological and astrophysical implications. For example, cosmic strings have been proposed to have a role in the galaxy formation as a source of density perturbations [53]. Cosmic strings have also observable effects through gravitational lensing being the most known effect the formation of double images [45,54,55]. It is worth to point out that the space-time generated by a thin string (with a Dirac delta matter source) does not exert force on a test particle being locally flat, but the existence of a conical defect still generates double images. Perhaps one of the most fascinating aspects of cosmic strings is that, under certain conditions, they become superconducting as it was shown in the pioneering article of Witten [56]. This can have observable effects as such superconducting strings would act as sources of synchrotron radiation or high energy cosmic rays [57]. Moreover, it has also been proposed that superconducting strings moving in a magnetized plasma can be a mechanism for the production of gamma-ray bursts [58].
Due to the many important cosmological and astrophysical implications of cosmic strings, it is of great interest to find analytic non-singular solutions that can be derived from some fundamental action principle which leaves no arbitrariness in the choice of fields and their potentials. Indeed, as it has been explained before, in the case of the Einstein-Yang-Mills-Higgs action no such exact solutions are known. The known exact solutions are the thin string with Dirac delta matter source and the global string, both of them possess curvature singularities which would go against the cosmic censorship conjecture. An important point is how to choose a fundamental matter field. Here we will consider the (gauged) NLSM as it is an effective low energy description of the dynamics of Pions (as well as of their electromagnetic properties).
In this paper, we construct the first examples of analytic and singularity free cosmic tube solutions for the self-gravitating SU (2)-NLSM. At large distance the metric behave in a similar way to the one of a cosmic string boosted in the axis and has an angular defect related to the parameter of the theory, however near the axis the parameter of the solution can be chosen in such way that the solution is free of singularities and without angular defect. This means that these solutions are free of singularities everywhere and the angular defect depends on the distance from the axis. It is also worth pointing out that the matter field does not have a sharp boundary and the curvature reaches its maximum on a tube around the axis rather than on the axis itself. All these features make the new solutions similar to a global string with the big difference that they are regular instead of having a singularity at finite distance from the axis. The cosmic tubes found here are related to global strings in the same way as non-Abelian monopole are related to Dirac monopoles. In other words, the NLSM regularizes the global string keeping a similar behavior at large distances. Another relevant feature of the solutions presented here is that these possess non-trivial topological charge (the third homotopy class). It is important to mention that in [59,60] and [61], string solutions in the Skyrme model with mass term have been constructed. These configurations were obtained using a static ansatz that facilitates the resolution of the field equations numerically, but that leads to the strings possessing a zero topological density. Those strings, therefore, are classically unstable and are expected to decay into Pions. On the other hand, the temporal dependency in the U field that we introduce, as we will see before, allows us to construct analytical and topologically non-trivial solutions by circumventing Derrick's theorem (as we have previously pointed out) and without having to include the Skyrme term or additional potentials. Our solutions, having a topological charge, can not decay in the configurations found in [59,60] and [61]. Furthermore, these tubes can be promoted to full solutions of the Einstein-Maxwell-NLSM in which the NLSM field is minimally coupled both to the U (1) field as well as to gravity. These gauged solutions carry a persistent current even when the U (1) gauge field is zero and therefore becomes superconducting in the sense of [56]. In particular, the superconducting currents are tied to the topological charge so that they cannot be deformed continuously to zero (that is why they are persistent). It is worth to emphasize that the gravitating solitons constructed in the present paper only involve degrees of freedom arising from low energy QCD minimally coupled with General Relativity (GR) without the need of additional potentials.

Comparison with existing literature
Here it is interesting to compare our novel results with some of the results already existing in the literature. In particular, we will refer to the nice references [59][60][61][62][63][64][65][66][67][68][69]. Here below we list these four important new features of our solutions (for the sake of clarity, we will call these novel features A, B, C and D).
(A) First and foremost, our regular and topologically nontrivial superconducting tubes are also gravitating. Namely, we couple them to the Einstein equations since one of our main goals is to determine the peculiar features of the gravitational field generated by superconducting tubes. The reason is that, as it was already clarified in Witten's original reference, the spectacular observational effects of superconducting tubes appear in situations in which the gravitational field should not be neglected. Moreover (for the reasons explained here below in point B) we are willing to construct such gravitating superconducting tubes analytically rather than numerically. On the other hand, many of the relevant examples of superconductive tubes available in the literature are numerical and found in flat space-times (such as the examples in the very nice references [59][60][61][62][63]68]).
(B) The second main novel feature is that our regular gravitating solitons in a sector with non-vanishing topological charge are completely analytic: to the best of our knowledge, all the most relevant examples available in the literature are numerical. 1 An obvious question is: why should one insist so much on finding analytic solutions if these equations can be solved numerically? Indeed, very powerful numerical techniques are available in the literature to analyze gravitating superconducting tubes (as, for instance, in [59][60][61][62][63][64][65][66][67][68][69]). There are many sound reasons to strive for analytic solutions even when numerical techniques are available (especially when, as in the present paper, one is working in (3 + 1)-dimensions). Here below we will list two of them which are quite relevant for the present analysis. First of all, as general motivation, it could be enough to remind all the fundamental concepts that we have understood thanks to the availability of the Schwarzschild and Kerr solutions in GR and of the non-Abelian monopoles and instantons in Yang-Mills-Higgs theory. Much of what we now know about black hole physics in GR and instantons and monopoles in gauge theories arose from a careful study of the available analytic solutions. Hence, a systematic tool to construct analytic gravitating superconducting tubes can greatly enlarge our understanding of these very important topological defects.
Secondly, the analytic knowledge of these configurations greatly simplifies the gravitational lensing analysis (indeed, it is worth emphasizing that gravitational lensing is the main phenomenological tool with which one can discover such objects).
(C) The third novel feature of our gravitating superconducting tubes is that they are constructed entirely in the low energy limit of QCD minimally coupled both to GR and to Maxwell theory. 2 Indeed, at leading order in the large N limit, the SU (2) NLSM is the low energy limit of QCD describing the dynamics of Pions (and also of Baryons when the topological charge is non-vanishing as in our manuscript). Thus, in particular, our solutions are physically very different from the ones in the two nice references [62,63] (which are closely related to the electroweak sector of the standard model). Our solutions instead can be very relevant in situations in which Hadrons interact strongly both electromagnetically and gravitationally. Needless to say, there are many relevant situations from the astrophysical and cosmological point of view where both the electromagnetic and the gravitational interactions with Hadrons cannot be neglected (as it happens, for instance, in the physics of neutron stars). Hence, our solutions are not only of academic interest but they can play an important role in phenomenologically relevant situations.
(D) The fourth novel feature of our gravitating superconducting tubes is that the persistent character of the superconductive current is extremely transparent. In particular, we have taken great care to build an ansatz with the following properties. Firstly, we want a consistent ansatz with a nonvanishing topological charge (which, in our case, is interpreted as the Baryonic charge of the configuration) because configurations with non-vanishing topological charge cannot decay into the trivial configurations (neither, in fact, into configurations with different topological charges). Secondly, we want a consistent ansatz in which the electromagnetic current is protected by topology (for instance, the numerical solutions constructed in the nice references by Nitta-Shiki and Jackson have vanishing topological charge [59][60][61]). Thirdly, we need a time-dependent ansatz (in order to avoid the Derrick scaling argument: see [70,71] and [72].) which is, at the same time, compatible with a stationary metric and carries a non-trivial topological charge.
In other words, the electromagnetic current associated with our configurations is proportional to the topological density of the configurations themselves. The physical consequence of this fact is that our current cannot be deformed continuously to the trivial vanishing current as this would imply also a jump in the topological charge. However, there is no continuous deformation that can change the topologi- 2 The existence of superconducting strings in the low energy limit of QCD has been discussed in the nice works [67,68] and [69] cal charge. It is worth to emphasize that the energy scale of GUT symmetry breaking is 10 15 GeV while the non-linear σ -model considered in the present manuscript is motivated as a low-energy effective theory of QCD, is relevant at much lower energy scales. On the other hand, the superconducting nature of the gravitating solitons to be constructed here below may strongly affect the scattering of electromagnetic waves leading to potentially visible effects. The reason is that gravitational lensing does not take into account the very strong interactions of the electromagnetic field with the matter field. Unfortunately, we cannot be precise yet since the problem of the scattering of electromagnetic waves in these gravitating solitons background must be solved numerically. We hope to come back to this interesting point in a future publication. Therefore, our superconducting currents are entirely protected by topology and they cannot decay; these are genuine persistent (and, therefore, superconductive) currents. Moreover, in the original Witten's reference, superconducting currents can only appear when suitable inequalities on the parameters of the Higgs potential are satisfied. In our case instead we do not need any potential at all; it is just the topology of the low energy limit of QCD which keeps the currents alive forever without any extra ingredient.
The structure of the paper is the following: in Sects. 2 and 3, we present the model, give our ansatz and the corresponding field equations. In Sect. 4 we construct the exact regular cosmic tube solutions of the Einstein SU (2)-NLSM and show that they possess non-trivial topological charge. Also we study the geodesic equations and their main physical properties. In Sect. 5 it will be shown how to promote the found configurations to be solutions of the gauged Einstein NLSM system and how these are superconducting configurations. In Sect. 6 we discuss the flat limit of our configurations. In the last section, our conclusions are detailed.

The Einstein SU(2)-NLSM
The Einstein-NLSM theory is described by the action where R is the Ricci scalar and L μ are the Maurer-Cartan Here κ is the gravitational constant and the positive coupling K is fixed by experimental data. In our convention c =h = 1 and Greek indices run over the four dimensional spacetime with mostly plus signature. In this paper, we follow a standard convention that the Riemann curvature tensor, the Ricci tensor, and the Ricci scalar are given by respectively. In order to produce a correct physical interpretation of the topological solitons constructed here and compare with the cosmic string solutions already existing in the literature, in this work we will refer to configurations with vanishing cosmological constant, however the techniques used here are also effective in presence of a cosmological constant. We hope to come back on this interesting issue in a future publication.
The complete Einstein-NLSM equations read where G μν is the Einstein tensor and T μν the energymomentum tensor of the NLSM given by The winding number of the configurations reads When the topological density ρ B is integrated on a spacelike surface w B represents the Baryon number. We will only consider configurations in which ρ B = 0.
To construct analytical solutions in this theory we will use the generalized hedgehog ansatz (see [19][20][21][22][23][24][25][26][27][28][29][30][31][32][33][34][35]), which is defined as where t i ≡ iσ i , with σ i are the Pauli matrices. It should be noticed that the Maurer-Cartan forms are adapted for the process of solving the field equations, since they satisfy the SU (2) algebra. The choice of the ansatz in Eq. (8) is adapted to the commutation relations of the Maurer-Cartan form in such a way to simplify as much as possible the corresponding Einstein equations (as it will be discussed in a moment). A necessary (but, in general, not sufficient) condition in order to have non-vanishing topological charge is where α, and are the three scalar degrees of freedom appearing in the standard parametrization of the SU (2)valued scalar field defined in Eq. (8).

Ansatz and field equations
The ansatz in Eq. (8) is defined over the Weyl-Lewis-Papapetrou metric Since the metric determinant of the section spanned by t and z is negative definite (and e −2R is definite positive as we will show in the next sections), the space-time with this metric is always Lorentzian regardless of the metric components in the t − z. As it will be discussed in the next sections, in order to have an integer value of the topological charge, one must allow the whole range of real number as the domain of the coordinate X : and θ is an angular coordinate with range According to [34,35] we will consider a matter field in the form In Appendix A we have given a detailed explanation of how the above ansatz can be constructed in a systematic way and why it is so effective in the case of the gauged gravitating NLSM. Note that the ansatz defined here above satisfies the necessary condition in Eq. (9) in order to possess a nontrivial topological charge. The sufficient conditions will be discussed in the following sections. This ansatz is very useful for (at least) three reasons. The first one is because Eq. (13) implies the relations which simplify greatly the NLSM equations. The second reason is that Eq. (13) allows to avoid the Derrick's scale argument [12] as it is a time dependent ansatz which, however is compatible with a stationary metric. Thirdly, the three coupled field equations for the NLSM in the metric defined in Eq. (10) reduce to the single second order ODE for α(X ): It is important to note that this second order equation for α(X ) can be reduced to the following first order equation where E 0 is an integration constant. However, the compatibility with the Einstein equation requires that E 0 = 0. This situation is different from what happens in flat space-time, where the integration constant E 0 can be non-zero, see [34] and [35]. Quite remarkably, the Einstein equations with the energymomentum tensor corresponding to the NLSM configuration defined in Eq. (13) are reduced to only two solvable equations: where Note that once R and α are solved, G can be obtained directly. In fact, Eq. (17) is nothing but a flat linear Poisson equation in two dimensions in which the source term is known explicitly (as α(X ) and R(X ) have been determined in Eqs. (14) and (16)). Hence, Eq. (17) can be solved, for instance, using the method of Green's function. In the next sections we will construct the solution using a direct method. At this point, it is worth to emphasize that the function G depends explicitly on θ . As the coordinate θ plays the role of an angular coordinate going around the Hadronic tube, the fact that G depends on θ implies that the present family of gravitating solitons is not axi-symmetric. Indeed, our space-time admits two Killing vectors ∂ t and ∂ z so that there does not exist any closed curve generated by a spacelike Killing vector. Moreover, a direct computation reveals that the Killing equation associated to the present class of metrics in Eqs. (10) and (13) has no further solution other than the two solutions ∂ t and ∂ z : this means that (due to the presence of the function G(X, θ)) the metric is manifestly not axi-symmetric as, otherwise, ∂ θ would be a Killing field as well.
The physical role of the function G will be discussed in the analysis of the curvature invariants and of the geodesics.

Gravitating tubes
In this section we construct exact regular cosmic tube solutions of the Einstein-NLSM and we study its relevant physical properties.

Solving the system
Eqs. (14) and (16) can be solved analytically, and the expressions for α(X ) and R(X ) are given by where C 1 , C 2 and C 3 are integration constants. Replacing the above in Eq. (17) we obtain the final equation for G, At this point it is important to emphasize that considering only Eqs. (8), (10), (13), (18) and (19) the complete Einstein-NLSM system has been reduced to a single equation for the metric function G given in Eq. (20). The energy density (measured by a co-moving observer) of these configurations is given by Note that we should impose the constraint on the integration constants to avoid divergence of the energy density at X → ±∞.

An analytical solution
Using that sin 2 (x) = (1 − cos(2x))/2, Eq. (17) becomes The function G can be expressed by the sum of two functions G 1 (X ) and G 2 (X, θ) which satisfy The first equation is trivially solved by a double integral with respect to r , We can solve the equation for G 2 by separation of variables, obtaining for some real constant y, and ψ 1 (X ) satisfying Equation (28) can be solve in terms of Elliptic Functions using the method of variation of parameters. The solution of the homogeneous equation is On the other hand, a particular solution can be found through here W (ψ (1) h , ψ (2) h ) = −4q is the Wronskian, and Therefore, the solution of Eq. (28) is the sum of ψ h and ψ p , namely

Regularity
Obviously, the regularity of the metric must be discussed by analyzing the corresponding curvature invariants. Taking into account Eqs. (18) and (19), the Ricci scalar S, the Kretschmann scalar, the square of the Ricci tensor and of the Weyl tensor for the metric in Eq. (10) read and are all regular everywhere (in particular, S is regular at X → ±∞) if we impose the following condition on the integration constants: One may wonder whether the above condition in Eq. (30) is enough to ensure the regularity of the metric. Indeed, it is possible to compute explicitly the main fourteen curvature invariants that are usually considered in the literature to analyze, in four dimensions, the issue of regularity (see [73,74]). Such invariants are The regularity condition at X → ±∞ is the same as in Eq. (30). One may wonder if there exists a divergent higher differential invariant of our space-time, such as R ;μν··· R ;μν··· . A potential source of such divergences is the presence of an anisotropic pressure that is not C 1 [75]. The energymomentum tensor of the cosmic tube, however, is smooth everywhere so that such a divergence cannot occur. Thus, all the curvature invariants which can be built from the metric in Eqs. (10) and (19) are regular everywhere if the condition in Eq. (30) on the integration constants holds. Interestingly enough, the function G(X, θ) is not relevant at all as far as the regularity of the curvature invariants is concerned. Thus, in particular, the curvature invariants do not depend on θ . This implies that all the curvature invariants have a local regular maximum at finite distance from the axis of the tube (as it will be discussed in the next sections, in this radial coordinate r the axis of the tube is located at X → −∞ while the peaks of the curvature invariants are at X = 0 which is at finite distance from the axis). Since, as it has been already emphasized, the curvature invariants do not depend on θ , the regular peaks of the curvature invariants lie on a ring in the (X − θ ) plane. Such a ring (which is a tube from a three-dimensional perspective) is also the support of the superconducting currents associated with the present class of gravitating topologically non-trivial solitons.
One can also verify that the space-time is of Petrov type II. One more constraint on the integration constant C 2 will arise from the analysis of the geodesic in the following sections.

Periodicity of U and the topological charge
In this subsection we will discuss a simple but deep property of the SU (2) valued scalar field U which has a very important consequence. First of all, one can notice that when with n an integer, the topological charge is non-zero while, when q is an integer, the topological charge vanishes. This can be seen as follows. The topological density corresponding to the NLSM configuration in Eq. (13) reads Thus, as the range of θ is [0, 2π ], the integral of the above density is non-zero if and only if Eq. (31) holds. It is also worth to note here that the z coordinate goes along the axis of the tube and along the topological density. Thus, in a sense, the can be interpreted as the topological charge per unit of length of the tube, with L z = dz.
Here it is worth to emphasize that, when the condition in Eq. (31) holds, θ is a proper angular coordinate with range [0, 2π ]. First of all, the metric itself is periodic with period 2π as it depends on θ only through the function G(X, θ). The function G(X, θ) depends 3 on θ only through the factor cos(2qθ) in Eq. (17). Secondly, the energy-momentum tensor of the SU (2)-valued scalar field is also periodic in θ with the same period 2π when the condition in Eq. (31) holds.
However, the matrix U itself is not periodic in θ with the same period 2π as it reads On the other hand, one can be sure that the "lack of periodicity" can be compensated by an internal Isospin rotation since the energy-momentum tensor corresponding to the above NLSM configuration is periodic with the correct 2π −period: indeed, if such "lack of periodicity" could not be compensated by an internal rotation, then the energy-momentum tensor would not have the correct period. Consequently, θ is a proper angular coordinate with range [0, 2π ] when q is halfinteger. Thus, as it happens with the spin-from-Isospin effect, the internal symmetry group plays a fundamental role. In that case, the spin-from-Isospin effect is generated by the possibility to require "spherical symmetry up to an internal rotation". In the present case, the condition that U satisfies periodic boundary conditions up to an Isospin rotation is enough to ensure that T μν is periodic and w B = 0.

Coordinate transformation and the asymptotic behavior
To analyze the nature of the metric in Eq. (10), it is convenient to make the following coordinate transformation: It is obvious that r (X ) increases monotonically as X increases, since Moreover, to make r (X ) well-defined, it is necessary to impose that This assumption gives us the range of r given by With this coordinate transformation, the metric becomes, As we will show now, the function e −2R(r ) of the new cylindrical radial coordinate r both for r close to zero and for r → ∞ is proportional to r 2 : This implies that θ is an angular coordinate. Consequently, it is very important to determine the coefficients 0 and ∞ which determine the effective deficit angles seen from observers "very close to" and "very far from" the axis of the tube, respectively.
Note that we have the following two limits of the function R; so that we also have In the limit of X X + for a sufficiently large X + 0, we find that Here, we have defined a finite constant and in the last line, we used the assumption that C 2 +K κ|q| < 0. Thus, we obtain In a similar way, the limit of X X − for a sufficiently small X − 0 is found to be Here, we also have defined a finite constant and in the last line we have used the fact that X X − 0. Therefore, we get Now, we see that we can choose C 3 such that the angular deficit is 1 near the axis defined by r (X ) = 0 so that θ becomes a proper angular coordinate. The ratio of the values of g θθ at two infinities becomes, When we assume that G dies out at spatial infinity, as supported by the plot given in Fig. 2, the asymptotic form of the metric is found to be which describes a boosted cosmic string at spatial infinity.

Geodesics
Let λ be an affine parameter of geodesics in our space-time.
The geodesic equations are found to bë where and the dot denotes the derivative with respect to λ. From the first two equations, we obtain for some constants a and b. Then, Eqs. (47), (48), and (49) becomë On the hypersurface of constant z, the Eq. (47) becomes from which we obtain Then, Eq. (46) becomes For convenience, let's put Then, the remaining geodesic equations becomë Using Eq. (57) in a linear combination Ẋ × (54) +θ × (55) yields Thus, for some constant C 4 , we havė We should assume that to avoid infinity velocities at large X . The line element for the geodesic with constant z is Divided by the affine parameter t, it can be written as By an appropriate rescaling becomes −1 or 0 for time-like or null geodesic, respectively. Using Eq. (63), we have Thus, the line element of the geodesic with constant z gives a relation between constants as follows Using Eq. (63) in the geodesic equation of r , from Eq. (60) one finds thaẗ After solving this equation, one should plug the solution to the equation of motion for θ in Eq. (61) to find θ(t). This completes the solving process of the geodesic equations with constant z. In principle, this problem can be reduced to an effective one-dimensional Newtonian problem observing that, along the geodesics with constant z, one haṡ where G(X, θ) is defined in Eqs. (24), (26) and ( 29). The reason is that from Eq. (66) one can determine θ = θ (X ) along the geodesics with constant z. Once θ = θ (X ) has been determined, one can insert it into Eq. (63) obtaining a first order Newtonian-like equation of the forṁ However, Eq. (66) is a quite complicated first order nonautonomous differential equation for θ (X ) due to the explicit form of G(X, θ) in Eqs. (24), (26) and (29). The relevant issue of the behavior of geodesics in this family of gravitating solitons deserves a more detailed analysis on which we hope to come back in a future publication.
In the usual case of cosmic strings, the angle defect of gravitational lensing is independent of the initial distance of the particle from the source. But in our case, the distribution of source is smoothly spread out to X = ∞, so that the angle defect depends on X 0 , the initial location of geodesic motion. This would be one of the most distinguished properties of our solution.

Constraint on integration constant
Combining the regularity conditions on the integration constants in Eqs. (22), (30), (34) and (64) we obtain the following single inequality Thus, when C 2 satisfies the above condition and C 3 is chosen as in Eq. (43) all the curvature invariants of the metric are regular and the geodesics behave in a reasonable way. Note that the experimental value of the Pions coupling constant is such that 0 < K κ 1. On the other hand, the integration constant C 1 is a free parameter which fixes the location of the point about which the profile function α(X ) is symmetric.

Gauged gravitating tubes
In this section we promote our configurations to be solutions of the Einstein-Maxwell-NLSM system, and also show how these are superconducting gravitating tubes.

The Einstein-Maxwell-NLSM theory
The action of the U (1) gauged Einsten NLSM theory is The respective field equations are where the current J μ is given by andT μν is the energy-momentum tensor of the Maxwell theorȳ Note that in the action in Eq. (68) there is a quadratic term in A μ . The topological charge in this case [76] is given by where Note that the second term in Eq. (75), the Callan-Witten term, guarantees both the conservation and the gauge invariance of the topological charge.

Gravitating tubes coupled to the electromagnetic field
Let's consider the following Maxwell potential together with the metric defined in Eq. (10) and the matter field in Eq. (13). This is a very convenient ansatz [34,35], because a direct computation reveals that the three field equations for the gauged NLSM reduce (once again) to α − q 2 2 sin 2α = 0, so that α(X ) = 2 arctan(e q X+C 1 ).
On the other hand, the four Maxwell equations reduce to just one linear equation: Note that this linear equation for = ω s − 2u can be easily solved, at least numerically, since α(X ) is explicitly known and R(X ) can be also determined explicitly. In fact, there are two non-trivial Einstein equations that can be combined to obtain an uncoupled equation for R, namely R − 1 2 K κ(α 2 + q 2 sin 2 α) = 0, so that, as in the case without Maxwell, and G satisfies the following equation Hence, once again, the function G satisfy a flat twodimensional Poisson equation in which the source is explicitly known. Therefore, once the Maxwell equation in Eq. (78) has been solved, the function G can be determined explicitly using several methods (such as the Green functions). Resuming, using the ansatz for U and A μ in Eqs. (13), (76) and for the metric in Eq. (10), the Einstein-Maxwell-NLSM field equations reduce simply to Eqs. (78) and (80), where α(X ) and R(X ) are in Eqs. (77) and (79). It is worth to emphasize that this is a quite remarkable reduction of such a coupled system of non-linear PDEs in a sector with nontrivial topological charge as the NLSM (which is already a very complicated theory in itself) is coupled directly both to GR and to the U (1) gauge field A μ .
The energy density measured in an orthonormal frame is found to be where we used the tetrad given bŷ The topological charge density, including the Callan-Witten term, is The current is given by and the components of the electric and magnetic fields read The plots here below as well as the review of the Witten construction clearly show why these solutions of the Einsten-Maxwell-NLSM represent regular topologically non-trivial gravitating superconducting tubes. In order to plot the relevant physical functions we have set our parameters as From Fig. 1 and Fig. 2 we can see that even if the metric is not axi-symmetric due to the explicit θ -dependence of G(X, θ) in Eq. (80), all the curvature invariants are axisymmetric as they only depend on X . Thus, the plots of all the curvature invariants are very similar and they all show a smooth peak at finite distance from the origin (remember that, in the coordinate X , the origin is at X → −∞). As we expected, e −2R(r ) goes as r 2 when r → 0 as well as when r → ∞.
In Fig. 3, the energy density associated to a comoving observer has two parts. The first one only depends on "X " and, as the curvature invariants, has its smooth maximum at the same finite distance from the axis. The second part depends both on "X " and on "θ " and the corresponding peak is at the same distance from the axis as the peak of the first term and, in θ is localized around θ ∼ π . The peaks of the energy density and the peaks of the topological density coincide, as expected. From the plots in Fig. 4 one can see that J μ = 0 only at the position of the two peaks and tends to zero out. Here we have imposed the following boundary conditions,

Review of the Witten argument
Before discussing the superconducting nature of the present solutions, we will shortly review the results in [56] (which have been considerably generalized in many subsequent works: see for instance [77][78][79][80][81][82][83], and references therein).
The main motivation to introduce such topological objects is related, of course, to the spectacular observable effects that such objects could have (were they to exist: see the original  Fig. 2 The metric function G as a function of X and θ reference [56] as well as [84][85][86][87][88][89][90][91][92][93][94][95][96][97][98][99][100][101] and references therein). The second motivation is related to the fact that such remarkable objects can be constructed using quite reasonable ingredients. Many of the examples available in the literature do not use exclusively building blocks within the standard model. 4 For instance, extra U (1) gauge potential as well as Higgs-like scalar fields are often important ingredients while in [77][78][79][80][81][82][83], supersymmetry plays a fundamental role.
The starting point of [56] is the following Lagrangian 5 The above theory has an U (1) × U (1) gauge symmetry (the first corresponding to A μ and the second to S μ ). In order for the above theory to support superconducting strings, it is necessary to include an interaction potential between the two Higgs fields σ and ψ. The choice of [56] was The reasons behind this choice are the following: The first necessary ingredient is the breakdown of the gauge symmetry corresponding to S μ in order to ensure the existence of vortices. The Higgs field ψ in the core of the vortex field is usually assumed to only depend on the two spatial coordinates (say, r and θ ) transverse to the vortex axis (which is along the z-axis). Then, one must require that, in the vacuum, σ = 0. At this point, with a clever choice of the range of the parameters of the Higgs potential, one can achieve the following situation. Despite the fact that the (minimization of the) kinetic energy tends to suppress σ = 0 within the core, if one chooses m 2 to be positive, the potential energy will favor σ = 0 within the core. As it was shown in [56], this can indeed happen. In other words, there is an open region in the parameter space in which σ = 0 asymptotically but σ = 0 within the core of the vortex associated to ψ. This is a fundamental technical step since the superconducting currents (to be described in a moment) are sustained by the region in which σ = 0. If σ 0 (r, θ) minimizes the energy of the string, then the superconducting current is associated with the (slowly varying) phase of σ 0 (r, θ). One can achieve this by introducing the dependence on z and t in σ as follows: The expression of the current is which is made of two factors. The first factor − → ∂ + e − → A is responsible for the dynamics of the zero modes along the strings (associated to the phase ). However, such a factor Fig. 3 The Energy density T00 and topological density ρ B as functions of X and θ Fig. 4 From left to right, the vector density plots of the Electric field E, the magnetic field B and of the current J μ by itself would be "useless" as it needs to "rely on" something. This "something" is the factor σ . Thus, first of all, the first factor needs σ to be different from zero somewhere. For superconducting strings, as it has been already emphasized, the spatial region where σ = 0 is tube-shaped. However, this is not enough: the configurations in which σ = 0 within a tube-shaped region must be stable, otherwise the current would decay. 6 In the settings of [56,[77][78][79][80][81][82][83], the linear stability of the configurations where σ = 0 within the string (as well as σ approaching to zero outside) were established by direct methods (such as linear perturbation theory). Note also that the current defined above cannot have arbitrarily large values since σ has a maximum value determined by the Higgs potential. Once the stability of such tube-shaped regions (which are going to host the superconducting currents) has been established, one can ask: Under which circumstances the above current is superconducting?
One needs a mechanism that keeps such current perpetually alive even in the absence of an external gauge potential. At this point, topology comes into play. Since is only defined modulo 2π , the integral over a close loop of the above current will not vanish in general even when one turns off the electromagnetic field. Indeed, one can build a topological invariant associated to (the integral of) . Thus, such a cur-rent cannot relax when the topological invariant associated to is non-zero. Basic ingredients A very nice pedagogical construction has been achieved in [105] where, assuming that symmetry breaking is always generated by suitable Higgs or scalar potentials, the author reduced the basic ingredients to the skeleton. The author required that bulk theory should have an unbroken global non-Abelian symmetry (let us assume, just for concreteness, that such a group is SU (2) ) allowing the existence of string-like configurations. Moreover, on such string-like configurations, one has to break SU (2) down to a subgroup U (1). As it has been already emphasized, such requirements are usually taken care of by introducing suitable potentials for the scalars. In fact, the Einstein-Maxwell-NLSM has all the above ingredients already "built-in" and there is no need to introduce any potential: the typical interactions among the Maurer-Cartan forms associated to the Isospin degrees of freedom do the job.
Secondly, it would be very nice to use topology not only to ensure the persistent character of the current but also to guarantee the appearance of regions where σ = 0 (such that σ approaches zero outside). Higher topological charges can be useful, as we will discuss in the following subsection.

The superconductivity of the tubes
The main features of the above expressions in Eq. (82) for the current of these topologically non-trivial gauged crystals are the following.
(1) The current does not vanish even when the electromagnetic potential vanishes (u = 0).
(2) Such a "current left over" J (0)μ : (where has been defined in Eq. (13)) which survives even when the Maxwell field is turned off, is maximal where the energy density is maximal (namely, where sin 2 (α) sin 2 (qθ) = 1, which defines the positions of the peaks in the energy density as well as in the topological density) and vanishes rapidly far from the peaks.
(3) Such residual current J (0)μ cannot be turned off continuously. This can be seen as follows. There are three ways to "kill" J (0)μ . The first way is to deform α to an integer multiple of π (but this is impossible as such a deformation would change the topological charge). The second way is to deform qθ to an integer multiple of π (but also this deformation is impossible due to the conservation of the topological charge). The third way is to deform to a constant (but also this deformation cannot be achieved). Note also that is defined modulo 2π (as the SU (2) valued field U depends on cos and sin rather than on itself). This implies that the line integral of ∂ μ along a closed contour does not necessarily vanish (as it happens in the original Witten argument).
The above characteristics show that the above residual current is a persistent current that cannot vanish as it is topologically protected. This, by definition, implies that J (0)μ defined in Eq. (89) is a superconducting current supported by the present gauged tubes.

The flat limit
In this section we discuss the flat limit κ → 0 of our solution describing a gauged gravitating soliton. In this limit, the field equations in (14), (16) and (78) become The solutions to the second and third equations are found to be which corresponds to the flat limit of the solution in Eq. (79). When C 2 does not vanish, the flat limit of the coordinate transformation to the radial coordinate of a cylindrical system in Eq. (32) gives so that the space-time metric becomes To avoid a conical singularity, we should impose that We have taken the negative value of C 2 to satisfy the inequality (67). In this cylindrical coordinate system, the equations for α and u become 2u(r, θ)) sin 2 qθ (e |q|C 3 +C 1 r |q| + e −|q|C 3 −C 1 r −|q| ) 2 = 0.
The constraint on the parameters (67) restricts the range of the parameter q It should be noticed that |q| = 1 is forbidden in the flat spactime limit if C 2 = 0. The flat limit of the energy-density is found to be whereê 0 = ∂/∂t. It follows from Eq. (98) that the energy density (100) of the gauged soliton decays rapidly, whereas the slow radial decay of the scalar field in the Cohen-Kaplan solution results in the formation of a singularity at a certain distance r = r max from the core of the cosmic string [49]. It should be noticed that r max is not the location of the maximum of the energy density (that is close to the string core) but a maximal radius in a non-singular part of spacetime, including the core of the Cohen-Kaplan string. 7 Since we are seeking a regular solution for the topological solitons, it should be possible to expand u(r, θ) near the axis as for some constant m and some smooth function f (θ ). It follows from Eq. (100) that the energy-density of the system diverges at r = 0 unless m ≥ 1. Near the axis, Eq. (98) at leading order is found to be This equation admits a solution for an arbitrarily small r only when q = 0, which is forbidden by (99). For this reason, the flat limit of our solution for the gauged soliton is singular if we assume that C 2 = 0. Therefore, the genuine flat limit of the gravitating superconducting tubes described in the previous sections is obtained when the integration constant C 2 in Eq. (94) vanishes. In this way we recover the configurations of topologically non-trivial gauged solitons within the space-time metric given by which was studied in [35] (but with a difference, which we will see immediately). In order to shed further light on the subtleties of the flat limit, it is worth to remind that the analytic solutions in [35] represent gauged superconducting tubes in a box with finite volume described by the flat metric in which the three coordinates (r , θ and φ) have finite range in order to describe a situation with a finite amount of Baryonic charge within a finite spatial volume. Thus, at least locally, one can identify the coordinate X with the coordinate r in [35]. The ansatz in [35] gives rise to Eq. (90) for the profile α which can be easily reduced to the following first order equation where E 0 is an integration constant. As explained in details in [35], the role of the integration constant E 0 is very important since it allows to "squeeze" an arbitrary amount of Baryonic charge within the finite flat box 8 and it is fixed through the relation It is also important to remember that in the flat case studied in [35] the Maxwell equations are reduced to a single equation in the form with = 2L p u − 1. On the other hand, there is a subtle fingerprint of the Einstein equations which does not disappear in the flat limit. The reason is that the profile α corresponding to the (flat limit of the) gravitating hadronic tubes constructed in the previous sections corresponds to a very precise choice of the integration constant E 0 in Eq. (105). It is easy to see, comparing Eqs. (14) and (15) with the equation for α in the flat case in [35], that the flat limit of the present gravitating superconducting tubes corresponds to the case E 0 =0, namely The fact that E 0 vanishes has two important consequences. Firstly, it admits only one peak of the energy density along the direction of X as could be expected. Secondly, the solution for α necessarily describes a kink from −∞ to +∞ . The finite range of the coordinate θ is consistent with the fact that, when the flat superconducting tubes in [35] are promoted to gravitating superconducting tubes as described in this manuscript, the coordinate θ becomes an angular coordinate. Hence, the flat limit of the gravitating solitons turn out to be a single superconducting tube in a box that has been extended to infinity in the X -direction.
The analysis of the flat solutions in [35] also clarifies why the gravitating version of the flat gauged solitons presented in this manuscript is not axi-symmetric. The reason is already manifest in the (flat limit of the) energy-density: (where we have used the freedom to scale some integration constants to 1). Clearly, the above expression (which is the energy density of the flat limit of the present gravitating hadronic tubes) is not axi-symmetric. The obvious reason is the explicit dependence of the energy-density on θ (which, in the gravitating case, plays the role of angular coordinate). Thus, in a sense, it is natural that the gravitating version of the flat solitons in [35] should not be axi-symmetric. We see also that T00 decays rapidly in X as expected. Unlike what happens with the Cohen-Kaplan solution [52], our configurations have a maximum at a finite distance while in Cohen-Kaplan this is a delta, being the difference exponentially small. In fact, we can say that our solution approaches asymptotically to the Cohen-Kaplan solution.
The above arguments show that our solution exists for an arbitrarily small κ. As can be seen from the curvature invariants given in Sect. 4.3, the gravitational effects on the curvature would be small in the weak-field limit. It should be emphasized that, however, the most spectacular consequence of the gravitating tube is the large deviation of the worldlines of physical particles (such as photons and electrons) from the geodesics. In particular, photons will not follow null geodesics due to the effective mass term arising due to the minimal coupling with the NLSM. On the other hand, electrons will not follow time-like geodesics since they will be accelerated by the Lorentz force. Both effects will be especially strong close to the local maxima of the energy-density. Consequently, the observer at infinity will detect a sudden and large deviations of these particles from the geodesic motion in an almost flat space-time.

Conclusions
The first example of analytic and curvature singularity free cosmic tube solutions for the Einstein SU (2)-NLSM has been found. The metric at large distance from the axis looks similar to a boosted cosmic string. The matter distribution has no sharp boundary and the curvature is concentrated at a finite distance from the axis. The angular defect of the solutions depends on the distance from the axis and the parameters of the solutions can be chosen in such a way that it vanishes near the axis but not at large distance. These properties make the solution similar to a global string but with the fundamental difference that while the global string has a curvature singularity at a finite distance from the axis whereas the new solutions are singularity free.
Due to the non-Abelian symmetry group, the most natural way to impose a periodicity condition on the SU (2) field is up to an inner space rotation. This very natural boundary condition allows the solution to carry non-trivial topological charge.
One of the most remarkable aspects of these solutions is that they can also be promoted to solutions of the Einstein-Maxwell-NLSM theory, i.e. the SU (2) field is minimally coupled to both the U (1) field and gravity, and without neglecting the corresponding Maxwell equations "sourced" by the currents arising from the NLSM. The gauged solutions are characterized by the fact that they can carry a persistent current even when the Maxwell field is zero, which means that they are superconducting. Moreover, the superconducting current is also topologically protected.
It is worth pointing out that in cosmology one of the most important observational consequences of the existence of cosmic strings is gravitational lensing. The overwhelming majority of papers dealing with gravitational lensing assume axi-symmetry of the cosmic string. The analytic solution found here however is not axi-symmetric and therefore the geodesic equation becomes non-trivial. An interesting feature of the solution found here is that it is highly repulsive in the core of the tube, the matter distribution has no sharp boundary and therefore spreads to infinity, the deflection angle depends on the initial distance from the source. It is reasonable to suppose that these non-trivial features should have interesting observational consequences and will be an object of study in further investigations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .

Appendix A: How the ansatz solves the system
In this first Appendix, we provide all the technical details behind the ansatz that we have used in the main text to solve the system of coupled field equations corresponding to the Einstein-Maxwell-NLSM together with the technical steps needed to derive the field equations.
As it has been already mentioned in the main text, the action of the NLSM and the field equations read, respectively where In order to achieve a better understanding of the present framework and of why our ansatz works it is convenient to express the above action explicitly in terms of the most general parametrization of SU (2): U (x μ ) = (cos α) 1 + (sin α) n j t j , t j = iσ j , n 1 = sin cos , n 2 = sin sin , n 3 = cos . (A3) Here the three SU (2) Pionic degrees of freedom are encoded in the three scalar functions α = α(x μ ), = (x μ ) and = (x μ ). Consequently, the action of the NLSM in Eq. ( A1) can be expressed explicitly in terms of the three scalar functions α, and . A direct computation shows that the result is the following action (which, of course, is completely equivalent to the one in Eq. (A1)): It is very convenient also to express the field equations in Eq. (A2) in terms of the three scalar degrees of freedom. A direct computation shows that, when varying the action in Eq. (A4) with respect to α, and , the following field equations are obtained It is important to emphasize that the field equations in Eqs. . The next piece of information needed to build the ansatz is the topological charge density As we want to consider only topologically non-trivial configurations, we must demand that Now, the problem is to find a good ansatz which respects the above condition and simplify as much as possible the field equations. A close look at Eqs. (A5), (A6) and (A7) reveals that a good set of conditions is Moreover, one should not forget the Derrick's no-go theorem which prevents the existence of static solitonic solutions of the NLSM: thus, the ansatz should depend on time but in such a way to have a stationary energy momentum tensor (as expected from a superconductive tube). A choice that satisfies Eq. (A9) and, at the same time, circumvent Derrick's scale argument is where the space-time metric 9 and the coordinates r , θ , t and z have been defined in the main text in Eq. (10). Some other useful identities satisfied by the above ansatz are together with Eq. (A8). Such identities greatly simplify the field equations. A direct computation reveals that the three NLSM field equations Eqs. (A5), (A6) and (A7) are reduced to the following ODE for α: keeping alive, at the same time, the topological charge.
Remarkably enough, these very intriguing properties of the ansatz are not destroyed by the inclusion of the minimal coupling with Maxwell field. The coupling of the NLSM with the Maxwell theory is introduced replacing, in the action in Eq. (A1), the partial derivatives acting on the SU (2) -valued scalar field U with the following covariant derivative A straightforward computation shows that the above replacement in Eq. (A13) in the Lagrangian in Eq. (A1) is completely equivalent to the replacement here below (in terms of α, and ) in the action in Eq. (A4). It is worth to emphasize that D μ determines the "direction" of the electromagnetic current (as it will be discussed here below). Thus, in this case the field equations for the gauged NLSM are where, for notational simplicity, we will still define as Obviously, when the derivative is replaced with the Maxwell covariant derivative (as defined in Eq. (A13) or, equivalently, in Eq. (A14)), in the field equations of the gauged NLSM many new terms will appear, coupling the SU (2) degrees of freedom with the U (1) gauge potential. Thus, one may ask: Which is the best choice of the ansatz for the gauge potential A μ that keeps as much as possible the very nice properties of the ansatz of the SU (2)-valued scalar field in Eqs. (A8), (A9) and (A11), and allowed (as explained in the main text) the complete analytic solutions in the Einstein-NLSM case?
In order to achieve this goal, it is enough to demand The above conditions determine that the Maxwell potential A μ must be of the form From the above, the explicit expressions of the components L μ and J μ can be evaluated directly as i sin 2 α sin 2 (qθ) e −i(z+tω s ) sin α sin(qθ)(cos α − i sin α cos(qθ) −e i(z+tω s ) sin α sin(qθ)(cos α + i sin α cos(qθ) −i sin 2 α sin 2 (qθ) , iq sin α cos α sin(qθ) , From the expressions for L μ one can see that, despite the explicit presence of A μ in the U (1)-covariant derivative, the three field equations of the gauged NLSM still reduce to Eq. (A12). The reason is that all the potential terms which, in principle, could couple the SU (2)-valued scalar field U with A μ in the field equations in Eq. (A15) actually vanish due to the identities in Eqs. (A11), (A16) and (A17) satisfied by the choice of our ansatz (that is why we have chosen the ansatz in that way). Moreover, due to the presence of a quadratic term in A μ in the action in Eqs. (A1) and (A4), which couple A μ with the SU (2)-valued scalar field U (as it happens in the Ginzburg-Landau description of superconductors), even when A μ = 0 the current does not vanish (as it is clear from Eq. (A19)). Such a residual current (which survives even in the A μ = 0 limit) cannot be deformed continuously to zero: the reason is that the only way to "kill" would also kill the topological charge but, as it is  (17) and (80).

Appendix B: About the stability in the flat limit
A remark on the stability of the superconducting tubes in flat space-time constructed in [35] (which in the case E 0 = 0 is reduced to the flat limit of the gravitating tube) is in order.

Perturbations on the profile
In many situations, when the hedgehog property holds (so that the field equations reduce to a single equation for the profile) the simplest non-trivial perturbations are those of the profile which keep the structure of the hedgehog ansatz (see [109,110] and references therein). One of the key technical results of the present work is that such a property holds even when the minimal coupling with GR is considered. In particular, it holds in the flat case (as shown in [34,35]). In the present case the simplest non-trivial perturbations are of the following form: which do not change the Isospin degrees of freedom associated with the functions and . The NLSM equations at first order in ε are reduced to α − q 2 2 sin(2α) + (ξ − q 2 cos(2α)ξ )ε = 0.
We can see that, imposing the on-shell relations we get to the following equation for ξ(X ): ξ − q 2 cos(2α)ξ = 0, so that the system always has the following zero-mode: ξ (X ) = ∂ X α 0 (X ), with α 0 a solution of the NLSM equation.
Due to the fact that the integration constant E 0 in Eqs. (105) and (106) can always be chosen in such a way that ∂ X α 0 (X ) never vanishes (see [35]), the zero-mode ξ (X ) has no node so that it must be the perturbation with lowest energy. Thus, the present solutions are stable under the above potentially dangerous perturbations. This is a very non-trivial test. Now we discuss the relation between the stability and the existence of a topological charge. The on-shell Hamiltonian density H obtained evaluating the ansatz in Eqs. (A3) and (A10) in Eq. (A2) turns out to be H on-shell = K 2L 2 (α 2 + q 2 sin 2 (α)).
Consequently, the on-shell energy E of the system can be written as E on-shell = π 2 K L (α 2 ∓ 2q sin(α)α d X, which turns out to be a positive definite term plus a total derivative. The integral of this total derivative defines a new topological charge Q(n), namely Q(n) = 2q 2π 0 (sin α) α d X = 2q cos α(0) 1 − (−1) n , which depends on n. It is worth to note that the above topological charge Q(n) is different from the Baryonic charge as the Baryonic charge is non-vanishing for any integer n appearing in the boundary conditions (see [35]) while the above topological charge is non-vanishing only when n is odd. The corresponding BPS bound for the on-shell energy is, One can see that there exist configurations which saturate the bound: such configurations satisfy α ± q sin α = 0.
Interestingly enough, also in this case the saturation of the above BPS bound does imply the general field equations. Thus, the present analytic solutions have two topological labels: the Baryonic charge |w B | = |np| and Q(n). Obviously, the configurations which saturate the above bound are stable. This happens precisely when the integration constant E 0 in Eqs. (105) and (106) is such that the solution for α is also compatible with the Einstein equations (see Eq. (16)). The present stability argument also extends to the gravitating configurations constructed in the main text which have |n| = 1.

Electromagnetic perturbations
Another very useful approach to study the stability is to consider only electromagnetic perturbations. For simplicity (see Eqs. (71), (72) and (107)), let us consider the following type of electromagnetic perturbations: (u, 0, 0, −Lu) → (u + εc 1 , 0, 0, −Lu + εc 2 ), c i = c i (t, X, θ, φ), ε 1. In our case, to ensure that the perturbed Maxwell equation can be solved we need to demand that V eff > 0, with V eff = ( m 2π ) 2 − V . Since V takes its maximum value when sin 2 (α) sin 2 (qθ) = 1, we get to the following relation We can see that there is a maximum value of the size of the box ∼ 10 f m.
The complete stability analysis requires to study the most general perturbations of the solutions of the present work, but this is a very hard task even numerically as it involves a coupled system of linear PDEs. We would like to revisit the perturbative stability of these solutions in a future publication.
Note that the only change that needs to be done with respect to classic procedure of [18] is the use of the light-cone time (instead of the usual time-like coordinate in [18]). The reason is that in [18] the authors considered the following type of Isospin perturbations of the original static Skyrmions U 0 introduced by Skyrme himself: For the same reason, in the present case is more convenient to choose the light-cone time.
A perturbation of the type defined in Eqs. (B6), (B7 ), (B8) and (B9) is called Isospin perturbation because of the following reason: Therefore, the perturbations of the SU (2) valued field U defined in Eqs. (A3) and (A10) that we are considering in the present subsection are fundamentally different from the perturbations in Eq. (B1) since the latter only changes the profile α(X ) while they keep the functions and unchanged while the perturbations in Eqs. (B6), (B7), (B8) and (B9) do not change the profile while they only affect the functions and . Since we are studying linear perturbations, it makes sense to analyze these two different types of perturbations separately.
It is direct but rather cumbersome to check that the perturbations on the functions and that follows the criteria of [18], this is, in the case with t-dependence as well when they depend on the light-cone coordinate τ = t/L − φ, are in both cases zero-modes of the system, meaning, the energy of the system does not decrease.