In full Horndeski theories, we show that the static and spherically symmetric black hole (BH) solutions with a static scalar field φ whose kinetic term X is nonvanishing on the BH horizon are generically prone to ghost/Laplacian instabilities. We then search for asymptotically flat hairy BH solutions with a vanishing X on the horizon free from ghost/Laplacian instabilities. We show that models with regular coupling functions of φ and X result in no-hair Schwarzschild BHs in general. On the other hand, the presence of a coupling between the scalar field and the Gauss-Bonnet (GB) term R 2 GB , even with the coexistence of other regular coupling functions, leads to the realization of asymptotically flat hairy BH solutions without ghost/Laplacian instabilities. Finally, we find that hairy BH solutions in power-law F (R 2 GB ) gravity are plagued by ghost instabilities. These results imply that the GB coupling of the form ξ(φ)R 2 GB plays a prominent role for the existence of asymptotically flat hairy BH solutions free from ghost/Laplacian instabilities.
I. INTRODUCTIONGeneral Relativity (GR) has been tested by numerous experiments in the Solar System. While gravity can be well described by GR on the weak gravitational background in our local Universe [1], the dawn of gravitational-wave (GW) astronomy [2] and black hole (BH) shadow measurements [3] have started to allow us to probe the physics of extremely compact objects like BHs and neutron stars [4][5][6][7]. On the other hand, we also know that the Universe recently entered a phase of accelerated expansion [8,9]. While the cosmological constant is the simplest candidate for the source of this acceleration, i.e., dark energy, the theoretical value of vacuum energy mimicking the cosmological constant is enormously larger than the observed dark energy scale [10]. The cosmological constant has also been plagued by tensions of today's Hubble constant H 0 constrained from high-and low-redshift measurements [11,12]. These facts led to the question for the validity of GR on large distances relevant to today's cosmic acceleration. A simple and robust alternative to GR is provided by scalar-tensor theories possessing a scalar degree of freedom coupled to gravity [13].The most general class of scalar-tensor theories with second-order Euler-Lagrange equations of motion is called Horndeski theories [14-17]. 1 There have been numerous attempts for the theoretical construction of dark energy models compatible with observations [25][26][27][28][29][30]. Although such a new scalar degree of freedom potentially manifests itself in the Solar System, fifth forces mediated by the scalar field can be screened [31][32][33][34] by Vainshtein [35] or chameleon [36] mechanisms around a compact body on the weak gravitational background. In the vicinity of a BH, on the other hand, a nonvanishing charge of the scalar field, i.e., scalar hair, gives rise to a nontrivial field profile affecting the background geometry. This offers an interesting possibility for probing the possible deviation f...