Horizon curvature and spacetime structure influences on black hole scalarization

Black hole spontaneous scalarization has been attracting more and more attention as it circumvents the well-known no-hair theorems. In this work, we study the scalarization in Einstein-scalar-Gauss-Bonnet theory with a probe scalar field in a black hole background with different curvatures. We first probe the signal of black hole scalarization with positive curvature in different spacetimes. The scalar field in AdS spacetime could be formed easier than that in flat case. Then, we investigate the scalar field around AdS black holes with negative and zero curvatures. Comparing with negative and zero cases, the scalar field near AdS black hole with positive curvature could be much easier to emerge. And in negative curvature case, the scalar field is the most difficult to be bounded near the horizon.

Experimental progress on gravitational waves [1][2][3] and the shadow of the M 87 black hole [4] further demonstrates the great success of Einstein's general relativity (GR). Yet it is unabated that the GR theory should be generalized, and in the generalized theories extra fields or higher curvature terms are always involved in the action [5][6][7]. So physicists have proposed various modified gravitational theories which indeed provide richer framework and significantly help us further understand GR as well as our universe.
Among them, the scalar-tensor theories which introduces a scalar field into the action attract lots of attention [8]. When the scalar field backreacts to the background metric, one could expect that hairy black hole solutions would be generated. One of the first hairy black hole solution in an asymptotically flat spacetime was discussed in [9] but soon the solution was argued to be unstable because the scalar field is divergent on the event horizon [10]. However, such irregular behaviour of the scalar field on the horizon was avoided in asymptotically AdS/dS spacetime with a presence of a cosmological constant in the gravity theory. It was found that the resulting hairy black hole solutions have a regular scalar field behaviour and all the possible divergence could be hidden behind the horizon [11][12][13][14][15][16][17][18][19][20][21][22].
The effects of higher-order curvature terms become more significant as we are exploring the strong field regime of gravity via detections of gravitational waves and black hole shadows. It is known that the inclusion of such terms probably bring in the well-known ghost problem [23], and Gauss-Bonnet (GB) corrections is a counter-case which is ghost-free. But it becomes a topological term in four-dimensional spacetime and has no dynamics in field equations when it is minimally coupled with Einstein-Hilbert action. A way to make this term meaningful in four-dimensional spacetimes is to consider its coupling with a scalar field [24]. As a special scalar-tensor theory with higher derivatives [25,26], the Einsteinscalar-Gauss-Bonnet (ESGB) gravity recently has attracted a lot of attention. Especially, the introduction of this coupling could lead to hairy black holes. Various black hole solutions and compact objects in the four-dimensional ESGB theories were studied in the literature [27][28][29][30][31][32].
Recently, the spontaneous scalarization with particular coupling function in ESGB theory is widely investigated. In this setup, besides GR solutions with a trivial scalar field configuration, the scalarized hairy solutions for black holes and stars could also exist, which evades the no-hair theorems [33][34][35]. It was shown in [36] that below a certain mass the Schwarzschild black hole background may become unstable in regions of strong curvature, and then when the scalar field backreacts to the metric, a scalarized hairy black hole emerges and it is physically favorable. A natural extension of these results was to consider the case of nonzero black hole charge. A large set of coupling functions between the Gauss-Bonnet invariant and the scalar field was considered in order to understand better the behaviour of the scalarized solutions [37]- [39] .
Various generalizations of the spontaneous scalarization procedure were followed in the literature. The scalarization due to a coupling of a scalar field to Ricci scalar was studied in [40] and in the presence of Chern-Simons invariant was studied in [41]. The generalized study of scalarized black hole solutions and compact objects in asymptotical flat spacetime of ESGB theories was discussed in [42]- [52]. The spontaneous scalarization of asymptotically AdS/dS black holes in ESGB theory with a negative/positive cosmological constant was extended in [53][54][55][56][57]. Especially, the connections of asymptotically AdS black holes scalarization with holographic phase transitions in the dual boundary theory was studied in [58,59]. Recently the spontaneous scalarization in f(R) gravity theories in the presence of a scalar field minimally coupled to gravity with a self-interacting potential was discussed in [60]. Also in ESGB theory the spin-induced black hole spontaneous scalarization, which is the outcome of linear tachyonic instability triggered by rapid rotation, was explored in [61][62][63][64][65][66].
The scalarization in asymptotical flat, AdS and dS black holes in ESGB gravity theory was studied, however detailed information on how different spacetime structures influence the scalarization process is not disclosed. An interesting question is how the cosmological constant leaves its imprint on the black hole scalarization and in which kind of spacetime the scalar field could be formed more easily. Besides, the studies of spontaneous scalarization in the existed literatures were focused on the black hole with spherical horizon, but as known in AdS spacetime, one can have AdS black holes with different topologies and it is interesting to study how these topologies of the black hole horizon affect the scalarization procedure 1 . Thus, the aim of this work is to probe the signal of scalarization in various black holes with different horizon curvature and try to understand its effects on the scalarization processes. To this end, we first consider the test scalar field perturbation and analyze its effective potential in various black hole backgrounds. Then we study the tachyonic instabilities of the scalar field near the horizon by solving the perturbation equation. We find that the scalar hair around the black hole with spherical horizon is easier to form in AdS spacetime than in flat case. While in dS case, the regular condition near cosmological horizon leads to divergence at large distance, so that the test field approximation breaks down. Moreover, our study on the effect of horizon curvatures indicates that the scalar hair around AdS black hole with toroidal horizon could be easier to be formed than that with hyperbolic horizon, but more difficult than that with spherical horizon. Both the effective potential analysis and the test field solution conclude the above observations. The work is organized as follows. We write down the general equation of motion of a test scalar field in the background of four dimensional black holes with different topologies in ESGB gravity theory in section II. In section III, we study the signal of scalarization of a black hole with positive curvature in different spacetimes. Then we investigate the signal of a black hole scalarization with negative curvature and zero curvature in section IV and section V respectively. The last section is devoted to conclusions and discussions.

BONNET THEORIES
We consider four-dimensional ESGB theories with the action where Λ is the cosmological constant, λ is the coupling constant between Gauss-Bonnet term and scalar field with the dimension of length and ∇ µ denotes the covariant derivative. The scalar field coupling function f (Φ) only depends on the neutral scalar field Φ, and the Gauss-Bonnet term is given by where R, R µν and R µναβ are the Ricci scalar, Ricci and Riemann tensors, respectively. The Einstein equation and Klein-Gordon equation can be derived from the above action as where the energy-momentum tensor Γ µν is given by with Ψ µ = λ 2 df (Φ) dΦ ∇ µ Φ. It is obvious that different forms of the coupling function f (Φ) shall give different properties of the ESGB theory. As addressed in [36], f (Φ) could satisfy conditions df (Φ) dΦ | Φ=0 = 0 and d 2 f (Φ) dΦ 2 | Φ=0 = b 2 > 0, where b is a constant, in order to admit Schwarzschild black hole as a background solution and further explore the black hole scalarization in the ESGB theory. Moreover, one usually assumes that the scalar field vanishes at infinity and normalizes the constant b to unity. Thus, in this work, we shall choose the form of the coupling function as When the scalar field vanishes, i.e., Φ = 0, the gravity theory admits black hole solutions, whose metric takes the general form, where M is the integral constant related the black hole mass, and k determines the curvature of two dimensional geometry near the horizon with the line element d − → x 2 k,2 given by where k = 1, 0 and −1 correspond to black hole horizon with spherical, planar and hyperbolic topology, respectively. Then the corresponding GB term can be easily computed as In order to explore the (in)stability of the background black holes, we consider a small fluctuation in the form . We consider that Y (x 1 , x 2 ) is spherical harmonics Y ,m for positive curvature so that A = ( + 1); For zero curvature Y (x 1 , x 2 ) is plane wave case as e iK(x+y) so that A = K 2 ; and for negative curvature Y (x 1 , x 2 ) is spherical harmonics Y − 1 2 ±iζ,m and A = 1 4 + ζ 2 [68,69]. Then we use the tortoise coordinate dr * = dr g(r) , and bring the Klein-Gordon equation into a Schrodinger-like form by reforming φ(r) = ϕ(r) r ∂ 2 ϕ(r) ∂r 2 * where the effective potential takes the form Noted that in the Schrodinger-like form (9) and (10) of the perturbation, we have only kept the leading term of the expansion df (Φ) ) in the potential function. We then analyze the signal of scalarization of the system by computing the effective potential with the lowest mode ( = K = ζ = 0). For simplification, we could focus on the static perturbation with ω = 0, thereby we shall solve the radial perturbed equation for φ(r), and see how the scalar field with the lowest mode emerges as the Gauss-Bonnet coupling increases. To control the variables of the black holes and compare the results, we shall set r h = 1 in various cases, namely, we fix the size of the black holes such that the black hole mass is solved via g(r = r h ) = g(r h = 1) = 0 in the following discussion. We would numerically employ the spectral method to solve the differential equations.
Since the existence of above black hole with different curvatures is dependent on the type of asymptotic background, i.e., the sign of Λ, we would discuss the signal of scalarization of various black holes in terms of the curvature. Moreover, we shall rewrite Λ = ±3/L 2 where L is the curvature AdS radius. As it was discussed in [59,70] because the scalar field is coupled to GB invariant, the effective mass of the scalar field receives a contribution from this coupling. Then for various values of the coupling, the effective mass of the scalar field becomes tachyonic outside the background black hole and the fluctuations are unstable destabilizing in this way the background metric.

VATURE
In this section, we focus on the background solution with positive curvature, i.e., k = 1 in (6), in which case there exist asymptotical flat (Λ = 0), AdS (Λ < 0) and dS (Λ > 0) black holes as possible backgrounds. Consequently, we shall separately study the signal of scalarization in all these asymptotical backgrounds, and then do comparisons on the emergence of non-trivial scalar field which is introduced by the Gauss-Bonnet coupling in these different spacetime backgrounds.

A. Asymptotical flat black hole
For Λ = 0, the asymptotical flat black hole (AFBH) with spherical horizon (k = 1) is just the Schwarzschild black hole. In this case, the redshift function in the field equation (11) is ). In this case, the potential (10) for the lowest mode = 0 of the massless scalar field (m = 0) is where we have inserted M = 1 in the second equality. A sufficient condition for the existence of an unstable scalar field mode is [36,71] which leads to λ 2 > 5 12 M 2 = 5 12 . This indicates that the tachyonic instability emerges when the coupling constant satisfies λ 2 > 5 12 . Especially, near the horizon r → 1, we can reduce from equation (11) that the scalar field φ(1) and its derivative φ (1) are both regular with a relation It is obvious that the contribution of the GB coupling could control the near horizon behaviour of the scalar field. While at infinity r → +∞, the scalar field behaves as We would set φ ∞ = 0 as the Dirichlet boundary condition so that the scalar field converges to zero as r → ∞ which is well behaved. Thus, we can consider the scalar field as a probed field to study the signal of the scalarization. Also φ − could be treated as the scalar charge. In Fig. 1, we show the effective potential of the scalar field as the function of radial coordinate. When λ = 0 without GB coupling, there is a potential barrier. As the GB coupling is turned on and increased, the barrier is suppressed, and a negative potential well form near the horizon which could trigger an instability and then lead to the development of scalar hair. We then study the behaviour of scalar field by solving the perturbation equation (11) and expect to see the signal of scalarization in the probe limit. By scanning the coupling parameter, it is found that there exists a threshold value λ t 0.60234 2 . When λ is smaller than the threshold value, the scalar field is zero everywhere which is trivial. When the coupling is larger than λ t , the scalar hair would emerge and increase as the coupling increases. The profile of the scalar field with various values of λ is shown in the left plot of Fig. 2. Moreover, in the middle plot of Fig. 2, we present the value of the scalar field near the event horizon φ h as the function of λ and corresponding scalar charge φ − at infinity is shown in the right figure. The black line on the horizontal axis is the trivial scalar field. As the coupling parameter increases larger than λ t , there emerges the bifurcation dividing the trivial solution, the non-trivial scalar field formes and grows dramatically. This phenomenon indicates that due to the instability brought by the GB coupling, the scalar field may emerge and grow 2 Noted that previously we discussed that for λ 2 > 5 12 M 2 , i.e. λ 0.645, the background black holes would be unstable. Here the threshold value λ t 0.60234 is a bit smaller. This is reasonable because (13) is only a sufficient condition, which was also verified in [36]. near the horizon of AFBH background, such that a new branch of scalarized hairy black hole develops. It is noticed that the scalarized hairy black hole solutions with backreaction were numerically explored in [37], and it was found that the entropy of the scalarized black hole is larger than of AFBH solution which indicates the hairy solution is physically favorable.

B. Asymptotical AdS black hole
In the AdS case with Λ = −3/L 2 < 0, the redshift function in the solution (6) with spherical horizon is g(r) = 1 − M r + r 2 L 2 (M = 1 + 1 L 2 as r h = 1) and the GB term As is known, the AdS background could possess a classical instability for the real scalar field, but provided its effective mass is above the Breitenlohner-Freedman (BF) bound, the perturbation do not break the AdS behavior [72]. In our case, the effective mass derived from (4) is given by m 2 ef f = m 2 − λ 2 2 R 2 GB , and so the regular tachyonic instability could exist as m 2 ef f > m 2 BF = −9/4L 2 . Specially, near infinity, the Klein-Gordon equation becomes which gives the boundary behavior of the scalar field as Here and m 2 e L 2 = m 2 L 2 − 12 λ 2 L 2 is the asymptotical behavior of m 2 ef f L 2 . It is obvious that m 2 e > m 2 BF = −9/4L 2 could also ensure the scalar field at infinity to be convergent to be zero. Thus, in AdS case, for the parameters satisfied the above analysis, we are safe to solve the probed perturbative equation (11) to see how exactly the GB coupling and AdS curvature radius affect the fate of the scalar field.
For a proper comparison in the flat and AdS cases, we set m 2 e L 2 = −2 and then the asymptotic behaviour becomes φ(r → ∞) = φ − /r + φ + /r 2 , such that the scalar charge φ − term has the same decay behavior as that in flat case. In the numeric, we set φ + = 0 as the boundary condition. It is noted that in general, the asymptotic behaviour of the scalar charge term could be different in flat and AdS cases. Thus, to compare the signal of scalarization, we have to carefully control the variables and then we require the decay behavior of the scalar charge term is coincident. It means that we shall probe the scalar field with the same scalar charge but not the same scalar mass, which is reasonable because non-vanishing scalar charge implies the occurrence of scalarization and modifies the general relativity black hole [34]. Actually, though the mass of scalar field m 2 L 2 = m 2 e L 2 − 12 λ 2 L 2 in AdS case seems to be very different from m = 0 in flat case, it is obvious that the mass in AdS case m 2 = −2/L 2 + 12λ 2 /L 4 go to zero as L → ∞, which is the flat limit. Our numerical result will also support the analysis, We shall show soon that all properties for very large L in AdS case recover those in flat case as expected.
To proceed, we first show the potential function (10) in Fig. 3. We see that the potential near the horizon could be negative for certain parameter regime and the negative potential may lead to instability of the scalar field and a small scalar field perturbation may finally destabilizes the AdS black hole background. In the left plot without the GB coupling, the effective potential never becomes negative unless the scalar field is charged as addressed in [73,74]. In the middle plot with λ = 0.3, the potential well becomes deeper as L decreases, meaning that there exists a critical L * that the potential is negative enough for the scalar field emerging near the horizon. In the right plot with fixed L = 0.5, as λ increases, the potential becomes deeper. This behaviour is similar as that in AFBH case, however, by careful comparison, we find that for the same GB coupling, the well in the right plot is more deeper than that in AFBH case (see Fig. 1). This implies that comparing to the flat spacetime, the scalar field in AdS spacetime would be easier bounded near the horizon, saying with smaller GB coupling. Considering the scalar field near the horizon and the boundary condition (17) by choosing φ + = 0, we scan the parametric space (λ, L) to find the solution of equation (11). We find a critical line for the threshold of scalar field as shown on the left panel of Fig. 4. Below the critical line, the scalar field is zero which is trivial, while above the line the scalar field behaves as a non-trivial solution. We can extract two properties from the left panel of Fig. 4. Firstly, the threshold value of λ first increases as L becomes larger, and then tends to be a constant as L → ∞. The constant is around λ t 0.60234 which matches the aforementioned value in AFBH case (see the blue dot). This is reasonable because as L → ∞, the cosmological constant goes to zero which is the case we studied in the previous subsection. Secondly, comparing to AFBH case, the threshold value of λ in AdS case is smaller, which implies that weaker GB coupling could introduce spontaneous scalarization for spherical AdS black hole. This behaviour supports our previous analysis on the effective potential of Fig. 3. In the middle panel of Fig. 4, we present the scalar field near the event horizon φ h as a function of λ. For larger L, the value of λ for the emergence of non-zero φ h is larger, which is consistent with the observed on the left of Fig. 4. The scalar hair near the horizon emerges at the threshold coupling, and the value at the horizon increases as λ increases and then tends to be constant. The scalar charge of the black hole φ − shown in the right figure of Fig. 4 supports our results. However, when λ is small, φ h decreases as L increases, but when λ is large enough, φ h increases as L increases. Starting from the probe limit, we can obtain certain signal of spontaneous scalarization of the AdS background brought by GB coupling with the form f (φ) = 1 2 (1 − e −φ 2 ). Numerical methods are called for to handle the fully backreacted system and get the hairy solution, as was explored for various different f (φ) in [53].

C. Asymptotical dS black hole
In the dS case with Λ = 3/L 2 > 0, only the black hole with spherical horizon exists, so the redshift function in the solution (6) is g(r) = 1 − M r − r 2 L 2 . As r h = 1, the mass is M = 1 − 1 L 2 , so that L > 1 is required for the positive mass. Moreover, the cosmological horizon is located at r c = 1 2 (−1 + √ 4L 2 − 3) > r h = 1, which further requires L > √ 3 in our setup. Similar as the above subsections, we first discuss the validity of the probe approximation considering the scalar field as a test field in the dS spacetime. Following the analysis of [54], we consider the Klein-Gordon equation Eq.(11) in pure dS background with g(r) = 1 − r 2 L 2 for which the cosmological horizon is located at r c = L. By setting φ(r) = ϕ(r) r and x = r rc , the Klein-Gordon equation is reduced as where u = 3χ−1 2 , and where we definem 2 ef f = m 2 − 9λ 2 r 4 c in the second equality. We see that the Gauss-Bonnet coupling contributes to χ as a modification of that in [54].
The solution of (19) is where P u and Q u are Legendre functions and s is an arbitrary constant. It is obvious that the solution diverges at r = 0, and also the second term diverges at r = r c because of Q u r rc diverges. So the regularity at r = r c requires s = 0, then the solution near r c is while the solution at r r c behaves as The precondition of scalarization is that a scalar field must have a tachyonic behaviour somewhere in between the BH and cosmological horizon, which requiresm 2 ef f < 0 and χ > 0. But this implies that the scalar field could diverge at r r c . Therefore, the regularity of the scalar field at the cosmological horizon makes the tachyonic test field divergent at large distance, which breaks the probe approximation. Even though the probe limit may be valid between the event horizon and cosmological horizon, we can not check the signal of scalarization in probe limit with the same boundaries as working in flat and AdS cases. We could improve the numerical skills to involve the backreaction and do the comparison, and we hope to address this issue in the future.

IV. SIGNAL OF SCALARIZATION OF BLACK HOLES WITH NEGATIVE CUR-VATURE
In this section we will study the signal of scalarization if the background metric is a topological black hole (TBH). For the action the presence of a negative cosmological constant allows the existence of black holes with a topology R 2 ×Σ, where Σ is a two-dimensional manifold of constant curvature. The simplest solution of this kind, when Σ has negative constant curvature, reads where dσ 2 2 is the line element of Σ which is just (6) with k = −1, Λ = −3L −2 and M = 2Gµ. The line element dσ 2 2 is locally isomorphic to the hyperbolic manifold H 2 and Σ must be of the form where Γ is a freely acting discrete subgroup (i.e., without fixed points). The configurations (24) are asymptotically local AdS spacetimes and these black holes are known as TBHs [75]- [82]. It has been shown in [83] that the massless configurations where Σ has negative constant curvature are stable under gravitational perturbations and the stability of the topological black holes was discussed in [69,84]. In a series of papers [68,85,86] it was shown that there is evidence that a vacuum TBH goes over to a hairy configuration, the MTZ black hole [13], through a second-order phase transition. Therefore it is interesting to see if in the presence of the coupling between scalar field and the GB term, there is evidence of TBH scalarization.
In this case, the GB term (8) is evaluated as R 2 GB = 4 r 2 (g (r) 2 + (g(r) + 1)g (r)). To fulfill g(r h ) = 0, the black hole mass should satisfy M = r 3 h /L 2 − r h . The positivity of black hole mass forces 0 < L < r h , while the positivity of Hawking temperature requires L < √ 3r h . Thus for r h ≤ L < √ 3r h , the topological black hole have a positive Hawking temperature, but the black hole mass could be zero or negative, and the stability of the these massless and negative mass black hole has been investigated in [84]. It is noted that the topological black holes have interesting structures in their own right, here for security we focus on topological black hole with positive mass, such that we have 0 < L < 1 as we have set r h = 1. Then the effective potential is shown in Fig. 5. Without the GB coupling (see the left plot), the effective potential is always positive, so there is no signal of tachyonic instability to induce the scalarization, so that the TBH keeps stable. With λ = 0.3 (see the middle plot), the potential well becomes deeper as L decreases . When L is small enough, the potential well would be deep enough to allow the tachyonic instability to happen. In the right plot with fixed L = 0.5, as λ increases, the potential well becomes deeper. We note that the rule for negative curvature in AdS case extracted from Fig. 5 is qualitatively the same as that in the positive curvature case (see Fig. 3). But with the same parameters, the potential well for k = −1 is shallower than that in the case with k = 1. This means that the onset of spontaneous scalarization is easier to happen for AdS black hole with positive curvature than that with negative curvature. Again, to further fix the parameters for the threshold of the nontrivial scalar field, we need to call back the perturbation equation (11). Near the event horizon, the scalar field While in asymptotical infinity, the behavior is also (17). We shall also set φ + = 0 in the numeric.
The critical threshold parameters in (λ, L) space is shown on the left of Fig. 6 where in the green shadow region the scalar field in nontrivial. As L increases, the critical coupling also increases which is similar as that in the cases with k = 1. The curve tends to be divergent as L → 1 because there is no black hole solution when L ≥ 1 in our setup. By comparing with left of Fig. 6 and Fig. 4, we find that the scalar hair is more difficult to form in the negative curvature than in the positive case. This is consistent with the conclusion from the comparison of the potential function. We also present the scalar hair near the horizon and the scalar charge at infinity as the function of λ for different values of L in the middle and right of Fig. 6, respectively. With the increase of L, the threshold value λ becomes larger, which matches the observation on the left of Fig. 6. Moreover, as L becomes larger, the curves become lower and lower rapidly and the value of φ h (φ − ) decreases. Especially, when L goes closely to the unity, φ h becomes very small, which indicates that the scalarization of the background black hole becomes very difficult. This behaviour is consistent with that in positive curvature at small coupling (see the middle plot of Fig. 4) as well as that in zero curvature shown in next section.
We have observed that scalarizations in hyperbolic AdS black hole backgrounds depend on two factors, the coupling strength between the scalar field and the GB term and the cosmological constant through the radius L. For the chosen L, once the coupling parameter λ increases, the effective potential in Fig. 5 becomes deeper which can ignite tachyonic instability and start the hair formation more easily. When we have big enough L, in Fig. 6 we see that the required critical λ to start the scalarization becomes bigger. In the middle panel of Fig. 6, we find that for bigger L, the formed scalar hair is weaker even we have stronger coupling strength λ. The growth of L counteracts the effect of the λ and hinders the scalarization process.

V. SIGNAL OF SCALARIZATION OF BLACK HOLES WITH ZERO CURVA-TURE
AdS black holes with a flat space (k = 0) is described as toroidal black hole spacetimes. Thus the corresponding redshift function in the metric (6) is g(r) = − M r + r 2 L 2 (M = 1 L 2 as r h = 1) and the GB term is R 2 GB = 4 r 2 (g (r) 2 + g(r)g (r)). We plot the effective potential in Fig. 7. The potential function with λ = 0 does not show negative well which insures the stability of the original toroidal hole. When the coupling between the scalar field and GB term is turned on, for the chosen AdS radius L, the potential well near the horizon appears and becomes deeper as λ increases. For the same coupling strength λ, similar potential well behavior will appear when we decrease the AdS radius L. The deep potential well can bound the scalar field near the horizon and lead to the instability of the original toroidal black hole. This property is similar to those observed in the case with positive curvature (Fig. 3) and negative curvature (Fig. 5). Through careful comparisons, we find that for the same parameters, the potential well in zero curvature case is deeper than that in the negative curvature case, while shallower than that in the positive curvature background. This indicates that in ESGB theory, the horizon curvature influences the scalarization. For AdS black hole with toroidal horizon the scalar hair can be formed more easily than the hyperbolic case, but if compared with the spherical AdS black hole, the scalarization process in toroidal AdS black hole is still more difficult. To see the signal of the scalar field growing near the horizon, we solve the perturbation equation (11) with the boundary condition near the horizon The behavior is also (17) in asymptotical infinity, and we set φ + = 0 as usual in the numeric. The phase diagram in parameters (L, λ) space is shown in the left panel of Fig. 8 where below the blue line, the scalar field perturbation will trivially die out, but above this line the black hole scalarization will happen. The values on this line are critical threshold values to ignite the scalarization. When the AdS radius L → ∞, we see that the threshold coupling λ linearly approaches infinity, which means that the toroidal AdS black hole cannot be scalarized in such limits. In the limit L → ∞, the background returns to the flat spacetime and there is no black hole with the zero curvature.
With some fixed L, we study the scalar field near the horizon and the scalar charge at infinity as the function of λ, and the results are shown in the middle and right of Fig. 8 respectively. Note that for large GB coupling, the probe limit may be not reasonable and backreaction should be involved. But in some sense our result for small GB coupling is reliable and we see the signal of scalarization from the scalar field perturbation. In [59], we constructed the scalarized hairy black hole with planar horizon in AdS space by considering the full backreaction. There we have set L = 1 and the threshold value is λ 0.64 which is consistent with that in Fig. 8 via careful comparison.

VI. CLOSING REMARKS
In this paper, we discussed the signal of scalarization for black holes with different curvature topologies in various spacetimes. We considered the scalar field as a probe to the general relativity solution in Einstein-scalar-Gauss-Bonnet theory. Our discussion is based on the analysis of its effective potential and the onset of nontrivial solution to the perturbed scalar equation.
We studied the signal of black hole scalarization with positive curvature in asymptotical flat, AdS and dS spacetime. We first analyzed the validity of the probe approximation in various spacetimes. We found that in flat and AdS cases, a regular condition of test scalar field near event horizon could correspond to convergent (vanishing) condition at the asymptotical infinity, meaning the test field approximation could be valid; while in dS case, the regular condition near cosmological horizon leads to divergence at large distance, so that the test field approximation breaks down. We then compared the effect of asymptotical flat and AdS spacetime spacetime structures on the scalarization process. Without the GB coupling λ, the effective potential of the scalar field are alway positive and no well emerges. As we increased the coupling, the negative potential well would show and becomes deeper which may lead to tachyonic instability. With the same parameters, the potential well for AdS black hole is deeper than that for Schwarzschild black hole, meaning that the scalar hair around the spherical AdS black hole may be easier to form. We solved the perturbed scalar field equation in the parameter space (λ, L), and obtained the critical line above which the non trivial scalar hair can be found. Fig. 4 shows that the scalar cloud around AdS black hole could form easier than that around the Schwarzschild black hole. Moreover, in the limit L → ∞, the threshold coupling in AdS case approaches the value in Schwarzschild black hole from below as shown in the left of Fig. 4.
We also investigated the signal of AdS black hole scalarization process with different curvatures of the horizon. For AdS black hole with positive (k = 1), negative (k = −1) and zero (k = 0) curvatures, the properties of effective potential is qualitatively the same (see Fig. 3, Fig. 5 and Fig. 7, respectively), that's to say, no negative potential well appears with λ = 0 for instability. The potential well near horizon appears and becomes deep as λ (or L) increases (decreases), which could bound the scalar field near the horizon and lead to instability. However, the potential well for zero curvature is deeper than that for negative curvature but shallower than that for positive curvature. Notice that the effect of horizon curvature on the potential could be roughly analytically evaluated from the expression (10), i.e, V (r) = g(r) g (r) r + A r 2 + m 2 − λ 2 2 R 2 GB . With fixed L and λ, we can diagnose the potential as follows. (i) m 2 = −2/L 2 + 12λ 2 /L 4 for different horizon curvatures is the same in our setup. (ii) For the lowest mode, the topological term A satisfies A(k = 1) = A(k = 0) = 0 < A(k = −1) = 1/4 . (iii) Fixing horizon radius implies that the black hole mass satisfies M (k = 1) > M (k = 0) > M (k = −1), then it is easy to obtain that the group term T = g (r) r − λ 2 2 R 2 GB in the potential satisfies T (k = 1) < T (k = 0) < T (k = −1) < 0 for fixed r near horizon. (iv) The overall factor g(r) fulfills g(k = 1) > g(k = 0) > g(k = −1) > 0 for fixed r outside the horizon. Combining (i)-(iv), it is not difficult to reduce that the effect of horizon curvature on the negative potential well is V (k = 1) < V (k = 0) < V (k = −1) near horizon. It is noted that only the angular eigenvalues term A mentioned in point (ii) is truly the topological effect, and the terms mentioned in points (iii-iv) are actually geometrical effect. Thus, in this sense, the effect from the horizon curvature on the potential as well as the scalarization is reflected by the interplay between the topology and geometry.
The behavior of potential indicates that in ESGB theory, the scalar hair around AdS black hole with toroidal horizon may be easier to be formed than that with hyperbolic horizon but more difficult than that with spherical horizon. Moreover, we collected the critical threshold curves for different curvatures in Fig. 9 where the blue dot represents the threshold value of Schwarzschild black hole in the limit L → ∞. It is obvious that for the AdS black hole with fixed L, the scalar hair is the easiest to be formed around spherical horizon, then around the toroidal horizon and it is the most difficult to be bounded near hyperbolic horizon. We shall present two comments on Fig. 9. i) In the limit L → 0, the threshold coupling in all three cases is very small in which case the scalarization could easily realized. Moreover, we know from the effective potential that the GB coupling should be involved for the tachyonic instability and it further leads to the spontaneous scalarization. So we may argue that the threshold value would be tiny as L → 0, or the scalarization could cut off at certain tiny L. This direction with large Λ deserves further study. ii) In the limit L → ∞, the threshold coupling for zero or negative curvature are divergent, suggesting that the scalar field can not be bounded in these cases. This is understandable because as L → ∞, the spacetime become flat and only Schwarzschild black hole with positive curvature exists as background in probe limit. Moreover, it is interesting to involve the full backreaction of the scalar field into the geometry and verify our results.