Leading higher-derivative corrections to Kerr geometry

We compute the most general leading-order correction to Kerr solution when the Einstein-Hilbert action is supplemented with higher-derivative terms, including the possibility of dynamical couplings controlled by scalars. The model we present depends on five parameters and it contains, as particular cases, Einstein-dilaton-Gauss-Bonnet gravity, dynamical Chern-Simons gravity and the effective action coming from Heterotic Superstring theory. By solving the corrected field equations, we find the modified Kerr metric that describes rotating black holes in these theories. We express the solution as a series in the spin parameter $\chi$, and we show that including enough terms in the expansion we are able to describe black holes with large spin. For the computations in the text we use an expansion up to order $\chi^{14}$, which is accurate for $\chi<0.7$, but we provide as well a Mathematica notebook that computes the solution at any given order. We study several properties of the corrected black holes, such as geometry of the horizon, ergosphere, light rings and scalar hair. Some of the corrections violate parity, and we highlight in those cases plots of horizons and ergospheres without $\mathbb{Z}_{2}$ symmetry.


Introduction
General Relativity (GR) describes the gravitational interaction as the effect of spacetime curvature. Einstein's field equations, that rule the dynamics of the gravitational field, can be derived from the Einstein-Hilbert (EH) action which is essentially the simplest non-trivial covariant action one can write for the metric tensor. This beautiful theory has passed a large number of experimental tests -including the recent detection of gravitational waves coming from black hole and neutron star binaries [1][2][3][4][5][6]-and it is broadly accepted as the correct description of the gravitational interaction.

JHEP05(2019)189
However, there are good reasons to think that GR should be modified at high energies. One of these reasons is that GR is incompatible with quantum mechanics. Although we still lack a quantum theory of gravity, it is a common prediction of many quantum gravity candidates that the gravitational action (1.1) will be modified when the curvature is large enough. For instance, String Theory predicts the appearance of an infinite series of higher-derivative terms [7][8][9] correcting the Einstein-Hilbert action. The precise terms and the scale at which they appear depend on the scheme and on the compactification chosen. Nevertheless, whatever the modification of GR is, it should be possible to describe it following the rules of Effective Field Theory (EFT): we add to the action all the possible terms compatible with the symmetries of the theory and we group them following an increasing order of derivatives (or more generally, an increasing energy dimension). In the case of gravity, we would like to preserve diff. invariance and local Lorentz invariance, 1 and this means that the corrections take the form of a higher-curvature, or higher-derivative gravity [14]. A more general possibility -that we will also consider here -is to increase the degrees of freedom in the gravitational sector, by adding other fields that are not active at low energies [15].
Generically, the introduction of higher-derivative interactions means that Ricciflat metrics no longer solve the gravitational field equations. As a consequence, the Schwarzschild [16] and Kerr [17] metrics, that describe static and rotating black holes (BHs) in GR, are not solutions of the modified theories. One has to solve the modified field equations in order to determine the corrected black hole solutions, and it is an interesting task to look at the properties of these corrected geometries.
On general grounds, the higher-derivative corrections modify the gravitational interaction when the curvature is large, and they usually improve the UV behaviour of gravity [18]. The effect of the corrections will be drastic precisely in situations where GR fails, such as in the Big-Bang or black hole singularities, and it is expected that higher-derivative terms can resolve these divergencies [19][20][21][22][23][24][25][26][27]. However, the corrections can also significantly modify the properties of a black hole at the level of the horizon if its mass is small enough. For example, the divergence of Hawking temperature in the limit M → 0 in Einstein gravity (EG) black holes can be cured by higher-derivative interactions [28][29][30]. In this way, one learns about new high-energy phenomena that might be interpreted as the signature of a UV-complete theory of gravity.
Besides its intrinsic interest, there is another reason why studying higher-derivativecorrected black hole geometries is interesting: they can be used to obtain phenomenological implications of modified gravity. Thanks to the LIGO/VIRGO collaborations [31,32] and the Event Horizon Telescope [33], amongst other initiatives [34], it will be possible in the next years to test GR with an unprecedented accuracy, and to set bounds on possible modifications of this theory [35][36][37][38][39][40][41]. But in order to do so, we first need to derive observational signatures of modified gravity. In order to measure deviations from GR on astrophysical black holes, the corrections should appear at a scale of the order of few kilometers, which is roughly the radius of the horizon for those BHs. Although this seems to be an enor-

JHEP05(2019)189
In most of the literature, only few terms in this expansion are included, so the solutions are only valid for slowly-rotating black holes. However, astrophysical black holes -and in particular those created after the merging of a black hole binary [75] -can have relatively high spin. Moreover, some effects of rotation -such as the deformation of the black hole shadow [33,[76][77][78] -are barely observable when the spin is low, and other phenomena only happen for rapidly spinning black holes [79][80][81]. Although numerical solutions are not in principle limited to small values of the spin, analytic solutions are most useful for evident reasons. Hence, it would be interesting to provide analytic solutions valid for high-enough angular momentum. Finally, instead of having a large catalogue of alternative theories of gravity and their black hole solutions, it would be desirable to describe a minimal model that captures all the possible modifications of GR at a given order -probably, up to field redefinitions -and to characterize the black holes of that theory.
The preceding discussion motivates the three main objectives of the present work. First, to establish a general effective theory that can be used to study the leading-order higher-derivative corrections to Einstein gravity vacuum solutions. Second, to obtain the corrections to Kerr black hole in these theories, providing an analytic solution that is accurate for high enough values of the spin. And third, to study in detail some of the properties of these rotating black holes, such as the shape of horizon or the surface gravity, that have often been disregarded in the literature.
The paper is organized as follows. In section 2 we describe the models with higherderivative terms and non-minimally coupled scalars that we will use throughout the text, and we argue that they capture the most general corrections of this kind. In section 3 we describe the ansatz used in order to find the rotating black holes in the previous theories and we solve the equations performing a series expansion in the spin. For the computations in the text we use a series up to order O(χ 14 ), that is accurate up to χ ∼ 0.7. We provide as well a Mathematica notebook that computes the solution at any given order (it is contained in the supplementary material). In section 4 we study some properties of the modified Kerr black holes: horizon, ergosphere, photon rings and scalar hair, and we highlight the interesting geometry of black holes in parity-violating theories, that do not possess Z 2 symmetry. We summarize our findings in section 5, commenting on possible extensions and applications of the present work. Finally, there are several appendices with supplemental information.

Leading order effective theory
The most general diffeomorphism-invariant and locally Lorentz-invariant metric theory of gravity is given by an action of the form S = d 4 x |g|L (g µν , R µνρσ , ∇ α R µνρσ , ∇ α ∇ β R µνρσ , . . .) .
(2.1) This is, the most general Lagrangian for such theory will be an invariant formed from contractions and products of the metric, the Riemann tensor, and its derivatives. However, the theory above can be generalized by slightly relaxing some of the postulates. We may

JHEP05(2019)189
construct the Lagrangian using as well the dual Riemann tensor: These terms generically lead to violation of parity, hence the theory is not (locally) invariant under the full Lorentz group, but only under one of its connected components. However, we know that parity is not a symmetry of nature, so in principle there is no reason to discard terms constructed withR µναβ . In general, one expands this Lagrangian in terms containing increasing numbers of derivatives, being the first one the Einstein-Hilbert term R, with two derivatives. The rest of the terms can symbolically be written as Since this term contains 2n + p derivatives, it should be multiplied by a constant of dimensions of length 2n+p−2 with respect to the Einstein-Hilbert term. This is the length scale at which the higher-derivative terms modify the law of gravitation. When the curvature is much smaller than this length scale (||R µνρσ || −2 ), the effect of the higher-derivative terms can be treated as a perturbative correction, and terms with increasing number of derivatives become more and more irrelevant. Thus, it is an interesting exercise to obtain the most general theory that includes all the possible leading-order corrections. Here we summarize how we construct this theory, but we refer to the appendix A for the details. The first terms one may introduce in the action are quadratic in the curvature and hence they contain four derivatives. These terms would induce corrections in the metric tensor at order 2 , but in four dimensions it turns out that all of these terms either are topological or do not introduce corrections at all. Thus, the first corrections in a metric theory appear at order 4 and they are associated to six-derivative terms. As we show in appendix A, it turns out that, up to field redefinitions, there are only two inequivalent six-derivative curvature invariants, one of them parity-even and the other one parity-odd. However, one could consider a more general theory, allowing the coefficients of the higher-derivative terms to be dynamical i.e., controlled by scalars. This is actually a very natural possibility that is predicted, for instance, by String Theory [15]. In that case, some of the four-derivative terms do contribute to the equations and they also correct the metric at order 4 . For simplicity, we will restrict ourselves to massless scalars, but we will allow, in principle, to have an undetermined number of them. Within this large family of theories, it is possible to show that the most general leading correction to Einstein's theory is captured by the action where is the Gauss-Bonnet density and φ 1 , φ 2 are scalar fields. Besides the overall length scale , there are only five parameters: α 1 , α 2 , λ ev , λ odd and θ m . The parameter λ odd violates JHEP05(2019)189 parity, while the "mixing angle" θ m represents as well a sort of parity breaking phase. For θ m = 0, π (no mixing between scalars), φ 2 is actually a pseudoscalar and the quadratic sector is parity-invariant. For any other value (θ m = 0, π), parity is also violated by this sector. The theory (2.4) contains, as particular cases, some well-known models that have been frequently used in the literature. The case λ ev = λ odd = θ m = 0, α 2 = −α 1 = 1/8 corresponds to the prediction of String Theory, where the length scale of the corrections in that case is the string length 2 = 2 s ≡ α . As we show in the appendix B, the corresponding action can be obtained from direct compactification and truncation of the Heterotic superstring effective action at order α . In that case, φ 1 is identified with the dilaton, while φ 2 is the axion, which appears after dualization of the Kalb-Ramond 2-form. Another well-known possibility (which is also claimed to proceed from the low-energy limit of String Theory) is λ odd = λ ev = α 2 = 0, which corresponds to the Einsteindilaton-Gauss-Bonnet theory. Rotating black holes in EdGB gravity have been studied both numerically [65,66] and in the slowly-rotating limit [43,63,64]. The case θ m = π/2, which represents an extension of EdGB gravity, has also been considered [72] (note that this case only contains one dynamical scalar and violates parity). On the other hand, the case α 2 = 0 with the rest of couplings set to zero corresponds to dynamical Chern-Simons gravity, whose rotating black holes were studied in refs. [68][69][70] in the slowly-rotating approximation, while ref. [71] performs a non-perturbative numerical study. As for the cubic theories, the parityeven term (controlled by λ ev ) can be mapped (modulo field redefinitions) to the Einsteinian cubic gravity (ECG) term [58], for which static black hole solutions have been constructed non-perturbatively in the coupling [59,60]. Phenomenological signatures of static black holes in ECG have also been recently studied in [82,83], where a first bound on the coupling was provided, and the possibility to detect deviations from GR in gravitational lensing observations was discussed. Rotating black holes in ECG have not been studied so far. Lastly, to the best of our knowledge, the parity odd cubic term has never been used in the context of black hole solutions.
The theory (2.4) has been constructed following the sole requirement of diff. invariance, but there are some other constraints that could be imposed on physical grounds. For instance, if one wants to preserve parity, then one should set θ m = λ odd = 0. Nevertheless, we know that nature is not parity-invariant, so keeping these terms is not unreasonable. If one does not wish to include additional light degrees of freedom the scalars should be removed, which amounts to setting α 1 = α 2 = 0 (in that case the scalars are just not activated). On the other hand it is known that higher-derivative terms may break unitarity by introducing ghost modes -non normalizable states. In the case, for instance, of the Gauss-Bonnet term in (2.4) this problem does not exist since it produces second-order equations. The equations of the cubic terms do contain higher-order derivatives -namely of fourth order -, but the mass scale at which we expect the new modes to appear is This is simply telling us that Effective Field Theory works up to the scale ||R µνρσ || ∼ JHEP05(2019)189 −2 , which is something we already knew. Finally, it is also possible to study causality constraints [84]. In relation to this, the results in [42] impose a severe bound on the coupling constants λ ev 4 , λ odd 4 of the cubic terms. If one wants to observe any effects of higher-derivative corrections on astrophysical black holes, necessarily the corrections should appear at a scale of the order of few kilometers (otherwise the effect would be too small to be detected). Such large couplings are very unnatural, since the natural scale of (quantum) gravity should be Planck length. According to [42], these large couplings lead to violation of causality, that could only be restored by adding an infinite tower of higher-spin particles of mass ∼ −1 . Since, obviously, this is not observed, it was concluded that the couplings associated to the cubic terms should be of the order of Planck scale, hence those corrections would not be viable phenomenologically. 3 In any case, nothing prevents us from studying the effect of the cubic curvature terms on black holes, no matter the scale at which they appear. These corrections give us valuable information about the effects of modified gravity at high energies, and this is intrinsically interesting, even if those corrections are not viable on an observational basis.
If, for some reason, all the theories in the model (2.4) were discarded, then the leading correction to GR would be given by the quartic-curvature terms introduced in [74]. These terms modify the metric at order O( 6 ) hence they are subleading when the couplings in (2.4) are non-vanishing. Rotating black holes in those theories were recently studied in [73] up to order χ 4 in the spin. The methods that we present in this work could be applied to the quartic theories as well and could be used in order to extend some of the results in [73]. For instance, one might compute the solution for higher values of the spin or obtain the form of the horizon, as we do in section 4.1.

Equations of motion
Our goal is to compute the leading corrections to vacuum solutions of Einstein's theory. Thus, our starting point is a metric g (0) µν that satisfies vacuum Einstein's equations while the scalars φ 2 take a constant value that can be chosen to be zero without loss of generality. 4 But this field configuration is not a solution when we take into account the higher-derivative terms. First we note that the coupling between scalars and the curvature densities in the action (2.4) induce source terms in the scalar equations of motion, so that they will not be constant anymore. More precisely the first correction is of order 2 , and it satisfies (2.10) 3 Let us note that, according to [74], the conclusion might be different if one considers different UV completions of the effective theory from the one assumed in [42]. 4 The action 2.4 is invariant (up to a surface term) under constant shifts of the scalars.

JHEP05(2019)189
On the other hand, the modified Einstein equations, derived from the action (2.4), can be written as where we have passed all the corrections to the right-hand-side in the form of some energymomentum tensors, that read (2.12) and Since the scalars are of order O( 2 ), we can see that the leading correction to the metric associated to the scalar sector is of order O( 4 ), the same order at which cubic curvature terms come into play. Thus, we expand the metric as where g (4) µν is a perturbative correction. Now, taking into account that g (0) µν solves Einstein's equations, we get the value of the Einstein tensor to linear order in g (4) µν : where ∇ is the covariant derivative associated with the zeroth order metric, andĝ (4) µν is the trace-reversed metric perturbation Then,ĝ µν satisfies the equation

The corrected Kerr metric
After introducing the theory (2.4), here we present the rotating black hole ansatz that we will use in the rest of the text, and in section 3.1 we sketch how to solve the equations JHEP05(2019)189 of motion. From now on we set G = 1. Let us first consider Kerr's metric expressed in Boyer-Lindquist coordinates: Let us very briefly recall some of the properties of this metric.
• Being a solution of vacuum Einstein's equations, it is Ricci flat: R µν = 0.
• It is stationary and axisymmetric, with related Killing vectors ∂ t and ∂ φ respectively.
• It represents an asymptotically flat spacetime with total mass M and total angular momentum J = aM .
• When M > |a| the solution represents a black hole, whose (outer) horizon is placed at the largest radius r + where ∆ vanishes: Since Ricci flat metrics do not solve the modified Einstein's equations, the rotating black holes of the theory (2.4) will not be described by Kerr metric. The search for an appropriate metric ansatz that can be used to parametrize deviations from Kerr metric is a far from trivial problem that has been studied in the literature [85,86]. However, as long as the mass is much larger than the scale at which the higher-derivative terms appear, M , the deviation with respect General Relativity will be small -at least outside the horizon. In that case, we can build the rotating black hole solution of (2.4) as a perturbative correction over Kerr metric. Since we want to describe an stationary and axisymmetric solution, the corrected metric has to conserve the Killing vectors ∂ t and ∂ φ . On the other hand, we do not expect to "activate" additional components of the metric, so that the corrections appear in the already non-vanishing components. Taking into account these observations, we can write a general ansatz for the corrected Kerr metric where H 1,2,3,4,5 are functions of x = cos θ and ρ only, and they are assumed to be small |H i | 1. Note that we have introduced the coordinate ρ in order to distinguish it from the coordinate r in Kerr metric. We have also introduced the functions (3.5)

JHEP05(2019)189
However, the ansatz (3.4) is far too general, and it turns out that we can fix some of the functions H i by performing a change of coordinates. In particular, it can be shown that there exists a (infinitesimal) change of coordinates (ρ, x) → (ρ , x ) that preserves the form of the metric and for which H 5 = H 3 . Thus, we are free to choose H 3 = H 5 , and in that case, the metric reads (3.6) Note that we are choosing the coordinates x and ρ such that the form of the (ρ, x)metric is respected -up to a conformal factor -when the corrections are included. It is easy to see that this choice of coordinates has a crucial advantage: the horizon of the metric (3.6) will still be placed at the (first) point where ∆ vanishes: ρ + = M + √ M 2 − a 2 . If we were not careful enough choosing the coordinates, the description of the horizon could be very messy, and this is perhaps the reason why in previous studies the horizon of the corrected solutions has not been studied in depth.
We note that, whenever we consider the corrections, the coordinate ρ does not coincide asymptotically with the usual radial coordinate r. Advancing the results in next subsection, we get that the functions H i behave asymptotically as where h (k) i are constant coefficients. Then, we can see that the usual radial coordinate r that asymptotically measures the area of 2-spheres is related to ρ according to Using this coordinate, the asymptotic expansion of the metric (3.6) reads 4M a sin 2 θ r dtdφ (3.9) When we solve the equations, we see that we are free to fix the asymptotic values of the coefficients h correct normalization at infinity), and we want the parameters M and a to still represent the mass and the angular momentum per mass of the solution, so the asymptotic expansion should read (3.10) From this, we derive the asymptotic conditions that we have to impose on our solution: Apparently, the condition M h 1 = 0 is also required, but this is actually imposed by the field equations.

Solving the equations
Once we have found an appropriate ansatz for our metric, eq. (3.6), we have to solve the equations of the theory (2.4). The first step is to solve the equations for the scalars (2.9), (2.10), from where we obtain φ 1 and φ 2 at order O( 2 ). Using this result we determine the right-hand-side of (2.17), while in the left-hand-side we introduce the metric correction g (4) µν , which can be read from (3.6). In this way, we get a (complicated) system of equations for the functions H i , that we have to solve. Unfortunately, these equations (including the ones for the scalars) are very intricate and we are not able to obtain an exact solution. However, a possible strategy is to expand the solution in powers of the angular momentum a, assuming that it is a small parameter. In previous works [43,63,64,68,69,72], this method has been employed in order to obtain a few terms in the expansion, which yields an approximate solution for slowly rotating black holes. Nonetheless, if one includes enough terms in the expansion, the result should give a good approximation to the solution also for high values of the spin. One of the goals of this paper is precisely to provide a method that allows for the construction of the solution at arbitrarily high-orders in the spin. For simplicity, let us first introduce the dimensionless parameter that ranges from 0 to 1 in Kerr's solution, χ = 0 corresponding to static black holes and χ = 1 to extremal ones. 5 Then, we expand our unknown functions in a power series in χ 14)

JHEP05(2019)189
where we recall that all the functions depend on ρ and x. Then, the idea is to plug these expansions in (2.9), (2.10), (2.17), expand the equations in powers of χ, and solve them order by order. The equations satisfied by the n-th components are much simpler than the full equations, and we are indeed able to solve them analytically. These are second-order, linear, inhomogeneous, partial differential equations, so that the general solution can be expressed as the sum of a particular solution plus all the solutions of the homogeneous equation. In general, the "homogeneous part" of the solution represents infinitesimal changes of coordinates, and the physics is contained in the inhomogeneous part, which is the one sourced by the higher-derivative terms. So, we have to find the solution that captures the corrections but does not introduce unnecessary changes of coordinates. We observe that the appropriate solution can always be expressed as a polynomial in x and in 1/ρ. More precisely, we get 6 are constant coefficients and in each case the value of k max depends on n and p. When we solve the equations we also observe that all the terms in these series are determined except the constant ones: those with p = k = 0. However, those coefficients are fixed by the boundary conditions. In the case of the scalars, their value at infinity is arbitrary, so we can set it to zero for simplicity (this does not affect the rest of the solution) On the other hand, for the H i functions we take into account the relations (3.11) that we derived previously, which imply that In this way, the solution is completely determined. Since this process is systematic, we can easily program an algorithm that computes the series (3.14) at any (finite) order n. We provide a Mathematica notebook as supplementary material that does the job. Using this code, we have computed the solution up to order χ 14 . As we show in appendix D, this expansion provides a minimum accuracy of about 1% everywhere outside the horizon for χ = 0.7, and much higher for smaller χ. Thus, we have an analytic solution that works for relatively high values of χ, and we will exploit this fact in next section. Due to the length of the expressions, in appendix C we show the solution explicitly up to order χ 3 , but the full series up to order n = 14 is available in the Mathematica notebook. Before closing this section, we would like to clarify the following point. In the preceding scheme the corrections are expressed as a powers series in the spin, but we are taking JHEP05(2019)189 the zeroth-order solution to be the exact Kerr's metric, which is non-perturbative in the spin. Thus, for consistency sake, one should imagine that we also expand the zeroth-order solution in the spin up to the same order at which the corrections were computed. However, for evident reasons we do not do this explicitly. Thus, in the next section, we will write the formulas for several quantities as the result for Kerr's metric, exact in the spin, plus linear corrections, perturbative in the spin, but one should bear in mind that the zeroth-order result should also be expanded.

Properties of the corrected black hole
In this section we analyze some of the most relevant physical properties of the rotating black hole solutions we have found. We study the geometry of the horizon and of the ergosphere, light rings on the equatorial plane, and scalar hair.

Horizon
In order for the metric (3.6) to represent a black hole, we have to show that it contains an event horizon. We have argued that, with the choice of coordinates we have made, the horizon is defined by the equation ∆ = 0, whose roots are ρ = ρ ± , where (4.1) The largest root ρ + corresponds to the event horizon, while ρ = ρ − is in principle an inner horizon. 7 In this work we will only deal with the exterior solution ρ ≥ ρ + . Then, let us show that ρ = ρ + is indeed an event horizon. More precisely, we will show that it is a Killing horizon, i.e. a null hypersurface where the norm of a Killing vector vanishes. Let us first check that the hypersurface defined by ρ = ρ + is null. In order to do so, we consider the induced metric at some constant ρ, which is given by Then, we can see that when we evaluate at ρ = ρ + , the previous metric is singular, namely it has rank 2. Evaluating the determinant of the (t, φ)-metric at ρ + , we get where, for consistency with the perturbative approach, we have expanded linearly in the H i functions. When we expand the combination between brackets in powers of χ using the solution we found, we see that all the terms vanish. Thus, the determinant vanishes,

JHEP05(2019)189
which proves that this hypersurface is null. The next step is to show that there exists a Killing vector whose norm vanishes at ρ = ρ + . Such vector is a linear combination of the two Killing vectors ∂ t and ∂ φ : for some constant Ω H . One can check that the only possible choice of Ω H for which ξ is null at ρ + is which represents the angular velocity at the horizon. It is then clear that the norm of the vector ξ vanishes at ρ = ρ + , since where in the last step we have used (4.4). However, the crucial point here is whether Ω H , given by (4.6), is constant. A priori, this quantity could well depend on x, in whose case ξ would not be a Killing vector, and therefore ρ = ρ + would not be a Killing horizon. Nevertheless, expanding this quantity in powers of χ we do find that it is constant (see (4.8) below), a fact that provides a very strong check on the validity of our results. Thus, we have shown that ρ = ρ + is a Killing horizon, and hence it should correspond to the event horizon of the black hole.
We can now evaluate the angular velocity in order to study deviations with respect to Kerr's solution. A useful way to express it is the following, where the first term is the value in Kerr black hole and we made explicit the linear corrections related to the different terms in the action. It turns out that the parity breaking terms do not contribute to this quantity -nor to many others, as we will see. The dimensionless coefficients ∆Ω H depend on the spin, and the first terms in the χ-expansion read 8 The profile of these coefficients is shown in figure 1. This plot was done using the expansion up to order χ 15 , which provides an accurate result up to χ = 0.7. Interestingly, we observe that the correction related to α 1 increases the angular velocity, while the one related to JHEP05(2019)189 α 2 decreases it. The one associated to λ ev can have either effect, since the sign of λ ev is in principle arbitrary. We observe that the effect of these terms is larger for smaller masses: the quantity that controls how relevant the corrections are is 4 /M 4 times the corresponding coupling. They become of order 1 when M ∼ , which marks the limit of validity of the perturbative approach.
Surface gravity. At this stage, the natural step is to compute the surface gravity κ, defined by the relation that the Killing vector (4.5) must satisfy on the horizon. The computation of κ is not straightforward because the coordinates we are using are singular at the horizon. A possibility in order in order to circumvent this problem consists in working in Eddington-Finkelstein coordinates, that cover the horizon. However, there exist a number of alternative methods that can be used in order to obtain the surface gravity even if the coordinates are not well-behaved. Here we will follow a trick proposed in [87]. First, let us rewrite (4.12) as where we made use of the Killing property ∇ (µ ξ ν) = 0. Then, let us focus on the left hand side of the equation. The norm ξ 2 is a function of x and ρ, so that ∂ µ ξ 2 only has nonvanishing µ = x, ρ components. However, one can explicitly check that lim ρ→ρ + ∂ x ξ 2 = 0, hence the only non-vanishing component is µ = ρ, and it is given by 14) JHEP05 (2019)189 where, as usual, we are expanding linearly in the H i functions. On the other hand, since ξ is normal to the horizon, we must have ξ µ = C δ µ ρ for some constant C. Of course, this is not true in general: one should imagine that the previous formula holds only on the horizon, where the coordinate ρ is singular. The exact factor C is computed by taking the norm ξ 2 = C 2 g ρρ and evaluating at the horizon, so that we get (4.15) Then, we can plug (4.14) and (4.15) into (4.13) to find (4.16) Finally, evaluating this expression on the solution and expanding order by order in χ we find Again, we observe that parity-breaking terms do not modify this quantity. In addition, the fact that we obtain a constant surface gravity is another strong check of our solution, since this is a general property that any Killing horizon must satisfy. The profile of these coefficients as functions of χ is shown in figure 2, using an expansion up to order χ 14 . We see that both quadratic curvature terms controlled by α 1 and α 2 increase the surface gravity, with the difference that the α 2 correction vanishes for static black holes. On the other hand, the contribution from λ ev has a different sign depending on χ. For χ < 0.5 the surface gravity is greater than in Kerr black hole, while for χ > 0.5 it is lower, or viceversa, depending on the sign of λ ev . Another aspect that we can mention is that these contributions do not seem to be vanishing when χ → 1. This would imply that κ χ=1 = 0, hence the solutions with χ = 1 would not be extremal. In fact, there is a priori no reason to expect that the relation JHEP05(2019)189 Figure 2. Variation of the surface gravity ∆κ (i) due to every correction, as a function of χ. We can observe that contributions coming from curvature-squared terms always increase the temperature, since ∆κ (1) and ∆κ (2) are positive and also the coefficients multiplying them. The contribution from λ ev has different sign depending on χ.
between mass and angular momentum of extremal black holes is preserved when higherderivative corrections are taken into account. 9 Our results in figure 2 suggest precisely this; extremality will be reached for a value of χ slightly different from 1. However, let us note that the series expansion in χ breaks down for χ = 1, so the perturbative approach is not reliable in order to analyze the corrections to the extremal Kerr solution. One would need to compute directly the corrections to extremal rotating black holes in order to confirm whether the extremality condition is indeed corrected.
Horizon geometry. Let us finally study the size and shape of the horizon, which will be affected by the corrections. The induced metric at the horizon is First, we can find the area, which is given by the integral Computing the integral order by order in χ, we can write the area as In the case of charged black holes, it is known that higher-curvature corrections modify the relation between mass and charges, see e.g. [88,89]. where every contribution ∆A (i) depends on χ, and the first terms read In figure 3 we show the profile of these quantities as functions of χ, using the expansion up to order χ 14 . We observe that the quadratic corrections always reduce the area (except α 2 in the static case, that does not contribute). On the other hand, the cubic even correction reduces or increases the area depending on whether λ ev > 0 or λ ev < 0, respectively. So far, we have not observed the effect of the parity-breaking corrections -they do not contribute either to the area, the surface gravity or the angular velocity of the black hole. This is expected on general grounds since these corrections contain only odd powers of x, and it is easy to see that the contribution, for instance, to the area, must vanish. Nevertheless, these terms do change the geometry and they will affect the shape of the horizon. Indeed, these parity-breaking corrections break the Z 2 symmetry of the solution, i.e. the reflection symmetry on the equatorial plane x → −x. It is expected that this loss of symmetry is manifest in the form of the horizon.

JHEP05(2019)189
In order to visualize the event horizon, we perform an isometric embedding of it in 3-dimensional Euclidean space E 3 . In terms of Cartesian coordinates (x 1 , x 2 , x 3 ), we can parametrize the most general axisymmetric surface as where f (x) and g(x) are some functions that must be determined by imposing that the

JHEP05(2019)189
induced metric on the surface, given by coincides with (4.21). We get immediately that these functions are given by However, it can happen that the solution does not exist if the argument of the square root in the integral becomes negative. In that case, the horizon cannot be embedded completely in E 3 . It turns out that this only happens for quite large values of χ (around χ ∼ 0.9), and for the values we are considering here, the complete horizon can be embedded. As usual, we expand the expressions (4.29) and (4.30) linearly on H i and at the desired order in χ and we obtain explicit formulas for f and g that we do not reproduce here for a sake of clarity. Now we can use the result to visualize the horizon. In figure 4 we show the horizon for parity-preserving theories. We fix the mass to some constant value and χ = 0.65 and we compare the horizon of Kerr black hole with the one in the corrected solutions for different values of the couplings. In this way, we can observe clearly the change in size and in shape of the horizon. As we already noted, both α 1 and α 2 reduce the area, but it turns out that they deform the horizon in different ways: α 1 squashes it while α 2 squeezes it. We also show the deformation corresponding to the "stringy" prediction α 1 = α 2 . In that case we observe that the effect of both terms together is to make to horizon rounder than in Einstein gravity. As for the cubic even correction, it mainly changes the size of the black hole while its shape is almost unaffected.
In figure 5 we present the horizon in the parity-breaking theories (characterized by the two parameters θ m and λ odd ). In the top row we plot the horizon for a fixed choice of higher-order couplings and for various masses, keeping χ = 0.65 constant. The visualization is clearer in this way since these corrections do not change the area. In addition, we can see that for large M the horizon has almost the same form as in EG, but as we decrease the mass the corrections become relevant and it is deformed. We observe in this case that the Z 2 symmetry is manifestly broken. Due to exotic form of these horizons we include as well a 3D plot in which we can appreciate them better. Very recently other works have described black hole solutions that do not possess Z 2 symmetry [73,90]. However, to the best of our knowledge, these are the first plots of black hole horizons without Z 2 symmetry in purely gravitational theories.

Ergosphere
Another important surface of rotating black holes is the ergosphere, which marks the limit in which an object can remain static outside the black hole. When g tt < 0, there are no timelike trajectories with constant (ρ, x, φ), so the ergosphere is identified by the condition JHEP05(2019)189 g tt = 0, which for the metric (3.6) can be written as This equation determines the value of the "ergosphere radius" ρ erg . Unlike the horizon radius ρ + , that does not receive corrections due to the clever choice of coordinates, the ergosphere radius is modified with respect to its value in Kerr metric. We may express the JHEP05(2019)189 corrections to ρ erg as where the first term represents the result in Einstein gravity and we have to determine the value of the coefficients ∆ρ (i) . Plugging this into (4.31), we find these coefficients, whose first terms in the χ-expansion are shown in eq. (E.1). In this case, we do get a non-vanishing contribution from the parity-breaking terms, though this is not directly relevant, since ρ erg has no physical meaning by itself. However, an interesting property that we note by looking

JHEP05(2019)189
at (E.1) is that all the corrections to ρ erg vanish at x = ±1, corresponding to the north and south poles of the ergosphere. There is a nice interpretation of this fact: the ergosphere and the horizon overlap at the poles. Indeed, the horizon radius ρ + does not have corrections, and the zeroth-order value of the ergosphere radius ρ In order to study the geometry of the ergosphere, we can compute the induced metric for ρ = ρ erg (x) at a constant time t = t 0 , which reads Using the value of ρ erg that we have found yields a complicated expression that we omit here for clarity sake. The most useful way to visualize the geometric properties of the ergosphere is to find an isometric embedding of the previous metric in Euclidean space, as we have just done with the horizon. The embedding is shown in figure 6 for paritypreserving theories, and in figure 7 for parity-breaking ones. In the former case, we plot the ergosphere for a fixed mass and χ = 0.65, and for different values of the couplings, including the GR result. We observe that the corrections change the size and shape of the ergosphere. The quadratic terms α 1 and α 2 both reduce the area of the ergosphere, while the cubic even term reduces its size for λ ev > 0, and increases it for λ ev < 0. The characteristic conical singularity at the poles of the ergosphere is also considerable affected by some corrections. In particular, we see that α 2 and λ ev < 0 have the effect of making the cone less sharp. In the top row of figure 7 we show instead the embedding of the ergosphere for several values of the mass, while keeping the couplings and χ = 0.65 constant. This helps the visualization since parity-breaking interactions do not change the area of the ergosphere. As the mass decreases, the effect of the corrections becomes relevant and we observe, as in the case of the horizon, that the ergosphere does not possess Z 2 symmetry. This is more explicit for the cubic odd correction λ odd that deforms the ergosphere giving it a characteristic "trompo" shape. The effect of Z 2 symmetry breaking is less obvious for the θ m deformation, but nevertheless it can still be observed. To the best of our knowledge, these are the first examples of ergospheres without Z 2 symmetry.

Photon rings
Another aspect of the modified Kerr black holes we would like to explore is their geodesics. The analysis of geodesics is necessary in order to obtain some observable quantities, such as the form of the black hole shadow [37]. However, a detailed analysis of geodesics will require of an independent study due to their intricate character. 10 For that reason, here we JHEP05(2019)189 consider only a special type of geodesics that are particularly interesting: circular orbits (ρ = constant) for light rays at the equatorial plane, i.e. at x = 0, known as the photon rings or light rings of the black hole. However, an appropriate question that we must answer first is whether there are geodesics contained in the equatorial plane at all. In the case of Kerr metric, the reason of their existence is the reflection symmetry x → −x, but we have seen that in our black holes this symmetry does not exist if we include parity- M 4 λ odd = 0.6 (right). In the latter case we observe clearly that the Z 2 symmetry is broken and the ergosphere acquires a characteristic "trompo" shape. The effect is more subtle in the left picture, but the Z 2 symmetry is also broken. breaking terms. In fact, in those solutions there is no equatorial plane! Therefore, we should not expect the existence of geodesics contained in the plane x = 0 if we include those corrections. In order to understand this better, let us examine the geodesic equations:

JHEP05(2019)189
whereẋ µ = dx µ dλ and λ parametrizes the curve x µ (λ). Let us evaluate these equations for a trajectory withρ = 0 and x = 0, which represents a circular orbit. We find that the µ = x JHEP05(2019)189 component of (4.35) reads In order for the truncation x = 0 to be consistent, the left-hand-side should vanish independently of the value ofṫ andφ. This does not always happens, and the reason is precisely the presence of parity-breaking interactions, controlled by λ odd and sin θ m . Note that all the terms appearing in (4.36) are proportional to ∂ x H i | x=0 . When the theory preserves parity, the solution possesses Z 2 symmetry and the functions H i only contain even powers of x, so that ∂ x H i | x=0 = 0. On the contrary, the parity-breaking terms introduce odd powers of x in the H i functions -in particular terms linear in x -implying that ∂ x H i | x=0 = 0. Thus, in such theories setting x = 0 is not consistent and there are no orbits contained in the plane x = 0 (probably there are no orbits contained in a plane at all, besides the radial geodesics at the axes x = ±1). 11 For simplicity, from now on we set the parity-violating parameters to zero, λ odd = θ m = 0, so that we can study equatorial geodesics. However, we believe that studying the geodesics in those theories is an interesting problem that should be addressed elsewhere. Let us then focus on the remaining equations. When they are evaluated onρ = 0 and x = 0, the µ = t and µ = φ components of the geodesic equations (4.35) tell us thaṫ t = const andφ = const and, consequently, the angular velocity ω ≡ dφ/dt is also constant. On the other hand, the component µ = ρ gives an equation for ω: where the Christoffel symbols are shown in eq. (E.6). Finally, we take into account that for massless particles we have g µνẋ µẋν = 0, that gives the following equation Now, using the equations (4.37) and (4.38) we can solve for ρ and ω. We get two solutions that we can express as the result in Einstein gravity plus corrections: (1) (4.40) 11 In ref. [73], rotating black holes were studied in the presence of quartic-curvature corrections, including a parity-violating combination, and it was stated that this interaction does not have effects on equatorial geodesics. Apparently, the analysis of geodesics in that paper missed the fact that those geodesics are not permitted if the parity-violating term is activated. On the other hand, that analysis should be perfectly valid if the problematic term is removed.

JHEP05(2019)189
The "+" solution corresponds to the prograde photon ring (the photons rotate in the same direction as the black hole), while the "−" solution represents the retrograde photon ring. We reproduce here the values of the coefficients ∆ω (i) ± expanded up to order χ 7 in the spin while the coefficients ∆ρ (i) ph± are shown in eq. (E.9) of the appendix. However, ρ ph± is a meaningless quantity, since ρ does not have a direct interpretation as a radius. What we should really consider as the radius of the light rings is Since the light ring (more precisely, the photon sphere) determines the shape of the black hole shadow, this quantity give us information about the deformation of the shadow (near the equator) due to the corrections. On the other hand, ω ± is also an interesting quantity, since it is related to the time-scale of the response of the black hole when it is perturbed. In fact, there is a known quantitative relation between the orbital frequency of the light ring and the quasinormal frequencies of static black holes in the eikonal limit [92,93]. Although the relation probably does not extend to the rotating case, we do expect that ω ± captures qualitatively the (real) frequencies of the first quasinormal modes. Hence, we can use ω ± in order to perform a first estimation of the effects of the corrections on the black hole quasinormal frequencies.
In figure 8 we show the frequencies ω ± and the radius R ± for several values of the higher-order couplings and we compare them to the GR values. These plots were computed using an expansion up to order χ 14 of both quantities. We note some characteristic features for each correction. In the case of the quadratic correction controlled by α 1 we see that both ω + and |ω − | increase with respect to the Einstein gravity values. On the other hand, JHEP05(2019)189 Figure 8. Frequencies and radii of the light rings in parity-preserving theories. In blue we plot the quantity corresponding to the prograde orbit and in purple that corresponding to the retrograde one. In the left column we show the frequencies for different values of the couplings and compare them to GR. In the right column we plot the radii R ± for the same values of the couplings.
for α 2 corrections we observe that ω + decreases while |ω − | increases so that the difference between the two frequencies is reduced. As for the cubic correction, it increases or decreases ω + if λ ev > 0 or λ ev < 0 respectively. It has little effect on ω − , but interestingly the sign is different depending on the value of χ. However, in order to characterize deviations from GR it is more useful to look at the ratio of frequencies ω + /|ω − |, that we show for a few cases in figure 9. In GR, this quantity is completely determined by the spin parameter χ, but in these theories it also depends on the combination 4 /M 4 . Thus, if one is able to determine χ by other means, the ratio ω + /|ω − | can be used to constrain the higher-order couplings.

Scalar hair
So far, we have only focused on the geometry, but one of the most remarkable features of the solutions of (2.4) is that the scalar fields acquire a non-trivial profile. In fact, the coupling of the scalars to the quadratic curvature invariants prevent these from being constant whenever the invariants are non-vanishing. A slightly less trivial fact -though also well-known [15,[51][52][53]94]-is that the scalars actually get a charge that can be measured at infinity. More precisely, the scalar φ 1 gets a charge Q while φ 2 gets dipolar moment P , that can be identified by looking at the asymptotic behaviour 12 Using the solution in powers of χ that we have found, we obtain

(4.48)
Remarkably enough, it is possible to guess the general term of these series and to sum them. We find (4.50) 12 The reason for the negative sign in front of Q is that the charge is conventionally defined as where the integral is taken on spatial infinity.

JHEP05(2019)189
One can check that the series expansion of these expressions matches those in (4.47) and (4.48), so they are most likely correct, and they give the exact value of the charges as functions of the spin. In the case of the charge Q, we also check that it agrees with previous results [39,41,95]. Despite having non-vanishing scalar charge, we note however that the solution has no "hair", because the charge is completely fixed in terms of the mass and the spin. In other words, the charge cannot be arbitrary. The reason is that the previous value of the charge is the only one compatible with the requirement of regularity of the solution at the horizon. If we introduce, by hand, any other value of the scalar charge, the resulting solution would develop a singularity at the horizon.
As we mentioned in section 2, in the context of String Theory φ 1 is related to the dilaton, while φ 2 is the axion. In appendix B we show that the precise identification with the effective action of Heterotic Superstring Theory is Then, the dilaton charge D associated to a rotating black hole reads, at leading order in α , This can be expressed in a very appealing form as D = α πT , where T = κ/(2π) is the Hawking temperature of the black hole. It turns out that this intriguing connection between asymptotic charge and temperature (or surface gravity) is not a coincidence, but a general phenomenon that happens in EdGB theory with linear coupling [95].
The field φ 2 gets a dipolar moment instead of charge because it is sourced by the parity-violating Pontryagin density -φ 2 is essentially the scalar that appears in dynamical Chern-Simons gravity [67]. When the spin vanishes we get P = 0, and in fact, φ 2 = 0, so that this kind of scalar hair is not present in spherically symmetric solutions [96].
Besides the asymptotic behaviour, it is also interesting to study the profile of the scalar fields as a function of x. The field φ 2 is odd under the Z 2 transformation x → −x, while φ 1 is even only for θ m = nπ, n ∈ Z. For other values of θ m , φ 1 does not have a defined parity, which is a manifestation of the breaking of the Z 2 symmetry. For instance, when evaluated on the horizon, ρ = ρ + , the field φ 1 is given by (4.52) We only show here a few terms in the χ-expansion for definiteness, but using the solution up to order χ 14 we can determine accurately the profile of φ 1 on the horizon for high values of χ. In figure 10 we plot φ 1 as a colormap on the horizon for χ = 0.65, and 2 α 1 = 2 α 2 = 0.4M 2 . From left to right, the parity-breaking parameter θ m takes the values θ m = 0, π/4, π/2. For θ m = 0 the profile is Z 2 -symmetric and has a mild variation, taking a maximum value at the equator. When θ m = 0, we observe the deformation of JHEP05(2019)189 Figure 10. Profile of the scalar field φ 1 on the horizon. We show here the case for 2 α 1 = 2 α 2 = 0.4M 2 and for a parity-breaking phase θ m = 0, π/4, π/2, from left to right.
the horizon that we reported in section 4.1, plus a "polarization" of the scalar field, that develops a maximum at the north pole and a minimum at the south one. Interestingly enough, the scalar profile provides an intuitive picture of the deformation of the horizon. The northern "hemisphere" grows due to the θ m correction, while the southern one has a smaller size, and this coincides with the fact that the scalar field is "concentrated" on the northern hemisphere, producing a larger energy density there. Thus, the horizon is enlarged in the region that has a greater scalar energy density.

Conclusions
In this work, we have computed the modified Kerr black hole solution in the effective theory (2.4), which provides a general framework to study the leading-order deviations from GR associated to higher-derivative corrections. We expressed the solution as a power series in the spin parameter χ and we showed that including enough terms we get an accurate result even for large values of χ. In this text we have worked with an expansion up to order χ 14 , that provides a good approximation for χ ≤ 0.7, but with the software we supply it should be possible to compute the series to higher orders in χ and to get a solution valid for χ ∼ 1. Although the series expansion involves lengthy expressions, it has obvious advantages with respect to numerical solutions, since it allows for many analytic computations, as we have illustrated in section 4.
We have studied some of the most remarkable properties of these rotating black holes, with special emphasis on the horizon. We have shown that the corrections modify the shape of the horizon, and in particular, that parity-violating interactions break the Z 2 symmetry of Kerr's black hole. We observed the same phenomenon in the case of the ergosphere, and, as far as we know, figure 7 contains the first example of ergospheres without Z 2 symmetry.
In addition, we have computed some quantities that were disregarded in previous studies on rotating black holes in modified gravity. In particular, we have obtained the surface gravity of these black holes, from which one obtains the Hawking temperature according to T = κ 2π , in natural units. Thus, from the results in section 4.1 we conclude that the quadratic curvature terms with non-minimally coupled scalars always increase the temperature of black holes, for any value of the spin. On the other hand, the cubic curvature JHEP05(2019)189 term raises or lowers the temperature depending on the sign of the coupling λ ev and on the value of the spin χ. The modification of Hawking temperature may have important consequences for the evaporation process of black holes [30], and it would be interesting to extend these results by obtaining the value of the temperature non-perturbatively in the coupling and in the spin.
As a first step in analyzing the geodesics of the modified Kerr black holes, we studied the photon rings, i.e. circular light-like geodesics on the equatorial plane. Remarkably, we have found that for parity-breaking theories there are no such orbits: indeed, there are no orbits contained in the equatorial plane because there is no equator at all. Thus, we computed the photon rings for parity-even theories, characterizing the deviations from GR.
Finally, we also noticed the non-trivial scalar fields, and we were able to obtain exact formulas for the monopole and dipole charges. We also computed the profile of the scalar φ 1 on the horizon and we observed how the Z 2 symmetry is broken when the parity-violating phase θ m is activated.
Let us now comment on some possible extensions and future directions. As we already mentioned, it would be interesting to obtain the solution for even larger values of the angular momentum, since the effects of rotation are more drastic when the spin is close to the extremal value. It would also be more or less straightforward to extend the results of this paper to other theories that we did not consider here, particularly the quartic ones in [73,74]. Another possible extension would entail adding a mass term for the scalars in (2.4), though this would considerably increase the difficulty of finding an analytic solution.
We have studied some basic properties of the modified Kerr black holes, but the next natural step is to derive observational signatures of these spacetimes. Analyzing the geodesics of these black holes is an interesting task, as one would potentially observe effects coming from the loss of integrability or from the absence of Z 2 symmetry in parity-violating theories. Once the geodesics are determined, one could study gravitational lensing or the black hole shadow, similarly as done e.g. in [76,77]. However, the most sensitive quantity to the corrections -and that we expect to measure in the near-future thanks to gravitational wave detectors [97] -is the quasinormal mode spectrum of the black hole. Hence, determining the quasinormal modes and frequencies of the rotating black holes presented here is a very relevant task, for which one needs to perform perturbation theory. Analyzing electromagnetic and scalar perturbations for test fields (or for one of the scalars contained in the model (2.4)) should not involve outstanding complications. On the other hand, the study of gravitational perturbations presents a more challenging problem, since one would need to derive the analogous of the Teukolsky equation [98] for the modified Kerr black holes. We feel that this problem should be addressed in future work.
The observation of deviations from General Relativity in astrophysical black holes would represent a tremendous breakthrough that would revolutionize our current understanding of gravity. But even if this is not the case, studying the effects of higher-derivative corrections on black hole geometries provides us with a rich source of new physics, and allows us to learn about new phenomena that could be inherent to an underlying UV-complete theory of gravity.

JHEP05(2019)189
Acknowledgments We wish to thank Vitor Cardoso and Leo Stein for useful comments. The work of PAC is funded by Fundación la Caixa through a "la Caixa -Severo Ochoa" International predoctoral grant. AR was supported by a "Centro de Excelencia Internacional UAM/CSIC" FPI pre-doctoral grant and by a grant from the "Residencia de Estudiantes". Both authors were further supported by the MINECO/FEDER, UE grant FPA2015-66793-P, and from the "Centro de Excelencia Severo Ochoa" Program grant SEV-2016-0597.

A Higher-derivative gravity with dynamical couplings
In this appendix, we are going to motivate our choice of effective action (2.4). Since our goal is to parametrize the leading corrections to vacuum solutions, we will start writing down an action including all possible curvature invariants containing at most 2n derivatives, and then we will discuss which terms are going to induce corrections. By dimensional analysis, a term with 2n derivatives will be multiplied by a factor 2n−2 , where is some length scale that we will assume to be small as compared with the size of the black hole, i.e. GM . It is clear then that the effective action can be always written as where in S (2n) we will include the terms with 2n derivatives. Up to four-derivative terms, we can add the following terms to the Einstein-Hilbert action It turns out that, if the coefficients α i are constants, none of these terms will modify a vacuum solution of GR at O 2 . The reasons are the following: both X 4 and R µνρσR µνρσ are topological terms and therefore do not contribute to the equations of motion. The last two terms are quadratic in Ricci curvature, which means that their contributions to the equations of motion will vanish when evaluated on a GR vacuum solution. In other words, Ricci flat metrics are also solutions of EG plus four-derivative terms. However, we can think of adding dynamical couplings, i.e. promoting α i → α i f i φ 1 , . . . , φ N , where φ A A=1,...,N is a set of N massless scalars. 13 To this aim, we have to include also a kinetic term for them in the action (A.1) so that it becomes where M AB (φ) is the (symmetric) matrix that characterizes the non-linear σ-model. However, as we check a posteriori, the scalars will be excited by the higher-derivative terms at order 2 . Then, we only need to include terms that are at most quadratic in the scalars,

JHEP05(2019)189
which contribute to the gravitational equations at order 4 . Thus, we can expand M AB in a Taylor series and only the constant term will contribute at leading order. By means of a redefinition of the scalar fields, this constant term can always be taken to be the identity matrix: M AB | φ A =0 = δ AB . On the other hand, the generalized action for the four-derivative terms, that we denote again by S (4) , is where we already neglected the R µν R µν and R 2 terms, that do not induce corrections at leading order, and we have now added the term α 5 f 5 (φ) ∇ 2 R that was neglected in (A.2) because in the non-dynamical case it is just a total derivative. In the dynamical case, this term can be written (ignoring total derivatives) as α 5 ∇ 2 f 5 (φ) R and, it is possible to prove that it can always be eliminated, at leading order, by a field redefinition of the metric, so that we can set α 5 = 0.
Indeed it is possible to show that if K µν is a symmetric tensor and we consider a term 4 K µν R µν in the action, the contribution to the field equations is trivial since it can always be eliminated by a field redefinition: g µν → g µν − 4K µν , whereK µν = K µν − 1 2 g µν K α α . To show this, let us compute the contribution of this term to the field equations. Passing this contribution to the right-hand-side of the equations. it can be written as an effective energy-momentum tensor where the dots indicate other possible contributions that vanish when evaluated on the zeroth-order Ricci-flat metric. Then by comparison with (2.15), it is clear that the corrected Einstein equation is solved by g µν = g (0) µν + 4K µν + . . . , beingK µν trace-reversed with respect to K µν . Since the equation is integrable, it is equivalent to preforming a field redefinition, so this kind of terms do not really contain new physics. We can use this result to demonstrate that other type of terms such as φR or G µν ∇ µ φ∇ ν φ, that appear for instance in Horndeski theories [99], can be also removed by a field redefinition.
Let us now analyze the couplings f 1 (φ) and f 2 (φ). The first we can do is to expand the functions around φ i = 0 and neglect O φ 2 terms, which is equivalent to neglect O 6 corrections in the metric. Thus, f i = a i +b iB φ B +O φ 2 and, for the same reasons exposed above, the constant coefficients a i can be neglected. Finally, observe that we still have the freedom to perform a SO (N ) rotation of the scalars that leaves invariant the kinetic terms. Using this freedom, up to global factors that can be reabsorbed in a redefinition of α 1 and α 2 , we can always choose This implies that the theory contains at most two active scalars. In summary, for our

JHEP05(2019)189
purposes the action (A.3) reduces to Then, corrections to vacuum solutions due to these curvature-squared terms will be parametrized by three parameters: α 1 , α 2 and θ m . These terms will induce O 4 corrections in the metric of the solution, since the scalars will be of order O 2 . Therefore, these corrections are equally important to those coming from the six-derivative terms (with constant couplings), which will also induce O 4 corrections in the metric. Since our goal is to parametrize the leading corrections to vacuum solutions in the most general way possible, we shall also include them.
The most general parity-invariant action formed with curvature invariants with six derivatives is There are other six-derivative terms that could be added, such as ∇ α ∇ β R µανβ R µν and ∇ α R µνρσ ∇ α R µνρσ , but these can be reduced to a combination of the terms included in the action. In addition, not all the terms in the previous action are linearly independent. In four dimensions we have two constraints that can be expressed as The first of these constraints actually corresponds to the vanishing of the cubic Lovelock density, X 6 = 0. These relations allow us to express the terms proportional to λ 1 and λ 3 as a combination of the rest of the terms since (A.9) can be rewritten as Hence, we can always set λ 1 = λ 3 = 0. The remaining terms, except those controlled by λ 2 and λ 4 are at least quadratic in Ricci curvature and do not induce corrections on Ricci-flat metrics, so we can ignore them: λ 5 = λ 6 = λ 7 = λ 8 = λ 9 = λ 10 = 0. As already discussed, the term proportional to λ 4 can be eliminated by a field redefinition, since it is proportional JHEP05(2019)189 to Ricci curvature. Consequently, we will not take it into account from now on, so we set λ 4 = 0. Therefore, we are left with only one term out of the initial ten. However, as we did with the four-derivative terms, we can also add parity-breaking densities by using the dual Riemann tensor. One finds again that there is only one independent term, and then the action S (6) reads where we have renamed the parameter λ 2 for evident reasons. Finally, we combine (A.7) and (A.12) to get the action of the effective field theory considered in the main text (2.4) and that we repeat here for convenience S (2n) . (A.13)

B Compactification and truncation of the effective action of the Heterotic String
Let us consider the effective action of the Heterotic Superstring, at first-order in the α expansion, without gauge fields. The ten-dimensional action is given by 14 where α = 2 s , being s the string scale, G (10) is the ten-dimensional Newton's constant, and g s is the string coupling constant. The curvature-squared term 15 is needed in order to supersymmetrize the action at first order in α , which otherwise would not be supersymmetric due to the presence of the Chern-Simons terms in the definition of the 3-form field strength H (see [101] for more details). As a consequence, the Bianchi identity is no longer dH = 0 but it is corrected by where R a b = 1 2! R µν a b dx µ ∧ dx ν is the curvature 2-form. Now, let us perform the dimensional reduction of (B.1) on a six torus, truncating all the Kaluza-Klein degrees of freedom. We get exactly the same action but now with the 14 With respect to the conventions of [88,100], here we are using mostly plus signature gµν → −gµν and the definition of the Riemann tensor differs by a minus sign, i.e. Rµνρ σ → −Rµνρ σ . 15 The curvature-squared term in the Bergshoeff-de Roo scheme [101] is is the curvature of the torsionful spin-connection Ω (−) a b = ω a b − 1 2 Hµ a b dx µ . For our purposes, however, R (−)µνρσ R (−) µνρσ = RµνρσR µνρσ + . . . , where the dots are terms that can be ignored.

JHEP05(2019)189
indices µ, ν running from 0 to 4 where G is the four-dimensional Newton's constant, related to the ten-dimensional one by and we have also introduced e φ∞ = g s . Let us show that, ignoring terms whose contribution to the equations of motion is either zero or trivial, this action can be rewritten in a form such that it is manifestly a particular case of (2.4). First of all, let us rewrite the Bianchi identity (B.2) as Secondly, we have to dualize the 3-form into a (pseudo)scalar ϕ. Following the usual procedure, we introduce a Lagrange multiplier into the action (B.3), The relation between the 3-form field strength and the scalar is found by imposing that the variation of the action with respect to H vanishes, Now, we rewrite (B.6) in terms of ϕ, getting S = 1 16πG d 4 x |g| e −2(φ−φ∞) R + 4 (∂φ) 2 − e 4(φ−φ∞) 2 (∂ϕ) 2 + α 8 R µνρσ R µνρσ + α e 2(φ−φ∞) 8 ϕR µνρσR µνρσ + . . . .

(B.8)
Since this action is not written in the Einstein frame, let us rescale the metric g µν → e −2(φ−φ∞) g µν in order to eliminate the conformal factor. Expanding in (φ − φ ∞ ) and keeping only the leading terms, we get where we have dropped some terms that can be removed with a field redefinition. Finally, defining φ 1 = 2φ − 2φ ∞ and φ 2 = ϕ, and ignoring terms that do not contribute to the equations of motion at leading order, we can write the action in the following final form . (B.10)

JHEP05(2019)189
We have upgraded the Riemann squared term to the Gauss-Bonnet density X 4 since both give the same contribution at leading order (the Ricci 2 and R 2 terms do not contribute). This can also be done by means of a field redefinition. Then, the choice of parameters that gives us the corrections predicted by the simplest compactification of the effective action of the Heterotic Superstring is

C The solution
In this appendix, we show the metric functions H 1 , H 2 , H 3 , H 4 as well as the scalars of the solution, φ 1 and φ 2 up to order O χ 3 , (C.

JHEP05(2019)189
and what is the radius of convergence for a given domain (ρ, x) ∈ Ω. We can study the convergence of the four functions H i at the same time by introducing the "norm" H n , as H n := H 2 1,n + H 2 2,n + H 2 3,n + H 2 4,n .

(D.3)
Thus, every H i,n converges if and only if H n converges. Since we are only interested in the exterior region of the black hole, it is sufficient to look at the convergence for ρ ≥ ρ + , −1 < x < 1. Using the terms of the solution up to order n = 14, we observe that the sequence H n converges in the exterior region if the spin is small enough. We wish to determine the maximum value of χ for which the expansion up to order n = 14 -the one we use thorough the text -provides an accurate approximation to the full series. This value of course depends on the point of the spacetime. For instance, far from the black hole the few first terms in the expansion already provide a very precise result, even for χ ∼ 1. On the contrary, the convergence is worse at the horizon ρ = ρ + , and, specially, at the axes x = ±1. Thus, we should look at the convergence at those points. It is useful to define the relative differences, For instance, when we evaluate this for α 1 = 0.5, α 2 = 0.7, θ m = π/4, λ ev = λ odd = 1 and χ = 0.7, at the north pole of the horizon, ρ = ρ + , x = 1, we get the following sequence (starting at n = 0): 1.3, -0.79, 2.7, -0.74, 2.9, -0.61, -0.52, 0.25, 0.72, -0.27, -0.21, 0.14, 0.076, -0.039, . . . We see that the sequence starts converging from n = 8, and the difference between H 13 and H 14 is barely a 4%. Since the difference with the term n = 15 will be even smaller, we are confident that the series up to n = 14 provides a precision of the order of 1%, for χ = 0.7 when evaluated at the north pole of the horizon. In the rest of points, and for smaller values of χ, the accuracy is significantly greater. We illustrate this in In figure 11 , where we show the profile of H n , for several values of n, in the line x = 1, ρ ≥ ρ + . Also, if for the same values of the couplings we set χ = 0.65, we get d 13 = −0.81%, so the series up to order n = 14 is around five times more accurate than for χ = 0.7.
Finally, one could try to determine what is the maximum value of χ for which the series will converge all the way up to the horizon. In order to find the radius of convergence, one can apply for example the root test to the coefficients H Using the coefficients up to order n = 14 it is difficult to provide a definitive answer, but the results seem to be consistent with χ max ∼ 1. So, it could be possible to get close to the extremal limit adding enough terms in the series expansion, though the number of terms required to get a good approximation increases quite rapidly as we approach χ = 1.