Nonlinear instability and scalar clouds of spherical exotic compact objects in scalar-Gauss-Bonnet theory

In this work, we present a new type of scalar clouds supported by spherically symmetric horizonless compact objects in the scalar-Gauss–Bonnet theory. Unlike the previous spontaneous scalarization that is triggered by the tachyonic instability, our scalarization arises from a nonlinear instability that is non-spontaneous. We explore two types of boundary conditions for the scalar field at the surface of the compact objects and find an infinite countable set of scalar clouds characterized by the number of nodes for both cases. Our study demonstrates that boundary conditions have a significant impact on the formation of scalar clouds. Specifically, for the Dirichlet boundary condition, scalarization is more likely to occur for compact objects with medium radii and becomes harder for ultra-compact and large ones. Conversely, for the Robin boundary condition, scalarization is easier for more compact objects.


I. INTRODUCTION
Black holes (BHs) are one of the most intriguing predictions of general relativity (GR) and have attracted extensive attention in recent decades due to their significance.The detection of gravitational waves [1] and the imaging of BHs [2] provide compelling evidence of their existence.
However, several challenging issues, such as the information loss paradox [3] and the emergence of unphysical curvature singularities [4], remain unresolved.To address these issues, physicists have proposed several alternatives to BHs, collectively known as exotic compact objects (ECOs) [5][6][7][8][9][10][11][12][13], whose mass and compactness can closely match those of BHs, making them potential candidates for explaining astronomical phenomena attributed to BHs.Theoretically, BHs and ECOs may be seen as distinct "species" of compact objects, which could co-exist in our universe.Notably, some recent gravitational wave mergers have the possibility to involve ECOs [14].Nevertheless, our understanding of ECOs is significantly less than that of BHs.The exploration of ECOs' physical properties may help us to gain a deeper understanding of BHs and quantum features of gravity.
With advancements in astronomical observations, particularly in the gravitational waves channel, we expect to access the near-horizon (surface) structure of BHs (ECOs) soon and differentiate between them.For a review, see [15,16].
A generic feature of ECOs is the replacement of the classical absorbing horizon by a reflective surface.A simple way to model this feature is proposed in [13], in which spacetime geometry is only modified at some microscopic scale near the would-be classical horizon while exterior geometry is still described by the usual Kerr metric.Without causing confusion, we will refer to this model as horizonless compact objects in the following.In the framework of GR, one of the most important properties of BHs is given by the no-hair theorem, which states that asymptotically flat static BHs are uniquely characterized by only two physical parameters, their mass and angular momentum, and cannot carry scalar or other types of hair [17][18][19].Therefore, it is interesting to see how modification of the structure close to the would-be horizon affects this theorem and whether ECOs can carry additional hairs.
Interestingly, it is found that the no-hair theorem applies not only to BHs but also to ECOs in general.In GR, it is found that spherically symmetric horizonless compact objects with a perfectly reflecting surface cannot support minimally coupled scalar, vector or tensor clouds [20,21].The theorem also holds when the scalar field is non-minimally coupled to the Ricci scalar [22,23].
However, if the objects begin to rotate, the situation changes and the no-hair theorem no longer holds.It has been discovered that scalar clouds, massless or massive, can form for rapidly and ultra-spinning horizonless compact objects with a perfectly reflecting surface at certain discrete radii [24,25].Physically, the formation of scalar clouds in this case can be understood from the mechanism of superradiant instability.Namely, with the absorbing horizon replaced by a reflecting surface, spinning ECOs in general suffer from superradiant instability under scalar perturbations [13], which results in the formation of scalar clouds.A similar phenomenon has already been observed in the BH scenario [26,27], but there only massive scalar cloud/hair can be formed.
Another physical mechanism that may endow ECOs with scalar cloud is the tachyonic instability.This phenomenon is inspired by the spontaneous scalarization of BHs in certain modified gravity theories [28][29][30][31][32][33] (see also [34] for a recent review).The Einstein-scalar-Maxwell (EsM) theory and the scalar-Gauss-Bonnet (sGB) theory are two well-known models of spontaneous scalarization, in which the scalar field is non-minimally coupled to the Maxwell term and the Gauss-Bonnet term, respectively.In both models, it is found that the scalar field perturbation of spherically symmetric horizonless compact objects with perfectly reflecting surface acquires a negative effective mass square, triggering the tachyonic instability, which results in the formation of scalar clouds [35,36].
Recently, A new physical mechanism that could induce scalar hair on BHs in EsM or sGB has been discovered [37,38], named as nonlinear instability.Different from the tachyonic instability, this type of instability occurs only when the scalar field perturbation becomes relatively large and scalar hair formed are thus non-spontaneous.Motivated by these works, we investigate whether ECOs can also support scalar clouds through this mechanism.This is the main goal of this paper.
The paper is organized as follows.In the next section, we will give a brief review of our model.In Sec.III, we perform time evolution of scalar perturbations of spherically symmetric horizonless compact objects with Dirichlet boundary condition and show the occurrence of the nonlinear instability.In Sec.IV, we construct the scalar clouds with Dirichlet boundary condition induced by the nonlinear instability.In Sec.V, we consider Robin boundary condition to show the influences of different boundary conditions on the formation of scalar clouds.The last section is the Summary and Discussions.

II. THE MODEL
We consider the scalar-Gauss-Bonnet (sGB) theory with the action [28][29][30][31] where the scalar field ϕ is coupled to the Gauss-Bonnet term through the Gauss-Bonnet coupling constant λ and the coupling function f (ϕ).As our main goal in this work is to search non-spontaneous scalarization of ECOs, we choose the coupling function to satisfy two following conditions With the first condition, the theory admits GR vacuum solutions with vanishing scalar fieldthe Kerr metric.In this work, following the approach proposed in [13], we consider a spherically symmetric horizonless compact object whose exterior geometry can be well described by the Schwarzschild metric where g(r) = 1 − 2M r with M being the mass of the object.The object surface locates at r = r s which should be outside the would-be horizon r = r h = 2M .On this background, dynamics of the scalar field outside the object is governed by the modified Klein-Gordon equation with R GB = 48M 2 r 6 .From it, one can see that the scalar field acquires an effective mass square dϕ 2 (ϕ = 0)R GB , which vanishes identically with the second condition of Eq. ( 2).This thus excludes the occurrence of the tachyonic instability, and so the scalarization, if exists, will be non-spontaneous.
Taking into account the conditions (2) and also the possible existence of stable scalar clouds, we consider the coupling function to take an exponential form with κ being a parameter.This kind of coupling function has also been considered to discuss non-spontaneous scalarization of BHs in sGB [38].

III. NONLINEAR INSTABILITY
In this section, we will first study the wave dynamics of the scalar field on the background (3).
As will show later, with this kind of coupling function (5), the scalar field will experience a kind of nonlinear instability.
For simplicity, we assume the scalar field perturbation to be spherical symmetric, ϕ = ϕ(t, r).
In tortoise coordinate dx ≡ dr g(r) , the scalar field equation (4) becomes which can be solved numerically by adopting the method of line [39].To solve the equations, physical boundary conditions are needed.At the object surface r = r s , there are usually two boundary conditions considered, the Dirichlet and Robin boundary conditions.We will give more comments on the two boundary conditions in the last section.Let us first consider a Dirichlet boundary condition, so the physical boundary conditions we need to impose are that the scalar field vanishes at the object surface and is outgoing at infinity.
We consider the initial perturbation to be a time-symmetric Gaussian pulse with x c = 20M, σ = 2M and A the perturbation amplitude.There are four free parameters in the model, {M, r s , κ, λ}.We fix M = 1 and so all quantities are measured in units of it.By performing time evolution of the scalar field perturbation, we found that, depending on values of the parameters, nonlinear instability may be triggered.In Fig. 1, two typical examples are given.
From the figure, one can see that the occurrence of the instability depends on the amplitude of the perturbation: When the amplitude A is small, the perturbation will exhibit a decaying late-time tail; While A becomes larger, the nonlinear term of ϕ in the equation enters the game and the scalar field finally settles down to a equilibrium state indicating the formation of scalar cloud.This kind of instability, named as nonlinear instability, has already been observed in the BH scenario [38].

IV. SCALAR CLOUDS WITH DIRICHLET BOUNDARY CONDITION
In the last section, we have observed that the scalar field may experience a nonlinear instability, which is a signal of the formation of the scalar cloud.In this section, we will construct the scalar cloud explicitly.For simplicity, we consider the scalar cloud outside the object to be static and spherical, i.e., ϕ = ϕ(r).After substituting the metric (3) into (4), the scalar field equation becomes To solve the above equation, physical boundary conditions are needed.Assuming that a Dirichlet boundary condition at the object surface and the scalar field outside the object is regular and bounded, we have the boundary conditions .n labels the number of nodes that the solution possesses between the object surface and radial infinity.We fix κ = 1, r s = 2.5M (left panel) and r s = 5M (right panel).Now, the model contains four free parameters {M, r s , λ, κ}.We can directly integrate the scalar field equation ( 8) from the object surface to infinity for fixed values of these parameters.
The number of free parameters can be reduced by noting that the scalar field equation ( 8) possess two scaling symmetries and where a, b are arbitrary scaling parameters.With the two symmetries, it is convenient for numerical calculations to fix M = 1 and ϕ 1 ≡ ϕ (r s ) = 1, where ϕ 1 is the first radial derivative of the scalar field at the surface.So finally, there are left only three free parameters {r s , λ, κ}.For a given κ and r s , the solution is determined uniquely by λ.However, not every value of λ will produce a bounded solution that satisfies ϕ(∞) = 0. Ony for certain discrete values of λ can produce bounded solutions.
For a fixed κ and r s , we found an infinite countable set of the coupling constant, {λ(κ, r s ; n)} n=∞ n=0 , which can support the bounded scalar clouds.Here a larger integer n labels a larger λ.
In Fig. 2, we show configurations of the scalar clouds for {λ(κ = 1, r s ; n)} n=5 n=0 with r s = 2.5M and r s = 5M as examples.From the figure, one can see that the scalar clouds exhibit a standingwave-like profile with n indicating the number of nodes that the solution has between the object surface and radial infinity.Moreover, radial oscillations of the scalar clouds are concentrated in the region near the object surface.As the object radius r s increases, the oscillations move outward.
To survey the relation between r s and the coupling constant λ, we plot {λ(κ, r s ; n)} n=5 n=0 as a function of r s in Fig. 3. From the figure, one can see that, for each n, λ(κ, r s ; n) is a convex function of r s with a minimum value λ(κ, r s ; n) min .It is interesting to note that for any n, the minimum value appears at r s ≈ 2.5M , which is close to but behind the photon sphere of the object.This suggests that scalar clouds are most likely to form when the object is very compact and has radius about r s ≈ 2.5M .When λ < λ(κ, r s ; n = 0) min ≈ 21.45, no scalar cloud can be supported for any r s .When λ(κ, r s ; n = 0) min ≤ λ < λ(κ, r s ; n = 1) min , only fundamental solution with n = 0 exist.As λ increases, excited solutions with higher n emerge.Moreover, as r s approaches the would-be horizon r h = 2M , {λ(κ, r s ; n)} show a divergent behavior, which implies that ultra-compact objects are hard to support scalar clouds.

V. SCALAR CLOUDS WITH ROBIN BOUNDARY CONDITION
In the above, we have shown that compact objects can support scalar clouds with Dirichlet boundary condition.Another boundary condition usually considered is the Robin boundary condition [13] As in the case with Dirichlet boundary condition, the scalar field perturbations may also trigger nonlinear instability, as demonstrated in Fig. 4. Compared to the Dirichlet case in Fig. 1, we can see that, for small A, the perturbation experiences a longer ringdown phase before the late-time decaying.  .n labels the number of nodes that the solution possesses between the object surface and radial infinity.We fix κ = 1, r s = 2.5M (left panel) and r s = 5M (right panel).
Similarly, we can construct the scalar cloud explicitly.In this case, it is convenient for numerical calculations to set M = 1 and ϕ(r s ) = 1 with the two scaling symmetries (10) and (11).At infinity, we also have ϕ(∞) = 0.As the case with Dirichlet boundary condition, for given values of parameters {κ, r s }, there also exists an infinite countable set {λ(κ, r s ; n)} n=∞ n=0 which supports the scalar clouds.Samples are shown in Fig. 5. Similar to the case with Dirichlet boundary condition, the scalar clouds also exhibit a standing-wave-like profile with the integer number n labeling the number of nodes that they possess between the object surface and infinity.Compared to Fig. 2, we can see that radial oscillations of the scalar clouds with Robin boundary condition are more concentrated in the region near the surface.And as the object radius r s increases, the oscillations move inward instead.The relation between the object radius r s and the coupling constant λ is shown in Fig. 6.From the figure, one can observe one important feature different from that of the case with Dirichlet boundary condition: For each n, λ is now a monotonically increasing function of r s .This suggests that with the Robin boundary condition, scalarization becomes easier for more compact objects.

VI. SUMMARY AND DISCUSSIONS
In this work, we construct a novel type of scalar clouds for spherically symmetric horizonless compact objects in sGB theory.These scalar clouds are not formed due to tachyonic instability, but rather due to nonlinear instability, which requires a certain threshold of perturbation magnitude.
We find that there exists an infinite countable set of coupling constants λ(κ, r s ; n) n=∞ n=0 (where n is the number of nodes between the surface of horizonless compact objects and infinity) that can support the scalar clouds for a fixed radius r s , given other parameter values.Moreover, we show that boundary condition (either Dirichlet or Robin) significantly affects the formation of scalar clouds.
With Dirichlet boundary condition, λ(κ, r s ; n) is a convex function of r s , indicating that scalarization is most likely for compact objects with medium radii and becomes more difficult for ultra-compact and large ones.However, with Robin boundary condition, λ(κ, r s ; n) is a monotonically increasing function of r s , suggesting that scalarization is easier for more compact objects.In both cases, there exists a threshold coupling constant λ below which no scalar clouds can be supported.
As λ increases, fundamental solution with n = 0 and excited solutions with higher n emerge successively.
In the context of Einstein's gravity with minimally coupled scalar field, these boundary conditions are commonly associated with a perfectly reflecting surface [13].However, this interpretation is no longer valid in our theory due to the presence of a nonlinear term in ϕ, which prevents the scalar field perturbation equation ( 6) from taking a wave-like form at the surface.Only when the perturbation is sufficiently small can we ignore the nonlinear term and regain this physical meaning.For general perturbations, the physical meaning of these boundary conditions is unclear and requires further investigation.
The coupling function we consider takes an exponential form with the exponent proportional to ϕ 4 .We have also examined some other exponential forms, such as the one with exponent ϕ 6 , and found similar instability and scalar clouds.
There are several possible extensions of this work.We work in the "decoupling limit" where the influence of the scalar clouds on the background geometry is neglected.This limit has shown its ability to capture the qualitative features of the fully nonlinear dynamics [40,41].However, it would be interesting to go beyond this limit and develop complete hairy solutions for further investigation.This would require specifying the explicit equation of state of the object and doing analysis model by model.In this paper we only consider spherically symmetric ECOs.A natural extension is to study the rotating case, where more than one physical mechanism may affect the formation of scalar clouds.Another intriguing question is whether other types of ECOs can support such scalar clouds.

FIG. 1 .
FIG. 1.Time evolution of the scalar field perturbations on the spherically symmetric horizonless compact object with Dirichlet boundary condition.The perturbation is a Guassian pulse ϕ(t = 0, x) = Ae − (x−xc ) 2 2σ with x c = 20 and σ = 2. Signal is extracted at x = 20.We fix M = 1 and all quantities are measured in units of M .

FIG. 2 .
FIG. 2. Configurations of the scalar clouds with Dirichlet boundary condition for various coupling constant{λ(κ, r s ; n)} n=5 n=0.n labels the number of nodes that the solution possesses between the object surface and radial infinity.We fix κ = 1, r s = 2.5M (left panel) and r s = 5M (right panel).

FIG. 4 .
FIG. 4. Time evolution of the scalar field perturbations on the spherically symmetric compact object with Robin boundary condition.The perturbation is a Guassian pulse ϕ(t = 0, x) = Ae − (x−xc ) 2 2σ

FIG. 5 .
FIG. 5. Configurations of the scalar cloud with Robin boundary condition for various coupling constant{λ(κ, r s ; n)} n=5 n=0.n labels the number of nodes that the solution possesses between the object surface and radial infinity.We fix κ = 1, r s = 2.5M (left panel) and r s = 5M (right panel).