Causal structure of black holes in shiftsymmetric Horndeski theories
Abstract
In theories with derivative (self)interactions, the propagation of perturbations on nontrivial field configurations is determined by effective metrics. Generalized scalartensor theories belong in this class and this implies that the matter fields and gravitational perturbations do not necessarily experience the same causal structure. Motivated by this, we explore the causal structure of black holes as perceived by scalar fields in the Horndeski class. We consider linearized perturbations on a fixed background metric that describes a generic black hole. The effective metric that determines the propagation of these perturbations does not generally coincide with the background metric (to which matter fields couple minimally). Assuming that the metric and the scalar respect stationarity and that the surface gravity of the horizon is constant, we prove that Killing horizons of the background metric are always Killing horizons of the effective metric as well. Hence, scalar perturbations cannot escape the region that matter fields perceive as the interior of the black hole. This result does not depend on asymptotics but only on local considerations and does not make any reference to nohair theorems. We then demonstrate that, when one relaxes the stationarity assumption for the scalar, solutions where the horizons of the effective and the background metrics do not match can be found in the decoupling limit.
I Introduction
Black holes are among the simplest objects in the Universe and general relativity (GR) predicts that their defining feature, the event horizon, acts as a causal boundary for all fields. Black holes are the perfect probes of the strongfield regime, especially in the era of gravitational wave astronomy. This regime is also where alternative theories of gravity are expected to yield new predictions in the description of compact astrophysical objects. When new fundamental fields become a part of the picture, the simplicity of black holes in GR, as expressed by the nohair conjecture Israel (1967); Carter (1971); Robinson (1975), is often replaced by solutions with a more complex structure.
Even the simplest extensions of GR, scalartensor theories which postulate that a scalar field takes part in mediating the gravitational interaction, can introduce significant deviations Sotiriou (2015a, b). Consider the most general scalartensor theory that leads to second order equations upon direct variation, described by the Horndeski action Horndeski (1974); Deffayet et al. (2009); Kobayashi et al. (2011)
(1) 
where
(2)  
(3)  
(4)  
(5)  
(6) 
, , and . Nohair theorems do exist for broad subclasses of the Horndeski action under certain assumptions: when the scalar does not exhibit derivative selfcouplings, , , , and the configuration is stationary and asymptotically flat Chase (1970); Hawking (1972); Bekenstein (1995); Sotiriou and Faraoni (2012); when the theory respects shift symmetry, constant, , assuming staticity and spherical symmetry Hui and Nicolis (2013) or slow rotation Sotiriou and Zhou (2014a), and provided that the scalar does not couple to the GaussBonnet invariant, Sotiriou and Zhou (2014a, b). However, there is no known nohair theorem that covers the complete action. On the contrary, it is known to admit hairy solutions even for simple static, spherically symmetric and asymptotically flat configurations, e.g. Kanti et al. (1996); Yunes and Stein (2011); Sotiriou and Zhou (2014a, b); Silva et al. (2018); Doneva and Yazadjiev (2018); Antoniou et al. (2018). Moreover, in certain cases it is possible to obtain stationary hairy black hole spacetimes by relaxing the stationarity assumption for the scalar only Herdeiro and Radu (2014); Babichev and Charmousis (2014); Charmousis et al. (2014); Kobayashi and Tanahashi (2014).
The existence of hairy black hole solutions in theories where the scalar exhibits derivative selfinteractions might have implications for causality. Noncanonical kinetic terms allow for perturbations of the scalar field to propagate superluminally as long as the background field is nontrivial, i.e. has nonvanishing derivatives Babichev et al. (2008, 2017, 2018). Indeed, as a simple approximation, consider the Lagrangian (1) as describing a scalar field on a fixed curved background. Expanding to second order in scalar perturbations on top of a background , we get the second order Lagrangian for scalar perturbations as follows
(7) 
The explicit form of is presented in Appendix A for the special case where respects shift symmetry. The propagation of linear scalar perturbations is determined by effective metric . Hence, this metric, or its inverse , where , define causality for such excitations. The matter fields are assumed to couple minimally to and hence massless matter perturbations will follow null trajectories of that metric. For a generic background scalar field , does not coincide with . This suggests that the causal structure of black holes in certain classes of theories described by the Horndeski action should be particularly intriguing. If and are not conformally related then they have different null cones, so massless excitations that follow null geodesics of can be superluminal or subluminal.
One can take this a step further: for stationary metrics, event horizons will be Killing horizons, so one can ask under which conditions the Killing horizons of and coincide and what is the causal structure when they do not. These are precisely the questions that we explore below. We restrict ourselves to a setup where the background metric is fixed and consider only linear scalar perturbations, as in the discussion above. This has obvious limitations: it only probes the local counterpart of the causal structure as perceived by scalar excitations and ignores nonlinear effects. It also does not take into account the role of metric perturbations, which will generally also propagate along null geodesics of a different effective metric. Nonetheless, this setup makes calculation tractable and already provides very interesting insights into the causal structure as perceived by the scalar.
The causal structure of black holes in Horndeski theories has been studied before in Tanahashi and Ohashi (2017) and Minamitsuji (2015) (see Reall et al. (2014) and Izumi (2014) for similar approaches) using the method of characteristics. In Ref. Tanahashi and Ohashi (2017) both the metric and the scalar field were taken to be dynamical and it was claimed that if a surface is a Killing horizon of , it is also a characteristic surface for all degrees of freedom when the metric and the scalar are stationary. It was then suggested that this implies that effectively acts as a horizon for all excitations. This last statement entails some implicit assumptions not spelled out in Ref. Tanahashi and Ohashi (2017). Our approach follows a different path and it provides a rigorous proof about when Killing horizons for different excitations coincide. Though less general due to the decoupling approximation, it is also more transparent physically. Hence, it simplifies the interpretation of the results and highlights certain subtleties and assumptions in the characteristics approach.
The paper is organized as follows: in the next section, we provide a proof that the Killing horizons of the fixed background metric will also be Killing horizons of , provided that both the metric and the scalar respect stationarity. In Sec. III we relax the assumption of stationarity for the scalar and we show that, in the zero backreaction limit, it is rather straightforward to construct solutions where the horizons of these two metrics do not coincide. Though it is not clear if this feature will survive once backreaction is taken into account, our results certainly motivate further work in this direction. Sec. IV contains a discussion of our results and their potential implications.
Ii Nogo results for black holes with multiple horizons
As argued above and do not coincide for a generic background. In fact, interestingly, they do not coincide even for the trivial background, constant. However, we will prove below that the black hole Killing horizon perceived by the metric acts as a Killing horizon for metric as well. In order to do so, we assume that the black hole spacetime has a Killing vector field which is timelike outside the black hole region, is null and orthogonal to on the horizon and the background scalar field satisfies this symmetry, . Moreover, we assume those derivatives of ’s appearing in the expression for in Appendix A are finite.
Let us start with the easiest case, where the scalar background is trivial, constant. The effective metric for linear perturbations around this solution is given by (see Appendix A)
(8) 
where and are constants. Consider a spacetime where no extra matter field is present. Then, the variation of the action w.r.t. yields the vacuum Einstein equations . As a result, and are related by a conformal factor and their causal structure is the same.
In general, when the scalar field background is trivial, the metric satisfies GR field equations, and thus describes a GR black hole. The most general GR black hole solution is described by a KerrNewman metric. Though not obvious, we have verified by direct calculation that the Killing horizon of the KerrNewman black hole acts as a horizon for the metric as well.^{*}^{*}*Presumably this statement can be proven using the properties of the electromagnetic field, but we have not attempted that as it is covered by the more general proof given below. In the argument above, we have exploited the field equation for the metric to derive our result. In what follows, we provide a proof to show that the Killing horizon of is a Killing horizon of the effective metric without any reference to the field equation for (or to the equation for the background scalar field ).
is a composite metric constructed from the spacetime metric and the scalar field derivatives. Consequently, it satisfies stationarity, i.e. where is the Lie derivative along . The Killing horizon is a Killing horizon for the effective metric, iff
(9)  
(10) 
where means equality holds on and is a scalar. We ignore which corresponds to a degenerate perturbation metric, signalling the instability of background scalar configuration. Condition I states that the Killing horizon is a null surface w.r.t. the effective metric. Condition II requires that the Killing vector field be hypersruface orthogonal on  a necessary condition for the hypersurface to be a Killing horizon. Indeed, multiplying both sides of Eq. (10) with , we get
(11) 
is the normal to , thus this equation implies that is normal to .
In Appendix B, we prove both conditions. For condition I we only need to assume that the background metric and scalar configuration satisfy stationarity. To prove condition II without having to resort to global considerations, we need to make the additional assumption that the surface gravity of the horizon is constant. Overall, we have proved that the spacetime metric horizon is a Killing horizon for the perturbation metric. Hence scalar perturbations cannot escape the region that matter fields perceive as a black hole. Note that if the scalar perturbations are subluminal, it is conceivable that could have another horizon outside .
Clearly, condition II implies condition I. Multiplying both sides of the Eq. (11) with yields
(12) 
which means is null w.r.t. the perturbation metric on . We have chosen to consider it as separate condition to make better contact with the results of Ref. Tanahashi and Ohashi (2017). Condition I also identifies the characteristics for scalar perturbations, as is discussed in more detail in Appendix C. As shown there, if one applies the approach of Ref. Tanahashi and Ohashi (2017) to the decoupling limit and considers perturbation around a stationary background configuration, one does indeed obtain condition I.
The characteristics can be thought of as maximal speed of propagation surfaces and the boundary of the future of a given subset of the manifold is a characteristic. Nonetheless, characteristics are clearly not horizons in general (any null surface in Minkowski space is a characteristic for the wave equation) and condition I does not imply condition II without further assumptions. So, proving the former condition is not sufficient to claim that the Killing horizon of will also be a Killing horizon of given that the two metrics are distinct. One can argue around this technical obstruction. Condition I does show that scalar excitations can cross the Killing horizon of in a single direction only, which depends on the orientation of the characteristic. Hence, assuming that both and are asymptotically flat and identifying their asymptotic regions should suffice to claim that scalar perturbations are trapped by the Killing horizon of and cannot reach infinity, as claimed in Ref. Tanahashi and Ohashi (2017).
It is hard to imagine a physical solution that would not satisfy these assumptions and yet it is interesting that they are not sufficient to prove that is a Killing horizon of without further assuming that the surface gravity is constant. Carter has proven without resorting to field equations that Killing horizons in axisymmetric spacetimes that satisfy “– orthogonality” have constant surface gravity Carter (1973).^{†}^{†}†By “– orthogonality” one refers to the requirement that the – plane be orthogonal to a family of 2dimensional surfaces, which is a prerequisite for the line element to have as the only offdiagonal component. The is also trivially true in spherical symmetry. Indeed, in what follows we present a complementary proof that is valid in spherical symmetry and does not make use of this assumption.
It is also worth mentioning that, if one is willing to assume that the Killing vector is timelike w.r.t. in the exterior of , then there seems to be an alternative way to argue that condition II has to hold. The Killing flow cannot pierce , as the latter is a Killing horizon for . Hence, either has to be orthogonal to and null w.r.t. as well, or is has to reside in and be spacelike w.r.t. . However, the latter case is excluded by continuity if is timelike immediately outside .
Special case: Spherical symmetry
For a static spherically symmetric spacetime and scalar field configuration, we make use of the coordinates. In this coordinate system, both spacetime and effective metrics are diagonal and given by
(13)  
(14) 
where . The Killing vector is given by and the Killing horizon of the effective metric is where , or equivalently .
The Killing horizon of the effective metric cannot be located away from the metric Killing horizon at where . The argument goes as follows: away from the metric Killing horizon, all components of the spacetime metric and its inverse, in coordinates, are finite. This means that all components of the Riemann tensor and field derivatives (with upper or lower indices) are finite. is given by different combinations of the Riemann tensor and field derivatives contractions (see Appendix A), and thus it is finite and cannot satisfy the condition for the effective metric Killing horizon .
Iii Timedependent scalars and black holes with multiple horizons
The theorem in the previous section ruled out the possibility of having multiple horizons for most stationary configurations in a large class of shiftsymmetric Horndeski theories. However, as mentioned in the Introduction, static, spherically symmetric metric with timedependent scalar configurations are known to exist in shiftsymmetric Horndeski theories Babichev and Charmousis (2014); Charmousis et al. (2014); Kobayashi and Tanahashi (2014). In particular, if depends linearly on Killing time , i.e. in a suitable coordinate system and is a constant, then its derivatives respect staticity. Shift symmetry implies that is present in the equations only through its derivatives, and hence the mismatch between the symmetries of the metric and the scalar does not necessarily lead to an inconsistency. This motivates the study of causality and the potential existence of multiple horizons in solutions with stationary metrics and timedependent scalars. We will not attempt an exhaustive analysis here. Instead, we will just provide a simple example of a solution in which the Killing horizons of and do not match, as evidence that this issue deserves further investigation.
One might actually expect the solutions of Refs. Babichev and Charmousis (2014); Charmousis et al. (2014) to readily provide such an example. Unfortunately, this is not the case within the confines of the zerobackreaction limit that we are employing here. This becomes apparent when one briefly examines the hairy black hole solutions presented in Babichev and Charmousis (2014). The action used there is a combination of and
(15) 
which corresponds to and . Shift symmetry implies that the equation for the scalar field can be written as a current conservation equation
(16) 
where
(17) 
The field equation for is linear in . Hence, one considers scalar perturbations in a fixed hairy background, consistent with our assumptions in the previous sections, . If that background is a solution of vacuum Einstein’s equation with or without a cosmological constant, the and are conformally related and share their horizons.
Let us instead assume that is one of the static, spherically symmetric metrics of the hairy solution reported in Ref. Babichev and Charmousis (2014). The following ansätze were used
(18)  
(19) 
and all of the solutions found turn out to satisfy
(20) 
everywhere in the spacetime. However, Eq. (20) implies that everywhere and, hence, the effective metric for linear perturbations is degenerate throughout spacetime when backreaction is neglected. Clearly, in such solutions one cannot discuss in a meaningful way Killing horizons of in the zerobackreaction limit. The degeneracy of is clearly related to the choice of theory and the expression for the effective metric. Nonetheless, it does not actually imply that there is some general pathology of the solutions or the theory, as it appears only when one ignores the metric perturbations but still uses a background with a nontrivial scalar (i.e. when backreaction is neglected selectively at first order only.) See Refs. Ogawa et al. (2016); Babichev et al. (2017, 2018) for discussions on the linear stability of these solutions.
In order to provide an example of a solution with multiple horizons without taking into account backreaction, we turn to a different Lagrangian, known as kessence,
(21) 
where . Varying this action with respect to , we arrive at the equation of motion
(22) 
The causal structure of linear perturbations on top of a background field is governed by
(23) 
Inverting (23) yields
(24) 
which determines how scalar perturbations propagate through the spacetime.
Let us assume that describes a static spherically symmetric black hole spacetime with timetranslational Killing vector . The Killing horizon is located where becomes null according to the spacetime metric, i.e.
(25) 
Scalar perturbations can evade the metric horizon provided that remains timelike with respect to the effective metric on the Killing horizon. This happens when
(26)  
(27) 
on the metric horizon. The first condition expresses that the scalar field is time dependent and the second one means that the perturbations are superluminal. It is important to note that the above conditions are local, i.e. they need to be satisfied only in a neighborhood of the Killing horizon.
We now proceed to construct an explicit solution for the scalar field in a Schwarzschild background with mass , and demonstrate that the horizon does not coincide with that usual horizon. To do so, we make a specific choice of kessence theory (21), with
(28) 
where () corresponds to subluminal (superluminal) propagation. According to Eq. (26), time dependence in is necessary to have multiple horizons. We choose the following ansatz
(29) 
where is the ingoing EddingtonFinkelstein coordinate, and is the tortoise coordinate defined as . We further assume that our solution is continuously connected to constant as , and we expand as follows,
(30) 
We then solve the equations perturbatively in . At first order in , Eq. (22) reads
(31) 
whose solution, after imposing regularity on the horizon, is
(32) 
With this solution and ignoring the conformal factor which is irrelevant for the causal structure, the effective metric (24) yields
(33) 
We want to express metric (33) in EddingtonFinkelsteinlike coordinates, thus we perform the following coordinate transformation
(34) 
Defining , up to second order in the metric will read
(35) 
Starting from , outgoing radial null rays of the effective metric can escape to infinity. For , corresponding to the superluminal case, this region extends to both sides of the Killing horizon of the background metric .
Our solution clearly provides an example of an effective metric whose horizon does not coincide with that of the background metric, albeit is a much simplified setup where there is no backreaction and where the background is that of GR. It can be seen as having been obtained at the decoupling limit, where there is no backreaction at all orders in perturbation theory. The potential caveats are: (i) once backreaction is taken into account, the causal structure of the effective metric can change; (ii) it is not guaranteed that beyond decoupling, a static metric with a timedependent scalar actually exists. On the other hand, one can argue that if the scalar field configuration evolves very slowly in time compared to the characteristic timescale defined by the mass of the black hole, then treating the black hole as quasistationary might be a good approximation that justifies our treatment. In any case, our sole purpose here is to demonstrate that the possibility of having multiple horizons deserves further investigation when it comes to timedependent scalars.
Iv Discussion
The causal structure of black holes in Horndeski theories can deviate from the causal structure dictated by the spacetime metric. Due to derivative couplings of the scalar field and metric, the perturbation metrics are not conformal to the spacetime metric. Here we have focused for simplicity on the causal structure as seen by the linear perturbations of the scalar field on a fixed background metric.
The main result of this paper is the following: although null cones of the linear perturbation metric () and the spacetime metric () are generally different, we have shown that for stationary black holes and scalar fields, and provided that the surface gravity of the horizon is constant, a Killing horizon of is also a Killing horizon for . This means that the black hole region as perceived by matter fields, which are minimally coupled to , is a subset of or the same as the one defined by linear perturbations of the scalar field.
Our results agree with previous results in the literature that used the method of characteristics Tanahashi and Ohashi (2017) when there is overlap, but they go a bit further to prove that is actually a Killing horizon for . They also help interpret previous result in the intuitive physical picture of effective metrics and pin down certain subtleties of the characteristics approach. Remarkably, to rigorously complete the proof that is a Killing horizon for for a generic stationary configuration one needs the local assumption of constant surface gravity. Alternatively, one can argue that is a causal boundary of all excitations simply by being a characteristic, and hence is expected to be a Killing horizon for . This correspondence of assumptions begs the conjecture that Killing horizons have constant surface gravity for stationary, asymptotically flat black hole configurations in Horndeski theories. Though there is no known proof of that statement (or a theoryindependent proof), it is worth mentioning the following. Spherical black holes always possess this property and Carter has proven that axisymmetric black holes will also have this property if they satisfy “– orthogonality” Carter (1973). The latter is a necessary condition for having a line element with as the only offdiagonal component (the usual ansatz for stationary, axisymmetric solutions). Finally, the zeroth law of black hole thermodynamics is related to the fact that surface gravity is constant. Hence, finding black hole solutions that do not satisfy this property would indeed be remarkable.
It is important to clarify the complementarity between our proof and nohair theorems. In our general proof we do not impose that the background configuration satisfies any specific field equations and, hence, we did not need to make any reference to nohair theorems as it applies equally well to hairy solutions. As has been discussed above, for backgrounds with a trivial scalar configuration the proof becomes rather trivial. This might give the false impression that the physical significance of our result is limited. To clarify this point, we would like to emphasize that stationary hairy solutions do exist in shiftsymmetric Horndeski theories and that our proof does not rely on most of the assumptions employed in nohair theorems. In particular, we bring the following points to the reader’s attention:

Our proof uses only local properties of Killing horizons and makes no reference to the asymptotics. Hence, it applies to black holes with nonflat asymptotics or matter in their vicinity, which are not covered by nohair theorems.

Interactions between the scalar and other fields, including the metric, could provide a nontrivial configuration for the scalar field. However, as long as the new interaction terms do not spoil the structure of the perturbation metric our proof holds. That is, the interaction term only needs to respect shift symmetry upon linearization.
To elaborate on this last point, consider a theory that belongs to the class of theories given by the Lagrangian (1) and assume is a solution to that theory. Then add , where is a coupling constant, to the Lagrangian. The resulting theory is still in the class of shiftsymmetric Horndeski theories, since is a total derivative term in 4 dimensions. The new term in the Lagrangian adds to the scalar field’s equation of motion and does not vanish in black hole spacetimes. Consequently, the new theory must admit a hairy solution only. However, the new term in the scalar field equation of motion vanishes from the equation of motion for scalar linear perturbations on a fixed background. Hence, it does not change the structure of the effective metric and our proof applies to this class of theories as well.
Clearly, our central simplifying assumption to consider scalar perturbations on fixed metric backgrounds is a very strong limitation. It is pertinent to revisit the problem taking metric perturbations into account. The perturbative analysis presented in Ref. Babichev et al. (2018), though focused on the role of causality for stability rather than exploring the horizon structure, provides some useful background for extending our analysis.
We also discussed briefly the possibility of having different horizons in configurations where the scalar does not respect stationarity and we provided a simple example in the decoupling limit, in order to motivate further work. Another crucial piece of motivation to look further into this issue is the following. As has been pointed out in the context of Lorentzviolating theories, where solutions with nested horizons for different excitations are common Eling and Jacobson (2006); Barausse et al. (2011); Barausse and Sotiriou (2013), the region between two such horizons can be seen as an ergoregion. In that region, the slowest of the two excitations can carry negative Killing energy as it is behind its Killing horizon and at the same time the faster excitation can still escape. If the two modes interact, then there is a possibility of energy extraction Dubovsky and Sibiryakov (2006); Eling et al. (2007); Jacobson and Wall (2010). Assuming that such a process can take place Benkel et al. (2018), it could lead to thermodynamics conundrums Dubovsky and Sibiryakov (2006); Eling et al. (2007); Jacobson and Wall (2010), or perhaps just provide a natural decay channel from hairy to nonhairy solutions Benkel et al. (2018).
V Acknowledgements
We are grateful to Norihiro Tanahashi for enlightening discussions. T.P.S. acknowledges partial support from the STFC Consolidated Grant No. ST/P000703/1. M.S. is supported by the Royal Commission for the Exhibition of 1851. M.S. and T.P.S. would like to thank Instituto Superior Tecnico for its hospitality during the final stages of this work. We would also like to acknowledge networking support by the COST Action GWverse CA16104.
References
 Israel (1967) W. Israel, Phys. Rev. 164, 1776 (1967).
 Carter (1971) B. Carter, Phys. Rev. Lett. 26, 331 (1971).
 Robinson (1975) D. C. Robinson, Phys. Rev. Lett. 34, 905 (1975).
 Sotiriou (2015a) T. P. Sotiriou, Proceedings of the 7th Aegean Summer School : Beyond Einstein’s theory of gravity. Modifications of Einstein’s Theory of Gravity at Large Distances.: Paros, Greece, September 2328, 2013, Lect. Notes Phys. 892, 3 (2015a), arXiv:1404.2955 [grqc] .
 Sotiriou (2015b) T. P. Sotiriou, Class. Quant. Grav. 32, 214002 (2015b), arXiv:1505.00248 [grqc] .
 Horndeski (1974) G. W. Horndeski, Int. J. Theor. Phys. 10, 363 (1974).
 Deffayet et al. (2009) C. Deffayet, S. Deser, and G. EspositoFarèse, Phys. Rev. D80, 064015 (2009), arXiv:0906.1967 [grqc] .
 Kobayashi et al. (2011) T. Kobayashi, M. Yamaguchi, and J. Yokoyama, Prog. Theor. Phys. 126, 511 (2011), arXiv:1105.5723 [hepth] .
 Chase (1970) J. E. Chase, Communications in Mathematical Physics 19, 276 (1970).
 Hawking (1972) S. W. Hawking, Commun. Math. Phys. 25, 167 (1972).
 Bekenstein (1995) J. D. Bekenstein, Phys. Rev. D51, R6608 (1995).
 Sotiriou and Faraoni (2012) T. P. Sotiriou and V. Faraoni, Phys. Rev. Lett. 108, 081103 (2012), arXiv:1109.6324 [grqc] .
 Hui and Nicolis (2013) L. Hui and A. Nicolis, Phys. Rev. Lett. 110, 241104 (2013), arXiv:1202.1296 [hepth] .
 Sotiriou and Zhou (2014a) T. P. Sotiriou and S.Y. Zhou, Phys. Rev. Lett. 112, 251102 (2014a), arXiv:1312.3622 [grqc] .
 Sotiriou and Zhou (2014b) T. P. Sotiriou and S.Y. Zhou, Phys. Rev. D90, 124063 (2014b), arXiv:1408.1698 [grqc] .
 Kanti et al. (1996) P. Kanti, N. E. Mavromatos, J. Rizos, K. Tamvakis, and E. Winstanley, Phys. Rev. D54, 5049 (1996), arXiv:hepth/9511071 [hepth] .
 Yunes and Stein (2011) N. Yunes and L. C. Stein, Phys. Rev. D83, 104002 (2011), arXiv:1101.2921 [grqc] .
 Silva et al. (2018) H. O. Silva, J. Sakstein, L. Gualtieri, T. P. Sotiriou, and E. Berti, Phys. Rev. Lett. 120, 131104 (2018), arXiv:1711.02080 [grqc] .
 Doneva and Yazadjiev (2018) D. D. Doneva and S. S. Yazadjiev, Phys. Rev. Lett. 120, 131103 (2018), arXiv:1711.01187 [grqc] .
 Antoniou et al. (2018) G. Antoniou, A. Bakopoulos, and P. Kanti, Phys. Rev. Lett. 120, 131102 (2018), arXiv:1711.03390 [hepth] .
 Herdeiro and Radu (2014) C. A. R. Herdeiro and E. Radu, Phys. Rev. Lett. 112, 221101 (2014), arXiv:1403.2757 [grqc] .
 Babichev and Charmousis (2014) E. Babichev and C. Charmousis, JHEP 08, 106 (2014), arXiv:1312.3204 [grqc] .
 Charmousis et al. (2014) C. Charmousis, T. Kolyvaris, E. Papantonopoulos, and M. Tsoukalas, JHEP 07, 085 (2014), arXiv:1404.1024 [grqc] .
 Kobayashi and Tanahashi (2014) T. Kobayashi and N. Tanahashi, PTEP 2014, 073E02 (2014), arXiv:1403.4364 [grqc] .
 Babichev et al. (2008) E. Babichev, V. Mukhanov, and A. Vikman, JHEP 02, 101 (2008), arXiv:0708.0561 [hepth] .
 Babichev et al. (2017) E. Babichev, C. Charmousis, G. EspositoFarèse, and A. Lehébel, (2017), arXiv:1712.04398 [grqc] .
 Babichev et al. (2018) E. Babichev, C. Charmousis, G. EspositoFarèse, and A. Lehébel, (2018), arXiv:1803.11444 [grqc] .
 Tanahashi and Ohashi (2017) N. Tanahashi and S. Ohashi, Class. Quant. Grav. 34, 215003 (2017), arXiv:1704.02757 [hepth] .
 Minamitsuji (2015) M. Minamitsuji, Phys. Lett. B743, 272 (2015).
 Reall et al. (2014) H. Reall, N. Tanahashi, and B. Way, Class. Quant. Grav. 31, 205005 (2014), arXiv:1406.3379 [hepth] .
 Izumi (2014) K. Izumi, Phys. Rev. D90, 044037 (2014), arXiv:1406.0677 [grqc] .
 Carter (1973) B. Carter, in Proceedings, Ecole d’Eté de Physique Théorique: Les Astres Occlus: Les Houches, France, August, 1972 (1973) pp. 57–214.
 Ogawa et al. (2016) H. Ogawa, T. Kobayashi, and T. Suyama, Phys. Rev. D93, 064078 (2016), arXiv:1510.07400 [grqc] .
 Eling and Jacobson (2006) C. Eling and T. Jacobson, Class. Quant. Grav. 23, 5643 (2006), [Erratum: Class. Quant. Grav. 27, 049802 (2010)], arXiv:grqc/0604088 [grqc] .
 Barausse et al. (2011) E. Barausse, T. Jacobson, and T. P. Sotiriou, Phys. Rev. D 83, 124043 (2011), arXiv:1104.2889 [grqc] .
 Barausse and Sotiriou (2013) E. Barausse and T. P. Sotiriou, Class. Quant. Grav. 30, 244010 (2013), arXiv:1307.3359 [grqc] .
 Dubovsky and Sibiryakov (2006) S. L. Dubovsky and S. M. Sibiryakov, Phys. Lett. B638, 509 (2006), arXiv:hepth/0603158 [hepth] .
 Eling et al. (2007) C. Eling, B. Z. Foster, T. Jacobson, and A. C. Wall, Phys. Rev. D 75, 101502 (2007), arXiv:hepth/0702124 [hepth] .
 Jacobson and Wall (2010) T. Jacobson and A. C. Wall, Found. Phys. 40, 1076 (2010), arXiv:0804.2720 [hepth] .
 Benkel et al. (2018) R. Benkel, J. Bhattacharyya, J. Louko, D. Mattingly, and T. P. Sotiriou, (2018), arXiv:1803.01624 [grqc] .
 Gourgoulhon and Jaramillo (2006) E. Gourgoulhon and J. L. Jaramillo, Phys. Rept. 423, 159 (2006), arXiv:grqc/0503113 [grqc] .
Appendix A Effective metric in shiftsymmetric Horndeski
In this section we show the exact form of the effective metric defined in Eq. (7). We split
(36) 
where the subscript index refers to the contribution coming from , and use the simplifying notation and . We have explicitly
(37)  
(38)  
(39)  
(40) 
where
(41)  
(42)  
(43)  
(44)  
(45)  
(46)  
Appendix B Killing horizon of
Here, we prove Eqs. (9) and (10). For the proof of Eq. (9) we only use stationarity of the metric and the scalar field. In the proof of Eq. (10), we also assume that the Killing horizon has constant surface gravity. We use the notation introduced in Appendix A for simplicity, namely and .
b.1 Proof of condition I
The proof of condition I, , is rather straightforward, but it involves going through each term of in Appendix A to show it vanishes. For an interested reader, the following identities are helpful to reproduce the result:
where and are tangential spacelike directions on the horizon.
b.2 Proof of condition II
Let us introduce the following definition:
(47) 
i.e. two (co)vectors are equivalent if they differ by a multiple of .
Then, we can express Eq. (10) as
(48) 
Given the expression for (Appendix A), the following three relations
(49)  
(50)  
(51) 
where is any symmetric tensor satisfying , are sufficient for Eq. (48) to be satisfied.
In what follows, we prove the relations in Eqs. (4951). In the proof of (50) and (51) we assume that the surface gravity of the horizon is constant.
b.2.1 .
The scalar field satisfies stationarity, i.e. . Taking a derivative of this condition, we get
(52) 
We can use the expansion of on the horizon
(53) 
where is a null direction transverse to satisfying , is a spacelike direction orthogonal to and and is the surface gravity of . Equation (52) yields
(54) 
b.2.2 .
The Killing horizon has zero expansion, shear and vorticity. Thus applying Raychaudhuri equation on yields
(55) 
Moreover, for any tangent spacelike direction on
(56) 
Thus for a constant surface gravity of
(57) 
Eq. (55) implies that is tangent to and eq. (57) implies that this vector is orthogonal to all spacelike directions on . The only vector on satisfying these properties has to be a multiple of . Consequently,
(58) 
b.2.3 .
In order to prove Eq. (51), we consider a family of null hypersurfaces in an open neighborhood around with tangent null direction such that
(59) 
As a result,
(60) 
By the properties of null hypersurfaces,
(61) 
where is the null direction transverse to the null hypersurfaces satisfying , is the projection of onto 2dimensional spacelike submanifold of the null hypersurfaces (orthogonal to and ) and (Eq. (5.20) of Gourgoulhon and Jaramillo (2006)).
Taking the derivative of Eq. (61) and substituting in Eq. (60), we arrive at
(62)  
We show that each line in the above equation gives a contribution proportional to on .

Let us define . We prove by showing and where is any tangent spacelike direction on .
Note that and . Thus
(63) (64) which yields .
In order to prove , notice that
(65) As a result,
(66) From , it follows
(67) where is a scalar. Thus, we can express as
(68) where . Consequently,
(69) The index of has to be along a tangent direction to which implies . Moreover, and this proves .

If we define