Gravitating compact $Q$-ball and $Q$-shell solutions in the $\mathbb{C}P^N$ nonlinear sigma model

We study compact gravitating $Q$-ball, $Q$-shell solutions in a sigma model with the target space $\mathbb{C}P^N$. Models with odd integer $N$ and suitable potential can be parameterized by $N$-th complex scalar fields and they support compact solutions. A coupling with gravity allows for harboring of the Schwarzschild black holes for the $Q$-shell solutions. The energy of the solutions behaves as $E\sim |Q|^{5/6}$, where $Q$ stands for the $U(1)$ Noether charge, for both the gravitating and the black hole solutions.Notable difference from the solutions of the flat space is that upper bound of $|Q|$ appears when the coupling with gravity is stronger. The maximal value of $|Q|$ quickly reduces for larger coupling constant. It may give us a useful hint of how a star forms its shape with a certain finite number of particles.


I. INTRODUCTION
Q-balls are stationary soliton solutions of certain complex scalar field theories with self-interactions [1][2][3][4]. The U(1) invariance of the scalar field leads to the conserved charge Q which is identified with the electric charge of the constituents for theories coupled to electromagnetic field. Q-balls have attracted much attention in the studies of evolution of the early Universe [5,6]. In supersymmetric extensions of the standard model, Q-balls appear as the scalar superpartners of baryons or leptons and they form coherent states with baryon or lepton number. They may survive as a major ingredient of dark matter [7][8][9].
In order to get the standard Q-balls one usually considers solutions with constant angular velocity ω in the internal space of fields and with spherical symmetry in position space. An absolute value of the Noether charge |Q| has no upper bound while the angular velocity always has some limitations. The solutions exist when the charge Q is greater than the lower critical value QC. When concerning the absolute stability then another critical point emerges. Namely, a Q-ball is stable against decay into smaller Q-balls (single consituents) if its energy, above the critical value QS, is less then mQ, where m is the mass of a single constituent.
Some are other types of Q-balls associated with local symmetries have been extensively studied in [10][11][12]. Such Q-balls are unstable if they possess sufficiently large charge. The instability is caused by repulsion that originates in the gauge force. It means that for stable solutions the charge |Q| has upper bound.
Compactons are field configurations that exist on finite-size supports. Outside this support the field is identically zero. In the case of spherical configurations the support is delimited by outer radius rout. For instance, the signum-Gordon model i.e. the scalar field model with standard kinetic terms and V −shaped potential gives rise to such solutions [13,14]. Interestingly, when the scalar field is coupled with electromagnetism then the inner radius emerges, i.e., the scalar field vanishes also in the central region r < rin. Thus, the matter field exists in the region rin ≦ r ≦ rout. Such configurations * pawel.klimas@ufsc.br † sawadoph@rs.tus.ac.jp ‡ phyana0513@gmail.com of fields are called Q-shells. Such shell solutions have no restrictions on upper bound for |Q|. The authors claim that the energy of compact Q−balls scales as ∼ Q 5/6 and of Q-shells for large Q as ∼ Q 7/6 . It clearly indicates that the Q-balls are stable against the decay while the Q-shells may be unstable. Many results concerning compact boson stars and shells are presented in [15][16][17][18]. For the boson shell configurations, one possibility is the case that the gravitating boson shells surround a flat Minkowski-like interior region r < rin while the exterior region r > rout is the exterior of an Reissner-Nordström solution. Another and even more interesting possibility is the existence of the charged black hole in the interior region. The gravitating boson shells can harbor a black hole. Since the black hole is surrounded by a shell of scalar fields, such fields outside of the event horizon may be interpreted as a scalar hair. Such possibility has been considered as contradiction of the no-hair conjecture. The higher dimensional generalizations have been considered in [19,20].
Recently, a new model supporting both Q-balls and Qshells has been proposed by one of us [21]. The model is defined in 3+1 dimensions and it has the Lagrangian density where the coupling constants M , µ have dimensions of (length) −1 and (length) −2 , respectively. The potential V is chosen in the way that the model supports compact solutions. For models with standard kinetic terms (it is the case of model (1)) it is sharp at its minima and so it is called V-shaped potential. The model (1) is a direct extension of the CP 1 nonlinear sigma model on a model with target space CP N . The field X is called the principal variable and it successfully parameterizes the coset space SU (N + 1)/U (N ) ∼ CP N . There exist compact solutions for all odd numbers N ,i.e., N := 2n + 1, n = 0, 1, 2, · · · . The salient feature of the model from other models containing only Q-balls, is that it supports both Q-ball and Q-shell solutions and the existence of Q-shells does not require the electromagnetic field. For n = 0, 1 the solutions form Q-ball while for n ≧ 2 they always are Q-shell like. Again, there is no upper bound for the charge |Q| and also there is no limitation from above for ω. In the complex signum-Gordon model with local symmetry Q-balls exist due to presence of the gauge field whereas in the case of the CP N model they appear as the result of selfinteractions between scalar fields.
In this paper, we consider model containing complex scalar fields coupled to gravity and obtain the compact Q-ball and Q-shell solutions. The resulting self-gravitating regular solutions can be identified with boson stars [22]. We also study the systems composed of compact Q-shell solutions and the Schwarzschild-like black holes located in the center of the shell. The space-time inside r < rin and outside the shell r > rout is the Schwarzschild space-time. The systems with these property might be considered as possible examples of violation of the no-hair conjecture. A notable difference between these solutions and solutions in flat space-time is that the first ones possess upper bound for |Q| and the value of this bound drastically decreases with increasing of gravitational coupling constant. The maximal charge |Q| is attained for large n for sufficiently small gravitational constant which explain, to some extent, why the boson stars have definite size and mass. In this paper we only present results for the Schwarzschild space-time, however, the analysis can be extended on the case of electrically charged solutions with Reissner-Nordström metric.
The paper is organized as follows. In Section II we describe the model, coupled to gravity and discuss its parametrization. In Section III we present numerical results. Section IV is devoted to detailed analysis of stability of the Q-ball, Q-shell solutions. Conclusions and remarks are given in the last Section.

A. The action, the equations of motion
We consider the action of a self-gravitating complex variable X coupled to Einstein gravity where G is Newton's gravitational constant. The CP N space has a nice parametrization in terms of the principal variable X [23] (see [21,24] for explicit construction of the solutions), defined as It parametrizes the coset space CP N = SU (N + 1)/SU (N ) ⊗ U (1) with the subgroup SU (N ) ⊗ U (1) being invariant under the involutive automorphism (σ 2 = 1). It follows that X(gk) = X(g) for σ(k) = k, k ∈ SU (N ) ⊗ U (1). The CP N model possesses some symmetries. It is easy to see that the kinetic term of the Lagrangian has the symmetry X → AXB † , A, B ∈ SU (N + 1). Thus the potential is chosen so as to break the symmetry down to the diagonal one, i.e., X → AXA † . The existence of compact solutions require special class of potentials. An example of such a potential which we shall adopt in this paper has the form As we shall see later, the behavior of fields at the outer border of compacton (and also at spatial infinity) implies X → I.
We assume the (N +1)-dimensional representation in which the SU (N + 1) valued group element g is parametrized by the set of complex fields ui: which fixes the symmetry to SU (N ) ⊗ U (1). The principal variable takes the form where ϑ := √ 1 + u † · u. Thus the CP N Lagrangian of the model (2) takes the form where The variation of the action with respect to the metric leads to Einstein's equations and where the stress-energy tensor is of the form Next, varying the action with respect to fields ui and u * i one obtains the CP N field equations The obtained system of coupled equations (9) and (11) is quite complex so we shall adopt numerical techniques to solve it. It is convenient to introduce the the dimensionless variables For the solutions with spherical symmetry, we employ the standard Schwarzschild-like coordinates such that the line element reads We also restrict N to be odd, i.e., N := 2n + 1. The ansatz um(t, r, θ, ϕ) = 4π 2n + 1 f (r)Ynm(θ, ϕ)e iωt (15) allows for reduction of partial equations to the system of radial ordinary equations, where Ynm, −n ≤ m ≤ n are the standard spherical harmonics and f (r) is the profile function. Each 2n + 1 field u = (um) = (u−n, u−n+1, · · · , un−1, un) is associated with one of 2n + 1 spherical harmonics for given n. The relation n m=−n Y * nm (θ, ϕ)Ynm(θ, ϕ) = 2n+1 4π is very useful for obtaining an explicit form of many inner products. Substituting (13) and (15) into the Einstein field equations (9), we get their components Then from (16) and (17) we obtain equations for the metric tensor functions A and C where α is a dimensionless coupling constant defined as We shall solve numerically such obtained equations in dependence on the coupling constant α and the frequency ω. Variation of the parameter α is equivalent to variation of the model parameter µ (or of M in fully dimensional case). In this context, the parameter µ is a kind "susceptibility" of gravity to matter. The CP N equations (11) read Here the prime ′ means differentiation with respect to r. The dimensionless Lagrangiañ allows us to obtain the dimensionless Hamiltonian of the CP N model, which reads The symmetry SU (N + 1) is reduced to SU (N ) ⊗ U (1) for parametrization (5). It is the symmetry of the Lagrangian (22). It contains subgroup U (1) N given by the transformation where βi are some global parameters. The Noether charge associated with this global symmetry plays important role for stability of solutions. The Noether currents corresponding with the symmetry transformation (24) read Using the ansatz (15) we find following form of the Noether currents and J The conservation of currents is almost straightforward because non-vanishing components (26) and (27) depend only on r and θ. Thus The integral of (28) When J (m) ϕ decreases sufficiently quickly at the spatial boundary (and so it does not contributes to the integral) the Noether charge is conserved. The spatial components of the Noether currents do not contribute to the charges, however, they can be used to introduce some auxiliary integrals  With the help of Q m and q m the total energy E can be expressed in the form

B. The behavior of solutions at the boundary
In this section we shall discuss behavior of solutions at the boundary, which means that we mainly look at the origin r = 0 and the border(s) of the compacton. First we consider expansion at the origin and so the solution is represented by series After substituting these expressions into equations (19), (20), where f0, A0 are free parameters. For n = 1 we obtain with free parameters f1 and A0.
For n ≧ 2 we have no nontrivial solutions at the vicinity of the origin r = 0. It means that the solution is identically zero inside a ball with certain radius. Thus we shall study expansion at the sphere with certain finite radius Rin. Such a solution is known as the Q-shell solution because its radial The functions f (r), A(r) and m(r) are represented by series  After determining few first coefficients of expansion we get In the region r > Rout the metric functions reduce to the Schwarzschild solutions because there is no matter function in this region. Thus we have It is well-known that certain global black hole observables can be directly deduced from asymptotic expansion of the solutions [25]. For instance, a mass of black hole is defined in terms of M0.

III. THE NUMERICAL RESULTS
In this section we present numerical solutions describing gravitating compact Q-balls and Q-shells. We use shooting algorithm to solve differential equations (19), (20), (21) with different n and look at dependence of these solutions on the parameters ω and α . According to (34) and (35) the solutions with n = 0, 1 are regular at the origin. Thus, they are Q-ball type solutions and they can be designated as boson stars. In Fig.1 we plot the solution for the CP 1 case. The matter profile function f (r) has non-zero value at the origin and reaches zero at the compaction radius. We also plot the metric functions A(r), C(r). With increasing of the coupling constant α the solutions tend to shrink. The bottom right figure in Fig.1 shows the shooting parameter f0 versus the coupling constant α. We observe the existence of branches of the solutions, i.e., there are two solutions for the unique shooting parameters. The solution of the second branch (the dotted line) is stronger localized at the origin.
As discussed in the last section, the solutions can not be nontrivial at the origin for n ≧ 2. It leads to existence of Q-shells with the matter field localized in the radial segment r ∈ (Rin, Rout). First we consider the case of self-gravitating solutions, i.e., assume that there is no extra mass at the center i.e., M0 = 0. It means that the space-time in the interior of the Q-shell must be regular. In Fig.2 we present examples of typical results for CP 5 and CP 101 . In both cases the metric functions are regular at the origin. The space-time inside empty region of the shell, r ≦ Rin, is Minkowski-like with constant metric functions A(r) = A0 (the dashed lines) and C(r) = 1.0. Our results show that the shell solutions corresponding with higher n become larger and thinner and that the difference between self-gravitating Q-shells and Q-shells in flat space-time is greater for higher n. We also note that the size of compactons varies significantly with, even small, changes of parameter α.
Next we have considered the case of Q-shells with a massive body immersed in their center. In particular, an interesting possibility is having this body as a Schwarzschild-like black hole [15]. We shall report here on such a case. We choose the mass M0 such that the event horizon is localized in a central (empty) part of the shell and we solve the equations in the First, we assume that in the region between the event horizon and the inner boundary of the shell r ∈ [r h , Rin] hold Einstein's vacuum equations as the consequence of vanishing of matter fields in this region. A similar approach is used in the region outside the shell r ∈ [Rout, ∞), where all matter fields vanish. Next, we solve the equations in the region r ∈ (Rin, Rout) and then smoothly connect the metric functions with both vacuum solutions at inner and outer boundary of the shell. Note that the values of both radii of the shell are not a priori known because they are determined by shooting procedure. In Fig.3 we present typical results for harbor solutions. We also present figures which show behavior of the ratio of the radii r h /Rin of the event horizon r h and the inner boundary of the shell Rin. The ratio r h /Rin = 0 corresponds with the self-gravitating, i.e. M0 = 0, solutions. We have also looked at dependence of the solution on the radius of the event horizon (mass of the black hole). The radii of the shell solution Rin and Rout smoothly vary with increasing of r h . As r h approaches Rin (the ratio tends to 1) the shell becomes thinner. Also the change of the coupling constant α (see middle figures in Fig.3) affects the form of the solution. We observe that the larger coupling constant is the thinner the shell is. Moreover, we also see that shells again become thinner with increasing the rotational velocity ω. A notable difference is visible between the CP 5 (the left) and the CP 11 (the right). For the CP 5 case the event horizon can be localized very close to the inner border of the shell, however, without touching it exactly. On the other hand, for CP 11 the solution collapses before the event horizon goes very close to the inner boundary. Also, for larger n the distance between these two radii grows significantly. Therefore we conclude that the black hole can not possesses the hair and the matter field is just a harbor for the black hole. Fig.4 is quite illustrative and it shows three-dimensional plots of the harbor solutions where, in concordance with the metric, the event horizon is depicted by black spheres. According to the ansatz (15) the density of matter fields possesses spherical symmetry, i.e., |u † · u| = f (r) 2 . The distance between the horizon and the inner boundary grows with increasing of n. (Note that the reason why the coupling constant was chosen α = 0.001 for the case n = 50 was the fact that we were not able to find solutions above a critical value α = αcrit < 0.01.) The stability of the Q-balls was studied in many literatures, see for instance [1,2,[26][27][28]. An important relation for analysis of the (absolute) stability of classical solutions is the scaling relation of total energy with the Noether charge of Q-ball or Q-shell. This relation was examined analytically in [21] where authors study the model (2) in the flat spacetime. In the limit of ω → ∞ (small amplitude fields) the model reduces to the signum-Gordon model [13] which has exact solutions and so the relation between the energy and the Noether charge is expected to be E ∼ Q 5/6 . Unfortunately, for the curved space-time such approximation results in equations which cannot be solved using analytical methods. For this reason we study the relation between the energy and the Noether charge numerically and without using any approximations. In Fig.5 we present the relation between E −1/5 and Q −1/6 for self-gravitating solutions in the models with n = 0, 1, 2, 5, 10, 20, 50. The dots indicate the solutions with given frequency ω. The energy decreases with increasing of ω. We draw dashed lines which connect the same value of ω on lines of different n. According to graphics presented in Fig.5 the relation between E −1/5 and Q −1/6 is linear with very good accuracy so one can conclude that the energy scales with the Noether charge as E ∼ Q 5/6 . We have also examined the energy-charge relation for harbor type solutions. The results of this analysis are presented in Fig.6 where we plot E −1/5 in function of Q −1/6 for n = 2, 5, 50 (i.e.,CP 5 , CP 11 , CP 101 ) for several values of ω and r h . The dots with the same color   Table I. represent solutions with given ω and different r h . The energy decreases with increasing of the horizon radius. It means that dots of the same color are ordered from r h = 0 (lower left) to the maximum r = r h,max (upper right). The values of r h,max corresponding to this plot are listed in Table I. The gradients of the same color are always close to the gradient of the line of r h = 0 (see Fig.5). It means that the presence of a black hole is not so apparent for the scaling relation. As a result, a linear relation between quantities E −1/5 and Q −1/6 holds for both Q-ball and Q-shell solutions.
In [21], the authors examined the Q−ω relation. They numerically confirmed that it has the form Q ∼ ω −1/6 . It thus enjoys a condition of the classical stability, 1 Q dQ dω < 0, proposed in [1]. We have examined this relation for the gravitating solutions. The result is plotted in Fig.7. It follows that all solutions with different values of n are stable.
In our opinion, some further study is required for better understanding of a relation between the energy and the Noether charge. The reasons are as follows: First, most of the previous studies exhibit some limitations on the frequency ω. Our solutions were obtained starting with a moderate value ω = 1 and then gradually increasing its value. An interesting question is if its value can be taken arbitrary small. We shall study this question in more detail. Second, Q-ball solutions reported in the literature have usually no restrictions on maximal value of |Q|. In the gauged U (1) symmetry model, the solutions are unstable for large |Q|, which can be overcome taking fermions of the scalar condensates with opposite charge, i.e. −|Q|, in the interior and so the charge can reach values around ∼ 1 × 10 8 [11]. Also in the scalar electrodynamics, presented in [14], the authors found Q-balls with scaling relation E ∼ |Q| 7/6 for large Q. This result may indicate instability of solutions for large Q and their tendency to evaporate. In [21], the authors do not report on possible existence of upper limit for Q. (Instead, they discuss the scaling relation based on the signum-Gordon limit of the model ω → ∞.) Thus there is certain variety of contexts in which the relations between Q, ω and E are presented. These analyses were conducted in the flat space-time, i.e., the effects caused by gravity has always been ignored. Certainly it is worth to investigate deeper such relations for the systems where the gravity is taken into account. One of our motivation for such a study is checking if there can exist an upper limit for the charge |Q| of gravitating systems. In the next section, we shall report on this issue in detail.

IV. FURTHER ANALYSIS
According to results reported in the literature the frequency ω is restricted to a certain finite interval. Compact Q-balls/shells are rather exceptional in this context because their existence do not require such restrictions. In fact, for a V -shaped potential there is no upper bound for ω and the lower bound of ω is zero (for the signum-Gordon model) or given by certain positive constant (the case of model (1)). In both cases admissible frequencies ω belong to infinite intervals. For this reason we shall look in more detail at the region of Fig.5 corresponding with ω ≤ 2.0. The results of our analysis are shown in Fig.8. In the case scaling relation between E −1/5 and Q −1/6 the linearity holds strictly even without assuming that ω → ∞. A more interesting relation is that between Q −1/6 and ω. For small n, the relation is linear, however, as n grows, its exhibits distortion from the linearity and simultaneously the lower bound of ω is gradually shifted in direction of higher values of ω. The derivative dQ −1/6 dω decreases with increasing of n. In particular, for n = 50 the derivative approaches zero for small ω. This means that the condition for classical stability is not fulfilled any longer. It indicates the existence of a lower bound of ω for gravitating solutions. Now we would like present some arguments which would shed light on the role of number of fields N := 2n + 1. In Fig.9 and Fig.10, we plot the Noether charge Q in dependence on the number n for several coupling constants α and two values of the frequency ω = 1.0 and ω = 0.95. For better transparency we plot also the interpolating function which depends on a continuous variable (also called n). In the case of flat space-time, α = 0, the relation is almost linear, which means that the model with given N has a solution with corresponding particle number Q. Therefore, both Q and N indicate the particle number of the constituents. In right figures we present plots of the derivative dQ dn for the interpolating functions. Notable property of solutions in the curved spacetime is that the correspondence between N and Q is broken. It means that the increment of n by one causes in generality a moderate change of the charge Q. However, at critical value of α the charge Q has maximum (which is seen from vanishing of dQ dn at this point), and for larger n the charge Q decreases with increasing of n. The peak moves to the left and became higher with decreasing of ω. When ω grows it moves to the right and became lower. In Fig.11, we draw for some cases of ω for α = 10 −2 and one can easily confirm the behavior. We conclude that above the critical value of the coupling constant α there is a maximal value of the charge Qmax. Although this value became larger for smaller values of ω the system eventually tends to instability as discussed above. We conclude that the gravitating solutions above the critical coupling have also an upper bound for |Q|.
It may give us a hint how a boson star with certain number of particles forms its shape. The systems with large n can be achieved via a weak coupling constant α and also a small frequency ω.

V. SUMMARY
In this paper, we have considered the family of CP N nonlinear sigma models coupled with gravity which possess solutions with compact support. Such solutions, in form of compact Q-balls and Q-shells, were obtained solving numerically the system of coupled equations. The self-gravitating solutions of this kind are interpreted as boson stars. We have also considered the case of spherical Q-shells with the Schwarzschild-like black hole immersed in their center. The space-time in the region where hold the Einstein's vacuum equations (interior end exterior of the shell) is described by the Schwarzschild metric. Although such solutions are possible candidates which would contradict the no-hair conjecture, we have shown that in its strict sense, this is not the case of our solution. Notable difference from solutions in the flat space is the existence of upper bound for |Q| for coupling constants α bigger than certain critical value. We have checked that the value of this bound drastically decreases with increasing of gravitational coupling constant. For sufficiently small coupling constant, the maximal value of the |Q| is reached for large n, It might be a possible explanation why the boson stars have definite size and mass.
We restricted our analysis to the case of asymptotically flat Schwarzschild space-time. We managed to obtain the stable shell solutions for n ≧ 2 even though the matter fields are not coupled with electromagnetism. The the model (22) possesses symmetry SU (2) ⊗ U (1). In current research we just focused on the subgroup U (1) N of the full symmetry group. Further, the ansatz (15) reduces the symmetry to U (1) and the corresponding Noether charge plays central role in stability of our solutions.  The n /dn Figure 10. Same as the plots of Fig.9 but for ω = 0.95.
This paper constitutes an initial step for construction of gravitating Q-balls (-shells) with non-Abelian symmetry SU (2) ⊗ U (1). The following problems require solution: • Extensions of our model onto the model with finite cosmological constants, i.e. onto the de Sitter or anti-de Sitter space-time background. Such extension is almost straightforward and it will be reported in the subsequent paper.
• A coupling with electromagnetism which will allow us to treat the charged boson stars and field configurations with the Reissner-Nordström-type metric. The analysis of these models is in progress..
• A study of non-Abelian Noether charges. There are know Q-balls in scalar field models with non-Abelian symmetry [29,30]. The boson stars with a non-Abelian SU(2) symmetry are discussed in [31,32]. The existence of Q-balls with the symmetry SU (2) ⊗ U (1) is highly probable. It would be an interesting possibility because it would two types of Noether charges which have crucial role for the stabilization of nontopological solitons.
We shall report on these issues in near future. The n Charge Q Figure 11. The Noether charge Q in function of n. We plot some cases around ω = 1.0 with the value α = 10 −2 .