One-dimensional soliton system of gauged kink and Q-ball

In the present paper, we consider a $\left(1 + 1\right)$-dimensional gauge model consisting of two complex scalar fields interacting with each other through an Abelian gauge field. When the model's gauge coupling constants are set equal to zero, the model possesses non-gauged Q-ball and kink solutions that do not interact with each other. It is shown that at nonzero gauge coupling constants, the model possesses the soliton solution describing the system consisting of interacting Q-ball and kink components. The kink and Q-ball components of the kink-Q-ball system have opposite electric charges, so the total electric charge of the kink-Q-ball system vanishes. Properties of the kink-Q-ball system are researched analytically and numerically. In particular, it was found that the kink-Q-ball system possesses a nonzero electric field and is unstable with respect to small perturbations of fields.


Introduction
It is known that in the case of Maxwell electrodynamics, any one-dimensional or two-dimensional field configuration with a nonzero electric charge possesses infinite energy. The reason is simple: at large distances, the electric field of such a field configuration does not depend on coordinate in the one-dimensional case and behaves as r −1 in the two-dimensional case, so the energy of the electric field diverges linearly in the one-dimensional case and logarithmically in the two-dimensional case. That is why there are no electrically charged solitons in one and two dimensions; such solitons appear only in three dimensions (e.g. three-dimensional electrically charged dyon [1] or Q-ball [2,3,4,5]). It should be noted, however, that electrically charged two-dimensional vortices exist in both the Chern-Simons [6,7,8,9,10] and the Maxwell-Chern-Simons [11,12,13,14] gauge models. Furthermore, it was shown in [15,16] that Chern-Simons gauge models also possess one-dimensional domain walls. The domain walls have finite linear densities of magnetic flux and electric charge, so there is a linear momentum flow along the domain walls.
Nevertheless, even in the case of Maxwell gauge field models, there are electrically neutral low-dimensional soliton systems with a nonzero electric field in their interior areas. In particular, the one-dimensional soliton system consisting of electrically charged Q-ball and anti-Q-ball components has been considered in [17] and the two-dimensional soliton systems consisting of vortex and Q-ball components interacting through an Abelian gauge field has been described in [18].
In the present paper, we research the one-dimensional soliton system consisting of Q-ball and kink components, possessing opposite electric charges, so the system with a nonzero electric field is electrically neutral as a whole. Properties of this kink-Q-ball system are investigated by analytical and numerical methods. In particular, we found that in contrast to the non-gauged onedimensional Q-ball, the kink-Q-ball system does not go into the thin-wall regime.
There is an interesting problem concerning the stability of the kink-Q-ball system with respect to small perturbations of fields. Recall that the Abelian Higgs model possesses an electrically neutral kink solution [19,20]. Formally, this gauged kink solution is the usual kink of a self-interacting real scalar field up to gauge transformations. However, properties of these two kink solutions differ considerably, because the classical vacua of the corresponding field models have a different topol-2 ogy. While the real kink is topologically stable, the gauged kink has exactly one unstable mode. From the topological point of view, the gauged kink lies between two topologically distinct vacua of the Abelian Higgs model, so it is a sphaleron [19,20]. Note, however, that the gauged kink is a static field configuration modulo gauge transformation, whereas the kink-Q-ball system will depend on time in any gauge. Due to this fact, the kink-Q-ball system cannot be a sphaleron, so its classic stability requires separate consideration.
The paper has the following structure. In Sec. 2, we describe briefly the Lagrangian, the symmetries, the field equations, and the energy-momentum tensor of the Abelian gauge model under consideration. In Sec. 3, we research properties of the kink-Q-ball system. Using the Hamiltonian formalism and the Lagrange multipliers method, we establish the time dependence of the soliton system's fields. The important differential relation for the kink-Q-ball solution is derived and the system of nonlinear differential equations for ansatz functions is obtained. Then we establish some general properties of the kink-Q-ball system. In particular, we research asymptotic behavior of the system's fields at small and large distances, establish some important properties of the electromagnetic potential, and derive the virial relation for the soliton system. In Sec. 4, we study properties of the kink-Q-ball system in the three extreme regimes, namely, in the thick-wall regime and in the regimes of small and large gauge coupling constants. We also establish basic properties of the plane-wave field configuration of the model. In Sec. 5, we present and discuss the numerical results obtained. They include the dependences of the energy of the kink-Q-ball system on its phase frequency and Noether charge, along with numerical results for the ansatz functions, the energy density, the electric charge density, and the electric field strength.
Throughout the paper, we use the natural units = c = 1.

The gauge model
The Lagrangian density of the (1 + 1)-dimensional gauge model under consideration has the form where F µν = ∂ µ A ν −∂ ν A µ is the strength of the Abelian gauge field and φ, χ are complex scalar fields minimally interacting with the Abelian gauge field through the covariant derivatives: The self-interaction potentials of the scalar fields have the form Let us suppose that the self-interaction potential U (|χ|) has a global zero minimum at χ = 0 and admits the existence of usual non-gauged Q-balls. Then the parameters g χ and h χ are positive and satisfy the condition 3g 2 χ < 16h χ m 2 χ . Unlike the sixth-order potential U (|χ|), the fourth-order potential V (|φ|) reaches a zero minimum on the circle |φ| = η. The potential V (|φ|) allows the existence of the complex non-gauged kink solution where m φ = √ 2λη is the mass of the scalar φ-particle and δ is an arbitrary phase.
Besides the local gauge transformations: the Lagrangian (1) is also invariant under the two independent global gauge transformations: As a consequence, we have the two Noether currents: and the two separately conserved Noether charges: Q φ = j 0 φ dx and Q χ = j 0 χ dx. Note also that in addition to the local and global gauge transformations, the Lagrangian (1) is invariant under the discrete C, P , and T transformations.
The field equations of the model are written as where the electromagnetic current j ν is written in terms of two Noether currents (8) From Eq. (12) it follows that the electric charges eQ φ and qQ χ of the complex scalar fields φ and χ are conserved separately. This fact is a consequence of the neutrality of the Abelian gauge field A µ .

3
The symmetric energy-momentum tensor of the model can be obtained by using the well-known formula T µν = 2∂L/∂g µν − g µν L: Thus, we have the following expression for the energy density of the model where 3 The kink-Q-ball system and its properties It is known [23,21] that any nontopological soliton, in particular, a Q-ball, is an extremum of an energy functional at a fixed value of the corresponding Noether charge. Using this basic property of a Q-ball and taking into account that the self-interaction potential U (|χ|) admits the existence of Q-balls formed from the complex scalar field χ, we shall search for a soliton solution of model (1) that is an extremum of the energy functional E = ∞ −∞ Edx at a fixed value of the Noether charge Q χ = ∞ −∞ j 0 χ dx. According to the method of Lagrange multipliers, such a solution is an unconditional extremum of the functional where ω is the Lagrange multiplier. Let us determine the time dependence of the soliton solution. To do this, we shall use Eq. (15) and the Hamiltonian formalism. In the axial gauge A x = A 1 = 0, the Hamiltonian density of model (1) is written as We see that in the adopted gauge, the model is described in terms of the eight canonically conjugated fields: φ, π φ = (D 0 φ) * , φ * , π φ * = D 0 φ, χ, π χ = (D 0 χ) * , χ * , and π χ * = D 0 χ, while the time component A 0 of the gauge field is determined in terms of the canonically conjugated fields by Gauss's law (17) and so it is not an independent dynamic field. Although energy density (14) is not equal to Hamiltonian density (16): the integral of Eq. (18) over the one-dimensional space vanishes provided that field configurations of the model possess finite energy and satisfy Gauss's law (17).
In the adopted axial gauge A x = 0, field equations (10) and (11) can be recast in the Hamiltonian form: where we use the relation E = On the other hand, the first variation of functional (15) vanishes on the soliton solution: where the first variation of the Noether charge Q χ is expressed in terms of the canonically conjugated fields as follows: Combining Eqs. (19) - (22), we find that in the adopted gauge, only the time derivatives of the canonically conjugated fields χ, π χ , χ * , and π χ * are different from zero: while the time derivatives of φ, π φ , φ * , and π φ * are equal to zero. Recalling that π χ = (D 0 χ) * = ∂ t χ * − iqA 0 χ * and taking into account Eqs. (24) and (25), we conclude that the time derivative of A 0 also vanishes. It follows that only the scalar field χ of the soliton system has nontrivial time dependence: Let us return to Eq. (21). This equation holds for arbitrary variations of fields on the soliton solution, including those that transfer the soliton solution to an infinitesimally close one. It follows that the energy of the soliton system satisfies the important relation where the Lagrange multiplier ω is some function of the Noether charge Q χ . Since the energy E and the Noether charge Q χ of the soliton system are gauge-invariant, relation (28) is also gauge-invariant. Like the case of non-gauged nontopological solitons [21,22,23], relation (28) determines basic properties of the gauged kink-Qball system. In Eqs. (27), functions f (x) and s (x) are assumed to be some complex functions of the real argument x. Substituting Eqs. (27) into field equations (9) -(11), we can easily check that the real and imaginary parts of f (x) satisfy the same differential equation with real coefficients. Similarly, the real and imaginary parts of s (x) also satisfy the same differential equation with real coefficients. It follows that the functions f (x) and s (x) have the form: f (x) = exp (iα)f (x), s (x) = exp (iβ)s (x), wheref (x) ands (x) are real functions, whereas α and β are constant phases. However, these phases can be cancelled by global gauge transformations (7), so the functions f (x) and s (x) can be supposed to be real without loss of generality. The functions a 0 (x), f (x), and s (x) satisfy the system of ordinary nonlinear differential equations: which is obtained by substituting Eqs. (27) into field equations (9) - (11). The most important among the local quantities of the kink-Q-ball system are the electromagnetic current density and the energy density. Their expressions in terms of a 0 (x), f (x), and s (x) are written as The energy E = ∞ −∞ Edx of the kink-Q-ball system must be finite. Using this fact and Eq. (33), we obtain the boundary condition for a 0 (x), f (x), and s (x): Note that the finiteness of the electric field's energy 0 /2dx leads to one more boundary condition for a 0 (x): This condition, however, is equivalent to Eq. (34a) provided that a 0 (x) is regular as x → ±∞. Gauss's law (29) can be written as a ′′ 0 = −j 0 , where j 0 is electric charge density (32). Integrating this equation over x ∈ (−∞, ∞) and taking into account boundary conditions (35), we conclude that the total electric charge Q = ∞ −∞ j 0 dx of a field configuration with finite energy vanishes: where Q φ and Q φ are the Noether charges defined by Eqs. (8).
It can easily be checked that system (29) This invariance is a consequence of the C-invariance of the Lagrangian (1). Using Eqs. (32), (33), and (37), we find the behavior of the energy E and the Noether charges Q φ and Q χ under the transformation ω → −ω: We see that the energy of the kink-Q-ball system is an even function of ω, whereas the Noether charges Q φ and Q χ are odd functions of ω.
The P -invariance of the Lagrangian (1) leads to the invariance of system (29) -(31) under the space inversion x → −x. Due to the space homogeneity, the system (29) -(31) is also invariant under the coordinate shift x → x + x 0 . Furthermore, due to Eqs. (7), the system (29) -(31) is invariant under the two independent discrete transformations: f → −f and s → −s. These facts and symmetry properties of boundary conditions (34) lead to the conclusion that a 0 (x) and s (x) are even functions of x, while f (x) is an odd function of x. This is consistent with the fact that the non-gauged kink solution is an odd function of x, whereas the non-gauged Q-ball solution is an even function of x.
The asymptotic form of the soliton solution for small x is obtained by substitution of the power expansions for a 0 (x), f (x), and s (x) into Eqs. (29) -(31) and equating the resulting Taylor coefficients to zero. By acting in this way, we obtain: where the next-to-leading coefficients are expressed in terms of the three leading coefficients a 0 , f 1 , s 0 , and the model's parameters. For large |x|, system (29) -(31) is linearized and we obtain the asymptotic form of the soliton solution satisfying boundary conditions (34): Let us discuss the global behavior of the electromagnetic potential a 0 (x). Since the total electric charge Q = +∞ −∞ j 0 (x) dx of the kink-Q-ball system vanishes, the electric charge density j 0 (x) must vanish at some points of the x-axis. Because of the symmetry j 0 (−x) = j 0 (x), these points (nodes of j 0 (x)) are symmetric with respect to the origin x = 0. Next, according to Gauss's law a ′′ 0 (x) = −j 0 (x), the second derivative a ′′ 0 (x) vanishes at the nodes of j 0 (x). Thus the nodes of j 0 (x) are the inflection points of the electromagnetic potential a 0 (x). From Eq. (29) it follows that in an inflection point x i , the electromagnetic potential a 0 (x i ) can be expressed in terms of f (x i ) and s (x i ): Two conclusions follow from Eq. (43). Firstly, at an inflection point x i , the sign of a 0 (x i ) coincides with the sign of ω (we suppose that the gauge coupling constants are positive by definition): Secondly, at an inflection point x i , the following inequality holds: Next, from Eq. (32), we obtain the expression for the electric charge density at the origin: from which it follows that the sign of the curvature of a 0 (x) at x = 0 is opposite in sign to ω − qa 0 (0): An elementary graphical analysis made using Eqs. (43) -(47) leads us to the following conclusions about the behavior of a 0 (x): and where ±x i1 is the two symmetric inflection points closest to the origin x = 0. From Eqs. (46), (48), and (49) it follows that the sign of the electric charge density at the origin coincides with that of the phase frequency We can also make conclusions about the behavior of a 0 (x) for |x| > x i1 . In particular, a 0 (x) cannot vanish at any finite x. Indeed, let x n be a conjectural point in which a 0 (x) vanishes. Then from Eq. (29) we have the relation We see that the sign of the curvature of a 0 (x) at the point x n is opposite to the sign of ω. Let ω be positive. Then from Eq. (51) it follows that in some neighborhood of x n , the functions a 0 (x) and a ′′ 0 (x) are negative. But according to Eq. (44), there are no inflection points for negative a 0 (x), so a ′′ 0 (x) can never change the sign, a 0 (x) decreases indefinitely, and boundary condition (34a) cannot be satisfied. It follows that a 0 (x) cannot vanish at any finite x. The case of negative ω is treated similarly. Thus, we come to an important conclusion that the electromagnetic potential a 0 (x) cannot vanish at any finite x, and so the sign of the electromagnetic potential coincides with that of the phase frequency over the whole range of x: sign (a 0 (x)) = sign (ω) (52) 6 for all x. Of course, this conclusion is valid only for adopted gauge (27c). Let a 0 (x), f (x), and s (x) be a solution of system (29) -(31) that satisfy boundary conditions (34). When we perform the scale transformation x → λx of the argument of the solution, the Lagrangian L = ∞ −∞ Ldx becomes a simple function of the scale parameter λ. The function L (λ) must have an extremum at λ = 1, so the derivative dL/dλ vanishes at this point. Using this fact, we obtain the virial relation for the soliton system: where is the electric field's energy, is the gradient part of the soliton's energy, is the kinetic part of the soliton's energy, and is the potential part of the soliton's energy.
The energy E of the soliton system is the sum of terms (54) -(57). Using this fact and virial relation (54), we obtain the two representations for the soliton system's energy: Another representation for the soliton system's energy can be obtained by integrating the term a ′ Finally, using Eq. (59), we obtain the relation between the Noether charge Q χ , the electric field's energy E (E) , and the kinetic energy E (T ) : 4 Extreme regimes of the kink-Q-ball system In this section, we will first study the properties of the kink-Q-ball system in the thick-wall regime [24,25,26]. In this regime, the parameter ∆ = m 2 χ − ω 2 1/2 tends to zero, so the absolute value of phase frequency tends to m χ . From Eqs. (42) it follows that in the thick-wall regime, the functions s (x) and a 0 (x) are spread over the one-dimensional space, whereas the asymptotic behavior of f (x) remains unchanged. In the thick wall regime, functions s (x) and a 0 (x) uniformly decrease as ∆ and ∆ 2 , respectively, while the function f (x) tends to non-gauged kink solution (5). For this reason, we perform the following scale transformation of the fields and x-coordinate: while the field f (x) is taken equal to that of kink solution (5). To research the properties of the kink-Q-ball system in the thick-wall regime, we shall use functional (15) that is related to the energy functional through the Legendre transformation: Using scale transformation (61), we can determine the leading term of the dependence of the functional F (ω) on ω in the thick-wall regime: where E k = 4η 3 √ 2λ/3 = 4η 2 m φ /3 is the rest energy of the non-gauged kink and the dimensionless functional F does not depend on ω: In Eq. (62), higher-order terms in ∆ may be neglected in the thick-wall regime, so we obtain sequentially: where known properties of Legendre transformation are used. Using Eqs. (64) and (65), we obtain the energy of the kink-Q-ball system as a function of its Noether charge in the thick-wall regime: It has been found numerically that the kink-Q-ball system does not turn into the thin-wall regime as the magnitude of the phase frequency tends to its minimum 7 value. Such behavior can be qualitatively explained as follows. Gauss's law a ′′ 0 (x) = −j 0 (x) has an obvious mechanical analogy. It describes a one-dimensional motion of a unit mass particle along the coordinate a 0 in time x. The particle is subjected to the time-dependent force −j 0 (x). It starts to move at the time x = 0 from the point a 0 (0) > 0 (we suppose that ω > 0) with zero initial velocity a ′ 0 (0). Since j 0 (x) > 0 as x < x i (x i is the inflection point of a 0 (x)), the particle's coordinate a 0 is decreased with increasing of x. At the time x i , the particle is at the point a 0 (x i ) given by Eq. (43) and possesses the velocity a ′ 0 (x i ) = −Q + /2, where Q + = xi −xi j 0 (x) dx is the positive electric charge of the central area of the soliton system. The negative electric charge ∞ xi j 0 (x) dx of the side area x > x i corresponds to the impulse of force acting in the positive direction. According to boundary conditions (34a) and (35), this impulse of force must bring the particle to rest at the point a 0 = 0 as time x tends to infinity.
As the magnitude of the phase frequency ω approaches to the minimum value, the spatial size of the kink-Q-ball system increases, whereas the positive electric charge of the central area weakly depends on ω. As a result, the velocity of the particle at the inflection point x i remains approximately constant, whereas the decelerating force −j 0 (x) decreases because of the spreading of the impulse of force − ∞ xi j 0 (x) dx (which is equal to half the positive electric charge Q + of the central area) over the side area [x i , ∞). Thus the decelerating force j 0 (x) decreases and so it cannot stop the particle. Because of that, the particle reaches the point a 0 = 0 for a finite period of time, so boundary conditions (34a) and (35) are not satisfied and the kink-Q-ball system cannot exist.
Next, we consider the regime of small gauge coupling constants. When the gauge coupling constants e and q vanish, the gauge field A µ = (a 0 (x) , 0) is decoupled from the kink-Q-ball system, which thus becomes the set of non-gauged kink and Q-ball that do not interact with each other. In this connection, we want to ascertain the behavior of the gauge potential a 0 (x) as e = ̺q → 0, where ̺ is a positive constant. To do this, we use asymptotic expressions (40a) and (42c) for a 0 (x) that are valid for small and large values of |x|, respectively. We also suppose that Eqs. (40a) and (42c) describe qualitatively the behavior of a 0 at intermediate |x|. Eqs. (40a) and (42c) depend on the free parameters a 0 and a ∞ , respectively. These parameters can be determined by the condition of continuity of a 0 (x) and a ′ 0 (x) at some intermediate x.
As a result, a 0 and a ∞ become functions of the model's parameters, including the gauge coupling constants e and q. It can be shown that a 0 and a ∞ tend to the same nonzero limit as e = ̺q → 0. It follows that at any finite x, the gauge potential a 0 (x) tends asymptotically to a constant as e = ̺q → 0: where the limit value α depends on the model's parameters and ̺. This fact and linearization of Eqs. (30) and (31) lead us to the asymptotic forms of f (x) and s (x) at small gauge coupling constants: where f k (x) and s q (x) are the non-gauged kink and Q-ball solutions, respectively, whereas f 2 (x) and s 1 (x) are some regular functions that depend on the model's parameters (except for e and q) and ̺. Note that due to the relation e = ̺q, we use only one expansion parameter e. Substituting Eqs. (67) and (68) into Eq. (32), we find the asymptotic behavior of the Noether charges Q φ and Q χ as e = ̺q → 0: where Q 0 = 2ω +∞ −∞ s 2 q dx is the Noehter charge of the non-gauged Q-ball and the coefficient Q 1 depends on the model's parameters (except for e and q) and ̺. Similarly to Eq. (69), we obtain the asymptotic behavior of the soliton energy's components (54) -(57) and the total soliton energy: where E dx are the gradient, kinetic, and potential parts of the non-gauged soliton system's energy, respectively, is the total energy of the non-gauged soliton system, and the coefficients E Note that the non-gauged solutions f k and s q can be expressed in analytical form, as well as the corresponding energies and the Noether charges. The corresponding expressions for the one-dimensional non-gauged Q-ball are given in [17]. Thus the coefficients Q 0 , E 0 , and E 0 can also be expressed in analytical form.
Next let us consider the opposite regime in which both gauge coupling constants tend to infinity: e = 8 ̺q → ∞. We suppose that the behavior of the electromagnetic potential a 0 (x) in a neighborhood of x = 0 is regular, so the coefficient a 2 given by Eq. (41a) is either finite or tending to zero as e = ̺q → ∞. But from Eqs. (29) and (31) it follows that the electromagnetic potential is an odd function of q: a 0 (x, −q) = −a 0 (x, q). Thus, we conclude that in the leading order, a 2 ∝ q −1 , so from Eq. (41a) it follows that a 0 ∼ ωq −1 − a −3 q −3 , where a −3 is a positive constant. This fact suggests that the electromagnetic potential a 0 (x) has a similar asymptotic expansion in the inverse powers of e: where a −1 (x) and a −3 (x) are some regular functions depending on the model's parameters (except for e and q) and ̺. Using Eqs. (30), (31), and (71), we obtain the general form of asymptotic expansions for f (x) and s (x): where f 0 (x), f −2 (x), s 0 (x), and s −2 (x) are regular functions depending on the model's parameters (except for e and q) and ̺. Similar to Eqs. (67) and (68), we can use Eqs. (71) and (72) to obtain the asymptotic expansions for the soliton energy's components, the total energy, and the Noether charges: where we use the tilde to distinguish corresponding coefficients from those of Eqs. (70a) -(70d). We see that as e = ̺q → ∞, the gauge field a 0 (x) tends to zero, so the electric field's energy E (E) also vanishes in this regime. At the same time, the products ea 0 (x) and qa 0 (x) tend to nonzero limits a −1 (x) and ̺ −1 a −1 (x), respectively, so the gauge field a 0 (x) does not decouple from the kink-Q-ball system. Due to this, the limit solutions f 0 (x) and s 0 (x) are different from the corresponding non-gauged solutions f k (x) and s q (x), respectively. From Eqs. (73b) -(73f) it follows that the soliton energy's components E (G) , E (T ) , E (P ) , the total soliton energy E, and the Noether charges Q χ and Q φ also tend to some finite values as e = ̺q → ∞. It follows that the electric charges of the kink and the Q-ball increase indefinitely in this regime, despite the fact that the electric field's energy E (E) tends to zero. This is because the electric charges of the kink and the Q-ball tend to cancel each other at any spatial point as e = ̺q → ∞. Note that the behavior of the kink-Q-ball system in the extreme regimes e = ̺q → 0 and e = ̺q → ∞ was investigated by numerical methods. It was found that it is in accordance with Eqs. (69), (70), and (73). Finally, we consider the plane-wave solution of gauge model (1). In this case, the gauge field A µ and the scalar fields φ and χ spread over the one-dimensional space and fluctuate around their vacuum values. Since the scalar field φ has nonzero vacuum value |φ vac | = η, the classical vacuum of model (1) is not invariant under local gauge transformations (6), so the local gauge symmetry is spontaneously broken. For this reason, the research of the plane-wave solution is convenient to perform in the unitary gauge Im (φ (x, t)) = 0. In this gauge, the Higgs mechanism is realized explicitly, so we can read off the particle composition of model (1). In the neighborhood of the gauge vacuum φ vac = η, χ vac = 0, we have the complex scalar field χ with the mass m χ , the real scalar Higgs field φ H with the mass m φ = √ 2λη, and the massive gauge field A µ with the mass m A = √ 2eη. We want to find the spatially uniform solution of field equations (9) -(11) possessing the Noether charges Q φ and Q χ (recall that eQ φ + qQ χ = 0 for any finite energy field configuration) and to determine its energy. For this, we use field equations (9) -(11) in the unitary gauge. We suppose that with unlimited spreading, the amplitudes of the complex scalar field χ and the real scalar Higgs field φ H tend to zero, so we can neglect higher-order terms in the Lagrangian (1). On spatially uniform fields, the field equations for A µ and φ H become algebraic ones, whereas the field equation for χ determines the time dependence of χ. The results obtained are presented as series in inverse powers of the plane-wave solution's spatial size L: We see that as L → ∞, the amplitudes of the fields A 0 , φ H , and χ tend to zero, whereas the phase frequency 9 of χ tends to m χ . Note that the fields A 0 , φ H , and χ of the plane-wave solution tend to zero as L −1 , L −2 , and L −1/2 , respectively, so the Higgs field φ H tends to zero much more quickly than the complex scalar field χ. Substituting Eqs. (74) -(76) into Eqs. (8), we obtain the Noether charge densities j 0 φ and j 0 χ of the planewave solution: From Eqs. (77) and (78) it follows that the electric charge of the plane-wave solution vanishes: as it should be. Next, we calculate the energies of the φ and χ components of the plane-wave solution: and so the total energy of the plane-wave solution turns out to be equal to Let us discuss the results obtained. First of all, from Eqs. (80) -(82) it follows that the φ component does not contribute to the energy of the plane-wave solution as L → ∞. At the same time, from Eqs. (77) and (79) it follows that the electric charge of the φ component does not vanish and is opposite to that of the χ component, so the total electric charge of the plane-wave solution vanishes. Thus, the φ component contributes to the electric charge of the plane-wave solution, but does not contribute to its energy. This can be explained as follows. In the unitary gauge Im (φ) = 0, the real Higgs field φ H fluctuates around the real vacuum average η, so the initial scalar field φ is written as φ = η + φ H . Further, from Eqs. (8a) and (14) we obtain the electric charge and energy densities of the φ component of the plane-wave solution: and where Eqs. (74) and (75) have been used. We see that in the leading order in L −1 , the Higgs field φ H contributes neither to j 0 φ nor to E φ , whereas the vacuum average of the scalar field φ contributes in both cases. We also see that as L → ∞, the behavior of j 0 φ and E φ is determined only by the electromagnetic potential A 0 . At the same time, from Eq. (74) it follows that for the plane-wave solution, so the φ component does not contribute to the plane-wave solution's energy as L → ∞.

Numerical results
To study the kink-Q-ball system, we must solve system of differential equations (29) -(31) satisfying boundary conditions (34). This first boundary value problem can be solved only numerically. To solve the boundary value problem, we use the method of finite differences and subsequent Newtonian iterations realized in the Maple package [27]. To check the correctness of a numerical solution, we use Eqs. (28), (36), and (53).
Let us denote by ∆E the difference between the energies of the kink-Q-boll system and the non-gauged kink: Figures 1  and 2 present the dependence of the dimensionless energy difference ∆Ẽ = ∆E/m χ on the dimensionless phase frequencyω. The curves in these figures correspond to gauge coupling constantsẽ =q taking the values from the set 0.1, 0.2, 0.3, 0.4, and 0.5. Figure 1 presents the curves in the range from the minimum values ofω, which we managed to reach by numerical methods, to the valueω = 0.88. Figure 2 presents the same curves in the range fromω = 0.88 to the maximum possible valueω = 1. We use the two figures for a better representation of the dependences ∆Ẽ (ω). For the same values of gauge coupling constants, the dependences ∆Ẽ (ω) and Q χ (ω) are qualitatively similar, so the dependences Q χ (ω) are not given in the present paper. Let us discuss the main features of the curves in Figs. 1 and 2. First of all, we note that the energy of the kink-Q-ball system does not tend to infinity asω tends to its minimum values (that depend on the gauge coupling constants). Indeed, it was found numerically that the dependencesẼ (ω) and Q χ (ω) have a branching point atω min : where A, B, and C are positive constants. We were unable to find any solutions of the boundary value problem forω <ω min , so we conclude that the kink-Q-ball system does not turn into the thin-wall regime in which bothẼ and Q χ must tend to infinity. The behavior of the curves in the neighborhood of ω = 1 is also rather unusual. We see that for the gauge coupling constantsẽ =q from the set 0.2, 0.3, 0.4, and 0.5, the dependenceẼ (ω) consists of two separate curves. The left curve starts from the minimal phase frequencyω min (where it has branching point (85)) and continues until the maximal phase frequencyω r . The behavior of the left curve in the neighborhood ofω r is similar to that in the neighborhood ofω min : where D, F , and G are positive constants. The right curve starts from some phase frequencyω l <ω r and The curve has no singularity in the neighborhood ofω l . According to Seq. 5, the kink-Q-ball system goes into the thick-wall regime asω → 1. Indeed, it was found numerically that in the neighborhood ofω tk = 1, ∆Ẽ ∼ Q χ ∼ H (ω tk −ω) 1/2 in accordance with Eqs. (64) and (65). Figure 3 shows the dependences ∆Ẽ (ω) forẽ =q = 0.05 andẽ =q = 0.1 in the neighborhood ofω tk = 1. We see that with decreasing gauge coupling constants, the left and right curves are merged into one, so the de-  Figs. 1 and 2. The second derivative d 2Ẽ /dQ 2 χ changes the sign when passing through the cusps or discontinuities, so convex and concave sections of the curves change each other. Note that in Figs. 4 and 5, the energy of the kink-Q-ball system turns out to be more than the energy of the plane-wave solution with the same value of Q χ . This also turns out to be true for all other cases considered in the present paper. It follows that the kink-Q-ball system may transit into the plane-wave field configuration through quantum tunneling. Figure 6 shows the kink-Q-ball solution with the dimensionless phase frequencyω = 0.3. The energy density, the electric charge density, and the electric field strength corresponding to this kink-Q-ball solution are presented in Fig. 7. We see that in accordance with Seq. 3, a 0 (x) and s (x) are even functions of x, whereas We also see that f (x) and s (x) reach neighborhoods of their boundary values (34) faster than a 0 (x). From Fig. 7 it follows that the kink-Q-ball system possesses the symmetric energy and electric charge densities. It also possesses the nonzero electric field strength that is an odd function of the space coordinate. The distribution of the electric charge density is a central symmetric peak with a positive j 0 surrounded by two areas with a negative j 0 , so the total electric charge of the kink-Q-ball system vanishes. The central positive peak is due to the contribution of the field χ, whereas the two side negative areas are due to the contribution of the field φ. Note that the energy density E reaches a close neighborhood of zero faster than the electric charge density j 0 and the electric field strength E x . The reason is the similar behavior of the electromagnetic potential a 0 in Fig. 6. Indeed, from Eq. (32) it follows that the electric charge density of the field φ is −2a 0 e 2 f 2 . We see that the electromagnetic potential a 0 can induce a nonzero electric charge density even if the scalar field φ reaches a close neighborhood of the vacuum value |η|. As a result, a substantial part of the electric charge of the complex scalar field φ comes from the parts of the two side areas where |φ| ≈ |η|, χ ≈ 0, and E ≈ 0.
Let us now discuss the issue of stability of the kink-Q-ball system. As has already been pointed out, the energy of the kink-Q-ball system turned out to be more than the energy of the plane-wave solution with the same Q χ for all cases considered in the present paper. It follows that the kink-Q-ball system is unstable against transit into a plane-wave configuration through quan- It is known that the gauge model described by the first line of the Lagrangian (1) possesses the kink solution [19,20]. In the adopted gauge A x = 0, this kink solution is given by Eq. (5). The gauged kink has zero electric charge, so it possesses finite energy. However, unlike the kink of a self-interacting real scalar field [28,29], the gauged kink is not a topologically stable field configuration. Due to the topological structure of the vacuum of the Abelian Higgs model, the gauged kink is a sphaleron [19,20]. The existence of the sphaleron is due to the paths in the functional space that connect topologically distinct vacua of the Abelian Higgs model [30]. The sphaleron lies between two topologically distinct neighboring vacua and has exactly one unstable mode.
The gauged kink is a static solution modulo gauge transformations. However, in the case of the kink-Q-ball solution, we have a different situation. It can easily be shown that the kink-Q-ball solution will depend on time in any gauge, so it is not a static solution. It follows that the point in the functional space corresponding to the kink-Q-ball solution will vary with time in any gauge. This fact does not allow the kink-Q-ball solution to be a sphaleron, so the question about the unstable modes of the kink-Q-ball solution should be investigated separately.
To investigate the classic stability of the kink-Q-ball system, it is necessary to study the spectrum of the operator of second variational derivatives in the functional neighborhood of the kink-Q-ball solution. Wherein, the model's fields must fluctuate so that the Noether charges Q φ and Q χ remain fixed and the perturbed electromagnetic potential A 0 + δA 0 continues to satisfy Gauss's law. All these factors make it difficult to study the spectrum even using numerical methods. However, these difficulties can be avoided if we numerically solve field equations (9) -(11) in the temporal gauge A 0 = 0 and with a perturbed initial field configuration in the close neighborhood of the kink-Q-ball solution. Indeed, Gauss's law can be easily implemented in the temporal gauge at t = 0, whereupon Gauss's law will be automatically satisfied for t > 0. The perturbed initial field configuration must have the same Q φ and Q χ as the kink-Q-ball system; then the field equations guarantee that the field configuration will also have the same Q φ and Q χ at later times. Having perturbed and unperturbed kink-Q-ball solutions of the field equations, we can observe how field fluctuations behave as time increases. If any fluctuation of fields oscillates in a close neighborhood of the kink-Q-ball solution then the solution is classically stable. If there exists at least one fluctuation of fields that increases exponentially with time then the kink-Q-ball solution is classically unstable.
In the present paper, we research the stability of the kink-Q-ball system forẽ =q = 0.05, 0.1, 0.2, 0.3, 0.4, and 0.5. For each value of the gauge coupling constants, we took the step of changing ofω equal to 0.1. To solve field equations (9) -(11), we use the solver of partial differential equations realized in the Maple package [27]. We have found that there are at least two unstable modes for the all considered gauge coupling constants and phase frequencies. The first unstable mode corresponds to an initial symmetric perturbation of Imφ, whereas the second one corresponds to an initial antisymmetric perturbation of Imχ. It follows that the kink-Q-ball system is not a sphaleron since it has at least two unstable modes, whereas a sphaleron must have exactly one unstable mode.

Conclusion
In the present paper, we consider one-dimensional model (1) consisting of two self-interacting complex scalar fields interacting through an Abelian gauge field. It was shown that the model possesses the soliton solution consisting of a gauged kink and a gauged Q-ball. Since the finiteness of the energy of the one-dimensional soliton system leads to its electric neutrality, the gauged kink and the gauged Q-ball have opposite electric charges. Due to the neutrality of the Abelian gauge field, the opposite electric charges of the kink and Q-ball components are conserved separately. Despite the neutrality of the kink-Q-ball system, it possesses a nonzero electric field.
The kink-Q-ball system has rather unusual dependences of the energy and the Noether charge on the phase frequency. Indeed, it was found that the energy and the Noether charge of the kink-Q-ball system do not tend to infinity as the phase frequency tends to its minimum value, but instead have the branch point. It follows that there is no thin-wall regime for the kink-Qball system. We also found that when the magnitude of the phase frequency is in the neighborhood of m χ , the dependences of the energy and the Noether charge on the phase frequency consist of two separate branches provided that the model's gauge coupling constants are large enough. In all cases, however, the kink-Q-ball system goes into the thick-wall regime as the magnitude of the phase frequency tends to m χ .
In addition to the kink-Q-ball solution, the model also possesses a plane-wave solution. For all sets of the model parameters considered in the present paper, the energy of the kink-Q-ball solution turns out to be more than the energy of the plane-wave solution with the same value of the Noether charge. Due to the topological structure of the model's vacuum, the kink-Q-ball solution is not topologically stable, so it can transit into the plane-wave configuration through quantum tunneling.
It is known that the Abelian Higgs model possesses the gauge kink solution. The gauge kink is electrically neutral and has one unstable mode. From the viewpoint of topology, the gauge kink is a static (modulo gauge transformations) field configuration lying between the two topologically distinct adjacent vacua. Unlike this, the kink-Q-ball solution depends on time in any gauge, so it is not a static field configuration. Hence, the kink-Q-ball solution cannot be a sphaleron, and so the question of its classic stability requires separate consideration. We research the classic stability of the kink-Q-ball system by means of numerical solution of the field equations with initial field configurations perturbed in the close neighborhood of the kink-Q-ball solution. It was found that in all considered cases, the kink-Q-ball solution has at least two unstable modes, so it is even more unstable than the gauged kink.