Perturbations of stealth black holes in degenerate higher-order scalar-tensor theories

Among the scalar-tensor modified theories of gravity, degenerate higher-order scalar-tensor (DHOST) models could play a special role for dark energy while being consistent with current observations, notably those constraining the speed of gravitational waves. Schwarzschild–de Sitter black holes were shown to be exact solutions of a particular subclass of quadratic DHOST theories, while carrying a nontrivial scalar profile that linearly evolves in time and hence potentially providing exciting new phenomenological windows to explore this model. We investigate the physical perturbations about such black holes and find that the odd-parity tensor perturbations behave in a way indistinguishable to general relativity. On the other hand, the effective metric for the (even-parity) scalar perturbations is singular, indicating that those exact black hole solutions are infinitely strongly coupled and cannot be trusted within the regime of validity of the DHOSTeffective field theory. We show how this strong coupling result is generalizable to a whole class of solutions with arbitrary manifolds both for DHOST and Horndeski.


I. INTRODUCTION
The discovery of the cosmic accelerated expansion has motivated numerous studies on modifications of gravity in the infrared. The uniqueness of general relativity (GR) implies that any model of model gravity necessarily needs to involve additional degrees of freedom (d.o.f.) or a breaking of locality or Lorentz invariance. Including a scalar field and exploring scalar-tensor theories is therefore one of the most natural and minimalistic way to gravity beyond GR. Usually, scalar-tensor theories are considered as effective theories of more fundamental theories and provide a general framework to explain the observed cosmic acceleration phenomenologically. Under this consideration, many attempts have been made to construct the most generic consistent theory that propagates only one scalar d.o.f. while interacting with gravity. For example, Horndeski theories [1,2] are constructed as the most general scalar-tensor theory in four-dimensional spacetime yielding only second order field equations. In theories with multiple fields, higher-order field equations can still propagate a single d.o.f. if their Lagrangian is degenerate [3,4], and this argument was used to further extend scalartensor theories to a more general class of degenerate higherorder scalar-tensor (DHOST) theories in [5] (see [6][7][8][9][10][11][12] for related discussions and [13] for a review).
On the other hand, the direct detections of gravitational waves (GWs) from binary black hole and neutron star mergers made significant effects on our understanding of gravity. Particularly, the observation of GW170817 [14] together with its optical counterpart GRB170817A [15] constraints the speed difference between GWs and light (propagating on a cosmological background) down to 10 −15 [16], which drastically restricts the viable candidates of scalar-tensor theories [17][18][19][20][21], provided that such scalartensor theories are still valid up to the LIGO frequency [22].
The development of the GW astronomy also stimulates the studies on black hole solutions in scalar-tensor theories, among which black holes with nontrivial scalar profile are of particular interests. No-hair theorem has been proved for the shift-symmetric Horndeski theory [23] and for the shift-symmetric Gleyzes-Langlois-Piazza-Vernizzi (GLPV) theory [24], which is a subclass of DHOST theories. Such theorems state that if the coupling functions of the theories are regular, the static, spherically symmetric, and asymptotically flat black hole solutions with static scalar field must have the Schwarzschild metric and the constant scalar field. Hairy black holes are allowed if some of the conditions are violated. For example, by violating the regularity condition, there are asymptotically flat hairy black holes in the shift-symmetric Horndeski theory [25,26] and in the shift-symmetric GLPV theory [24]. However, these solutions usually present metrics different from GR black holes. Another example is the solution found in the shift-symmetric GLPV theories that do not have the canonical kinetic term [24], which also circumvents the no-hair theorem. Also see [27][28][29] for hairy black holes in Einstein-Scalar-Gauss-Bonnet theories.
Hairy black hole solutions can also be found if one allows for time-dependence. 1 For instance, within the context of shift-symmetric DHOST theories, a branch of hairy black holes is constructed by considering a linearly time-dependent scalar field profile φðt; rÞ ¼ qt þ ψðrÞ and a constant kinetic term X ¼ ∂ a φ∂ a φ. The linearly timedependent part of the scalar field can be thought as the background field that is responsible for the cosmic acceleration. Such solutions are investigated in [24,[33][34][35][36][37][38][39][40]. Moreover, these solutions can be "stealth," which means that the nontrivial scalar hair does not gravitate at the background level. Stealth black hole solutions were initially introduced in [41]. Recently, the stealth linearly time-dependent solutions have been studied in the context of quadratic DHOST theories [42,43]. Taking a bottom up approach, these studies identify the theories that possess stealth Schwarzschild-de Sitter (SdS) black hole solutions.
Motivated by this wealth of fruitful solutions, the stability of physical implications of those solutions was rapidly explored. If the scalar field couples directly to matter sources, one would expect binary systems to radiate scalar gravitational waves which would typically be in tension with observations unless a Vainshtein mechanism or other type of screening is implemented [44][45][46][47][48][49]. Moreover, in DHOST theories, even when the scalar mode may not a priori directly couple to external matter sources, since the physical propagating d.o.f. are mixed between the metric and the covariant scalar, implicit matter coupling typically appears in such theories which can then also lead to gravitational scalar radiation.
As for the stability of such solutions, for shift-symmetric Horndeski theories it was initially argued in [50] that they could be unstable against odd-parity perturbation; however, this statement was more recently revisited in [51,52]. For specific time slicings, it was argued in [50] that the positivity of either the kinetic or radial gradient energy would be violated in the vicinity of the horizon. However, that particular statement is gauge dependent and in [51,52] it was shown how there could exist a slicing for which the theory was stable everywhere. Even if the kinetic or the gradient term where to switch sign in all slicings, this would only signal the breakdown of the Horndeski effective field theory when either the kinetic or the gradient term becomes sufficiently small and the predictability of the Horndeski effective field theory therefore fails before the instability can even occur (see [53] for a related discussion).
In this study, we will explore the stability and validity (in the EFT sense) of the exact quadratic DHOST solution found in [42,43]. Odd perturbations have already been the subject of a very interesting analysis [40] (which appeared during the final stages of this work), and in the rest of this paper we shall investigate both the odd and even-parity perturbations about the stealth black hole solutions in the shift-symmetric quadratic DHOST theories. In agreement with [40], we show that such solutions are stable against odd-parity perturbations. Indeed, the odd-parity perturbations are exactly the same as that of GR black holes. The even-parity perturbations however differ due to the presence of the scalar field and while the diagonalized would-be tensor modes could a priori be stable, the diagonalized scalar mode sees a singular effective metric, indicating the stealth black hole solution in the shift-symmetric quadratic DHOST theories can unfortunately not be trusted.
To put the current findings in perspective, we highlight that the nature of this problem is more severe than that observed in [50] for Horndeski black holes. Indeed, fluctuations about the shift-symmetric Horndeski black hole solutions are well-behaved at sufficiently large distances and the theory only becomes unreliable close to the horizon (where either the kinetic or gradient term becomes small, before they would become negative). Such solutions can therefore still potentially provide relevant phenomenology away from the horizon. For the stealth black hole solutions in the shift-symmetric quadratic DHOST theories on the other hand, the effective scalar fluctuations are always everywhere and all the time degenerate and those solutions can therefore never and nowhere be trusted. These types of pathologies are very similar in nature to those observed about the exact static black solutions found in massive gravity [31], where it was shown that solutions that perfectly mimic GR black hole solutions, actually do so by effectively suppressing the graviton mass on those backgrounds and therefore making the additional d.o.f. present in massive gravity infinitely strongly coupled. The case of DHOST appears to be very analogous to the case of massive gravity: black hole solutions can only precisely be the same as in GR if the effect of the scalar field is effectively entirely suppressed, and therefore making this d.o.f. infinitely strongly coupled. Instead a perhaps more promising direction for black hole solutions in DHOST theories just like in massive gravity is the presence of a (small but nonetheless nonzero) time dependence of the metric. This time dependence would typically be governed by the graviton mass scale in massive gravity or the dark-energy scale in DHOST theories and could therefore imply a time-evolution only visible on timescales of the order of the age of the Universe. 2 The rest of the paper is organized as follows. In Sec. II, we first present the shift-symmetric quadratic DHOST theories and review the linear time-dependent black hole solutions. In Sec. III, we derive equations of both odd-and even-parity perturbations and highlight the presence of an (infinitely) strongly coupled d.o.f. The strong coupling issue is then generalized to a larger class of solutions with generic manifolds in DHOST theories including rotating black hole solutions in Sec. IV. The implications for a class of solutions in Horndeski is also highlighted. Section V is then conclusions and outlook.

II. STEALTH BLACK HOLES IN SHIFT-SYMMETRIC DHOST THEORIES
The Lagrangian of the shift-symmetric DHOST theory up to quadratic order is given by [5] In the wake of GW170817 [21], the requirement that the speed of GWs should be the same as light in a cosmological background imposes the following conditions (if one were to assume that the DHOST effective field theory remained sufficiently under control at LIGO frequency scale [22]), where here a prime denotes the derivative with respect to the argument X. In addition, it also requires A 3 ¼ 0 to prevent a rapid decay of GWs into the scalar field [55] (see also [56]). This condition as well as the condition A 1 ¼ 0 was not imposed in the recent analysis presented in [40]. Thus, the subclass of DHOST theories that will be considered in this paper is where we used A 4 ¼ 6F 02 =F. Note that we have also included external matter fields ψ i that only couple to the metric g μν . Even though at this level there is no direct coupling between φ and the external matter sources, we will still see in what follows that the physical scalar d.o.f. in this theory does directly couple (already at tree-level) to external sources. This is due to the nontrivial mixing of the physical d.o.f. in these types of degenerate field theories. For the back hole solutions, we assume the ansatz, where a bar refers to the background. Note that AðrÞ has nothing to do with the A i ðXÞ previously introduced in (2.1). We further require X to be a constant, Under the ingoing Eddington-Finkelstein coordinates (v, r) defined by dv ¼ dt þ dr=A, it can be shown that φ ≃ qv at the vicinity of the future event horizon [33,43]. Therefore, we will pick the branch with the "þ" sign, in which case φ is regular at the future event horizon. In this paper, we will mostly focus on SdS black holes, i.e., The theories also allow Schwarzschild black holes with linear time dependence, which can be obtained by sending Λ → 0 [43], and our analysis applies as well. Following the analysis presented in [42,43], we can see that the ansatz (2.6) solves the equation of motions in two following cases: where all terms are evaluated at X ¼ X 0 . For theories with A 3 ¼ 0, i.e., theories have no significant GWs to scalar field decay, the conditions of having black holes in those two cases are degenerate. We point out that the model considered at the moment is equivalent to Horndeski after field redefinitions; however, we perform the analysis on perturbations in the DHOST frame as it allows us to easily generalize it to a more generic class of DHOST theories in what follows.

III. BLACK HOLE PERTURBATIONS
We first start with the covariant equations of motion. Varying the action (2.5) with respect to the inverse metric g ab , we obtain the modified Einstein equation, where we define X a ≡ ∂ a X. Equation (3.1) involves terms with three derivatives acting on φ. Here T ab is the (conserved) stress-energy tensor associated with the matter fields ψ i . Conservation of T ab , i.e., D a T a b ¼ 0 sets the equation of motion for the matter fields ψ i . Varying with respect to φ yields the equation of motion for the scalar field To perturb the equations about the black hole solution, we write g ab ¼ḡ ab þ h ab , φ ¼φ þ δφ, and T ab ¼ 0 þ δT ab , where a subscript 0 or a bar refers to the background. 3 For convenience, we also keep X ¼ X 0 þ δX with δX ¼ 2φ a δφ a −φ aφb h ab where indices are raised and lowered with respect to the background SdS metricḡ ab . Working to first order in perturbations about the background, we then have the perturbed equations where we have used conditions (2.9) or (2.10) and To make further progress and properly separate out the relevant d.o.f., we consider the following linear combinations of the equations of motion: ab , leading to D a ½φ a ð−2F 0 ðF 00 respectively. Here we have defined the source tensor S ab as In the following, we shall decompose the perturbations based on their behaviors under parity transformations ðθ; ϕÞ → ðπ − θ; ϕ þ πÞ. This decomposition allows us to consider odd perturbation and even perturbation separately.

A. Odd sector
The odd-parity perturbation of the metric can be written as [57][58][59] h odd where h 0 , h 1 , h 2 are functions of (t, r), Y lm are the spherical Harmonics, a comma denotes the partial derivate, and Note that perturbation of the scalar field δφ is even under the parity transformation, hence does not couple with h odd ab and can be omitted in the odd perturbation equations. Moreover, explicate calculation shows that h odd ab φ a φ b ¼ 0 and hence δX vanishes in the odd sector. Together with the fact that δR odd ¼ 0, Eq. (3.3) simplifies to where F 0 plays the effective role of the Planck scale and we can therefore conclude that the odd-parity perturbation of the stealth SdS black holes considered here is identical to that in GR.
Relaxing the assumptions: The previous result relied on the assumption A 3 ¼ 0 (so as to prevent GWs from decaying into dark energy); however, we may wonder what the effects would be if some of those assumptions were relaxed. Actually, we find that even in the case where A 3 ≠ 0 and hence A 5 ≠ 0, the odd sector of GWs still behave identically as in GR. We can see this by perturbing the Lagrangian to quadratic order and using the fact thatB 3 ¼B 5 ¼ δ ð2Þ B 5 ¼ 0, then we can see that the presence of A 3 and A 5 leads to two extra terms, 9q 2 ΛĀ 3φ aφbḡcd h odd ac h odd bd andĀ 3 δ ð2Þ B 3 , which eventually cancel each other given the background solution. 4

B. Even sector
Next we turn to the even-parity perturbations. The metric perturbations can be written as [57][58][59] h even where again H 0 , H 1 , H 2 , H 0 , H 1 , K, and G are functions of (t, r), and ∇ θ;ϕ are covariant derivatives on the two-sphere of radius one. We now also have the scalar perturbation δφ ¼ Φðt; rÞY lm : ð3:14Þ In the following, we will sketch how to solve the even perturbation. An observation is that, given some initial conditions 5 for δφ and h even ab , one can directly solve for δX as a whole from Eq. (3.7) (for a specific matter source distribution set by δT ab ). The expression for δX can then be plugged into the effective source term S ab defined in Eq. (3.8), and this can then be used to solve for the remaining even-parity effective tensor mode in a very similar way as in GR. To see this work in practice, it is convenient to set a gauge and we do so differently 4 While this work was in progress, the interesting analysis of [40] appeared on the arXiv, also discussing the odd perturbation of stealth black holes but with further relaxing the assumptions to A 1 ≠ 0. Our results agree in the case of A 1 ¼ 0. 5 Note that not all components in h even ab are independent as most of them are related constraint equations as we shall see later. To solve for the system, we only need to set initial conditions for χ, _ χ, δX, and Φ. depending on whether we are dealing with the monopole, dipole, or higher multipoles.

Higher multipoles
For multipoles with l ≥ 2, we may fix the gauge by setting G ¼ K ¼ H 0 ¼ 0 and refer to Appendix A to see how coordinate transformations affect the even sector and check that this gauge can be chosen. The gauge fixing for monopole and dipole is different and will be treated separately below.
where Y stands for the spherical harmonics, commas denote partial derivatives, and the coefficients c i , d i , and f i are given in Appendix B. Then the equation of motion for χ can be written as where sðt; rÞ depends on δX and δT ab and is given by the same combination as the right-hand side of Eq. (3.17) with E L ab replaced by S ab . Using the following relations: S tt ; S rr ; S tr ∝ Y; ð3:19Þ the angular dependence fully drops out from the right-hand side of Eq. (3.18).
In particular, the last relation (3.22) can be seen as follows. For a, b ¼ θ, ϕ, terms in S ab proportional toφ aφb vanish, terms proportional toḡ ab lead to contribution proportional to r 2 Y, and terms involving covariant derivatives form the Laplacian operator in the two-sphere and therefore lead to contribution proportional to lðl þ 1ÞY.
We can now (in principle) solve (3.18) explicitly for χ and infer the other components in the metric perturbations by considering the following constraint equations: where J ¼ lðl þ 1Þ. Again, the above constraint equations are accompanied with some "source" terms on their righthand side, which are given by the same combination with E L ab replaced by S ab , hence proving the constraints for H 0 , H The above analysis shows that, for multipoles with l ≥ 2, there are two propagating d.o.f. in the even sector. With a trivial choice of P, Q, and F, we can get back to GR, in which case χ becomes the usual propagating d.o.f. in the even-parity sector. Therefore, we may think of the dynamical equations for the two d.o.f. in the even-parity sector as being Eqs. (3.7) and (3.18). We shall comment on this in what follows but first we look at the monopole and dipole.

Monopole
In the case of monopole, the contributions from H 0 , H 1 and G vanish identically. We therefore instead set the gauge K ¼ H 1 ¼ 0 (see Appendix A for confirmation that such a gauge can be fixed for the monopole). Then we find the following two constraints equations: which confirms the fact that there are no monopole tensor modes and the relevant dynamics of the physical scalar monopole is given once again by Eq. (3.7).

Dipole
In the case of dipole, h even ab only depends on K and G through the particular combination K − G, and thus we can which also confirms the fact that there are no dipole tensor modes and the relevant dynamics of the physical scalar dipole is given also once again by Eq. (3.7).

Dynamics of the scalar mode
Whether we were dealing with the monopole, the dipole, or the higher multipoles, we have shown that the relevant dynamics for the physical scalar mode is governed solely by Eq. (3.7). With this in mind, we shall therefore focus on that equation more closely and instead of any of the gauge choices we used previously, we shall now set a gauge so thatφ aφb h even ab ¼ 0 irrespectively of which multipole we are dealing with. We emphasize that this is only for convenience, but none of the results depends on that precise gauge choice. In this case, it is easy to see that Eq. (3.7) becomes solely an equation for δφ of the form 32Þ Given the background solution, we see that the scalar fluctuation δφ sees a singular effective metric g ab eff ∼φ aφb , which only ever has one nonvanishing eigenvalue. This implies that the physical (diagonalized) scalar fluctuations living on this exact black hole solution would be thus infinitely strongly coupled and the background solution cannot be trusted.
The strong coupling issue can be understood as follows. Consider, for instance, the EFT of scalar field π with operators that enter at an arbitrary scale Λ (that could potentially be as large as the Planck Scale), : ð3:33Þ In this model, just like in the theory we have considered in this paper, the effective metric of π becomes singular as B → 0. The strong coupling problem is always present when B → 0, but becomes manifest after performing the following redefinitions so that the kinetic term is canonical: and π → π=B 1=4 : ð3:34Þ In these new variables, the action (3.33) is ð3:35Þ As B → 0, it is now manifest that the toy model becomes infinitely strongly coupled and there is no sense in which perturbations theories can be trusted. This means that the background solution itself cannot be trusted for these types of models. This is precisely the issue that arises in the stealth black hole solutions considered here. While the example considered here is simply that of a scalar field, the argument applies to any nonlinear relativistic field theory as is the case for any gravitational theory. Being a gravitational theory, the DHOST models considered here involve an infinite number of nonlinear interactions. In the theory at hand, perturbing the action (2.5) around the black hole solution leads to nonlinear interactions including terms of the form P 00 0φ b δφ b δφ a δφ a , P 000 0φ b δφ b ðδφ a δφ a Þ 2 , etc., involving spatial derivatives acting on the scalar d.o.f. As we have seen in the toy model, a singular effective metric implies strong couplings of these nonlinear interactions and hence the breakdown of perturbation theory. One may avoid the strong coupling issue by setting all the derivatives of P, Q, F, and A i to vanish to all orders. This would then lead to a trivial scalar field theory which would not enjoy the stealth black hole solution considered here. Notice that it is the fundamental relativistic structure of the theory which enforces the existence of operators with spatial derivatives of the form (3.33). Had we been dealing with a nonrelativistic theory or what that manifestly breaks Lorentz invariance, the existence of spatial derivative in the nonlinear interactions could have been avoided and strong coupling would be prevented. However, for theories such as Horndeski, DHOST or other relativistic theories, singularity of the effective metric leads to strong coupling. We refer to [53,60] for further discussions on the issue of strong coupling in these types of theories.

IV. STRONG COUPLING ISSUES FOR GENERIC MANIFOLDS
Quadratic DHOST: Before concluding, it is worth pointing out that the strong coupling results derived in this paper hold beyond the SdS metric considered here and are actually generalizable to much more generic manifolds 6 and scenarios so long as X is constant on the background manifold, X ¼ X 0 ¼ const.
Indeed, consider the full quadratic DHOST theory (2.1) with generic shift-symmetric functions P, F, Q, and A i ðXÞ, then we can show that any background solution (denoted by the subscript 0 ) that satisfies the following properties on that particular background solution suffers from infinitely strong coupling and cannot be trusted, 8 > > > > > > > > < > > > > > > > > : This can be seen by considering ab are the field equation and the modified Einstein equation derived from the full quadratic DHOST theory (2.1), and E ð1Þ ¼ḡ ab E ð1Þ ab . Subject to condition (4.1), Eq. (4.2) is solely an equation for δφ with the highest order derivative proportional toφ aφb D a D b δφ. In particular, this implies that the rotating black hole solutions found in [61] suffer from the same issue, and fluctuations of the scalar d.o.f. about the rotating black hole found in [61] with finite stealth hair are infinitely strongly coupled (apart in the special case of spherical symmetry where the constraint A 3 ðX 0 Þ ¼ 0 is relaxed).
Note that those conditions do not impose to be dealing with theories where for instance A 3 , A 5 , and Q 0 vanish identically, and strong coupling would still be an issue even if say A 0 To avoid strong coupling, at least one of the constraints in (4.1) should be violated but note also that the conditions (4.1) are not the unique conditions under which the issue may arise, and violating one or several of the conditions in (4.1) does not necessarily ensure the absence of strong coupling issue. It is possible that strong coupling occurs on particular solutions even if the previous conditions are not satisfied, or that strong coupling arises instead for the tensor d.o.f. [11].
Further relaxing the assumptions: One may raise the question of what would occur if for instance A 3 ðX 0 Þ did not vanish precisely but was simply taken to be (extremely) small so as to prevent too much GW decay into dark energy on a particular solution of interest (of course if the EFT is not valid on those scales, the constraint on A 3 could potentially be relaxed further). If, for instance, A 3 ðX 0 Þ was considered to be small but nonvanishing, then the effective metric of the scalar d.o.f. about that solution could include four nonvanishing eigenvalues but the magnitude of those would be governed by the (extremely) small scale present in A 3 ðX 0 Þ and would also indicate strong coupling issues (low cutoff). For instance, if one were to consider perturbations about a spherically symmetric configuration, this would imply that the higher multipoles would not be suppressed as compared to lower multipoles. Second, since δφ does couple to the trace of external matter fields in generic DHOST theories, as can be seen from the righthand side of (3.32), any small test particle would lead to arbitrarily large emission of scalar waves. 7 Horndeski: This also applies to any Horndeski theories [1] that satisfies an equivalent set of conditions, independently of how symmetric (or not) the background manifold is. Consider a shift-symmetric Horndeski theory of the form with minimal coupling to external sources. Then for this theory, any solution on which X 0 is constant and for which G 0 f. about this would be solution is infinitely strongly coupled and the existence of such a solution could not be trusted. This result is independent of any details of the manifold considered and the symmetry of the solution. This is in addition to potential strong coupling issue that may occur for the tensor modes. 6 We wish to thank Hayato Motohashi for very useful discussion on this point. 7 Unless F 0 0 ¼ 0 in which case (2.9) or (2.10) would also imply Q 0 0 ¼ P 0 0 ¼ 0 and we would have A 4 ðX 0 Þ ¼ 0, then on that background P 0 would effectively play the role of a cosmological constant, the term proportional to Q 0 would be a total derivative and we would effectively just be dealing with GR and a scalar field minimally coupled to gravity.
In particular, we may point out that for appropriate choices of functions A i , the DHOST theory considered in (2.1) reduces to a special case of Horndeski [1], for which the static black hole solutions and their stability were explored in [62] and the strong coupling results remain valid in that particular subclass of Horndeski.
Indeed, following the analysis performed in [62] (which applied for static solutions q ¼ 0), with G 3;X ðX 0 Þ ¼ K ;X ðX 0 Þ ¼ 0, G 4 to be constant and imposing G 5 to vanish we find that the effective metric for the dynamical evend.o.f. is singular, see Appendix C, confirming a strong coupling issue for that particular limit of the Horndeski black Hole solutions. In that case, this strong coupling issue appears to be closely linked to the requirement that X be a constant at the background level.

V. OUTLOOKS
In this paper, we investigated the perturbation of linearly time-dependent stealth SdS black holes in shift-symmetric quadratic DHOST theories. We focus on the subclass of DHOST theories described by action (2.5), i.e., those theories that (1) predict unitary GW speed and (2) have no significant decay of GWs into the scalar fluctuations. The linearly time-dependent stealth SdS black holes exist if the functions in action (2.5) satisfy condition (2.9) or (2.10). As usual, we decomposed the perturbation based on their parity and derived the perturbation equations, respectively.
By deriving the perturbation equations, we find that the odd-parity perturbation is the same as that of GR black holes. Actually, this is the case even if the DHOST theories involve a nontrivial A 3 (in which case the stealth SdS black holes also exist). Since the background geometry is exactly SdS, the scalar perturbation does not couple with the oddparity metric perturbation. We also find that the even-parity perturbations are different from that in GR in general. The metric perturbation could be sourced by the perturbation of external matter field in a different way due to the presence of the scalar field. More concerning, we find that the scalar fluctuation sees a singular effective metric and hence suffers from a strong coupling problem. The black hole solution considered is therefore beyond the regime of validity of the DHOST effective field theory and cannot be trusted.
Finally, we point out that the issue of strong coupling derived in this paper is very generic to a large class of DHOST and Horndeski solutions. In particular, those issues apply to other rotating black hole solutions with scalar hair found in the literature and in some subclasses of Horndeski theories. We show that under a set of conditions DHOST and Horndeski solutions suffer the same scalar strong coupling issue irrespectively of the specific manifold and symmetry of the system.

ACKNOWLEDGMENTS
We would like to thank Christos Charmousis, Hayato Motohashi, and Huan Yang for helpful discussion and suggestions. C. d. R. would like thank the Perimeter Institute for Theoretical Physics for its hospitality during part of this work and for support from the Simons Emmy Noether program.

APPENDIX A: EVEN-PARITY COORDINATE TRANSFORMATIONS
In this appendix, we shall see the effect of an even-parity coordinate transformation so as to motivate our gauge chosen in the study of the even sector of Sec. III B. Recalling that the even-parity metric perturbations can be written as and the scalar perturbation as we now consider an infinitesimal coordinate transformation x a →x a ¼ x a þ ξ a with ξ a ¼ Tðt; rÞ; Rðt; rÞ; Θðt; rÞ∂ θ ; Θðt; rÞ∂ ϕ sin 2 θ Y lm ðθ; ϕÞ: