Physics Letters B 804 (2020) 135372

Contents lists available at ScienceDirect

Physics Letters B

www.elsevier.com/locate/physletb

Scalarization of horizonless reflecting : Neutral scalar fields non-minimally coupled to Maxwell fields

Yan Peng

School of Mathematical Sciences, Qufu Normal University, Qufu, Shandong 273165, China a r t i c l e i n f o a b s t r a c t

Article history: We analyze condensation behaviors of neutral scalar fields outside horizonless reflecting stars in the Received 25 January 2020 Einstein-Maxwell-scalar . It was known that minimally coupled neutral scalar fields cannot Received in revised form 26 February 2020 exist outside horizonless reflecting stars. In this work, we consider non-minimal couplings between Accepted 16 March 2020 scalar fields and Maxwell fields, which is included to aim to trigger formations of scalar hairs. We Available online 24 March 2020 analytically demonstrate that there is no hair theorem for small coupling parameters below a bound. For Editor: M. Trodden large coupling parameters above the bound, we numerically obtain regular scalar hairy configurations supported by horizonless reflecting stars. © 2020 The Author(s). Published by Elsevier B.V. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/). Funded by SCOAP3.

1. Introduction the context of scalar-tensor theories [51]. At present, lots of spon- taneous scalarization models were constructed in the background The recent observation of gravitational waves may provide a of black holes [52–60]. way to test the nature of astrophysical black holes [1–3]. In clas- As mentioned above, some candidate quantum-gravity models sical general relativity, one famous property of black holes is no suggested that quantum effects may prevent the formation of hori- hair theorem, which states that a nontrivial static scalar field can- zons and a reflecting wall may lay above the would-be horizon po- not exist in the exterior region of asymptotically flat black holes, sition [14–21]. Interestingly, it was found that neutral scalar field see references [4–11] and reviews [12,13]. However, some can- hairs cannot exist outside such horizonless reflecting stars (even didate quantum-gravity models suggested that, due to quantum the is charged) [26]. When considering scalar-Gauss-Bonnet effects [14–18], the classical absorbing horizon should be replaced couplings, we showed that the coupling can lead to the formation by a reflecting surface [19–25]. Interestingly, no hair behaviors also of neutral scalar field hairs in the background of horizonless re- appear for such horizonless reflecting stars [26–38]. In particular, flecting stars [61]. In fact, spontaneous scalarization is even for static scalar fields non-minimally coupled to the Ricci cur- a very universal property, which also can be induced by another vature, no hair theorem still holds in backgrounds of black holes type of non-minimal couplings between scalar fields and Maxwell and horizonless reflecting stars [39–44]. fields [62–67]. As a further step, it is very interesting to exam- Intriguingly, for static scalar fields non-minimally coupled to ine whether scalar-Maxwell couplings can trigger condensations of the Gauss-Bonnet invariant, scalar field hairs can exist outside neutral scalar fields outside horizonless reflecting stars. asymptotically flat black holes [45–48]. Spinning hairy black holes This work is organized as follows. We start by introducing a were also numerically obtained in the scalar-Gauss-Bonnet the- model with a neutral scalar field coupled to the Maxwell field in ory [49]. In addition, analytical formula of the scalar-Gauss-Bonnet the charged horizonless reflecting star . For small cou- coupling parameter was explored in [50]. These scalarization mod- pling parameters, the neutral scalar field cannot exist. In contrast, 2 for large coupling parameters, we get numerical solutions of scalar els are constructed by introducing an additional term f (ψ)RGB, 2 hairy horizonless reflecting stars. Main conclusions are presented where f (ψ) is a function of the scalar field ψ and RGB is the Gauss-Bonnet invariant. In the scalar-Gauss-Bonnet gravity, under in the last section. scalar perturbations, the bald black hole is thermodynamically un- stable and it may evolve into a hairy black hole [46,47]. This 2. Investigations on the coupling parameter between scalar fields intriguing mechanism of hair formations is usually called sponta- and Maxwell fields neous scalarization, which was found long ago for neutron stars in We take the Lagrange density with scalar fields non-minimally coupled to Maxwell fields in the asymptotically flat background. It E-mail address: [email protected]. is defined by the following expression [62–64] https://doi.org/10.1016/j.physletb.2020.135372 0370-2693/© 2020 The Author(s). Published by Elsevier B.V. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/). Funded by SCOAP3. 2 Y. Peng / Physics Letters B 804 (2020) 135372

ν 2 2 2 L = R −∇ ∇ν − μ ψ + f ()I. (1) 2M Q N(rpeak) = 1 − + > 0. (10) r r2 Here R is the scalar curvature.  is the static neutral scalar field peak peak ρσ with μ. f () is a function coupled to I = Fρσ F . In the With (9) and (10), we get the relation 2 linearized regime, there is I =−Q and the general coupling r4 α Q 2 l(l + 1) function can be expressed as f () = 1 − α2, where α is the − − μ2 0. (11) r4 r2 model parameter describing coupling strength [62–64]. In the limit peak peak → of α 0, it returns to the usual Einstein-Maxwell-scalar gravity. According to the inequality (11), we deduce a bound on the The scalar field differential equation is coupling parameter  f I 2 4 2 ∇ν∇ − 2 +  = μ r + l(l + 1)r 2r4 + l(l + 1)r2 ν μ  0. (2) peak peak μ s s 2 α 2 2 Q Q  The charged static spherically symmetric background is μ2(M + M2 − Q 2)4 + l(l + 1)(M + M2 − Q 2)2 2 > . (12) dr Q 2 ds2 =−N(r)dt2 + + r2(dθ 2 + sin2θdφ2). (3) N(r) If compact reflecting stars are surrounded with static neutral In the weak-field limit, the metric function N(r) is scalar hairs, the parameter α should be above the bound (12). In other words, we obtain a no hair theorem for small coupling pa- 2 2M Q rameters N(r) = 1 − + (4)   r r2 μ2(M + M2 − Q 2)4 + l(l + 1)(M + M2 − Q 2)2 with M and Q representing the star mass and star charge respec- α . (13) Q 2 tively. We point out that this background metric is valid on the condition αψ 2  1 It implies that neutral massive static exterior scalar fields usually We take the scalar field decomposition cannot exist in cases of large field mass, large spherical harmonic  index, large star mass or small star charge. In particular, for α = 0, imφ (r,θ,φ)= e Slm(θ)Rlm(r). (5) the no hair condition (13) always holds, which means that charged lm stars cannot support minimally coupled neutral scalar hairs. In the following, we numerically show that scalar hairs can be induced For simplicity, we label Rlm(r) as ψ(r). With relations (2), (3), by large non-minimal coupling parameters satisfying (12). 2 (5) and I =−Q , we derive the ordinary differential equation [33, r4 68,69] 3. Neutral scalar field hairs non-minimally coupled to Maxwell fields  2 2  2 N  α Q l(l + 1) μ ψ + ( + )ψ + ( − − )ψ = 0. (6) r N r4 N r2 N N We numerically solve the equation (6)together with boundary + conditions (7). Besides parameters rs, M, Q , l and α, we also need Here l is the spherical harmonic index and l(l 1) is the charac-  initial values of ψ(rs) and ψ (rs) to integrate the equation. The re- teristic eigenvalue of the angular scalar eigenfunction Slm(θ). flecting condition of (7)gives the value ψ(rs) = 0. According to We label rs as the radial coordinate of the star surface. Since  the symmetry ψ → kψ of equation (6), we firstly set ψ (r ) = 1 we focus on the without a horizon, the star sur- s without loss of generality. Since the equation (6)also satisfies the face is outside the gravitational radius, which can be expressed  symmetry r → γ r, μ → μ/γ , M → γ M, Q → γ Q , we use di- 2 2 as rs > M + M − Q . At the star surface, we take scalar re- mensionless parameters μrs, μM, μQ , l and α to describe the flecting surface boundary conditions ψ(rs) = 0. In the far region, system. For given values of μrs, μM, μQ and l, using standard the physical massive static scalar fields asymptotically behave as shooting methods, we search for the proper α with the vanishing →∞ ∼ 1 −μr ψ(r ) r e . So the scalar field satisfies bound-state condi- condition ψ(∞) = 0. The equation of motion of the scalar field is tions linear with respect to ψ and the solution is scale free. After getting  numerical solutions, we modify the boundary condition ψ (rs) = 1 = ∞ = 2 2 ψ(rs) 0,ψ() 0. (7) so that αψ  1 holds. In this work, we take very small αψ sat- − isfying αψ 2 < 10 7. According to boundary conditions (7), one concludes that the In the case of μrs = 2.7, μM = 1.5, μQ = 1.0 and l = 0, we function ψ(r) must possess (at least) one extremum point r = r peak choose various α to try to get the physical solution with ψ(∞) = = between the star surface r rs and spatial infinity. It can be a pos- 0. As shown by red curves in Fig. 1, if we choose α = 322, the itive maximum extremum point or a negative minimum extremum solution diverges quickly to be ∞. For green curves in Fig. 1, if we →− point. With the symmetry ψ ψ of equation (6), without loss choose a little larger value α = 323, then the solution decreases of generality, we can only study the case of positive maximum ex- to be −∞ in the larger r region. It turns out that α = 322 and tremum points. Then the scalar field around the extremum point α = 323 are not related to the physical scalar field solution with is characterized by [26] decaying behaviors at infinity. In fact, general mathematical solutions of equation (6)behave   ψ(r )>0,ψ(r ) = 0,ψ(r ) 0. (8) ≈ · 1 −μr + · 1 μr →∞ peak peak peak as ψ A r e B r e with r . The red line of Fig. 1 corresponds to B > 0 and the green line of Fig. 1 represents the (6) and (8)yield the relation case of B < 0. As the value B should change continuously with = α Q 2 l(l + 1) μ2 α, indicating the existence of a critical α corresponding to B 0. − − 0 for r= r . (9) For this critical α, physical scalar fields asymptotically decay as 4 2 peak r N r N N ∝ 1 −μr = = = ψ r e at infinity. With μrs 2.7, μM 1.5, μQ 1.0 and Since we concentrate on horizonless stars, the extremum point l = 0, we numerically obtain a discrete value α ≈ 322.083016, is outside the gravitational radius satisfying which corresponds to the solution satisfying ψ(∞) = 0. We plot Y. Peng / Physics Letters B 804 (2020) 135372 3

Table 1 The parameter α(l) with μM = 1.5, μQ = 1.0and various μrs .

μrs 2.62 2.66 2.70 2.74 2.78 α(l = 0) 258.102979 298.304665 322.083016 343.522473 364.215671 α(l = 1) 280.963138 322.681432 347.375787 369.638470 391.118639 α(l = 2) 326.150845 370.872697 397.388732 421.291276 444.340681

4. Conclusions

We investigated formations of neutral scalar field hairs outside asymptotically flat spherical horizonless reflecting stars. We con- sidered scalar fields non-minimally coupled to Maxwell fields. We showed that the coupling parameter plays an important role in scalar condensations. We analytically got a bound on the coupling parameter expressed in the form   μ2(M + M2 − Q 2)4 + l(l + 1)(M + M2 − Q 2)2 α , Q 2 where α is the coupling parameter, μ is the scalar field mass, M is the star mass, Q is the star charge and l is the spherical harmonic index. For α below this bound, no neutral scalar hair theorem holds. In contrast, for large above this bound, with shooting Fig. 1. We show behaviors of ψ(r) in cases of μrs = 2.7, μM = 1.5, μQ = 1.0, l = 0 α and different values of α. The red line corresponds to α = 322 and the green line methods, we obtained regular scalar hairy configurations with a is with α = 323. horizonless reflecting star in the center. We also examined effects of star radii, star mass, star charge and spherical harmonic index on condensations of neutral scalar fields.

Declaration of competing interest

The authors declare that they have no known competing finan- cial interests or personal relationships that could have appeared to influence the work reported in this paper.

Acknowledgements

This work was supported by the Shandong Provincial Natural Science Foundation of China under Grant No. ZR2018QA008. This work was also supported by a grant from Qufu Normal University of China under Grant No. xkjjc201906.

References

Fig. 2. We show the function ψ(r) in the case of μr = 2.7, μM = 1.5, μQ = 1.0, s [1] B.P. Abbott, et al., Virgo and LIGO Scientific Collaborations, Phys. Rev. Lett. 116 l = 0and α = 322.083016. (2016) 061102. [2] B.P. Abbott, et al., Virgo and LIGO Scientific Collaborations, Phys. Rev. Lett. 116 (2016) 241103. the physical solution with blue curves in Fig. 2, which asymptoti- [3] B.P. Abbott, et al., Virgo and LIGO Scientific Collaborations, Phys. Rev. X 6 (2016) 041015. cally approaches zero in the far region. For higher modes l ≥ 1, we [4] J.D. Bekenstein, Transcendence of the law of baryon-number conservation in showed physical solutions with discrete α in Fig. 3. Similarly, in black hole physics, Phys. Rev. Lett. 28 (1972) 452. other cases of black holes, scalar hairy configurations are also char- [5] J.E. Chase, Event horizons in static scalar-vacuum space-times, Commun. Math. acterized by discrete coupling parameters in the linearized regime Phys. 19 (1970) 276. [6] C. Teitelboim, Nonmeasurability of the baryon number of a black-hole, Lett. [47]. Nuovo Cimento 3 (1972) 326. Now we study how parameters μrs, μM, μQ and l can affect [7] R. Ruffini, J.A. Wheeler, Introducing the black hole, Phys. Today 24 (1971) 30. the discrete coupling parameter α, which corresponds to the de- [8] J.D. Bekenstein, Novel “no-scalar-hair” theorem for black holes, Phys. Rev. D caying scalar field. In Table 1, for μM = 1.5, μQ = 1.0 and fixed l, 51 (12) (1995) R6608. [9] D. Núñez, H. Quevedo, D. Sudarsky, Black holes have no short hair, Phys. Rev. we show effects of r on discrete . It can be seen that a larger μ s α Lett. 76 (1996) 571. radius μrs corresponds to a larger discrete α. With μrs = 2.7, [10] S. Hod, Hairy black holes and null circular geodesics, Phys. Rev. D 84 (2011) μQ = 1 and fixed l, according to data in Table 2, the discrete α 124030. decreases as we choose a larger star mass μM. In Table 3, we see [11] Yan Peng, Hair mass bound in the black hole with non-zero cosmological con- stants, Phys. Rev. D 98 (2018) 104041. that the discrete decreases as a function of the star charge Q . α μ [12] J.D. Bekenstein, Black hole hair: 25-years after, arXiv:gr-qc /9605059. Results in Table 3 also implies that α becomes smaller when we [13] Carlos A.R. Herdeiro, Eugen Radu, Asymptotically flat black holes with scalar μQ = Q = hair: a review, Int. J. Mod. Phys. D 24 (09) (2015) 1542014. choose a small Q M μM , which is qualitatively the same as [14] Pawel O. Mazur, Emil Mottola, Gravitational condensate stars: an alternative to cases of black holes expressed by analytical formula (17) in [64]. black holes, arXiv:gr-qc /0109035. From data in Tables 1, 2 and 3, we see that larger spherical har- [15] Cecilia B.M.H. Chirenti, Luciano Rezzolla, How to tell a from a black monic index l leads to a larger discrete coupling parameter α. hole, Class. 24 (2007) 4191–4206. 4 Y. Peng / Physics Letters B 804 (2020) 135372

Fig. 3. We Plot the function ψ(r) in the case of μrs = 2.7, μM = 1.5and μQ = 1.0. The left panel corresponds to l = 1and α = 347.375787. The right panel is with l = 2 and α = 397.388732. The two physical solutions with different l behave very similarly to each other.

Table 2 The parameter α(l) with μrs = 2.7, μQ = 1.0and various μM. μM 1.48 1.49 1.50 1.51 1.52 α(l = 0) 329.865858 326.322158 322.083016 316.754944 309.433789 α(l = 1) 355.357693 351.725147 347.375787 341.903046 334.372292 α(l = 2) 405.760364 401.953940 397.388732 391.632426 383.690319

Table 3 The parameter α(l) with μrs = 2.7, μM = 1.5and various μQ . μQ 0.92 0.96 1.00 1.04 1.08 α(l = 0) 351.506957 339.345528 322.083016 304.458791 287.540951 α(l = 1) 380.541786 366.508043 347.375787 328.019191 309.520394 α(l = 2) 437.965671 420.222598 397.388732 374.603527 352.976403

[16] Kostas Skenderis, Marika Taylor, The fuzzball proposal for black holes, Phys. [33] S. Hod, Charged reflecting stars supporting charged massive scalar field config- Rep. 467 (2008) 117–171. urations, Eur. Phys. J. C 78 (2017) 173. [17] Vitor Cardoso, Luís C.B. Crispino, Caio F.B. Macedo, Hirotada Okawa, Paolo Pani, [34] Yan Peng, Static scalar field condensation in regular asymptotically AdS reflect- Light rings as observational evidence for event horizons: long-lived modes, ing star backgrounds, Phys. Lett. B 782 (2018) 717–722. ergoregions and nonlinear instabilities of ultracompact objects, Phys. Rev. D [35] Yan Peng, On instabilities of scalar hairy regular compact reflecting stars, J. 90 (4) (2014) 044069. High Energy Phys. 10 (2018) 185. [18] Mehdi Saravani, Niayesh Afshordi, Robert B. Mann, Empty black holes, firewalls, [36] Yan Peng, Hair formation in the background of noncommutative reflecting and the origin of Bekenstein-Hawking entropy, Int. J. Mod. Phys. D 23 (13) stars, Nucl. Phys. B 938 (2019) 143–153. (2015) 1443007. [37] M. Khodaei, H. Mohseni Sadjadi, No skyrmion hair for stationary spherically [19] Vitor Cardoso, Paolo Pani, Testing the nature of dark compact objects: a status symmetric reflecting stars, Phys. Lett. B 797 (2019) 134922. report, Living Rev. Relativ. 22 (1) (2019) 4. [38] Bartłomiej Kiczek, Marek Rogatko, Ultra-compact spherically symmetric dark [20] Carlos Barceló, Raúl Carballo-Rubio, Luis J. Garay, Gravitational wave echoes matter charged star objects, J. Cosmol. Astropart. Phys. 1909 (09) (2019) 049. from macroscopic quantum gravity effects, J. High Energy Phys. 1705 (2017) [39] Avraham E. Mayo, Jacob D. Bekenstein, No hair for spherical black holes: 054. charged and nonminimally coupled scalar field with selfinteraction, Phys. Rev. [21] S. Hod, Stationary bound-state scalar configurations supported by rapidly- D 54 (1996) 5059–5069. spinning exotic compact objects, Phys. Lett. B 770 (2017) 186. [40] S. Hod, No nonminimally coupled massless scalar hair for spherically symmet- [22] Jahed Abedi, Hannah Dykaar, Niayesh Afshordi, Echoes from the abyss: ten- ric neutral black holes, Phys. Lett. B 771 (2017) 521–523. tative evidence for Planck-scale structure at black hole horizons, Phys. Rev. D [41] S. Hod, No hair for spherically symmetric neutral black holes: nonminimally 96 (8) (2017) 082004. coupled massive scalar fields, Phys. Rev. D 96 (12) (2017) 124037. [23] Bob Holdom, Jing Ren, Not quite a black hole, Phys. Rev. D 95 (8) (2017) [42] S. Hod, No nonminimally coupled massless scalar hair for spherically symmet- 084034. ric neutral reflecting stars, Phys. Rev. D 96 (2017) 024019. [24] Elisa Maggio, Paolo Pani, Valeria Ferrari, Exotic compact objects and how to [43] S. Hod, No hair for spherically symmetric neutral reflecting stars: nonminimally quench their ergoregion instability, Phys. Rev. D 96 (10) (2017) 104047. coupled massive scalar fields, Phys. Lett. B 773 (2017) 208–212. [25] Paolo Pani, Emanuele Berti, Vitor Cardoso, Yanbei Chen, Richard Norte, Grav- [44] Yan Peng, No scalar hair theorem for neutral Neumann stars: static massive itational wave signatures of the absence of an . I. Nonradial scalar fields nonminimally coupled to gravity, Nucl. Phys. B 947 (2019) 114730. oscillations of a thin-shell gravastar, Phys. Rev. D 80 (2009) 124047. [45] Thomas P. Sotiriou, Shuang-Yong Zhou, Black hole hair in generalized scalar- [26] S. Hod, No-scalar-hair theorem for spherically symmetric reflecting stars, Phys. tensor gravity, Phys. Rev. Lett. 112 (2014) 251102. Rev. D 94 (2016) 104073. [46] Daniela D. Doneva, Stoytcho S. Yazadjiev, New Gauss-Bonnet black holes with [27] Yan Peng, No hair theorem for massless scalar fields outside asymptotically flat curvature induced scalarization in the extended scalar-tensor theories, Phys. horizonless reflecting compact stars, Eur. Phys. J. C 79 (10) (2019) 850. Rev. Lett. 120 (2018) 131103. [28] Srijit Bhattacharjee, Sudipta Sarkar, No-hair theorems for a static and stationary [47] Hector O. Silva, Jeremy Sakstein, Leonardo Gualtieri, Thomas P. Sotiriou, reflecting star, Phys. Rev. D 95 (2017) 084027. Emanuele Berti, Spontaneous scalarization of black holes and compact stars [29] S. Hod, Charged massive scalar field configurations supported by a spherically from a Gauss-Bonnet coupling, Phys. Rev. Lett. 120 (2018) 131104. symmetric charged reflecting shell, Phys. Lett. B 763 (2016) 275. [48] G. Antoniou, A. Bakopoulos, P. Kanti, Evasion of no-hair theorems and novel [30] S. Hod, Marginally bound resonances of charged massive scalar fields in the black-hole solutions in Gauss-Bonnet theories, Phys. Rev. Lett. 120 (2018) background of a charged reflecting shell, Phys. Lett. B 768 (2017) 97–102. 131102. [31] Yan Peng, Bin Wang, Yunqi Liu, Scalar field condensation behaviors around re- [49] Pedro V.P. Cunha, Carlos A.R. Herdeiro, Eugen Radu, Spontaneously scalarised flecting shells in Anti-de Sitter , Eur. Phys. J. C 78 (8) (2018) 680. Kerr black holes, Phys. Rev. Lett. 123 (2019) 011101. [32] Yan Peng, Scalar field configurations supported by charged compact reflecting [50] S. Hod, Spontaneous scalarization of Gauss-Bonnet black holes: analytic treat- stars in a curved spacetime, Phys. Lett. B 780 (2018) 144–148. ment in the linearized regime, Phys. Rev. D 100 (2019) 064039. Y. Peng / Physics Letters B 804 (2020) 135372 5

[51] T. Damour, G. Esposito-Farese, Nonperturbative strong-field effects in tensor- [61] Yan Peng, Scalarization of compact stars in the scalar-Gauss-Bonnet gravity, J. scalar theories of gravitation, Phys. Rev. Lett. 70 (1993) 2220. High Energy Phys. 1912 (2019) 064. [52] Yves Brihaye, Carlos Herdeiro, Eugen Radu, The scalarised Schwarzschild-NUT [62] Carlos A.R. Herdeiro, Eugen Radu, Nicolas Sanchis-Gual, José A. Font, Sponta- spacetime, Phys. Lett. B 788 (2019) 295–301. neous scalarization of charged black holes, Phys. Rev. Lett. 121 (2018) 101102. [53] Yves Brihaye, Betti Hartmann, Charged scalar-tensor solitons and black holes [63] Pedro G.S. Fernandes, Carlos A.R. Herdeiro, Alexandre M. Pombo, Eugen Radu, with (approximate) Anti-de Sitter asymptotics, J. High Energy Phys. 1901 Nicolas Sanchis-Gual, Spontaneous scalarisation of charged black holes: cou- (2019) 142. pling dependence and dynamical features, Class. Quantum Gravity 36 (13) [54] Carlos A.R. Herdeiro, Eugen Radu, Black hole scalarisation from the breakdown (2019) 134002. of scale-invariance, Phys. Rev. D 99 (2019) 084039. [64] S. Hod, Spontaneous scalarization of charged Reissner-Nordström black holes: [55] Daniela D. Doneva, Stella Kiorpelidi, Petya G. Nedkova, Eleftherios Pa- analytic treatment along the existence line, Phys. Lett. B 798 (2019) 135025. pantonopoulos, Stoytcho S. Yazadjiev, Charged Gauss-Bonnet black holes with [65] Yun Soo Myung, De-Cheng Zou, Instability of Reissner-Nordström black hole in curvature induced scalarization in the extended scalar-tensor theories, Phys. Einstein-Maxwell-scalar theory, Eur. Phys. J. C 79 (3) (2019) 273. Rev. D 98 (2018) 104056. [66] Yun Soo Myung, De-Cheng Zou, Stability of scalarized charged black holes in [56] Hayato Motohashi, Shinji Mukohyama, Shape dependence of spontaneous the Einstein-Maxwell-Scalar theory, Eur. Phys. J. C 79 (8) (2019) 641. scalarization, Phys. Rev. D 99 (2019) 044030. [67] De-Cheng Zou, Yun Soo Myung, Scalar hairy black holes in Einstein-Maxwell- [57] Masato Minamitsuji, Taishi Ikeda, Scalarized black holes in the presence of the conformally coupled scalar theory, arXiv:1911.08062 [gr-qc]. coupling to Gauss-Bonnet gravity, Phys. Rev. D 99 (2019) 044017. [68] S. Hod, Onset of superradiant instabilities in rotating spacetimes of exotic com- [58] De-Cheng Zou, Yun Soo Myung, Scalarized charged black holes with scalar pact objects, J. High Energy Phys. 1706 (2017) 132. mass term, arXiv:1909 .11859 [gr-qc]. [69] S. Hod, Ultra-spinning exotic compact objects supporting static massless scalar [59] Caio F.B. Macedo, Jeremy Sakstein, Emanuele Berti, Leonardo Gualtieri, Hec- field configurations, Phys. Lett. B 774 (2017) 582. tor O. Silva, Thomas P. Sotiriou, Self-interactions and spontaneous black hole scalarization, Phys. Rev. D 99 (2019) 104041. [60] Jose Luis Blázquez-Salcedo, Daniela D. Doneva, Jutta Kunz, Stoytcho S. Yazad- jiev, Radial perturbations of the scalarized EGB black holes, Phys. Rev. D 98 (2018) 084011.