Radially and azimuthally excited states of a soliton system of vortex and Q-ball

In the present paper, we continue to study the two-dimensional soliton system that is composed of vortex and Q-ball components interacting with each other through an Abelian gauge field. This vortex-Q-ball system is electrically neutral as a whole, nevertheless it possesses a nonzero electric field. Moreover, the vortex-Q-ball system has a quantized magnetic flux and a nonzero angular momentum, and combines properties of topological and nontopological solitons. We investigate radially and azimuthally excited states of the vortex-Q-ball system along with the unexcited vortex-Q-ball system at different values of gauge coupling constants. We also ascertain the behaviour of the vortex-Q-ball system in several extreme regimes, including thin-wall and thick-wall regimes.


Introduction
It is known that (1 + 1)-dimensional gauge models and (2 + 1)-dimensional gauge models without the Chern-Simons term do not allow the existence of electrically charged solitons because any electrically charged compact object will have infinite energy in these models. The reason for this is simple: at large distances, the electric field of a one-dimensional object does not depend on the distance and that of a two-dimensional object is inversely proportional to the distance. As a result, with increasing distance, the energy of the electric field diverges linearly in the one-dimensional case and logarithmically in the two-dimensional case.
Thus, in (1+1) and (2+1)-dimensional pure Maxwell gauge models, solitons should be electrically neutral. The neutrality, however, does not mean the absence of electric field. In Refs. [19,20], one and two-dimensional soliton systems composed of topological and nontopological components were described. The components interact with each other through an Abelian gauge field and possess opposite electric charges, so the soliton systems are neutral as a whole. Despite electrical neutrality, these soliton systems possess a nonzero electric field that tends to zero exponentially at spatial infinity, resulting in finite electrostatic energy.
The characteristic feature of nontopological solitons is the presence of radially and azimuthally excited states [21][22][23][24][25]. The Q-ball components of compound soliton systems may also be radially or azimuthally excited. The corresponding excited compound soliton systems will have some new features compared to unexcited systems. In the present paper, we study radially and azimuthally excited states of the vortex-Q-ball soliton system described in [19]. We also study the unexcited vortex-Q-ball system using different values of gauge coupling constants.
The paper is structured as follows. In Sect. 2, we describe the Lagrangian, the symmetries, the field equations, and the energy-momentum tensor of the gauge model under consideration. In Sect. 3, we list some properties of the vortex-Q-ball system; among them, the basic differential relation, the asymptotic behaviour of fields at small and large distances, some properties of the gauge potential, the virial relation, and the Laue condition for the vortex-Q-ball system. In Sect. 4, we study properties of the vortex-Q-ball system at extreme values of parameters. In Sect. 5, we present and discuss the numerical results for the unexcited vortex-Q-ball system at different values of gauge coupling constants, the radially excited vortex-Q-ball system, and the azimuthally excited vortex-Q-ball system. In all three cases, we present dependences of the vortex-Q-ball energy on the phase frequency and on the Noether charge along with radial dependences of the vortex-Q-ball ansatz functions.
Throughout the paper, we use the natural units = c = 1.

The gauge model
The Lagrangian density of the (2 + 1)-dimensional gauge model under consideration has the form The model describes the two complex scalar fields φ and χ that minimally interact with the Abelian gauge field A µ through the covariant derivatives The scalar fields φ and χ are self-interacting ones. The self-interaction of φ and χ is described by the fourthand sixth-order potentials, respectively where |φ| 2 = φ * φ and |χ| 2 = χ * χ. In Eqs. (3a) and (3b), λ, g, and h are positive self-interaction constants, m is the mass of the scalar χ-particle, and v is the vacuum average of the amplitude of the complex scalar field φ. From Eq. (3a) it follows that the potential V (|φ|) possesses the continuous family of minima lying on the circle |φ| = v. At the same time, we suppose that the potential U (|χ|) has a global isolated minimum at χ = 0.
For this to hold, the parameters of U (|χ|) must satisfy the inequality 3g 2 < 16m 2 h. The invariance of the Lagrangian density (1) under local gauge transformations (4) and the electrical neutrality of the Abelian gauge field A µ lead to the invariance of the model under the two independent global gauge transformations The invariance of the Lagrangian density under global transformations (5) results in the existence of the two conserved Noether currents The field equations for the model have the form where the electromagnetic current j ν is In Sect. 3, we shall need the form of the symmetric energy-momentum tensor T µν of the model In particular, we shall need the expression for the energy density E = T 00 of a field configuration of the model where E i = F 0i and B = F 21 are the electric field strength and the magnetic field strength, respectively.
3 The vortex-Q-ball soliton system and some of its properties When the gauge coupling constant q is equal to zero, model (1) possesses both the ANO vortex solution [26,27] formed of the complex scalar field φ and the gauge field A µ and the two-dimensional nongauged Q-ball solution formed of the complex scalar field χ. The ANO vortex and the two-dimensional Q-ball are electrically neutral, thus they do not interact with each other. The situation changes drastically if the gauge coupling constant q is different from zero. It was shown in Ref. [19] that in this case, the soliton system consists of interacting vortex and Q-ball components. This two-dimensional vortex-Q-ball system is electrically neutral as a whole, because its vortex and Q-ball components have opposite electrical charges. Nevertheless, the vortex-Q-ball system possesses a nonzero radial electric field in its interior.
Using the Hamilton formalism and Lagrange's method of multipliers [28,29], it was shown in Ref. [19] that there exists a gauge in which only the complex scalar field χ has nontrivial time dependence ∝ exp [−iωt], while the complex scalar field φ and the Abelian gauge field A µ do not depend on time. It was also shown that the vortex-Q-ball system satisfies the important differential relation where E = Ed 2 x and Q χ = j 0 χ d 2 x are the energy and the Noether charge of the vortex-Q-ball system, respectively, and ω is the phase frequency of the complex scalar field χ. Note that in Eq. (13), the phase frequency ω is treated as a function of the Noether charge Q χ . Eq. (13) is a consequence of the fact that the vortex-Q-ball system is an extremum of the energy functional E = Ed 2 x at a fixed value of the Noether charge Q χ = j 0 χ d 2 x. To describe the vortex-Q-ball system, we shall use the following ansatz where K and N are integers, ij are the components of the two-dimensional antisymmetric tensor ( 12 = 1) and n j are those of the two-dimensional radial unit vector n = (cos(θ), sin(θ)). The ansatz functions a 0 (r), a (r), F (r), and σ (r) satisfy the system of nonlinear differential equations The energy density of the vortex-Q-ball system can also be expressed in terms of the ansatz functions The regularity of the vortex-Q-ball solution at r = 0 and the finiteness of the energy E = 2π ∞ 0 E (r) rdr lead to the boundary conditions for the ansatz functions where δ K0 is the Kronecker symbol. Substituting the power expansions for the ansatz functions in Eqs. (15)- (18) and taking into account boundary conditions (20), we obtain the power expansion at the origin of the vortex-Q-ball solution with N = 1. The power expansion of the ansatz function a 0 (r) has the form where and The power expansions of the ansatz functions a (r) and F (r) are written as where the next-to-leading order coefficients are Finally, the power expansion of the ansatz function σ (r) is where the next-to-leading order coefficient is and the term τ K is To obtain the asymptotic form of the vortex-Q-ball solution as r → ∞, we linearize Eqs. (15)- (18) and use boundary conditions (20). As a result, we obtain the where m A = √ 2ev and m φ = √ 2λv are the masses of the gauge boson and the scalar φ-particle, respectively, and ∆ ω = √ m 2 − ω 2 is the mass parameter that defines the asymptotic behaviour of the scalar field χ. Note that the asymptotic forms (31)- (34) are valid only if the mass of any of the three particles (gauge boson, φ-particle, and χ-particle) does not exceed the sum of the mass of the two remaining particles, with the mass parameter ∆ ω playing the role of the χ-particle's mass. Only if this condition is met, can Eqs. (15)-(18) be linearised.
From Eqs. (21)- (30), it follows that the behaviour of the vortex-Q-ball solution in the neighbourhood of r = 0 is determined by four independent parameters: a 1 , b 2 , c 1 , and d K . At the same time, Eqs. (31)-(34) tell us that the behaviour of the vortex-Q-ball solution at spatial infinity is also determined by four independent parameters: a ∞ , b ∞ , c ∞ , and d ∞ . The coincidence of the numbers of parameters that define the behaviour of the solution at the origin and at infinity makes the existence of a solution of the boundary value problem in Eqs. (15)- (18) and (20) possible.
Having boundary conditions (20), we can obtain the constraint on the Noether charges Q φ and Q χ of the vortex-Q-ball system. To do this, we rewrite Eq. (15) (Gauss's law) in compact form where j 0 (r) is the electric charge density expressed in terms of the ansatz functions Then, we integrate both sides of Eq. (35) with respect to r from zero to infinity. Using boundary conditions (20) and asymptotic expressions (21)- (30) and (31)- (34) of the vortex-Q-ball solution, it is easily shown that the integral of the left-hand side of Eq. (35) vanishes. At the same time, the integral of the right-hand side of Eq. (35) is proportional to the electric charge of the vortex-Q-ball system. It follows that the electric charge of the vortex-Q-ball system vanishes. Combining this fact and Eq. (10), we obtain the constraint on the Noether charges Q φ and Q χ of the vortex-Q-ball system where the Noether charges are expressed in terms of the ansatz functions Gauss's law (15) allows us to ascertain some global properties of the time component of the gauge potential A 0 (r) = a 0 (r)/(er). To do this, we rewrite Eq. (15) in the form (rΩ (r)) = 2r q 2 σ (r) where the function According to Eq. (20a), the function Ω (r) → ω as r → ∞. Let the phase frequency ω be positive, then from Eq. (39) it follows that 0 < Ω (r) < ω. Indeed, if Ω (r) < 0 (Ω (r) > ω) at some r =r, then Ω (r) will be negative (positive) at r >r, and thus the boundary condition Ω (r) → r→∞ ω cannot be satisfied. If the phase frequency ω is negative, then ω < Ω (r) < 0. Thus, we have the global conditions on Ω (r) which can be rewritten in terms of the gauge potential A 0 (r) = a 0 (r)/(er) as Eq. (20b) tells us that the boundary conditions for the ansatz function a (r) are the same as those for the ANO vortex. It follows that the magnetic flux of the vortex-Q-ball system is quantized as it is for the ANO vortex where B(r) = −a (r)/(er) is the magnetic field strength.
In particular, the magnetic flux of the vortex-Q-ball system does not depend on the gauge coupling constant q that determines the strength of interaction between the gauge field A µ and the complex scalar field χ.
Having the symmetric energy-momentum tensor (11), we can form the angular momentum tensor Use of Eqs. (11), (14), and (44) results in the angular momentum density J = 1 2 ij J 0ij = J 012 expressed in terms of the ansatz functions where E r (r) = − (a 0 (r) / (er)) is the radial electric field strength. Next, we integrate the term −rBE r = −e −2 a (a 0 /r) by parts, taking into account boundary conditions (20) and using Gauss's law (15) to eliminate a 0 . As a result, the expression for the angular momentum J = 2π ∞ 0 J (r) rdr of the vortex-Q-ball system takes the form where the last line in Eq. (46) follows from Eqs. (38a) and (38b). Using Eq. (37), we can rewrite Eq. (46) in two equivalent forms Let suppose that the gauge coupling constants e and q are multiples of some minimal gauge coupling constant, then the ratio q/e is a rational number. Eq. (47) tells us that in this case, the angular momentum of the vortex-Q-ball system vanishes when the ratio K/N is equal to q/e. Conversely, the ratio q/e is a rational number if there exists a vortex-Q-ball system possessing zero angular momentum. Note that if K/N = q/e then field configuration (14) is invariant under an axial rotation modulo the corresponding gauge transformation, as it should be for a gauged soliton system possessing zero angular momentum.
Using Eqs. (15)- (19), (38), and (46), one can easily ascertain properties of the vortex-Q-ball system under the change of sign of the phase frequency and under the change of signs of the winding numbers N and K From Eq. (19) it follows that the energy of the vortex-Q-ball system can be presented as the sum of five terms where is the energy of the electric field, is the energy of the magnetic field, is the potential part of the energy, is the gradient part of the energy, and is the kinetic part of the energy. In Ref. [19], it was shown that the parts of the energy of the vortex-Q-ball system satisfy the virial relation The energy of the vortex-Q-ball system can be written in several equivalent forms. For this, we integrate the energy density of the electric field (A 0 /(er)) 2 by parts using Gauss's law (15) and taking into account the boundary conditions (20a). As a result, we obtain the following expression for the energy of the vortex-Qball system Combining Eqs. (50) and (57), we obtain the expression for the Noether charge Q χ in terms of E (E) and E (T ) Next, Eqs. (56) and (58) result in the linear relation between the parts of the energy of the vortex-Q-ball system Eqs. (56) and (58) form a system of two linear equations. Using this system, we can express the pairs of variables (E (E) , E (B) ), (E (E) , E (P ) ), (E (T ) , E (B) ), and (E (T ) , E (P ) ) in terms of corresponding remaining variables. Substituting these expressions in Eq. (50), we obtain four more expressions for the energy of the vortex-Q-ball system The spatial components of the energy-momentum tensor are expressed in terms of the ansatz functions as follows: where the radial functions and are the distribution of shear force and pressure, respectively. The conservation of the energy-momentum tensor ∂ i T ik = 0 results in the differential relation between the shear force s (r) and the pressure p (r) To obtain the Laue condition [30,31] for the pressure distribution we multiply Eq. (64) by r 2 and integrate by parts over r from zero to infinity. It can be shown that the Laue condition (65) is equivalent to virial relation (56).

Extreme regimes of the vortex-Q-ball system
In this section, we ascertain the properties of the vortex-Q-ball system in several extreme regimes. First, we consider the behaviour of the vortex-Q-ball system at extreme values of gauge coupling constants, then we discuss the vortex-Q-ball system at extreme values of phase frequency ω (thick-wall and thin-wall regimes). The two scalar fields φ and χ of the vortex-Q-ball system interact with the Abelian gauge field A µ . The intensity of this interaction is determined by the two gauge coupling constants e and q. Note that in the natural units = c = 1, the electric charges of the scalar φ and χ-particles also equal e and q, respectively. We suppose that the group associated with the gauge symmetry of model (1) is the compact Abelian group U (1). In this case, the electric charges (gauge coupling constants) e and q are commensurable and thus the ratio τ = q/e is a rational number. This means that some minimal elementary electric charge exists and that all others charges are multiples of this elementary charge.
Let us consider the extreme regime in which both e and q tend to zero, while the ratio τ = q/e remains finite. It can be ascertained both analytically and numerically that as q = eτ → 0, the time component of gauge potential where α 0 (r) is some function of r that remains finite as q = eτ → 0. It follows from Eq. (66) that a 0 (r)/(er) uniformly tends to zero in the limit q = eτ → 0. Hence, the energy of the electric field (51) also tends to zero (E (E) ∝ e 2 ). Next, after the change of radial variable r = ρ/m A = ρ/( √ 2ev), Eq. (16) will not explicitly depend on the infinitesimal parameters e and q. Instead, it will depend only on the ratio τ = q/e, which is finite. On the other hand, after the change of radial variable ) will concentrate in a small neighbourhood of ρ = 0. Hence, in the limit q = τ e → 0, the ansatz functions F (ρ) andσ(ρ) in Eq. (16) can be replaced by their limiting values of 1 and 0, respectively. After that, Eq. (16) takes the simple form The solution of Eq. (67) satisfying boundary condition (20b) can be written as where K 1 (ρ) is the modified Bessel function of the second kind. From Eq. (68) it follows that in terms of the initial variable r, the ansatz function a(r) spreads out over the interval r ∈ [0, ∞) as q = τ e → 0. Specifically, for any finite r =r, a (r) → 0 in this limit. At the same time, the ansatz functions F (r) and σ(r) remain concentrated in a finite neighbourhood of the origin as q = τ e → 0 because Eq. (17) does not explicitly depend on the gauge coupling constants, and Eq. (18) depends on them only through the combination τ = q/e, which is finite. These facts and Eq. (66) lead us to conclude that the gauge field A µ decouples from the Q-ball component of the vortex-Q-ball system, and thus the Q-ball component tends to the nongauged Q-ball solution as q = τ e → 0. Hence, the energy of the Q-ball component tends to the finite energy of the two-dimensional nongauged Q-ball in the limit q = τ e → 0. The situation is different for the vortex component, which is topologically nontrivial. The topological nontriviality prevents decoupling of the gauge field from the complex scalar field φ as q = τ e → 0. Indeed, it follows from Eqs. (52) and (68) that the energy of the magnetic field remains finite and tends to πv 2 N 2 as q = τ e → 0. Moreover, it follows from Eq. (43) that the magnetic flux diverges as e −1 in the limit q = τ e → 0. Nevertheless, the scalar field φ of the vortex component tends to the nongauged global vortex solution [32] because in the limit q = τ e → 0, the ansatz function a(r) uniformly tends to zero in the region where the ansatz function F (r) differs appreciably from the limiting value of 1.
According to Derrick's theorem [33], the energy of the global vortex is infinite, hence the energy of the vortex component of the vortex-Q-ball system tends to infinity in the limit q = τ e → 0. Indeed, let us integrate the term (N + a) 2 v 2 F 2 r −2 in gradient energy (54) over the region r m φ , where the ansatz function F (r) is close to 1. Using Eq. (68) for the ansatz function a(r), we obtain the expression 2π 2 N 2 v 2 ln (m φ /m A ) ∝ ln e −1 , which diverges logarithmically as e → 0. Now we consider another extreme regime in which both gauge coupling constants e and q tend to infinity: q = τ e → ∞. In this regime, the gauge field A µ ceases to be dynamic. Instead, it is expressed in terms of the complex scalar fields φ and χ in the whole space, except for the infinitesimal neighbourhood of the origin. Indeed, the Lagrangian density (1) is written in terms of ansatz functions (14) as We see that as q = τ e → ∞, the first two terms in Eq. (69) tend to zero (provided that the corresponding derivatives in Eq. (69) are finite) and thus can be neglected. In this case, the field equations for the ansatz functions a 0 (r) and a(r) become purely algebraic (terms with derivatives can be neglected in Eqs. (15) and (16)), so a 0 (r) and a(r) are expressed in terms of F (r), σ(r), and the model's parameters For K = 0, the ansatz functions a 0 (r) and a(r) uniformly tend to their limiting values (70) on the interval r ∈ [0, ∞). At the same time, Eq. (70) cannot be used in the infinitesimal neighbourhood of the origin (where r m −1 A ∝ e −1 ) for nonzero K. Indeed, it follows from Eq. (70b) that a(0) = −N if K = 0 (σ(0) = 0 in this case), but this contradicts boundary condition (20b).
It was found that the behaviour of the vortex-Qball system differs substantially for zero and nonzero K. From Eq. (70b) it follows that for K = 0, the ansatz function a(r) varies from 0 to the vicinity of −N on the interval ∆r A ∼ e 0 , which does not depend on e. Hence, we obtain the sequence of relations: where Eq. (24) is used. It follows that for zero K, the energy of the magnetic field tends to zero as q = τ e → ∞.
On the other hand, for nonzero K, the ansatz function a(r) varies from 0 to the vicinity of −N on the interval ∆r A ∼ m −1 A ∝ e −1 , which shrinks to the origin as e → ∞. Thus we have the chain of relations: rdr ∝ e 0 . We see that unlike the previous case, the energy of the magnetic field tends to a finite value as q = τ e → ∞. This is because the magnetic field strength B increases indefinitely (∝ e) in the infinitesimal neighbourhood (∆r A ∝ e −1 ) of the origin. At the same time, magnetic flux (43) of the vortex-Q-ball system tends to zero as e −1 for both K = 0 and K = 0.
From Eq. (70a) it follows that for all values of K, the time component a 0 /(er) of the gauge potential A µ uniformly tends to zero (∝ e −1 ) in the limit of large e. Furthermore, the major part of variation of the ansatz function a 0 (r) happens on the interval that does not depend on e. Thus it follows that the energy of the electric field (51) also tends to zero (E (E) ∝ e −2 ) in the limit q = τ e → ∞.
Substituting Eqs. (70a) and (70b) in Eq. (36), we find that the electromagnetic current vanishes as q = τ e → ∞ For K = 0, Eq. (71) is valid for all r ∈ [0, ∞), while for nonzero K, it is valid only on the interval (∆r A , ∞), where ∆r A ∼ m −1 A ∝ e −1 . Next, let us consider the thick-wall regime of the vortex-Q-ball system. In this extreme regime, the modulus of the phase frequency tends to the mass of the scalar χ-particle: |ω| → m. In this case, the Q-ball component of the vortex-Q-ball system spreads over the two-dimensional space, so the amplitude σ of the complex scalar field χ uniformly tends to zero. The time component A 0 = a 0 /(er) of the gauge potential also uniformly tends to zero in the thick-wall regime. From Eqs. (16) and (17) it follows that the influence of the Qball component on the vortex component of the vortex-Q-ball system can be neglected in the thick-wall regime. Hence, the vortex component of the vortex-Q-ball system tends to the ANO vortex solution as |ω| → m. At the same time, the energy and the Noether charge of the Q-ball component tend to finite values in the thickwall regime as they do for the corresponding nongauged two-dimensional Q-ball [34,35].
To ascertain the behaviour of the vortex-Q-ball system in the thick-wall regime, we rescale the ansatz functions A 0 (r) = a 0 (r)/(er), σ(r), and the radial variable r as follows: where r = ∆ −1 ωr and ∆ ω = m 2 − ω 2 1/2 . Next we introduce a new functional F that is related to the energy functional E = Ed 2 x by the Legendre transformation It can be shown that −F (ω) is equal to the Lagrangian L = 2π ∞ 0 Lrdr on field configurations that satisfy Gauss's law (15). Next, using Eq. (72), we write the functional F (ω) as the sum of the three terms In Eq. (74), F 0 , F 2 , and F 4 are written as the integrals of the corresponding densities: and Note that in Eq. (75), the prime means the differentiation with respect to the radial variable r, while in Eqs. (76) and (77), it means the differentiation with respect to the rescaled radial variabler = ∆ ω r. Note also that as ∆ ω → 0, we have replaced the ansatz functions a (r) = a ∆ −1 ωr ≡ā (r) and F (r) = F ∆ −1 ωr ≡ F (r) by their limiting values −N and 1, respectively.
Eq. (76) does not contain the derivative ofĀ 0 , and thusĀ 0 can be expressed in terms of the ansatz function σ and the model's parameters We see that F 2 is explicitly dependent on the small parameter ∆ ω . Hence, Eq. (74) can be rewritten as where F 0 is the same as in Eq. (74),F 2 = 2π F 2r dr, F 4 = 2π F 4r dr, and the corresponding densities arē We see that none of F 0 ,F 2 , andF 4 depend on the phase frequency ω.
Using the known properties of the Legendre transformation, we obtain the Noether charge Q χ and the energy E of the vortex-Q-ball system as functions of the phase frequency ω and (83) and (84) it follows that the energy and the Noether charge of the vortex-Q-ball system tend to the finite values E tk = F 0 + 2ω 2 tkF 2 and Q χtk = 2ω tkF2 , respectively, as ω → ω tk = ±m. Moreover, Eqs. (83) and (84) make it possible to obtain the dependence of the energy on the Noether charge, which is valid for both signs of Q χ in the thick-wall regime where ∆Q χ = |Q χtk | − |Q χ |.
Although the influence of the Q-ball component on the vortex component is negligible in the thick-wall regime, but the reverse is not true. Indeed, by varying the functionalF 2 = 2π F 2r dr inσ, we obtain the differential equation forσ, which is valid in the thick-wall regimē We see that Eq. (86) depends on the parameters N and v of the vortex component, hence, the influence of the vortex component on the Q-ball component can not be neglected in the thick-wall regime. Eq. (81) tells us that the functionalF 2 = 2π F 2r dr depends on the gauge coupling constants e and q only through the combination τ = q/e. Hence, the limiting value Q χtk = 2ω tkF2 depends on e and q only through the combination τ = q/e. At the same time, we see from Eq. (75) that the functional F 0 = 2π F 0 rdr (the energy of the ANO vortex for given e, λ, v, and N ) depends solely on the gauge coupling constant e. Hence, the limiting energy E tk = F 0 + 2m 2F 2 of the vortex-Qball system depends on e and q separately. In particular, it diverges logarithmically (∝ ln(e −1 )) when e tends to zero and τ remains fixed.
Finally, we consider the thin-wall regime of the vortex-Q-ball system. This extreme regime occurs when the absolute value |ω| of the phase frequency tends to the minimum possible value ω tn . As a result, the energy E = 2π ∞ 0 Erdr and the Noether charge Q χ = 2π ∞ 0 j 0 χ rdr increase indefinitely as ω → ω tn . In the thin-wall regime, the field configuration of the vortex-Q-ball system can be divided into three regions: the central transitional region, the basic interior region, and the exterior transitional region. The spatial size of the internal region increases indefinitely as the vortex-Q-ball system approaches the thin-wall limit. The characteristic feature of the thin-wall regime is that the ansatz functions σ(r), Ω(r) = ω −τ a 0 (r)/r, and F (r) tends to constant values in the internal region as ω → ω tn . Equating the derivatives in Eqs. (15), (17), and (18) to zero (Eq. (15) should be rewritten in terms of Ω(r) as it is in Eq. (39)), we obtain a system of three algebraic equations (ω − Ω) 2 + τ 2 λv 2 1 − F 2 = 0.
One more equation can be obtained from Eq. (59). Indeed, in the thin-wall regime, the energy of the magnetic field tends to zero, while the potential and kinetic parts of energy increase indefinitely. Hence, the term E (B) can be neglected in the thin-wall regime and Eq. (59) takes the form E (T ) − E (P ) = 0. This fact and Eqs. (53) and (55) result in the fourth algebraic equation Before we go any further, let us ascertain the behaviour of the ansatz function a(r) in the internal region of the vortex-Q-ball system in the thin-wall regime. From Eq. (16) it follows that in this case, the behaviour of a(r) is described by the differential equation where and F tn and σ tn are the limiting values of the corresponding ansatz functions in the thin-wall regime. The appropriate solution of Eq. (91) is where K 1 √ βr is the modified Bessel function of the second kind and C is a positive constant. We see that in the internal region of the vortex-Q-ball system, a(r) tends to the constant value Note that the limiting value (94) does not equal the boundary value −N of a(r) at spatial infinity. It follows from Eq. (94) that the behaviour of a(r) is determined by the thin-wall background values F tn and σ tn of F (r) and σ(r), respectively. At the same time, Eqs. (17) and (18) tell us that in the thin-wall regime, the backward influence of a(r) on F (r) and σ(r) can be neglected in the internal region of the vortex-Q-ball system. Indeed, the ansatz functions a 0 (r) and Ω(r) are connected by the relation a 0 = τ −1 r (ω − Ω), and the constancy of Ω in the internal region of the vortex-Q-ball system results in the linear growth of a 0 in r there. Conversely, the ansatz function a(r) is bounded in the internal region of the vortex-Q-ball system and the ratio a(r)/a 0 (r) tends to zero there, because the size of the internal region increases indefinitely in the thin-wall regime.
Thus, we have the system of four algebraic equations (87)-(90), which are valid in the thin-wall regime. These equations allow us to determine the limiting thin-wall values Ω tn , σ tn , and F tn of the corresponding ansatz functions and the limiting thin-wall value ω tn of the phase frequency. The solution of system (87)-(90) cannot be obtained analytically, in general. However, this solution can be obtained as series in the parameter τ = q/e, provided it is small enough F tn = 1 + 9 32 where σ tn0 = √ 3 2 g h and ω tn0 = m 1 − 3 16 are the limiting thin-wall values of σ(r) and ω for the nongauged Q-ball. Note that the system (87)-(90) does not depend on the integers K and N , hence the thinwall values Ω tn , σ tn , F tn , and ω tn also does not depend on K and N . Furthermore, the system (87)-(90), and consequently, Ω tn , σ tn , F tn , and ω tn depends on the gauge coupling constants e and q only through the combination τ = q/e. Using Eqs. (95a)-(95d), we can obtain the thin-wall values of the energy, Noether charge, and angular momentum density where are the corresponding densities for the nongauged Qball. It follows from Eqs. (95a), (97a), and (97b) that the ratio E tn j 0 χtn = ω tn0 1 + 3 8 as it should be in the thin-wall regime.

Numerical results
The system of differential equations (15)-(18) with boundary conditions (20) is a mixed boundary value problem on the semi-infinite interval r ∈ [0, ∞). It is obviously that this problem can be solved only by numerical methods. To solve this problem, we use the boundary value problem solver provided in the Maple package [36]. The correctness of numerical results is controlled with the help of Eq. (13) and the Laue condition (65). The mixed boundary value problem (15)- (18), and (20) depends on the eight parameters: m, g, h, λ, v, e, q, and ω. Without loss of generality, the number of the parameters can be reduced, if we rescale the radial variable r and the ansatz function σ as follows: After rescaling, the vortex-Q-ball system will be described solely by the six dimensionless parameters:h = m 2 g −2 h,λ = gv 4 m −4 λ,ṽ = g 1/2 m −1 v,ẽ = g −1/2 e, q = g −1/2 q, andω = m −1 ω, whereas the two dimensionless parametersm andg will be equal to 1. In most numerical calculations, we use the following values for the nongauged dimensionless parameters:h = 0.2, λ = 0.125,ṽ = 2. The dimensionless gauge coupling constantsẽ andq are taken to be equal (τ =q/ẽ = 1), and may vary in some interval. In addition to these parameters, the mixed boundary value problem (15)- (18), and (20) also depends on the two integers: N (topological winding number of the vortex component) and K (nontopological winding number of the Q-ball component). We shall consider the vortex-Q-ball systems with the vortex winding number N = 1, while the Q-ball winding number K may take a range of values. Figure 1 shows the dependence of the dimensionless energyẼ = gm −2 E of the vortex-Q-ball system on the dimensionless phase frequencyω = m −1 ω. The presented curves correspond to the vortex-Q-ball system at six different values of the gauge coupling constants and to the nongauged two-dimensional Q-ball. We see that with decreasingω, the system passes into the thin-wall regime in which the energy and the Noether charge increase indefinitely. Whenω → 1, the vortex-Q-ball system passes into the thick-wall regime. In this regime, the behaviour of the vortex-Q-ball system differs con-  Fig. 2 Dependence of the Noether chargeQ χ of the vortex-Q-ball system on the dimensionless phase frequencyω in the vicinity ofω = 1. The parameters of the vortex-Q-ball system are the same as in Fig. 1 siderably from that of the nongauged two-dimensional Q-ball. In particular, the energy of the Q-ball decreases monotonically asω → 1, whereas the behaviour of the energy of the vortex-Q-ball system is more complicated, as follows from the subplot in Fig. 1. We see that for the vortex-Q-ball system, the curvesẼ(ω) are s-shaped in the vicinity ofω = 1. This fact is due to the nontrivial interaction between the vortex and Q-ball components of the soliton system. Note that in Ref. [19], we have not managed to obtain, by numerical methods, the sshaped parts of theẼ(ω) andQ χ (ω) curves.
It follows from the subplot in Fig. 1 that the limit valueẼ(1) is not a constant for the vortex-Q-ball system, but increases with a decrease ofẽ. It was found numerically that for small values ofẽ, the valueẼ(1) ≈ a + b ln(ẽ −1 ), where a and b depend on the gauge coupling constants through the ratio τ =q/ẽ in accordance with the conclusion of Sect. 4.
Beyond the neighbourhood ofω = 1, the rescaled Noether chargeQ χ = gm −1 Q χ depends onω similar to the dimensionless energyẼ = gm −2 E in Fig. 1. However, the behaviour of the curvesQ χ (ω) andẼ(ω) is different in the neighbourhood ofω = 1. Indeed, Figs. 2 and 3 show the curvesQ χ (ω) in neighbourhoods of ω = 1 for different values of the gauge coupling constants. The curves in Fig. 2 correspond to the same values of the gauge coupling constant as in Fig. 1, while those in Fig. 3 correspond to larger values of the gauge coupling constants. We see that unlike the curvesẼ(ω)  Figures 2 and 3 show that the form of the curves Q χ (ω) changed with the increase in the gauge coupling constants. In Fig. 2, the curvesQ χ (ω) are s-shaped similar to the curvesẼ(ω) in the subplot in Fig. 1. At the same time, it follows from Fig. 3 that by increasing the gauge coupling constants, the curvesQ χ (ω) cease to be s-shaped, and the turning points of s-shapedQ χ (ω) curves turn into a single inflection point of monotonically increasingQ χ (ω) curves. Figure 4 presents the energyẼ of the vortex-Q-ball system as a function of its Noether chargeQ χ for the same values of the gauge coupling constants as in Fig. 1. We see that all the curvesẼ(Q χ ) that correspond to the vortex-Q-ball system have one cuspidal point, whereas there is no cuspidal point on the curveẼ(Q χ ) for the two-dimensional nongauged Q-ball. The next characteristic feature is that theQ χ -coordinates for the rightmost points of the upper branches of the curvesẼ(Q χ ) coincide. Of course, this is a consequence of the fact that the valueQ χ (1) depends onẽ andq only through the ratio τ =q/ẽ, which is the same for all vortex-Qball curves in Fig. 4.
In Fig. 4, the dashed straight line corresponds to the plane-wave field configuration of the complex scalar field χ. We see that the energy of the two-dimensional nongauged Q-ball is less than that of the plane-wave field configuration, except for the point of contact, at which they are equal. Hence, the two-dimensional nongauged Q-ball is stable against decay into massive scalar χ-bosons. It follows from Fig. 4 that except for the neighbourhood of cuspidal points, the energy of the vortex-Q-ball system is less than that of the plane-wave field configuration. Hence, the Q-ball component of the vortex-Q-ball system is stable against decay into massive scalar χ-bosons provided thatẼ <Q χ .
To ascertain the possibility of decay of the Q-ball component of the vortex-Q-ball system, we introduce the value ∆Ẽ =Ẽ −Q χ −Ẽ v , whereẼ v is the energy of the vortex component of the system. We defineẼ v as the energy of the ANO vortex at a given value of the gauge coupling constantẽ. The curvesẼ v (Q χ ) are presented in the subplot in Fig. 4. It is obvious that decay of the Q-ball component into scalar χ-bosons is possible only if ∆Ẽ is positive. It follows from the subplot in Fig. 4 that for all gauge coupling constants, ∆Ẽ is positive on the upper branches of curvesẼ(Q χ ). Furthermore, ∆Ẽ is also positive on the lower branches ofẼ(Q χ ) curves in neighbourhoods of their cuspidal points. Hence, the vortex-Q-ball system is unstable in the area presented in the subplot in Fig. 4. The instability, however, can be either classical (the presence of one or more unstable modes in the functional neighbourhood of the soliton system) or quantum-mechanical (the possibility of quantum tunneling of the soliton system to another state). It was shown in Refs. [34,37] that the appearance of a cusp on the energy-Noether charge curve indicates the onset of a mode of instability. Hence, the Q-ball components of the soliton systems lying on Next we present the ansatz functions of the vortex-Q-ball system for different values of gauge coupling constants. Figure 5 presents the ansatz functions a 0 (r)/r andσ(r), and Figure 6 presents the ansatz functions a(r) and F (r). It follows from Fig. 5 that a 0 (r) /r|r =0 increases monotonically and a 0 (r) /r tends to limiting form (70a) asẽ andq increase. In particular, a 0 (r) /r|r =0 reaches the maximum limit valueω/τ = 0.7 asẽ =q → ∞, which is consistent with Eqs. (42) and (70a). The ansatz functionσ(r) also tends to a limiting form as the gauge coupling constants increase indefinitely. Figure 6 shows that similar to a 0 (r)/r, the ansatz function a(r) is sensitive to the magnitude ofẽ andq. In particular, it tends to zero at any finiter asẽ =q → 0 and tends to limiting form (70b) asẽ =q → ∞. The dependence of F (r) on the the gauge coupling constants is not as strong as that of a(r). It follows from Fig. 6 that F (r) tends to a limiting form asẽ =q → ∞ and to the global vortex solution asẽ =q → 0.
Like two and three-dimensional nongauged Q-balls, radially excited states in the vortex-Q-ball system exist. In Fig. 7, we can see the curvesẼ(ω) for the unexcited (n = 0) and the first six (n = 1, . . . , 6) radially excited states of the vortex-Q-ball system. As in the previous (unexcited) case, the radially excited system passes into the thin-wall regime asω →ω tn and into the thick-wall regime asω → 1. Like curves in Fig. 1, the curvesẼ(ω) of radially excited vortex-Q-ball systems are s-shaped in the vicinity ofω = 1. Next, we see that at fixedω, the energy of the the vortex-Q-ball system increases with an increase in n. It was found numerically that at fixed ω, the energy and Noether charge increase quadratically in n starting with n = 2: where a and b are positive constants, whereas the signs of constants c and d coincide with that ofω. The behaviour ofQ χ (ω) curves is similar to that ofẼ(ω) curves in Fig. 7. In particular, the valuesQ χ (1) are different for different n. Figure 8 presents the curvesẼ(Q χ ) for the unexcited and the first six radially excited states of the vortex-Q-ball system. All curvesẼ(Q χ ) have cuspidal points in whichẼ andQ χ reach minimum values. Numerically, we found that starting with n = 2, the minimum values ofQ χ in the cusp points are well described by a quadratic dependence: Q χc = a c + b c n 2 . It follows that at a givenQ χ , the number of radially excited states of the vortex-Q-ball system does not exceed where x denotes the floor of x (the greatest integer less than or equal to x). We see that at largeQ χ , the number of radially excited states rises ∝Q It follows from Fig. 8 that at a givenQ χ , the lowerbranch energyẼ increases with n. Hence, radially excited vortex-Q-ball states lying on lower branches of E(Q χ ) curves are unstable with respect to the transition into less excited states. Furthermore, we can see from the subplot in Fig. 8 that the Q-ball components of radially excited vortex-Q-ball states are unstable with respect to decay into scalar χ-bosons in neighbourhoods of cuspidal points (the upper branches and small parts of the lower branches ofẼ(Q χ ) curves).
Figures 9 and 10 present ansatz functions for the first 10 radially excited states of the vortex-Q-ball system. We see from Fig. 9 that oscillations of theσ(r) curves are inharmonious. In particular, amplitudes of peaks and valleys ofσ(r) curves are decreased with the increase in the number of oscillations, whereas the distances between adjacent peaks and valleys ofσ(r) curves are increased. As seen in Fig. 10, the ansatz functionsã 0 (r)/r and a(r) also oscillate. In particular, the positions of local maxima ofã 0 (r)/r and a(r) approximately coincide with those of peaks and valleys of σ(r), whereas the positions of local minima ofã 0 (r)/r and a(r) are approximately the same as those of nodes (zeros) ofσ(r).
It was shown in Ref. [25] that the number of radially excited states of the three-dimensional gauged Q-ball is finite. At the same time, we found no indication of the finiteness of the number of radially excited states for the vortex-Q-ball system. Of course, the reason for this difference is that the three-dimensional gauged Q-ball possesses an electrical charge, whereas Ansatz functions a, 3a 0 /r, and F for the first 10 radially excited states of the vortex-Q-ball system. The parameters of the vortex-Q-ball system are the same as in Fig. 9 the two-dimensional vortex-Q-ball system is electrically neutral. As a result, the electric charge density of the three-dimensional gauged Q-ball is always either positive or negative, whereas that for the vortex-Q-ball system is alternating. This difference is apparent in the behaviour of the ansatz function Ω(r) defined in Eq. (40). Indeed, it is shown in Ref. [3] that in the case of the three-dimensional gauged Q-ball, Ω(r) is a bounded (0 < Ω(r) < ω) and monotonically increasing function of r. For the vortex-Q-ball system, the ansatz function Ω(r) is also bounded in the interval (0, ω) as it is seen from Eq. (41). However, Ω(r) need not be monotonic in this case. Indeed, the oscillating behaviour of a 0 (r)/r shown in Fig. 10 results in the oscillating behaviour of Ω(r) = ω−τã 0 (r)/r. The monotonic increase of Ω(r) leads to the restriction of the number of radially excited states for the three-dimensional gauged Q-ball, because in this case, Ω(r) reaches its maximum value ω for a finite number of oscillations of σ(r). In contrast, the nonmonotonic oscillating behaviour of Ω(r) makes it possible for vortex-Q-ball systems with an arbitrarily large number of oscillations (and hence nodes) of σ(r) to exist; thus, there are no restrictions on the number of radially excited states of the vortex-Q-ball system. Now we turn to a description of the vortex-Q-ball's states with nonzero integer-valued parameter K. Figure 11 shows the curvesẼ(ω) for the first few states with nonzero K. It follows from Fig. 11 that any twõ E(ω) curves whose parameters K and K satisfy the condition K + K = 2 tend to the same limit in the thick-wall regime whenω → 1. A similar statement  Fig. 11 Dependence of the energyẼ of the vortex-Q-ball system on the phase frequencyω for the first few azimuthally excited states. The curves correspond to parametersm = 1, g = 1,h = 0.2,λ = 0.125,ṽ = 2,ẽ =q = 0.1, and N = 1 is valid for theQ χ (ω) curves. This facts can be explained as follows. In the thick-wall regime, the limiting energy of the vortex-Q-ball system is written as E tk = F 0 + 2m 2F 2 , where F 0 is the energy of the ANO vortex with a given e, and 2m 2F 2 is the energy of the Q-ball component. It follows from Eq. (81) thatF 2 is expressed solely in terms of the rescaled ansatz functionσ. Next, the rescaled ansatz functionσ satisfies differential equation (86). We see that Eq. (86) depends on the integer-valued parameters K and N only through the combination κ = (K−τ N ) 2 . It follows that the parameters K and K , which satisfy the relation K + K = 2τ N , lead to the same value of κ. Because both K and K are integers, the value of 2τ N is also an integer, thus the parameter τ = e/q is an integer or half-integer. If 2τ N is a noninteger, thenẼ(ω) curves corresponding to different K will not tend to the same limit asω → 1. In our case, the parameters N and τ are equal to one, so 2τ N = 2. Thus, differential equation (86) will have the same form for any two vortex-Q-ball systems whose parameters satisfy the condition K + K = 2.
Eq. (86) is valid for allr except those in the interval 0 ≤r ∆ ω m −1 A,φ whose width tends to zero in the thick-wall regime. This is because the rescaled ansatz functionsā(r) ≡ a(∆ −1 ωr ) andF (r) ≡ F (∆ −1 ωr ) are different from their limiting values in this infinitesimal interval ofr. Next, the ansatz functionsσ K (r) and  Fig. 12 Dependence of the energyẼ of the vortex-Q-ball system on the Noether chargeQ χ for the first few azimuthally excited states. The straight dashed line corresponds to the linear dependenceẼ =Q χ . The subplot presents the dependence of the difference ∆Ẽ =Ẽ −Q χ −Ẽ v onQ χ in neighbourhoods of cuspidal points σ K (r) (the dependence ofσ on K is shown explicitly) satisfy the same boundary conditions:σ(0) = 0 and σ(∞) = 0. The only exception is in the case K = 0 for which the left boundary condition isσ (0) = 0. However, it was found numerically that in this case, σ(0) → 0 as ∆ ω → 0, and the left boundary condition for K = 0 becomes essentially the same as that for the complementary case K = 2. It follows that the ansatz functionsσ K (r) andσ K (r) that satisfy the condition K + K = 2 tend to the same limit in the thick-wall regime and so do the corresponding functionalsF 2K andF 2K . Hence, the energies (Noether charges) of the two vortex-Q-ball systems with K + K = 2 tend to the same value in the thick-wall regime. Note that there is no complementary vortex-Q-ball system for the system with K = 1. Indeed, in this case the parameters K and K are the same: K = 1, K = 1 ⇒ K + K = 2.
It follows from Eq. (47) that the two vortex-Q-ball systems with the same Noether charges and with integer parameters K and K such that K + K = 2τ N possess the opposite angular momenta. Hence, the angular momenta of two vortex-Q-ball systems with K + K = 2 tend to the opposite values in the thickwall regime, whereas the angular momentum of the state with K = 1 vanishes. Figure 12 presents theẼ(Q χ ) curves for the same K as in Fig. 11. We see that just as in Figs. 4 and 8, the curvesẼ(Q χ ) have cuspidal points. The only exception is theẼ(Q χ ) curve for K = 1; in this case, the absence of a cuspidal point follows from the monotonicity of the correspondingẼ(ω) curve in Fig. 11. The conservation of the Noether charge and the angular momentum leads to the conclusion that the parts ofẼ(Q χ ) curves lying below the lineẼ =Ẽ v +Q χ correspond to stable vortex-Q-ball systems. In the subplot in Fig. 12, we can see the parts ofẼ(Q χ ) curves that lie above the lineẼ = E v +Q χ and thus correspond to unstable vortex-Q-ball systems.
It follows from Fig. 12 that the curvesẼ K (Q χ ) and E K (Q χ ) with K + K = 2τ N = 2 approximately coincide starting with (K , K ) = (4, −2). Hence, the energies of the two vortex-Q-ball systems that correspond to coincident curvesẼ K (Q χ ) andẼ K (Q χ ) are approximately the same at given Q χ . Note that the parameter K − τ N has opposite values for the curves E K (Q χ ) andẼ K (Q χ ) with K + K = 2τ N . Hence, the angular momentaJ K (Q χ ) andJ K (Q χ ) are opposite as it follows from Eq. (47). It follows from Eqs.  The reason for this is the presence of the centrifugal term −r −2 (K + τ a (r)) 2 σ (r) in Eq. (18). The contribution of this term is proportional to r −2 (K − τ N ) 2 when a(r) is in the vicinity of the limiting value of −N . The growth of |K − τ N | leads to the growth of the factor κ = (K − τ N ) 2 in the centrifugal term. The growth of κ must be compensated, otherwise, the solution of mixed boundary value problem (15)- (18) and (20) will not exist. Such compensation is achieved through the shift of the ansatz function σ(r) towards higher values of r. The shift results in reducing the factor r −2 in the centrifugal term which compensates for the growth of the factor κ. With increasing r, it is convenient to pass to the new ansatz function ∆a(r) = N + a(r). Then the system (15)-(18) will depend on the parameters K and N only through the combination K − τ N . Figures 13 and 14 show the ansatz functions for the same K as in Fig. 11. We see that for K = 0, the ansatz functionsσ(r) have a slightly asymmetric bell-shaped form. Numerically, we found that the radial positions of maxima of theσ(r) curves increase approximately linearly with |K − τ N |. It follows from Fig. 13 that with an increase in |K − τ N |, the ansatz functions a 0 (r)/(ẽr) also form maxima whose radial positions are approximately the same as those of the correspondingσ(r) curves. Next, we see from Fig. 14 that the forms of a(r) curves strongly depend on K − τ N . In particular, with an increase in |K − τ N |, the minimum of a(r) appears for positive K and the maximum (together with the left adjacent minimum) of a(r) appears for negative K. The radial positions of these extremes of a(r) approximately coincide with those of the maxima of the correspondingσ(r). Such a difference in the behaviour of the a(r) curves with positive and negative K is explained by the fact that the driving term −2eq(K − τ N )σ(r) 2 in Eq. (16) has the opposite sign in these cases. Finally, it follows from Fig. 14 that the forms of the F (r) curves weakly depend on K.

Conclusion
In the present paper, we continued the study of the vortex-Q-ball systems started in Ref. [19]. We investigated properties of the unexcited vortex-Q-ball systems at different values of gauge coupling constants and properties of radially excited vortex-Q-ball systems. We also investigated properties of azimuthally excited vortex-Q-ball systems. The vortex-Q-ball system is composed of a vortex (topological soliton) and a Q-ball (nontopological soliton), thus it combines properties of both topological and nontopological solitons. In particular, the vortex-Q-ball possesses a quantized magnetic flux and satisfies basic relation (13) for nontopological solitons.
The electromagnetic interaction between the vortex and Q-ball components results in significant changes in their properties. Firstly, the vortex and Q-ball components of the vortex-Q-ball system acquire opposite electrical charges, so the system remains electrically neutral as a whole. At the same time, neither ANO vortex nor two-dimensional Q-ball possess electrical charge because in gauge models without the Chern-Simons term, any two-dimensional electrically charged object will have infinite energy. Secondly, the vortex-Q-ball system possesses nonzero angular momentum despite the fact that both the ANO vortex and the two-dimensional axially symmetric Q-ball have zero angular momenta. Thirdly, the interaction between the vortex and Q-ball components leads to substantial change in the E(ω) and E(Q χ ) curves for the vortex-Q-ball system in comparison with that of the two-dimensional nongauged Qball. In particular, in the case of the vortex-Q-ball system, some of E(ω) curves are s-shaped in the vicinity of ω = m and the E(Q χ ) curves have cuspidal points, whereas none of these features hold for the twodimensional nongauged Q-ball.
There exist several extreme regimes for the vortex-Q-ball system. We investigated four of them: the thickwall regime in which the phase frequency ω tends to the maximum value, the thin-wall regime in which ω tends to the minimum value, and regimes of small and large gauge coupling constants. In particular, we found that the limiting thick-wall value of the Noether charge depends on the gauge coupling constants e and q only through the ratio τ = q/e. We also found that the limiting thick-wall energies of two vortex-Q-ball systems whose azimuthal parameters K and K satisfy the condition K + K = 2τ N are the same. As for extreme values of gauge coupling constants, the gauge field A µ is decoupled from the vortex and Q-ball components as q = τ e → 0 and is expressed in terms of the scalar fields φ and χ as q = τ e → ∞. The latter means that the gauge field A µ ceases to be a dynamic object as q = τ e → ∞.
Note that the very possibility of the existence of radially excited vortex-Q-ball states results from the nonlinear character of mixed boundary value problem (15)- (18) and (20). Indeed, a linear homogeneous boundary value problem is the SturmLiouville problem. Solutions with different numbers of nodes correspond to different eigenvalues of the SturmLiouville problem. It follows that these solutions satisfy different differential equations because the eigenvalue is a parameter of differential equation. In contrast, a nonlinear boundary value problem with fixed parameters may have more than one solution, as it is in the case of the vortex-Q-ball system. In this case, the second and subsequent solutions of the mixed nonlinear boundary value problem correspond to the radially excited states of the vortex-Q-ball system.