Classical behaviour of Q-balls in the Wick-Cutkosky model

In this paper, we continue discussing Q-balls in the Wick--Cutkosky model. Despite Q-balls in this model are composed of two scalar fields, they turn out to be very useful and illustrative for examining various important properties of Q-balls. In particular, in the present paper we study in detail (analytically and numerically) the problem of classical stability of Q-balls, including the nonlinear evolution of classically unstable Q-balls, as well as the behaviour of Q-balls in external fields in the non-relativistic limit.


Introduction
Among various models admitting solutions of the Q-ball [1,2] type, 1 there exists a class of models with two different scalar fields (the complex one and the real one), which also admit the existence of the Q-ball-like solutions. The most known example of such a Q-ball is presented in the well-known paper [6]. One can also recall a simplified version of the model of [6], in which the potential of the real scalar field is neglected [7] (see also [8]). In paper [9] the simplest example of the two-field Q-ball was presented: instead of the quartic interaction of the scalar fields of [7], the triple Yukawa interaction between the scalar fields, resulting in the well-known Wick-Cutkosky model [10,11], was considered. A remarkable feature of the model is that the main characteristics of Q-balls can be obtained analytically, including the form of the energy-charge dependence. The latter simplifies the subsequent analysis of the model.
Despite Q-balls in this model are composed of two scalar fields, the Wick-Cutkosky model turns out to be a simple and useful toy model for examining various properties of Q-balls and for testing the methods, which can be applied to Q-balls in other theories (including theories with the standard one-field Q-balls). A part of the analysis was performed in paper [9], in the present paper we focus on the problem of classical stability of Q-balls and on the behaviour of Q-balls in external fields.
It should be noted that the Schrödinger-Poisson system, which appears, for example, when one considers the Newtonian limit for the Bose stars made of scalar fields [12][13][14], is very similar to the Wick-Cutkosky model. Thus, we think that the results obtained for the case of the Wick-Cutkosky model can be useful for examining the Bose stars too.
First, we will examine the classical stability of Q-balls in the Wick-Cutkosky model. The criterion of classical stability of two-field Q-balls was established in paper [6]. The proof presented in that paper is based on examining the properties of the energy functional of the system, while keeping the charge fixed. An alternative proof along the lines of the Vakhitov-Kolokolov method [15,16], which is based on the use of only the linearized equations of motion, can be found in [17]. The classical stability criterion states that a Q-ball is stable with respect to small perturbations (classically stable) if dQ dω < 0 for this Q-ball and if a certain operator (we will denote this operator as L + ), which arises in the linearized equations of motion for perturbations, has only one negative eigenvalue. Usually, the operator L + indeed has only one negative eigenvalue. However, it is not a rule in the general case (except for the standard one-field Q-balls in (1+1)-dimensional space-time), so in any particular case one should check that the number of negative eigenvalues of the operator L + does not exceed one in order to be sure that Q-balls with dQ dω < 0 are classically stable. Below we will show explicitly that the latter condition holds for Q-balls in the Wick-Cutkosky model.
Second, for the case of classically unstable Q-balls, we will find explicitly the exponentially growing (instability) mode which is a solution to the corresponding linearized equations of motion. Next, we will examine the nonlinear evolution of classically unstable Q-balls by simulating numerically the evolution of the perturbed Q-ball. We will show that, depending on the characteristics of the initial unstable Q-ball and on the contribution of the instability mode into the initial perturbation, unstable Q-balls in the Wick-Cutkosky model can evolve in three completely different ways (analogous evolution of classically unstable Q-balls may occur in other models admitting both classically stable and classically unstable Q-balls): the first way ends up at a classically stable Q-ball with the charge different from the one of the initial unstable Q-ball, the second way results in spreading of the Q-ball into spherical waves, the third way corresponds to the spherically symmetric collapse of the Q-ball. The latter process is analogous to the nonlinear self-similar evolution observed in [18], that forces the particles to fall into the Bose star center.
And third, we will briefly examine the behaviour of classically stable Q-balls in external fields. From the classical point of view, the most interesting cases are those of the longrange interactions. Thus, the most suitable example seems to be the one that is provided by the electromagnetic interaction. In particular, classical motion of the Q-ball with logarithmic scalar field potential [19] in external electromagnetic fields was thoroughly examined in [20,21]. Meanwhile, it is clear that if a Q-ball interacts with some sort of the field, it turns out to be the source of this field. For the case of electromagnetic fields, we get U(1) gauged Qballs [22,23] instead of ordinary Q-balls (in particular, for the logarithmic scalar field potential we get the model of [24]). Although in some cases the backreaction of the U(1) gauge field on the scalar field can be neglected for U(1) gauged Q-balls [25], the existence of the gauge field can considerably modify the Q-ball properties. In particular, the question about the classical stability of U(1) gauged Q-balls has not been solved yet. Moreover, it is possible (although it is not proved) that all U(1) gauged Q-balls are unstable with respect to non-spherically symmetric perturbations [17]. In such a case, the interaction of a U(1) gauged Q-ball with the external electromagnetic field may eventually destroy the Q-ball.
From this point of view, the Wick-Cutkosky model is very attractive. On the one hand, one may expect that there exist classically stable Q-balls. On the other hand, the real scalar field in this model is massless, which ensures the existence of long-range interactions. Moreover, the back-reaction of the real scalar field can not be neglected -it is the real scalar field that provides the attraction for the complex scalar field and ensures the Q-ball existence. The latter makes the analysis even more interesting. We will show that, as expected, the Q-ball as a whole indeed obeys the Newton law, at least in the case of small speeds and accelerations. Meanwhile, the accelerated motion of the Q-ball modifies its form in accordance with its speed and acceleration. Now let us turn to the description of the Wick-Cutkosky model.

The model
The action of the Wick-Cutkosky model [10,11] has the form whereφ → 0 for | x| → ∞, m 2 > 0, h > 0. For the Q-ball, we take the standard ansatz Here r = √ x 2 and f (r) is a real function without nodes (without loss of generality, we can set f (r) > 0 for any r). The functions f (r) and g(r) satisfy the boundary conditions From the very beginning, it is convenient to pass to the new variables for which the equations of motion, following from (1) with (2), take the simple form with the boundary conditions One can see that there is no dependence on the frequency ω in these equations. The numerical solution to equations (5), (6), which was found in [9], is presented in Fig. 1.
The charge and the energy of the Q-ball take the form [9] where Due to the symmetry ω → −ω =⇒ Q → −Q, E → E, without loss of generality from here on we will consider the case Q ≥ 0, which corresponds to 0 ≤ ω < m, see It is illustrative to look at the dimensionless energy densityρ(R) of the Q-ball, which is defined as The plots of the Q-ball dimensionless energy density for different values of ω are presented in Fig. 3. We see that there exist areas with negative energy density in the Q-ball center for ω m < −G(0)−1 1−G(0) ≈ 0.565. Analogous areas with negative energy density were observed in [26] for ordinary one-field Q-balls in a theory with the scalar field potential V (χ * χ) = m 2 χ * χ − κ (χ * χ) 2 , κ > 0. 2

Classical (in)stability
Now we turn to examining the classical stability of Q-balls in the Wick-Cutkosky model. We will discuss two aspects of this problem: stability with respect to small perturbations and nonlinear evolution of the classically unstable Q-balls. Let us start with the first aspect.

Stability with respect to small perturbations
As was mentioned in Introduction, in order to be sure that Q-balls with dQ dω < 0 are classically stable, it is necessary to check that the operator L + has only one negative eigenvalue [17]. The form of the operator L + in the Wick-Cutkosky model can be easily found using the results of paper [17]: For examining this operator, it is convenient to introduce the variables In the notations X, F (R) and G(R), where R = X 2 (see (4)), the corresponding eigenvalue problem for the operator L + takes the form where There is no dependence on ω in equation (13), which simplifies the subsequent analysis considerably. The eigenfunction Ψ N ( X) can be represented in spherical coordinates as which results in with the asymptotic ξ l n (R) ∼ R l , η l n (R) ∼ R l for R → 0. Now let us define the operator L l X as leading to Let Ψ l+1 0 (R) be the eigenfunction of the operator L l+1 X , corresponding to the minimal eigenvalue Λ l+1 Since the operator L l X is Hermitian, the function Ψ l+1 0 (R) can be decomposed into eigenfunctions of the operator L l X forming a complete orthonormal set: 3 3 Of course, there exists a continuous spectrum for Λ N > 0, so it is more appropriate to use . Moreover, the function Ψ l+1 0 (R) can belong to the continuous spectrum, in such a case it should be normalized to delta-function, not to unity. However, the continuous spectrum can be easily transformed into a discrete spectrum by putting the system into a "box" of a finite size, whereas technically it is simpler to work with the discrete spectrum.

Now we can write for Ψ
Since by the definition Λ l n ≥ Λ l 0 for n ≥ 1, then Now recall that there exists an exact solution to equations (15), (16). It has the form where C is a constant. This solution corresponds to the three translational modes of the Q-ball. It is expected that this solution has the smallest eigenvalue for l = 1. Thus, it is sufficient to check that eigenfunction (21) has the lowest eigenvalue among the eigenfunctions with l = 1 and to look for the spherically symmetric (l = 0) eigenfunctions with negative eigenvalues. 4 For problem (15), (16), this was done numerically, revealing that solution (21)  Thus, since the operator L + in the Wick-Cutkosky model has only one negative eigenvalue, we can be sure that Q-balls with dQ dω < 0, i.e., Q-balls with |ω| > m √ 2 (see (8)), are classically stable.

Exponentially growing mode
In order to find explicitly the instability mode for Q-balls with dQ dω > 0, let us consider the linearized equations of motion for perturbations against the Q-ball. As was shown in [17], the instability mode can be only of the form where u( x), v( x) and ϕ( x) are real functions. The corresponding linearized equations of motion take the form [17] − ∆u + (m 2 − ω 2 + hg)u + hf ϕ = −γ 2 u + 2ωγv, It is clear that the fields u( x), v( x) and ϕ( x) can be represented in spherical coordinates and decomposed as where u sp (r), v sp (r) and ϕ sp (r) are spherically symmetric, whereas u nsp ( x), v nsp ( x) and ϕ nsp ( x) are non-spherically symmetric (i.e., include all the modes with the nonzero momenta l in spherical coordinates). The equations of motion for the spherically symmetric parts of the fields and for the non-spherically symmetric parts of the fields decouple. In particular, for u nsp ( x), v nsp ( x) and ϕ nsp ( x) we get where L v = −∆ + (m 2 − ω 2 + hg). Now let us multiply equation (28) by u nsp ϕ nsp and integrate the result over the space. Then, let us multiply equation (29) by v nsp and integrate the result over the space. Combining the resulting expressions, we can obtain Recall that there is only one negative eigenvalue of the operator L + , the corresponding eigenfunction is spherically symmetric. Thus, the first integral in (30) is nonnegative. Next, recall that L v f = 0, which is just the equation for the Q-ball profile f (r). Since the function f (r) has no nodes, it corresponds to the lowest eigenvalue of the operator L v , which is equal to zero. Thus, the second integral in (30) is also nonnegative. Consequently, relation (30) cannot be fulfilled for γ 2 > 0, which means that instability mode can be only spherically symmetric. Usually, there exists only one exponentially growing mode in the spectrum of excitations against a Q-ball. 5 In order to find this instability mode explicitly, it is convenient to pass to variables (4): where Ω = ω m ,γ = γ m . Solution to equations (31)-(33) was found numerically, the result is presented in Fig. 5.

Nonlinear evolution of the classically unstable Q-balls
Now we turn to examining the nonlinear evolution of classically unstable Q-balls in the Wick-Cutkosky model, i.e., Q-balls with |ω| < m √ 2 . We will consider only the spherically symmetric evolution. In this case, it is convenient to use the variables In these variables, the nonlinear equations of motion take the form Here −1 < Ω < 1. It is clear that it is sufficient to examine only the case 0 ≤ Ω < 1. For the stationary Q-ball solution, we have Then we slightly perturb the stationary solution for the fields as δF (R), δḞ (R), δG(R), δĠ(R) are small perturbations generated randomly in the Fourier space with the white noise spectrum. It turns out that for the classically unstable Q-balls, the only relevant part of the perturbation is the exponentially growing instability mode, i.e., only the contribution of this mode defines the subsequent evolution of the dynamical system. Thus, in such a case we can take δF Then we numerically evolve these perturbed configurations forward in time. In the simulations, we use the stable second-order iterated Crank-Nicolson scheme [27] according to its realization presented in [28]. Outgoing spherical waves that leave the interaction region are removed with the help of the Kreiss-Oliger filter [29]. Actually, here we use exactly the same numerical method as the one that was used in [17] for examining the classical stability of U(1) gauged Q-balls; the detailed description of the method can be found in Appendix C of [17].
The results of the numerical simulations show that, as expected, classically stable Q-balls (Q-balls with Ω > 1 √ 2 ) survive and drop some charge outside its core by means of spherical waves. The case of classically unstable Q-balls (Q-balls with 0 ≤ Ω < 1 √ 2 ) turns out to be more involved: the perturbations eventually destroy these Q-balls, but they can evolve in three completely different ways, depending on the form of the perturbation of the initial Q-ball and its parameters. If Re δF (0) > 0, then all such Q-balls collapse. If Re δF (0) < 0, then such Q-balls evolve in two different ways depending on the frequency of the initial Q-ball. Namely, if the frequency of the initial Q-ball is such that 0 < Ω in 0.38 (here Ω in is the frequency of the initial classically unstable Q-ball, see equations (39)-(40)), then such Q-balls simply spread into spherical waves; 6 whereas initial Q-balls with 0.38 Ω in < 1 √ 2 drop some charge and turn into the classically stable Q-balls. Below we will discuss all these ways of the unstable Q-ball evolution in more detail.
First, let us consider the case in which the contribution of the field u of the instability mode into the perturbation is positive in the center of the initial Q-ball. In such a case, all classically unstable Q-balls collapse at finite time T * . 7 While approaching the collapse time T * , the profile of the complex scalar field Y (T, R) forms a central singularity |Y | ∼ 1 R 4 , see Fig. 6. It is interesting to note that the Q-ball collapse can be described by means of the self-similar coordinate z = R √ T * −T . As can be seen from the left plot of Fig. 7, the function (T * −T ) 2 |Y (T, R)| approaches some attractor function which depends on the self-similar coordinate z only. An   analogous self-similar evolution was described in [18] for the case of collapsing Bose stars, though with a different form of singularity (in the case of Bose stars, the form of singularity is ∼ 1 R [18]). This function is not equal to zero at z = 0, which means that |Y (T, 0)| ∼ 1 (T * −T ) 2 . More precisely, |Y (T, 0)| ≈ 6 (T * −T ) 2 , which corresponds to a solution of equations (36), (37) with the terms 1 R ∂ 2 R (RY ), 1 R ∂ 2 R (RΦ) and 1 1−Ω 2 Y neglected. In this approximation, G ≡ −Y . The latter is explicitly demonstrated on the right plot of Fig. 7.
It is interesting to look at the energy density of the collapsing system. The dimensionless energy densityρ(T, R) is defined as The energy density also becomes singular as T approaches T * , the behaviour of the energy density of all collapsing Q-balls is similar regardless the value of the energy density at R = 0, see the typical example in Fig. 8. One can see that lim energy density of the Q-ball with Ω in = 0.6 at R = 0 is positive (see Fig. 3). This behaviour of the energy density is analogous to the one of the collapsing Q-ball in [26]. Meanwhile, since the total energy E (see (41)) is conserved over time (it was checked numerically that it is conserved with a good accuracy during the collapse, which is an additional test of the correctness of the numerical simulation), it turns out that for any R s > 0 almost all energy of the initial Q-ball is concentrated within the sphere of radius R s at a fixed time T s < T * . One can choose this sphere to be of the Schwarzschild radius of the Q-ball or even smaller. In this connection, it is interesting to check whether or not the collapse of the classically unstable Q-balls leads to a production of non-rotating black holes along the lines of [30][31][32]. However, in order to check this possibility, it is necessary to perform numerical simulations taking into account the existence of the gravitational field. This problem calls for further detailed investigation. Now we turn to the case in which the contribution of the field u of the instability mode into the perturbation is negative in the center of the initial Q-ball, i.e., Re δF (0) < 0. As was noted above, in such a case there can be two different ways of the Q-ball evolution. The most interesting case is the one with 0.38 Ω in < 1 √ 2 , in which classically unstable Q-balls transform into classically stable Q-balls. In Table 1, the frequencies of the initial and resulting Q-balls are presented. It turns out that the frequency of the resulting Q-ball can be roughly approximated by Ω out ≈ √ 2 − Ω in (the fit of the numerical data results in Ω out ≈ 1.37 − 0.946 Ω in ), see Fig. 9. An example of such an evolution for the Q-ball with Ω in = 0.6 is presented in Fig. 10. An important point is that the charge of the resulting Q-ball is smaller (for some Q-balls it is much smaller) than the charge of the corresponding initial Q-ball, which contradicts the naive  Table 1: The frequencies of the initial Q-balls Ω in = ω in m and the frequencies of the resulting Q-balls Ω out = ωout m .  Table 1) are presented in Fig. 11. As for the case 0 < Ω in 0.38, all such Q-balls simply spread into spherical waves. An example of such a spreading is presented in Fig. 12.  Table 1) on the E(Q) diagram.

Q-balls in external fields
Finally, let us discuss how the classically stable Q-balls behave in the external long-range field which can be produced, say, by other Q-balls in the Wick-Cutkosky model. The equations of motion for the system in the presence of an external real scalar field, which is denoted by U( x), take the formχ − ∆χ + m 2 χ + hχφ + hχU = 0, where˙= ∂ t . Of course, in the general case this system of equations can be solved only numerically. Thus, in order to get some analytical results, we suppose that U( x) ≪ max|φ( x)|, hU( x) ≪ m 2 . Let us take the following ansatz for the scalar fields: where u(t, x), v(t, x) and ρ(t, x) are real functions and is the Q-ball solution. Here R(t) denotes a position of the Q-ball center at the moment of time t. In what follows, we will consider the case |Ṙ i | ≪ 1, |R i |/ω ≪ 1. Then, the fields u(t, x), v(t, x) and ρ(t, x) can be assumed to be rather small to ensure the validity of the linear approximation for these fields. Now we substitute (44), (45) into equations of motion (42), (43) and retain the terms linear in u, v, ρ (omitting the corrections containingṘ i ,R i in the terms with u, v, ρ) and the terms ∼R i , ∼Ṙ iṘj . The terms containing more than two derivatives in time t (like ∼R iṘj ) are also omitted. We get the system of coupled equations where ∂g( x− R(t)) ∂x i and the operator L + is defined by (11). For convenience, we denote x i = x i and R i = R i .
The form of equations (46), (47) suggests that in the leading order a solution for the fields u, ρ and v contains the terms ∼R i , ∼Ṙ iṘj . Thus, the time derivatives of these fields contain the terms with more than two derivatives in time t. Such terms were omitted in the derivation of equations (46), (47), thus, in the leading order we can also omit the termsu,v,ü,v,ρ. It is clear that the radiation which is produced by the accelerated motion of the Q-ball, turns out to be neglected too. The equation for the field v decouples in this case and we can simply set v ≡ 0. Thus, we arrive at the system of equations

Conclusion
In the present paper, we discussed Q-balls in the Wick-Cutkosky model. Because of its simplicity, the model turns out to be a useful toy model for examining various properties of Q-balls, which can be inherent to Q-balls in other theories. In particular, we performed a detailed study of the Q-ball stability, including the stability with respect to small perturbations and spherically symmetric nonlinear evolution of the classically unstable Q-balls. The results of the analysis demonstrate that, depending on the characteristics of the Q-ball and on the form of the initial perturbation, the nonlinear evolution of the classically unstable Q-balls can lead to completely different outcomes. Namely, such Q-balls can spread into spherical waves, collapse or evolve into a classically stable Q-ball. Analogous evolution of classically unstable Q-balls may occur in other models admitting classically unstable Q-ball solutions.
In addition, we examined the behaviour of classically stable Q-balls in external fields. The key feature of the analysis is that this long-range external field is of the same nature as the one that forms the Q-ball itself. It is shown that, at least in the non-relativistic limit, the Q-ball as a whole obeys the Newton law. Meanwhile, the accelerated motion of the Q-ball modifies its form in accordance with its speed and acceleration.
We hope that the results presented in this paper can be useful for the further study of the Q-ball behaviour in different models.
Substituting the latter result into (65), finally we get