More about Q-ball with elliptical orbit

Q-balls formed from the Affleck-Dine field have rich cosmological implications and have been extensively studied from both theoretical and simulational approaches. From the theoretical point of view, the exact solution of the Q-ball was obtained and it shows a circular orbit in the complex plane of the field value. In practice, however, it is reported that the Q-ball that appears after the Affleck-Dine mechanism has an $elliptical$ orbit, which carries larger energy per unit $U(1)$ charge than the well-known solution with a circular orbit. We call them"elliptical"Q-balls. In this paper, we report the first detailed investigation of the properties of the elliptical Q-balls by $3$D lattice simulation. The simulation results indicate that the elliptical Q-ball has an almost spherical spatial profile with no nodes, and we observed a highly elliptic orbit that cannot be described through small perturbations around the ground state Q-ball. Higher ellipticity leads to more excitation of the energy, whose relation is also derived as a dispersion relation. Finally, we derive two types of approximate solutions by extending the Gaussian approximation and considering the time-averaged equation of motion and we also show the consistency with the simulation results.


I. INTRODUCTION
Q-ball is a non-topological soliton in a complex scalar field theory with a global U (1) symmetry [1]. The solution corresponds to the spherically symmetric field configuration minimizing the energy of the system with the fixed U (1) charge. Such a solution exists for the potential V (φ) of the complex scalar field φ if V (φ)/|φ| 2 has the minimum at φ = 0. One of the characteristics of the solution is that the configuration of φ has time dependence and the orbit of φ in the complex plane is a true circle.
The formation of the Q-ball has been studied analytically and numerically in some extensions of the standard model. For example, in the supersymmetric extension of the standard model, some scalar condensates with the baryon number form Q-balls after Affleck-Dine baryogenesis [2][3][4][5][6][7][8][9][10][11][12][13][14]. There, the initial configuration of the complex scalar field is set to be the coherently oscillating state in the complex plane. Due to the instability of such a state, the initial fluctuation of the scalar field grows and eventually develops into the Q-balls.
One non-trivial result from the numerical analysis is the discovery of excited states of Q-balls [14][15][16]. 1 The excited state has larger energy compared to the energyminimizing solution with the same charge. The configuration of the complex scalar in the excited state draws an elliptical orbit, while the orbit of a genuine Q-balls is a true circle. We call those excited Q-balls as the "elliptical" Q-balls and genuine ones as the ground state Q-balls. It is numerically shown that the elliptical Qballs are semi-stable at least [14] and they are considered as the "transients", which appear during the formation of the circular Q-balls.
Even though the formation of the elliptical Q-balls was pointed out, their properties were not pursued sufficiently. For instance, in the numerical simulation, there is no detailed investigation on how the energy excitation depends on the ellipticity. The theoretical explanation of the elliptical Q-ball solution is also not present.
It should be noted that the elliptical Q-ball is different from the excited state solution that is discussed in e.g. Refs. [17,18]. The radial excitation with nodes [17,18] and the (non-spherical) angular excitation [17] of the Qball solution are studied as possible forms of the excitation. However, the Q-balls produced in the lattice simulation seem to have no nodes and are almost spherically symmetric (See Sec. III for more details.). Furthermore, perturbations around the ground state Q-ball [1,19,20] do not explain the dispersion relation of the elliptical Qball.
In this paper, we investigate the detailed properties of the elliptical Q-ball for a scalar potential motivated by the AD mechanism in gravity mediation where the properties of the ground Q-ball size and energy are given by the simple formulas. We perform 3D lattice simulations and present the spatial profiles and the dispersion relation of the elliptical Q-ball. Then, we will derive some approximate solutions of the elliptical Q-ball, which we will show are in rough agreement with the lattice simulation results. We also discuss the stability of the elliptical Q-ball by considering perturbations to the approximate solution.
The paper is organized as follows. In Sec. II, we review the ground state Q-ball solution and its properties in gravity mediation model. In Sec. III, we describe the setup of our numerical simulation. In Sec. IV, we show the numerical results on the formation of the elliptical Qballs and their detailed properties. In Sec. V, we derive the approximate elliptical Q-ball solutions and discuss the consistency with the simulations. The final section is devoted to our conclusion.

II. GROUND STATE Q-BALL
In this section, let us review the ground state Q-ball solution. In particular, we focus on the Q-ball in the gravity mediation model.
Let us consider the theory of a complex scalar field φ with the Lagrangian density, Here, V (φ) is a potential with the U (1) symmetry, which is invariant under the phase rotation of φ. The Noether charge of this U (1) symmetry is given by The energy is represented as To figure out the ground state Q-ball solution, let us consider the energy-minimizing condition for a fixed charge. Using the method of Lagrange multipliers, one must find an extremum of where ω is a Lagrange multiplier. From the first term of the integrand in the second line, the time-dependent part of φ is minimized for the choice Note that this formula indicates that the orbit in the complex plane is a true circle. Then, the Eq. (4) is reduced to As the extremum condition, we then obtain the equation of motion for Φ(x), where r is the radial coordinate and we assume that Φ is spherically symmetric in the position space. From the finiteness of the energy and the continuity of the solution, we impose the boundary condition, The condition [1] for the existence of the solution is given by This condition claims that V (Φ) must be "shallower" than the quadratic. In the following, the potential is specified to the gravity mediation type, where M * is the renormalization scale defining the mass parameter, m, of the complex scalar field. The constant K is for example determined as K = −0.01 ∼ −0.1 by the spectrum of the minimal supersymmetric standard model. Throughout this paper, we take, We have confirmed that the qualitative results in the following are the same for K = −0.01. In gravity mediation scenarios, the equation of motion is numerically solved and it is shown that the solution is well approximated by the Gaussian spatial profile [9], where R is the constant which corresponds to the radius of the Q-ball. By substituting this into the Eq. (8), one obtains the following relation, Then, the charge and the energy is obtained as We stress that one important feature of the ground state Q-ball in gravity mediation is the dispersion relation, E mQ. As we will see later, the Q-balls formed in the realistic situation have elliptical orbits and no longer follow this relation.
Hereafter, we redefine the parameters as

III. SIMULATION SETUP
In this section, we perform simulations in the gravity mediation model and show the formation and properties of the elliptical Q-balls. While there is no rigid theory as in the previous section on the elliptical Q-ball solution, we will give the approximate theory in the next section, which turns out to roughly agree with the simulation results.
In our analysis, we use the package PyCOOL [21], which we modified for our analysis. The PyCOOL is an object-oriented Python program that uses GPU(s) for simulating the scalar field dynamics in the early universe and performs the 3D simulations.
We assumed the matter-domination after the inflation 3 and set the initial scale factor unity: a(t = t i ) = 1. As the initial condition for the homogeneous AD field, we take where ini represents the initial ellipticity of the orbit in the complex plane. We also give the initial Gaussian fluctuation, though the eventual phenomena are known to be irrelevant to this initial condition on the fluctuation, except for the formation time of the objects. In the next section, we will show the result of our simulation where the elliptical Q-balls are formed.

IV. NUMERICAL RESULTS
As a result of our simulation, we confirmed that the non-linear lumps are formed at t form ∼ 8 × 10 2 , which are stable, that is, do not decay into the scalar waves, until the final time. We present the result in Fig. 1, where we set ini = 0.2. From the figure, we can see that the non-linear spherical objects are formed. We will identify these objects as the excited Q-balls by fitting the approximate semi-analytic profiles as we will discuss later. First of all, let us confirm that the produced Qball has an elliptical orbit. In Fig. 2, we present the orbit at the center of the lump, which is plotted around the final time t f = 5000. We can see that it is indeed an ellipse that has a rotating axis (left figure). We also find that the orbit itself slightly rotates around the origin when we observe the orbit for a longer time (right figure), whose behavior will be explained in the next section. The normalization of the figures is the same as Eq. (18).
In Fig. 3, we also present the charge density profile of the lump with the maximal charge. The horizontal axes denote the spatial distance from the Q-ball center. The vertical axis is the charge density. The red, green, and blue curves show the charge density profile for the x, y, and z coordinates. It should be noted that the charge density profile of the elliptical Q-ball is almost spherically symmetric. Such spherical profiles of the Q-balls are also seen in Fig. 1. We also note that there are no nodes for the profile. Those properties suggest that the elliptical Q-ball is different from the excited states discussed in Refs. [17,18]. 4 Next, let us see the dispersion relation of the elliptical Q-ball. In Fig. 4, we plot the relation between E/Q and of the Q-balls with the charge of 10 Q 100, where is defined by the short axis of the ellipse divided by the long axis of the ellipse. The difference in the color of the dots denotes that in the initial condition. We vary the initial ellipticity ini from 0.2 to 1.
We confirm that E/Q of the elliptical Q-balls are larger than those of the ground state Q-balls. In other words, it is shown that the energy of the elliptical Q-balls is larger than that of the ground state Q-balls for a fixed charge. This time we have run the simulations with N grid = 64 3 in order to save the time of computation. We confirmed for several samples that there is no essential difference from the case of N grid = 128 3 in the plot.

V. APPROXIMATE SOLUTION FOR ELLIPTICAL Q-BALL
In this section, we derive the approximate solution for the elliptical Q-ball. We consider two different approaches to obtain such a solution. Those approximations lead to the dispersion relation of the elliptical Qballs which roughly fits the simulation results.
elliptical Q-ball is not the same as the excited state discussed in Refs. [17,18], while the slight deviation from the spherical symmetry might be explained through the perturbative angular excitation of the Q-balls.

A. Extended Gaussian Approximation
The ground state Q-ball solution in gravity mediation is well fitted by the Gaussian approximation. The elliptical Q-ball solution is also approximated by considering the extension of the Gaussian profile as explained below.
Let us take the Gaussian ansatz, where ϕ(t) is a complex function of t. We decompose ϕ(t) into the radial direction Φ(t) and the phase direction θ(t), Here, we do not specify the time-dependence of Φ(t) and θ(t). We note that the time dependence of ϕ(t) in the ordinary Gaussian ansatz for the ground state Q-ball solution is just with constant ω. From the ansatz in Eq. (22), the energy of the system is rewritten as Performing the integration for the spatial coordinate, we obtain the energy formula, From the U (1) charge conservation, the following relation is also obtained, Using the above relation in Eq. (26), the energy of the system is represented only by Φ(t), where U (Φ) is an effective potential, and we defineẼ = 2 √ 2E/π 3/2 R 3 andQ = 2 √ 2Q/π 3/2 R 3 . This system is analogous to that of Keplerian motion of the planet in the gravitational potential. The lowest energy state described in this system is apparently (Fig. 5), where Φ 0 is a minimum of the effective potential and we used R −2 ∼ |K| 1. In this case,θ = Q/Φ 2 0 = constant and ϕ describes the circular orbit in the complex plane. This solution is nothing but the ground-state Qball solution. Then, the energy per unit charge is given by which is consistent with the result of a conventional formulation of the Q-ball.
To obtain the excited state solution, let us consider the caseΦ = 0 in the system. Because Φ obtains the kinetic energy, Φ oscillates around the potential minimum Φ 0 (Fig. 6).
In such a condition, the radial direction of the complex field φ oscillates between Φ min and Φ max , which means that the orbit of φ is elliptical. Due to the non-vanishing kinetic energy of Φ, the energy of this Q-ball is larger than circular (ground state) Q-ball: Since Φ min and Φ max are related to the ellipticity parameter as = Φ min /Φ max , the ratio E/Q is calculated as a function of . Neglecting the contribution of the order |K|, we obtaiñ Therefore, While we have neglected the K dependence in the above formulas for simplicity, for the actual computation using our approximate analysis we are solving the field equation including the K dependence. However, it will turn out that the K dependence in the dispersion relation, for instance, is negligible for K 0.1, so that the above formula is actually allowed in practice. We will discuss on this point at the end of this subsection.
The dispersion relation in Eq. (33) is an extension from that of circular Q-ball. While elliptical Q-ball solution coincides with the circular Q-ball solution for = 1, it carries larger energy per unit charge than the circular one for = 1. Therefore, in this sense, we can regard the elliptical Q-ball as a kind of "excited state" of the Q-ball.
Furthermore, we also show the orbit of φ in the complex plane in Fig. 7. There, we take Q = 50, R = 2/|K|. We can see that the axis of the ellipse is rotating in time. This behavior seems to agree with the simulation results in Fig. 2.

Comparing with Simulation Results
To confirm the availability of the extended Gaussian approximation, let us compare the approximation with the simulation results. In Fig. 8, we show the fitted lines for the orbit and the charge density. The gray dashed lines are the results by the extended Gaussian approximation with Q = 36.5, R = 2/|K|. 5 The other color lines are the same as the orbit and the charge density in Fig. 2 and Fig. 3. We confirm the rough agreement between the approximation and the simulation results. More quantitatively, we can obtain the parameters Q, Φ max , of the elliptical Q-ball from the simulations. On the other hand, we solve the field equation by inputting the parameters Q, R, Φ max and obtain in our approximate analysis. Then, we can compare from the simulation with from the approximate analysis, in order to examine the validity of the latter. We present the results for several samples in Tab. I. The table shows that our approximate theory is valid for these samples with about 10% precision level. We especially note that the slight deviation from the spherical symmetry does not seem to be that important in the description of the properties of the elliptical Q-ball, since the approximation R = 1/ 2|K| is roughly valid. For the better approximation, we may also need to take this deviation in the profile into account. We will leave the detailed analysis to the future work.   The data from the simulation is the same as that in Fig. 4. Here, we take R = 2/|K|. for given is also calculated in Fig. 9. The gray dashed line is the result for Q = 50 (and R = 2/|K|), which fits well with the simulation results. Note that we are solving the field equation including the K dependence to draw the gray dashed line, as we mentioned earlier.
However, the K dependence is almost negligible for K 0.1, so that we can use the simple formula Eq. (33) in practice: See the K dependence of the dispersion relation in Fig. 10. 6

B. Complete Ellipse Approximation
As another way, let us consider the time-averaged equation of motion like the I-ball/oscillon solution. As shown in Fig. 7, the orbit of the elliptical Q-ball is almost closed during one oscillation. Focusing on such adiabatic motion of the elliptical Q-ball, we take the complete elliptical orbit ansatz 7 , where is the ellipticity parameter in the defined range of 0 < ≤ 1. For = 1, the ansatz is nothing but the one for the ground state Q-ball in Eq. (6). Under this ansatz, the total charge Q is given by Apparently, the charge is time-invariant under the above ansatz.
The time-averaged equation of motion for Φ [22][23][24] is whereV (Φ) is the time-averaged potential, Here, 7 The elliptical orbit is not completely closed and slightly rotating as discussed in the previous section.
FIG . 11. Fitting in the complete ellipse approximation. The gray dashed line is the orbit and the charge density profile in the left and the right figures, respectively. The solid lines are the same as the one in Fig. 2 (left) and Fig. 3. Here, we take = 0.2, ω = 1.149.
FIG. 12. The dispersion relations obtained from the lattice simulation and the complete ellipse approximation. The gray dashed line is E/Q for Q = 50. The data from the simulation is the same as the ones in Fig. 4.
Using the above formulas, we obtain the dispersion relation, from the time-averaged energy, This result is similar to the results from the extended Gaussian ansatz in Eq. (33).

Comparing with Simulation Results
Let us first determine the variables ω and from the orbit in Fig. 2. From the figure, the ellipticity parameter is given as 0.2. The parameter ω is determined to fit the value of the field. 8 In Fig. 11 (the left figure), we present the fitted line with = 0.2, ω = 1.149 as the gray dashed line. The blue line is the same as the line in Fig. 2.
From the above solution, the charge density profile profile is calculated in Fig. 11. The gray dashed line is the result with the same parameters, = 0.2, ω = 1.149. The blue line is the same as the line in Fig. 3. Here we can also see the rough agreement between this approximate solution and the simulation result.  We can perform more quantitative examination of the agreement as in the previous subsection. First, we can obtain the parameters , Φ max , and Q from the simulation results, and we input the parameters , Φ max 9 and derive Q using the complete ellipse approximation, which we compare with the charge Q from the simulation. The results for several samples are shown in Tab. II. The table shows that the estimation of Q using the complete ellipse method is valid for these samples at least with about 10% precision level.
The dispersion relation E/Q for is compared with the simulation result in Fig. 12. There, the gray dashed line is the result for Q = 50, which is consistent with the simulations as well. Therefore, we find that this approximation is also useful for the discussion on the properties of the elliptical Q-ball.

C. Stability of Elliptical Q-ball
Before closing this section, let us briefly discuss the stability of the elliptical Q-balls. We know that the circular Q-ball = 1 is absolutely stable because it is the most stable state which the complex scalar can take. In the case of the elliptical Q-balls, however, they might move to the more stable state, that is, less elliptical Qballs. In Ref. [14], the fission of the elliptical Q-ball to the multiple Q-balls is observed, which are less elliptical than the original one. Therefore, the elliptical Q-balls are considered as the "transients", which appear during the formation of the circular Q-balls. The mechanism of how the elliptical Q-ball makes such transition is an open question. Here, we discuss that the growing perturbation around the elliptical Q-ball may be the seed for the transition to more stable Q-balls.
Since we now know that the simple approximations roughly describe the elliptical Q-ball, here we use the extended Gaussian approximation, and consider a little perturbation δΦ around the elliptical Q-ball solution: where Φ is the background elliptical Q-ball solution which satisfies the equation of motion. Then, the equation of motion for the perturbation δΦ is whose solution will give how the perturbation grows. In Fig. 13, we present the growth of the perturbation for several values of elliptical parameter . Here, we take Q = 50, R = 2/|K| for demonstration. The green, orange, blue lines denote the fluctuation |δΦ| for = 0.2, 0.5, 0.8 after 300 oscillation of Φ, respectively. We also show in Fig. 14 the relation between the elliptic-ity and the time of δΦ = 1 10 for some k from the approximation (τ ). The figures imply that δΦ has instability and grows exponentially in time, which is nothing but the broad resonance which is originated by the higher-order self-coupling of the complex scalar, and the perturbation growth is efficient for highly-elliptical case 1. For a certain time scale, the fluctuation δΦ dominates over the background mode Φ and the original Q-ball solution is no longer maintained.
This time scale is expected to be related to the magnitude of the background oscillation: If the amplitude of the oscillation of Φ is larger, the instability becomes stronger and the above time scale becomes shorter. This means that highly-elliptical Q-ball has a short lifetime. This feature leads us to expect that the perturbation growth may be the seed of the fission of the elliptical Q-ball, since it is known that the fission is efficient for highly-elliptical case 1 [14]. 11

VI. CONCLUSIONS AND DISCUSSION
In this paper, we have studied the detailed properties of the elliptical Q-ball by the simulation and the semianalysis. The elliptical Q-ball is the localized clump state with finite charge, which has the elliptical orbit in the complex field plane in Fig. 2. 12 This is different from the ground state Q-balls with the circular orbit. Those elliptical Q-balls correspond to the excited Q-ball which has larger energy per unit charge compared with the ground state Q-ball of the same charge. Although such an excited state of the Q-ball is discussed in many works of literature, the elliptical Q-ball formed in the numerical simulation is not described in detail. To investigate the detailed properties, we have performed the numerical lattice simulation in 3 + 1 space-time dimension. We have first studied the spatial profiles of the elliptical Q-balls and showed that they have no nodes and some non-spherical Q-balls are formed. Furthermore, the (highly) elliptical orbit in the complex field plane shows that the elliptical Q-ball cannot be described by small perturbation against the ground state Q-ball. We have also calculated the dispersion relation of the elliptical Q-balls, which depends on the ellipticity of the Q-balls as E/Q (1 + 2 )/2 . We note that this relation is not well understood through the previous studies. We have derived the approximate solution for the elliptical Q-balls that explains the numerical results, including the ellipticity dependence and 10 The value itself is not important here because we are interested in relative behavior depending on . 11 Direct comparison of the transition time given by our approximation with that from the simulation results seems to be difficult because we need to extract information on the value of perturbation when the elliptical Q-ball forms. Therefore, we only compare the qualitative properties in this paper. 12 More precisely, the elliptical orbit is not closed, with its axes rotating.
the rotating property of the elliptical orbit, which leads to our conclusion that the elliptical Q-ball is an almost spherical (nearly Gaussian) object that has the exciting energy depending on the ellipticity as E/Q (1+ 2 )/2 , supported by the simulation and our approximation approaches. We have also discussed the instability of the elliptical Q-ball against the small perturbation.
Several comments are in order. Throughout the paper, we have studied the elliptical Q-ball in gravity mediation. The gravity mediation model is useful since it is the simplest model in the sense that the properties including the Q-ball size and energy are given by the simple formulas. However, the existence and the properties of the elliptical Q-balls in the other models are non-trivial matter. For instance, in the gauge mediation models [5], it has not been reported by the simulations whether the excited states are formed or not as transients. Furthermore, the ansatz for the Q-ball solution and its dispersion relation in the gauge mediation models are different from those in the gravity mediation models and they should be examined if the similar discussion as in this paper is valid.
There is also another type of Q-ball called New type Q-ball, which is realized when the both gravity and gauge mediation effects exist and gravity mediation dominates the potential. The formulas for the Q-ball properties are similar to the case when there exist only the gravity mediation effects, but the scales are different, e.g. soft mass and gravitino mass, which leads to the difference in the cosmological scenario (For instance, while the gravity mediation type Q-ball is unstable against the decay into the baryons, New-type Q-ball is stable and considered as a dark matter candidate).
In this work, we assumed the matter-dominated era during the Q-ball formation. Several scalar field baryogenesis mechanism [25][26][27][28], however, occurs during the radiation dominated era. We expect that the elliptical Q-balls will be formed even during radiation dominated era (for generic potentials), since the reason that the elliptical Q-balls are formed as the excited transients is the energy gap between the initial coherent condensate and the final ground state Q-balls, which generally exists, especially independent of the way of cosmic expansion.
Let us also comment on the applications of our approximate theory to other non-topological solitons, e.g. oscillon/I-ball [22][23][24][29][30][31][32][33][34] and boson star [35,36]. For instance, Refs. [37,38], suggest that oscillon/I-ball and boson star can be understood as Q-ball solutions associated with the particle number U (1) charge, which might make it possible to derive some extended solutions of such non-topological solitons, using our extended theory of the Q-balls. These will also be explored in our forthcoming paper.