Higgs phase in a gauge $\mathbf{U}(1)$ non-linear $\mathbf{CP}^1$-model. Two species of BPS vortices and their zero modes

In this paper zero modes of fluctuation are dissected around the two species of BPS vortices existing in the critical Higgs phase, where the scalar and vector meson masses are equal, of a gauged $\mathbb{U}(1)$ nonlinear $\mathbb{CP}^1$-model. If $2\pi n$, $n\in \mathbb{Z}$, is the quantized magnetic flux of the two species of BPS vortex solutions, $2n$ linearly independent vortex zero modes for each species are found and described. The existence of two species of moduli spaces of dimension $2n$ of these stringy topological defects is thus locally shown.


Introduction
CP N models in two-dimensional Euclidean space-time were introduced by Golo and Perelomov, see [1], with the goal of discussing self-dual 2D-instantons in a complex/Khäler framework. These systems in Minkowski R 1,3 space-time are a variant of the O(N ) non-linear sigma model, which is constructed in turn as a restriction of the linear sigma model of Gell-Mann and Levy [2] to the dynamics of their Goldstone bosons. Several generalizations of the non-linear O(N )-sigma model appeared as low energy effective theories in hadronic physics related to current algebras before the surge of non-Abelian gauge theories in particle physics. In this work we shall focus in a model belonging to the family of gauged non-linear Sigma models, see e.g. [3], which might be treated also in interaction with gravity, see [4], granting to this kind of systems a rôle in Cosmology.
In particular, we shall concentrate in the simplest model studied in Reference [5] where we added to the action of the gauge CP 1 -model a potential energy density depending only on the fields and guaranteeing that both the scalar and vector particles are massive. Moreover, we showed that in this model there are two species of BPS vortices living respectively in the south and north charts. Because our vortex solutions have cylindrical symmetry they can be interpreted as cosmic strings conveying by themselves cosmological implications if the model is defined in curved space-time, see inter alia the monographic textbook [6] and References quoted therein.
Besides of letting the system to live in the Higgs phase, our choice of the potential energy density is guided by the existence of a self-dual structure: first-order static field equations together with a topological bound. This structure arises at the critical point where the scalar and vector particle masses are equal and it is the requirement to build generalizations of these gauge non-linear sigma models with extended supersymmetry. In fact, in the disguise of 2D instantons Nitta and Vinci discovered in [7] the same topological defects in the framework of two-dimensional N = (2, 2) SUSY sigma model that we described within a purely bosonic context in [5]. The promotion of the field theoretical model at the stake to a N = 2 supersymmetric status is possible through an standard procedure using the Kähler structure of CP 1 .
Our goal in this essay, however, is to analyze thoroughly the BPS vortex zero modes of the two species. Despite of being in a Higgs phase our model exhibits zero gap fluctuations around the BPS vortices. In Reference2 [8,9] we developed this task in a fully detailed manner for the BPS vortices in the Abelian Higgs model. The information that we acquired in this work was the key to the computations achieved in [10] of the one-loop BPS vortex string tension shifts in the AHM when these objects are immersed in the quantum world. We plan thus to establish the analogies and differences between the zero modes of fluctuation of the standard self-dual vortices and those of the two species of BPS vortices in our gauged CP 1 model. Because zero modes provide a bridge between the classical and quantum domains we hope that the results obtained here will be useful in future research about quantum properties of the two species of BPS CP 1 vortices.
2 A gauge U(1) massive non-linear CP 1 -sigma model with two self-dual vortex species We shall address a system of three scalar fields Φ = (Φ 1 , Φ 2 , Φ 3 ) which take values on a S 2 -sphere target space: We also denote x µ · x µ = g µν x µ x ν and choose the metric tensor in R 1,3 as: g µν = diag(1, −1, −1, −1). Stereographic projections respectively from the north and south poles of the sphere to the plane lead to the south and north charts that form a minimal atlas in S 2 intersecting at the equator. In these two planes/charts we define the complex scalar fields distinguished by the superscripts S and N . Use of complex coordinates characterizes the target S 2sphere as the complex manifold CP 1 , i.e., the complex projective line. The information coming from one of these charts is translated into the other one via the transition function φ S = ρ 2 /φ N * that allows the global definition of the CP 1 by prescribing how the two charts are related. We remark that our choice of transition function reverses the orientation of the north chart with respect to that of the south chart. The reason for choosing this option is to deal later with scalar fields coupled to the gauge field with identical electric charges in both charts. The massless Lagrangian L 0 describing the dynamics of the CP 1 model written in terms of the south chart fields, reads: The global U(1)-symmetry with respect to changes in the field phase φ S → e iχ φ S may be promoted to local invariance following the standard procedure: a gauge field A µ enters the system and supplements the local U(1) transformation φ S → e iχ(x µ ) φ S with the gauge transformation A µ → A µ +∂ µ χ. A potential energy density yielding spontaneous breaking of the gauge symmetry in the Higgs mechanism mode will be also allowed. All this leads to the Lagrangian of a gauged massive Abelian non-linear CP 1 model with dynamics governed by the Lagrangian: The covariant derivatives and the electromagnetic tensor are defined as usual: In order to write the Lagrangian (1) in the north chart we need to know how the covariant derivatives and the potential energy density read there: Thus, in the other chart the Lagrangian becomes: a formula showing, together with (1), that this model is globally defined on CP 1 . For the sake of concision we shall write the Lagrangians (1) and (2) in the unified form where φ stands for either φ S or φ N if either the south or the north charts are alternatively considered in the description of the model. In both charts the metric factor in the target manifold which appears multiplying the kinetic terms respectively in (1) and in (2) is written in the form: Notice that the Abelian Higgs model Lagrangian follows the generic form (3) when the metric factor g(|φ| 2 ) is the unity and there is only one chart in the non-compact target space C. We recall that the celebrated Abrikosov-Nielsen-Olesen vortex filaments arise in the AHM as cylindrically symmetric solutions of the field equations grown from planar topological solitons endowed with a quantized magnetic flux. The ANO vortices are static configurations, henceforth their investigation requires the temporal gauge A 0 = 0. Axial symmetry is realized in two steps: first, one searches for field configurations living in the x 3 = 0 plane: φ = φ(x 1 , x 2 ) and A j = A j (x 1 , x 2 ), j = 1, 2, a restriction that is consistent only in the axial gauge A 3 = 0. Second, in the restricted second-order field equations one looks for solutions with this characteristics and appropriate boundary conditions, which force the quantization of the magnetic flux. Repeated infinitely in the third dimension the ANO magnetic flux tubes are obtained from these planar vortex solutions having finite energy per unit length. If the parameters of the AHM are such that the system lives at the critical point between Type I and Type II superconductivity, there exists a first-order PDE (Bogomolny) system of field equations such that the vortex solutions, usually called self-dual or BPS vortices, saturate a topological bound, i.e., their energy per length unit is proportional to the magnetic flux.
Closely following the developments in [5], we summarize the investigation of BPS vortex solutions starting from the Lagrangian (3), which covers the two species of BPS vortices in our gauged massive Abelian CP 1 model. From this Lagrangian one easily derives the energy per unit of length functional for static and planar configurations: We then focus on the set of all the static planar configurations with finite energy per unit length: Membership to C space demands the following asymptotic behavior: where M is the set of zeroes of the non negative function U (|φ| 2 ) ≥ 0. Moreover, we shall restrict ourselves to a subclass of models governed by the Lagrangian (3) for which the energy per unit length functional (4) admits a Bogomolny arrangement. It was proved in Reference [5] that V [φ, A] is the sum of positive (perfect squares) terms plus a (positive) topological quantity: when the scalar field potential term U (|φ| 2 ) is of the form having into account that g(|φ|) > 0. Notice that G(0) = 0, G (z) = g(z) because of the definition in (7) and we choose 0 < α 2 < G(∞) in order to leave room for spontaneous symmetry breaking of the U(1) symmetry. Solutions of the first-order PDE system are thus minima of the energy per unit length and therefore also solve the second-order field equations. Moreover, the boundary conditions (5) ensure that the vector field A j r→∞ iφ * ∂ j φ is asymptotically purely vorticial, and, via Stoke's theorem, that the total magnetic flux is a topological quantity given by the degree of the map from the circle at infinity S 1 ∞ to the manifold M provided by the asymptotic behaviour of the scalar field. BPS or self-dual vortices are static and x 3 -independent field configurations that solve the PDE system (8) subjected to the boundary conditions (5). Thus, the BPS vortex solutions saturate the Bogomolny bound V [φ BP S , A BP S ] = α 2 πn, where n is an integer, becoming the configurations of minimum energy per unit of length in each disconnected sector of the configuration space. The BPS magnetic flux lines are accordingly stable minima of the functional V (φ, A) which do not exhibit fluctuation modes of negative energy. Saturation of the Bogomolny bound, however, is compatible with flat directions in configuration space C of neutral equilibrium. In [5] we applied the index theorem unveiling the existence of 2n zero modes of fluctuations around the new species of vortices. This means that the quanta of magnetic flux are, at least locally, free of moving independently from each other through the moduli space of vortices.
Following the same strategy as in Reference [8] on the research for vortex zero modes in the AHM we shall first investigate the planar PDE system (8) by using polar coordinates in the x 1 -x 2 plane, i.e., we shall explore circularly symmetric, or, cylindrically symmetric seen in three dimensions, BPS vortices. The remnant gauge freedom allows us to set the radial gauge condition A r = 0 and we shall plug the ansatz φ(r, θ) = f n (r) e inθ , rA θ (r) = n β n (r) , n ∈ Z (10) into the PDE system (9). We end with the non linear ODE system determining the radial profiles f n and β n of the scalar and vector field components of the self-dual vortices. From the boundary conditions (5) we derive the asymptotic behavior of these radial profiles : where G −1 denotes the inverse function to G. Because the vacuum orbit M = {v : |v| 2 = G −1 (α 2 )} is a circle in the target manifold, the scalar field at infinity φ(∞, θ) = G −1 (α 2 ) · e inθ provides a map from the boundary of the plane at infinity to M which belongs to the n-th homotopy class of the first homotopy group of the circle: Π 1 (S 1 ) = Z. The winding number n of this map is also encoded in the fact that the magnetic flux is classically quantized: We do not only require the asymptotic behavior (12) to the vortex solutions of (11), but we demand also regularity at the origin, which implies β n (0) = 0 in (10). But near the (0, 0, x 3 ) axis, the ODE system (11) can be solved analytically because if β n (0) = 0 the equation for the radial profile becomes: dfn dr ≈ ± n r f n (r), which is solved by where d n is an integration constant. We obtain next β n (r) near the origin by solving the other linearized equation in (11): n r To derive this linearization we used the radial profile (13) near the origin and the power series expansion of the G[f 2 n (r)] function around r = 0: G[f 2 n (r)] = g(0)d 2 n r 2n + o(r 2n+1 ) up to the r 2n+3 -th order. Solving for β n (r) the linear ODE (14) is an easy task: The vortex radial profiles written in (13) and (15) clearly show that f n (0) = β n (0) = 0, i.e., both the scalar and the vector vortex fields vanish at the vortex center (x 1 , x 2 ) = (0, 0) as demanded. The behavior in the intermediate region varies with the choice of potential U (|φ| 2 ), henceforth with the metric g(|φ| 2 ) of the target manifold. Since the ODE system (11) is in general not analytically solvable, a numerical scheme is usually used to obtain solutions in the whole radial range compatible with the prescribed behavior near r = 0 and close to r ∞.
We come back now to the non-linear CP 1 -sigma model characterized at the beginning of this Section. The compatibility relation (7) between metric in the target space and potential fixes the scalar field self-interaction in the south chart to be: If α ∈ (0, 2ρ) a spontaneously broken symmetry scenario arises from the structure of the potential (16). The family of degenerate vacua, which in the south chart coordinates reads constitutes a parallel circle on the S 2 -sphere target manifold. The mass of the vector particle is the coefficient of the quadratic term A µ A µ in the Lagrangian after choosing one point in the vacuum circle and expanding the scalar field near this point, Simili modo, expanding the purely scalar field contribution to the Lagrangian in the unitary gauge where the spurious Goldstone boson is set to zero, ϕ S 2 = 0, up to the quadratic order in the Higgs field φ S 1 (x µ ) one finds that the scalar particle mass is also: At least in the south chart our model lives: 1) In a Higgs phase, all the particles have mass. 2) At the critical point between Type I and Type II superconductivity, scalar and vector particles have the same mass.
Coordinates in the north chart are related to those in the south chart by means of the (orientation reversing) transition function: φ N = ρ 2 φ * S . In the norh chart the potential (16) thus reads: The vacuum orbit seen in the north chart is the parallel circle: Because |v N | 2 = ρ 4 · |v S | −2 one realizes that it is the same parallel as (17) but seen from the south pole. Moreover, the masses of the scalar and vector particles in the north chart are the same as in the south chart and the model is globally defined in target space.
The vacuum orbit parallel M splits the sphere into two spherical caps, each cap entirely belonging to either the south or the north chart. This fact is behind the existence of two species of self-dual vortices. Each species takes values only on one of these two caps sharing only the parallel M as common boundary. We recall also that stereographic coordinates describing the vortex scalar field Φ vanish at the vortex center. The two species of self-dual vortices in the gauged CP 1 -model are classified as follows: 1. The south vortex species: If the S 2 -valued Φ field describing the vortex solution points downwards at the vortex center, Φ 3 (0) = −1, only the south chart enters the game. The first-order ODE system (11) determining the south vortex species reads where the obvious notation φ S (r, θ) = f S n (r)e inθ has been used.
2. The north vortex species: If the vortex configuration in Φ-space valued at the origin points upwards, Φ 3 (0) = 1, only the north chart plays a rôle. The self-duality ODE system solved by the north vortex species corresponds to derived from (11) by taking the north chart fields and the appropriate function G.
Although the ODE systems (22) and (23) characterize two different species of vortex solutions, they exhibit related features. Notice that the potentials describing the scalar field interactions in each chart are connected by the relation U N (|φ| 2 ; α 2 ) = U S (|φ| 2 ; 4ρ 2 − α 2 ) where the unified notation has been again used. This identity implies that the functional form of the first-order ODE systems (22) and (23) is symmetric under the interchange of the parameters α 2 and 4ρ 2 − α 2 ; simply look at the right members of the ODE's on the right in (22) and (23). These two parameters are inversely related, the greater the first one is, the lesser the second one becomes since α 2 ∈ (0, 4ρ 2 ). The values of α 2 and 4ρ 2 − α 2 have a direct impact in the magnitudes of the derivatives of the function f and β as we can check from (22) and (23). The asymmetry between the values of α 2 and 4ρ 2 − α 2 causes the difference between the core sizes of the vortices belonging to different species.
Recall that the vortex core is characterized by small values of the scalar field. In Figure 1, below we have displayed the field profiles of these two species of vortices with vorticity (or winding number) n = 3 for α = 1 and ρ = 1, tantamount to 4ρ 2 − α 2 = 3.
In the first and second graphics in Figure 1 the scalar and vector field profiles of the 3-vortex belonging to the south species are respectively depicted by means of arrow plots. The shadowed discs are centered at the zeroes of the scalar and vector fields. In the third and fourth graphics analogous arrow plots for the north species of the 3-vortex are shown. The components of the arrows specify the real and imaginary parts of the complex scalar fields. We observe in these figures that the winding number of the two species of 3-vortex profiles is three: note that in any quadrant every arrow rotates a 3 2 π angle. On these figures we have plotted density representations of the modulus of the scalar field. The darker the colour is the lesser the value of the modulus is in such a way that the shadowed region in the scalar field graphics represents the vortex core. Notice that for our choice of α and ρ the south species vortices are thicker than their north species partners. The same routine has been followed for the representation of the vortex vector field profiles of the two species, see Figures 1-(2) and 1-(4).

Zero mode fluctuations around cylindrically symmetric BPS vortices of the two species
We shall devote this Section to investigate the zero modes (eigenmodes with null eigenvalue) which exist in the spectrum of cylindrically symmetric vortex fluctuation operators of the two species. In [5] the existence has been shown of 2n linearly independent zero modes of fluctuation around BPS vortices living alternatively in one of the two charts by using the Callias-Bott-Weinberg index theorem, see [11,12,13]. These eigenmodes are the key ingredient to describe the adiabatic motion of vortices in their moduli space, see e.g. [14], and therefore they are the cornerstone for studying the effective vortex dynamics in low energy scenarios.

Analytical description of the zero modes of fluctuation of cylindrically symmetric BPS vortices
In order to distinguish between the vortex solution and its fluctuations we shall denote the scalar and the vector field profiles corresponding to a self-dual vortex solution of vorticity n as while the scalar and vector field fluctuations of the self-dual vortex solution will be denoted as ϕ( x) = ϕ 1 ( x) + iϕ 2 ( x) and a( x) = (a 1 ( x), a 2 ( x)) respectively. Fluctuations of this type built around the BPS vortex are zero modes of fluctuation if the perturbed fields are still solutions of the equations (8) at first order in . This requirement compels the fluctuation fields a j ( x), j = 1, 2, and ϕ a ( x), a = 1, 2 to comply with the first order partial differential equations − ∂ 2 a 1 + ∂ 1 a 2 + g(|ψ| 2 )(ψ 1 ϕ 1 + ψ 2 ϕ 2 ) = 0 In order to discard pure gauge fluctuations we select the generalized background gauge in the fluctuation space Assembling the field fluctuations in a real four column vector the first-order PDE system of the four equations (24)-(25) may be re-written in the matrix form: Here, D denotes the first-order partial differential matrix deformation operator and the zero mode fluctuations ξ( x) are required also to be square-integrable vectors of a Hilbert space H We are interested in the zero mode fluctuations of cylindrically symmetric n-vortex solutions, and thus we plug the ansatz ψ 1 ( x, n) = f n (r) cos(nθ) , ψ 2 ( x, n) = f n (r) sin(nθ) V 1 ( x, n) = − nβ n (r) r sin θ , V 2 ( x, n) = nβ n (r) r cos θ into the PDE system (24)-(25). The outcome is the equivalent PDE system: where we have used the polar representation for the vector fluctuation field: a 1 = a r cos θ − a θ sin θ and a 2 = a r sin θ + a θ cos θ. In the system (28) a discrete symmetry between the vectors in the kernel of the first-order operator D is hidden: if ξ = (a r , a θ , ϕ 1 , ϕ 2 ) is a zero mode (a solution of the system (28)) then ξ ⊥ = (a θ , −a r , ϕ 2 , −ϕ 2 ) is another zero mode orthogonal with respect to ξ. This symmetry reduces the search for the 2n vortex zero modes of this field theory model to only n of them. We now investigate the analytical properties of these eigenfunctions. Given the cylindrical symmetry of the vortex solutions, we propose an angular dependence of the zero mode wave unctions of the form where k ∈ Z in order to deal with single valued eigenfunctions: ξ(r, θ, n, k) = ξ(r, θ + 2π, n, k). Plugging the ansatz (29) into the system of four first-order PDE's (28) reduces the problem to solve the following system of two coupled first-order ODE's: for the radial form factors s nk (r) and t nk (r) of the vortex zero mode fluctuations defined in (29). Notice that the norm (27) in the Hilbert space of normalizable fluctuations is simplified to ξ( x, n, k) 2 = 2π dr r s 2 nk (r) + g(f 2 n (r)) t 2 nk (r) .
We can go further in our analytical calculations by solving for the radial function t nk (r) in the first equation of (30) and inserting the result in the second equation of (30). This manoeuvre leads to the second-order linear differential equation for the radial form factor s nk (r): In order to analyze the structure of zero mode vortex fluctuations we must study thus the properties of L 2 -integrable solutions of (33). We start this task by analyzing the behavior of the solutions of (33) near the origin and subsequently we shall examine its asymptotic behavior far away from the vortex center: • Regularity of the function s nk (r) at the origin. The Frobenius method may be applied to the linear differential equations (33) at the regular singular point r = 0. We expand the function s nk (r) as a power series where the power s is chosen as the minimum value such that c (n,k) 0 = 0, i.e., h nk (r) is regular and does not vanish at the origin. The norm (31) is now rewritten in terms of the power series h nk (r): by using the relation (32) and (34). Plugging the power series expansion (34) into (33) we obtain the following recurrence relation between the coefficients c (n,k) j where we have used the power series expansion (13) and (15) of the n-vortex radial form factors f n (r) and β n (r) together with the expansion of the metric factor evaluated in the vortex solution near its center: Notice that in (36) the recurrences are cut at order 2n+1, we will see shortly that there is no need of taking into account more terms to ensure regularity at the origin and, henceforth, L 2 -integrability, accounting for only the dominant terms near the vortex center.
For j = 0, 1 the recurrence (36) is simply: Since from hypothesis c = 0. Near the origin, the first summand in the integrand of (35), recall that s = n − k − 1, is therefore: Poles at the origin in the integrand are skipped if a condition which implies that the integer number k is bounded by the vorticity n.
The two term recurrence relations for the next group of indices j = 3, . . . , 2n − 1 and the characteristic exponent s = n − k − 1 becomes: Starting from c (n,k) 1 = 0 it is easily checked that (39) implies c (n,k) 2i+1 = 0 for all the odd indices j = 2i + 1 in the range 3 < 2i + 1 ≤ 2n − 1. The recurrence (39) for even indices, however, j = 2i reads: = · · · = c (n,k) 2k = 0 vanish. c (n,k) 2 is zero because the factor i(k − i + 1) appearing in the left hand side of (40) is non null, while i − 1 present in the right-hand side is zero, for i = 1. If 2 ≤ i ≤ k a similar situation happens: all the right side members in (40) are zero because the coefficients are zero but the left hand sides must be also zero, restricting the values of the coefficients up to c (n,k) 2k to be zero. The first non-null coefficient after c (n,k) 0 is c (n,k) 2k+2 because k − i + 1 = 0 in this case. The first two terms of the s nk (r)-power series expansion near r = 0 are thus: where c (n,k) 0 and c (n,k) 2k+2 are arbitrary non-null constants. Together with the bound (38) this means that it is enough to identify the even coefficients up to c (n,k) 2n+1 in order to describe the zero mode wave functions near the origin, a fact that justify the truncation assumed in the recurrence relations (36). We finally pass to analyze the second summand in the integrand of (35) near the origin: seeing that it is regular at the origin if and only if 2k + 1 ≥ 0, i.e., k ≥ 0. Therefore regularity at the origin restricts the values of k to the first n natural numbers k = 0, 1, 2, . . . , n − 1 such that there are at most n zero modes, or rather 2n, if the orthogonal zero modes to these null potential eigenfunctions are accounted for.
• Asymptotic behavior of the function s nk (r). For large values of r the modulus of the scalar complex field tends to a constant value f (r) → f ∞ > 0 that belongs to the vacuum circle M, whereas the radial profile of the vector field tends to one: β(r) → 1. Bearing this asymptotic behavior in mind, we see that at large r the ODE equation (33) reduces to the modified Bessel differential equation: The general solution of this second-order ODE is well known: where I N [x] and K N [x] are modified Bessel functions respectively of the first and second kind. It is crystal clear that we must choose C 1 = 0 In formula (41) in order to obtain zero mode eigenfunctions with an exponential decaying tail which satisfy the L 2 -integrability condition. 2k+2 in order to obtain a solution with the adequate asymptotic behavior. In conclusion, there exists n zero modes ξ( x; k) of the generic form (29) whose radial profiles s nr (r) and t nr (r) are solutions of the linear first-order ODE system (30). Moreover, all these zero modes characterized by the wave number k are linearly independent.
Together with their corresponding orthogonal partners ξ ⊥ ( x; k), these whole set of 2n zero modes form a basis in the tangent space to the moduli space of BPS vortices at the cylindrically symmetric vortex carrying magnetic flux n and centered at the origin.
Sturm-Liouville theory is enough to ensure the existence of these null eigenfunctions in the intermediate range between a neighborhood of the origin and other one close to the infinite point. Nevertheless, there is no way of analytically finding the vortex solutions at intermediate range, It is possible, however, to gather good information about the BPS vortex zero mode profiles by using numerical methods. In this sense it is better than to directly attack equation (33) for s nk (r) simply to solve by numerical procedures the simpler equation in terms of the function h nk (r). Plugging in (33) we end with the second-order linear ODE −rg[f 2 n (r)] which will be our starting point to generate the zero mode fluctuation by means of the numerical scheme by some variant of a shooting procedure using the known solution near the origin as initial condition.

Deformations of BPS cylindrically symmetric vortices of the two species by their zero mode fluctuations
The perturbed fields up to first-order in the fluctuations ψ( x; n, k) = ψ( x; n) + ϕ( x, k) and V ( x; n, k) = V ( x; n) + a( x, k) are deformed vortex solutions of the general BPS equation (8). Recall that ψ( x; n) and V ( x; n) stand respectively for the cylindrically symmetric self-dual vortex scalar and vector field profiles obtained through the ansatz (10) and the solution of the BPS equations (11). ϕ( x, k) and a( x, k) are the vortex zero mode fluctuations determined by using the ansatz (29) to solve (32), a procedure simplified by assuming (42) to end with the solution of (43) as explained in the previous Section. The main rôle in the description of these perturbed solutions is played by the scalar field profile. In cylindrical coordinates adapted to the symmetry of the vortex flux lines the perturbed vortex complex field profile reads: Sufficiently close to the origin this profile is analytically known in terms of the integration constants d n and c (n,k) 2k+2 plus the value of the metric at the origin g(0): ψ( x; n, k) ≈ r k e ikθ d n r n−k e i(n−k)θ − (2k + 2)c (n,k) 2k+2 Because if r k = 0 and r n−k = (2k + 2)c (n,k) 2k+2 g(0)d 2 n e −i(n−k)θ the perturbed vortex scalar field profile (44) exhibits one zero of multiplicity k situated at the origin. We observe that the multiplicity of the zero at the origin decrease from n to k with respect to that of the unperturbed vortex scalar field. n − k new zeroes appear located at the vertices of a regular k-polygon: n−k and θ j = 2j + 1 n − k π , j = 0, 1, 2, · · · , n − k − 1 .
slightly displaced from the single zero of multiplicity n of the cylindrically symmetric vortex profiles placed at the vortex center. In sum, under the k = n − 1 zero mode perturbation one quantum of magnetic flux moves away from the origin along the x-axis while the remaining n − 1 quanta stay at the origin. Under the next zero mode perturbation k = n−2 two quanta move respectively along the half-axis forming angles π/2 and 3π/2 with the x 1 -axis, the other n − 2 quanta staying at the origin. Under the generic n − k zero mode n − k quanta of magnetic flux depart from the origin along the half-axis forming respectively angles θ j with the x 1 -axis, the other k quanta remaining at the origin. Perturbations of BPS cylindrically symmetric n = 3-vortices in the CP 1 -sigma model undergoing zero mode fluctuations are illustrated in Figures 2 and 3 where the parameters ρ = 1 and α = 1 are chosen. On BPS vortices with n = 3 quanta of magnetic flux there exist six zero modes of fluctuation ξ( x; k) and ξ ⊥ ( x; k), characterized by the "polarizations" k = 0, 1, 2. We recall that in the gauged U(1) CP 1 -model there are two species of BPS cylindrically symmetric vortices taking values respectively in the south and north charts of the target manifold. In Figure 2 the graphical information relative to the zero mode fluctuations ξ S ( x, 3, k) corresponding to the south species of BPS 3-vortices is collected. In the first row the k = 2 zero mode fluctuation ξ S ( x, 3, 2) of the generic form (29) is described. The scalar and vector field profiles ϕ S ( x, 3, 2) and a S ( x, 3, 2) of this null eigenmode ξ S ( x, 3, 2) are respectively depicted in the first and the third graphics by means of a vectorial plot. In the second and fourth pictures of the same row it is plotted the way in which the scalar and vector fields of the circularly symmetric BPS 3-vortex configuration are deformed by the zero mode fluctuation ξ S ( x, 3, 2). One of the three single quanta of magnetic flux superimposed at the origin in the unperturbed solution is displaced along the x 1 -axis while the remainder ones are untouched, see Figure 2 (first row, second plot). The same pattern is shown in the second row where the main features of the k = 2 zero mode fluctuation ξ S ( x, 3, 1) around a cylindrically symmetric BPS 3-vortex belonging to the south chart are graphically described. In the first and third plots of this second row we see vectorial graphs of the scalar and vector fluctuations. Perturbations where two quanta of magnetic flux are ejected from the vortex center are shown in the second and fourth graphs. Finally, the third row in Figure 2 includes the plots corresponding to the k = 0 zero modes ξ S ( x, 3, 0). It is remarkable to notice that the three single vortices initially situated at the origin are expelled in the directions determined by the vertices of a equilateral triangle, see Figure 2 (third row, second plot).
In Figure 3 the same graphical information is collected and offered for the zero mode fluctuations ξ N ( x, 3, k) of the BPS cylindrically symmetric 3-vortex belonging to the north chart species as well as the associated perturbed fields ψ N ( x, 3, k) and V N ( x, 3, k). The zero mode structure in this chart is almost identical to the structure shown in Figure 2 corresponding to the south species of vortices. The only important difference is the fact that the north species of vortices described in Figure 3 exhibits an smaller core than the south chart ones displayed in Figure 2, in agreement with the same discrepancy unveiled in Section 2 between the south and north vortex cores. Figure 2: Graphical representations of the scalar and vector components of the south class 3-vortex zero mode fluctuations ξ S ( x, 3, k) (displayed in the first and third columns respectively) and the perturbed scalar and vector fields ψ S ( x, 3, k) and V S ( x, 3, k) (displayed in the second and fourth columns respectively) for the values k = 2 (first row), k = 1 (second row) and k = 0 (third row).
Despite the similarity between the structures described in Figures 2 and 3 the differences offer a relevant information about how the two species of vortices behave when they suffer a zero mode fluctuation. For the sake of comparison, we have employed the same magnitude of the perturbation parameter in the graphics representing perturbed fields in Figure 2 and 3. It is clear that the effect of the zero mode fluctuation in the north class of vortices is smaller than in the south class. In the context of adiabatic vortex dynamics we can claim that the thick (south class) self-dual n-vortices are split in its single constituents faster than the thin (north class) self-dual n-vortices by zero mode fluctuations, as we can see by comparing Figures 2 and 3. Figure 3: Graphical representations of the scalar and vector components of the north class 3-vortex zero mode fluctuations ξ N ( x, 3, k) (displayed in the first and third columns respectively) and the perturbed scalar and vector fields ψ N ( x, 3, k) and V N ( x, 3, k) (displayed in the second and fourth columns respectively) for the values k = 2 (first row), k = 1 (second row) and k = 0 (third row).

Outlook
In this work we have thoroughly described the zero modes of fluctuation around cylindrically symmetric BPS vortices of the two species existing in a gauge U(1) non-linear CP 1 model that were discovered in Reference [5]. Besides their rôle in the scrutiny of low energy dynamics as achieved e.g. in [14], zero modes have an strong impact in the evaluation of one-loop shifts to classical masses of BPS topological defects. It is conceivable to apply the method developed in [15] to the kinks discovered in Reference [16] in order to improve the results obtained in [17] about the one loop correction to the S 2 -kink masses. In a similar vein it seems highly plausible that the calculations performed on the one-loop string tension shifts of the BPS vortices in the Abelian Higgs model in Reference [10] may be repeated successfully for the two species of BPS vortices in the U(1) gauge non-linear CP 1 model profiting from the results described in this paper. It is also tempting to extend this work to the CP 2 and CP N generalizations of the model treated here. In this case, besides Abelian vortices, semi-local topological defects also appear and one might follow the developments in [18] in this non-linear context.