Analyticity of Quasinormal Modes in the Kerr and Kerr–de Sitter Spacetimes

We prove that quasinormal modes (or resonant states) for linear wave equations in the subextremal Kerr and Kerr–de Sitter spacetimes are real analytic. The main novelty of this paper is the observation that the bicharacteristic flow associated to the linear wave equations for quasinormal modes with respect to a suitable Killing vector field has a stable radial point source/sink structure rather than merely a generalized normal source/sink structure. The analyticity then follows by a recent result in the microlocal analysis of radial points by Galkowski and Zworski. The results can then be recast with respect to the standard Killing vector field.


Introduction
When studying linear and nonlinear wave equations on black hole spacetimes, such as the Kerr spacetime and Kerr-de Sitter spacetime, quasinormal modes play a prominent role.Indeed, for linear equations, within certain limitations corresponding to trapped null-geodesics, solutions have an asymptotic expansion at timelike infinity in quasinormal modes.Such expansions, or the corresponding decay or lack thereof statements, have a long history which in the mathematics literature goes back to Sá Barreto and Zworski [SBZ97], Bony and Häfner [BH08], Dyatlov [Dya11,Dya12], Vasy [Vas13], Shlapentokh-Rothman [SR15], Hintz and Vasy [HV15] and Gajic and Warnick [GW20].In the physics literature the importance of these has been clear even longer, going back to Regge and Wheeler [RW57], Vishveshwara [Vis70], Zerilli [Zer70], Whiting [Whi89], Kodama, Ishibashi and Seto [KIS00] and others.For nonlinear equations the non-decaying quasinormal modes become an obstacle to solvability; for equations with gauge freedom, such as Einstein's equation, it is non-decaying modes that are not 'pure gauge' that play an analogous role [HV18].
Quasinormal modes are solutions of the homogeneous wave equation which are eigenfunctions of the covariant derivative along appropriate Killing vector fields.A key consideration for applications is that for similar covariant eigenfunctions as the forcing (right hand side of the wave equation), there should be a satisfactory Fredholm theory.In this case the covariant eigenvalues (resonances) form a discrete set, and the corresponding eigenspaces are finite dimensional.Since Fredholm theory is global, this necessitates working relative to Killing vector fields with suitable global behavior.
Recently Galkowski and Zworski [GZ21] showed that quasinormal modes for non-rotating black holes are real analytic at the horizon; indeed they obtained a substantially stronger microlocal result.In this paper we generalize their result to the case of rotating black holes whose importance is underlined by their ubiquity.Our proof relies crucially on the ability to locally transform the rotating black hole quasimode problem to the non-rotating one, and thus being able to apply the result of [GZ21].This transformation is facilitated by locally considering analogues of quasinormal modes with respect to a different Killing vector field that is lightlike on the horizon; this change is very simple in the Kerr and Kerr-de Sitter case as we discuss below, but in fact works in general for non-degenerate Killing horizons under an additional condition, as is also described below.While these modes are with respect to a different Killing vector field, we can in fact relate these to the quasinormal modes with respect to the original globally well-behaved vector field to obtain the real analyticity result.Indeed, a key feature of the Kerr-de Sitter setting is the presence of two horizons, and the well-behaved Killing vector fields with respect to each of these horizons, while globally well-defined, are ill-behaved at the other horizon.Thus, it is of central importance for our approach to be able to work locally near a horizon to obtain the analyticity conclusions.
It is conceivable that the results of [GZ21] could be proven in the more general setting of [Vas13], which would imply analyticity in the case of rotating black holes.However, such a proof would be significantly more technically involved than the proof we give here.The proof of analytic hypoellipticity in [GZ21] relies on a microlocal normal form of Haber [Hab14], which in turn relies on the relevant Lagrangian (in our case the conormal bundle of the horizon) being radial with respect to the Hamilton vector field.For rotating black holes, there is more intricate internal dynamics and one could therefore not directly apply the results of Haber.Indeed, our reduction to [GZ21] can be considered as a way of 'straightening' the dynamics and thus bringing it to a model form.
Furthermore, it could perhaps be possible to prove analytic hypoellipticity of Keldysh-type operators more explicitly using ODE theory and separation of variables and thereby avoid referring to [GZ21].This was the approach by Lebeau and Zworski in [LZ19] (see also the work by Zuily in [Zui17]), where an explicit class of Keldysh type operators where studied, and certain values for the subprincipal symbol had to be excluded.However, the main purpose of this paper is to reduce the hypoanalyticity of quasinormal modes in Kerr(-de Sitter) spacetimes to the irrotational, i.e.Keldysh-type, case.Applying the ODE approach directly in the rotational case seems cumbersome and it is not clear to us whether it would work.
In the rest of the introduction we describe the precise results in the rotating black hole setting, as well as the generalization to non-degenerate Killing horizons.Then in Section 2 we discuss geometric aspects of these Killing horizons.In Section 3 we then prove our general local result.In Section 4 we use these local results to obtain a global result for joint modes of two Killing vector fields on Kerr and Kerr-de Sitter spacetimes.Finally, in Section 5 we show how these results imply the real analyticity of the quasinormal modes on Kerr and Kerr-de Sitter spacetimes, with modes taken with respect to the standard Killing vector field.
1.1.Kerr and Kerr-de Sitter spacetime.Fix three parameters a ∈ R and m, Λ ≥ 0, such that the polynomial has four distinct real roots r − < r C < r e < r c if Λ > 0 and two distinct real roots r C < r e if Λ = 0.The latter condition is equivalent to |a| < m.
Assuming Λ > 0, the domain of outer communication in the sub-extermal Kerrde Sitter spacetime is given in Boyer-Lindquist coordinates by the real analytic spacetime where The domain of outer communication in the subextremal Kerr spacetime is defined analogously, by passing to the limit Λ = 0. We set r c = ∞ if Λ = 0.The Boyer-Lindquist coordinates could be thought of as spherical coordinates around the black hole, where the black hole is centered at r = 0.Even though they are not defined at the north and south poles θ = 0 and π, it is straightforward to check that the metric (2) extends real analytically to Since the metric extends real analytically, so does linear wave equations with a principal symbol given by the metric in these coordinates.We refer to [HV18, Sec.3], for a more thorough discussion of the geometry of Kerr-de Sitter spacetimes.
These coordinates are singular at the roots of µ(r).In order to define quasinormal modes, we need to extend this metric real analytically over the future event horizon and the future cosmological horizon, corresponding to the roots r = r e and r = r c , respectively.This can be done, for instance, by the following coordinate change: where Φ and Ψ satisfy In the case Λ > 0, we let f : (r e − δ, r c + δ) → R, for a small δ > 0, be a real analytic function such that f (r e ) = −1 and f (r c ) = 1.In the case Λ = 0, there is no cosmological horizon, so we instead assume that lim r→∞ f (r) = 1.
The metric (2) extends real analytically to the manifold φ * ,θ , and is given by (3) We will throughout the paper assume that f is chosen as in [PV, Rmk.1.1], so that the hypersurfaces are spacelike, for all c ∈ R, and that δ > 0 is small enough so that the hypersurfaces are spacelike.The two real analytic lightlike hypersurfaces are called the future event horizon and future cosmological horizon, respectively.Note that the real analytic Killing vector fields ∂ t and ∂ φ , in Boyer-Lindquist coordinates, extend to real analytic Killing vector fields ∂ t * and ∂ φ * on (M * , g * ).
We will consider wave equations on complex tensors.Fixing r, s ∈ N 0 , let T s r U denote the complex (r, s)-tensors on an open subset U ⊂ M * and let ∇ denote the Levi-Civita connection acting on T s r U. We let C ∞ (T s r U) and C ω (T s r U) denote the smooth and real analytic complex tensor fields, respectively.Let P be a wave operator, i.e. is a linear differential operator with principal symbol given by the dual metric, i.e.P = −g αβ ∇ α ∇ β + lower order terms.More precisely, there are complex tensor fields We consider solutions to wave equations P u = f , where the coefficients A and B are invariant under the Killing vector fields ∂ t * and ∂ φ * .This is a natural assumption for geometric wave equations, where A and B are typically given by curvature expressions.Our first main result is the following: Theorem 1.1.Let (M * , g * ) be the subextremal Kerr(-de Sitter) spacetime, extended real analytically over the future event horizon (and future cosmological horizon if Λ > 0).Assume that • A and B are real analytic in M * , . Smooth tensor fields satisfying (ii) and (iii) in Theorem 1.1 and P u = 0 are called quasinormal modes.For functions, these conditions are equivalent to assuming that which is perhaps the more common way to express quasinormal modes.
Combining Theorem 1.1 with the Fredholm theory developed by the second author in [Vas13] and [Vas] (see also [VZ00,Vas21]) and by both authors in [PV], we deduce our second main result, where we consider quasinormal modes only with respect to the Killing vector symmetry where r 0 ∈ (r e , r c ) is the unique point such that µ ′ (r 0 ) = 0, as opposed to modes with respect to both ∂ t * and ∂ φ * separately (as in Theorem 1.1).Concretely, this means that quasinormal modes are supposed to satisfy For solutions to linear scalar wave equations on any subextremal Kerr-de Sitter spacetime, there is an asymptotic expansion in these quasinormal modes up to an exponentially decaying term [PV,Thm. 1.5].This extends the result of [Vas13], by removing restrictions on the angular momentum.
For the Fredholm theory to go through in the case Λ = 0, we will need the induced operator on the modes to be a scattering operator with self-adjoint (i.e.real) scattering principal symbol near spatial infinity in the sense of Melrose [Mel94].Let us use the convention that if Λ = 0, then r 0 = ∞, giving the standard notion of quasinormal modes on the Kerr spacetime.This amounts to making appropriate decay assumptions on A and B: Theorem 1.2.Let (M * , g * ) be the subextremal Kerr(-de Sitter) spacetime, extended real analytically over the future event horizon (and future cosmological horizon if Λ > 0).Assume that iii) in case Λ = 0 we also assume that Im σ ≥ 0 and for some ǫ > 0, then u ∈ C ω (T s r M * ).Here we used the notation S ′ for tempered distributions and the notation F ∈ O ∞ (r α ) for a complex tensor field F if and only if for all k ∈ N 0 , there is a constant C k > 0, such that Remark 1.3.In the case when Λ = 0, one could weaken the assumptions on u, A and B at spatial infinity in various ways and still get a Fredholm problem following the arguments of [Vas].Indeed, the natural condition on u| t * =0 in the case σ ∈ R\{0} is formulated microlocally in terms of variable order Sobolev spaces, c.f. [Vas, Prop.5.28].Moreover, the threshold growth r 1 2 could be adjusted depending on A and B, to allow for more general coefficients, see [Vas,Sec. 5.4.8].We restrict for simplicity to this setting.
The restriction Im σ ≥ 0 for Kerr spacetimes is due to the lack of a directly applicable Fredholm theory for the Fourier conjugated (in −t * ) operators in this case, though alternatives are still available for studying these resonances.For functions, the condition (ii) in Theorem 1.2 is equivalent to assuming that which should be compared with equation (4) above.
In the special case when a = 0, the Kerr(-de Sitter) spacetime simplifies to the Schwarzschild(-de Sitter) spacetime.In this case, Theorem 1.1 and Theorem 1.2 can be immediately deduced from the framework developed by Galkowski-Zworski in [GZ21] as follows: Wave equations for modes with respect to (5) reduce in the coordinate system (t * , r, φ * , θ) to a Keldysh type operator, exactly of the type studied in [GZ21].Galkowski-Zworski prove in [GZ21, Thm.1] (generalizing [Zui17, Thm.1.3]) the analytic hypoellipticity of such operators, thus proving the real analyticity of quasinormal modes when a = 0.In fact, if a = 0, the argument goes through without assuming that the coefficients A and B are invariant under ∂ φ * .Due to the rotation in the Kerr(-de Sitter) spacetime when a = 0, this argument does not go through immediately.The key to be able to apply the analytic hypoellipticity theory by Galkowski-Zworski to the case a = 0 is the main new idea of this paper and is described in the next subsection.
1.2.Non-degenerate Killing horizons.By checking the formula (3) for the extended metric g * , one observes that the Killing vector field where r 0 ∈ (r e , r c ) is the unique point such that µ ′ (r 0 ) = 0, is lightlike at the horizons if and only if a = 0.This turns out to be exactly why the modes with respect to (6) satisfy the useful Keldysh type equation if and only if a = 0.In the Kerr(-de Sitter) spacetime, the Killing vector fields and are lightlike at the horizons H + e and H + c (if Λ > 0), respectively.We show that the mode solutions with respect to these Killing vector fields satisfy equations which are almost of Keldysh type.More precisely, the bicharacteristics associated to the mode equation will have the radial point structure assumed by Galkowski-Zworski in their analytic hypoellipticity result [GZ21, Thm.2].Now, if u satisfies the assumption of Theorem 1.1, then and similarly with r e replaced by r c .This shows that u is a mode solution with respect to both Killing vector fields ( 7) and ( 8).The analytic hypoellipticity result by Galkowski-Zworski thus shows that u is real analytic near the horizons H + e and H + c (if Λ > 0).This is the main step in the proof of Theorem 1.1, the rest follows by standard propagation of real analyticity for wave equations and analytic hypoellipticity of elliptic equations (c.f.[Mar02, Chapter 4]).
In fact, this method is not specific to the Kerr(-de Sitter) spacetime, but turns out to work for any Killing horizon in any real analytic spacetime, assuming the surface gravity of the Killing horizon is nowhere vanishing.Assume therefore that (M, g) is a real analytic spacetime, i.e. a time-oriented Lorentzian manifold, of dimension n + 1 ≥ 2, with sign convention (−, +, . . ., +) and with a real analytic lightlike hypersurface H ⊂ M .We assume in particular that the metric g is real analytic.
Definition 1.4.A real analytic Killing vector field W on M , such that W | H is lightlike and tangent to H, is called a horizon Killing vector field with respect to H.
For each Killing horizon H and horizon Killing vector field W , it is straightforward to check that for a real analytic function κ : Definition 1.5.Given a Killing horizon H and a horizon Killing vector field W , the surface gravity is the real analytic function κ defined in (10).
The key assumption to prove real analyticity of quasinormal modes is that the surface gravity κ is nowhere vanishing.All horizons in Kerr(-de Sitter) spacetimes have surface gravity proportional to µ ′ at the horizons, where µ was defined in (1) (c.f.Step 1 in the proof of Theorem 1.1).This is the reason our result only applies to subextremal Kerr(-de Sitter) spacetimes, since subextremality makes sure that µ ′ does not vanish at the roots of µ, i.e. at the horizons.
Remark 1.6.We note in Lemma A.1 that if for all X ∈ T H and H is connected, then the surface gravity κ is constant.In practice, the condition (11) is often satisfied.Indeed, it is for example satisfied if the spacetime satisfies the Einstein equation with a cosmological constant of any sign or if the spacetime satisfies the dominant energy condition (c.f.[Pet21a, Rmk.1.16]).
In case κ is constant, we get a dichotomy of non-degenerate Killing horizons, where κ = 0, and degenerate Killing horizons, where κ = 0.
As in the previous subsection, we fix r, s ∈ N 0 and consider linear wave operators on complex (r, s)-tensors T s r M and write with complex tensor fields Our third main result in this paper is the following theorem: • W is a real analytic horizon Killing vector field, • the surface gravity κ is nowhere vanishing, • A and B are real analytic and Note that all assumptions in Theorem 1.7 are local.As explained above, we will apply Theorem 1.7 with which will prove the main step in Theorem 1.1 and Theorem 1.2, namely the real analyticity near the horizons.
Our methods require the existence of a horizon Killing vector field.This allows to reduce the wave equation for the modes to the useful (almost) Keldysh form.Surprisingly, a horizon Killing vector field is quite often guaranteed to exist in vacuum spacetimes with horizons.Proving the existence of a horizon Killing vector field has been the central tool in various black hole uniqueness results for the subextremal Kerr spacetime.This line of argument was pioneered by Hawking, who showed that stationary real analytic vacuum black holes with a non-degenerate event horizon necessarily admit a horizon Killing vector field [Haw72,HE73].This result was later generalized to higher dimensional analytic vacuum black holes by Hollands-Ishibashi-Wald [HIW07] and Moncrief-Isenberg [MI08].
There is an analogous problem for compact (also called cosmological) Cauchy horizons in vacuum spacetimes.A conjecture by Moncrief and Isenberg [MI83] states that any compact Cauchy horizon in a vacuum spacetime admits a horizon Killing vector field.The existence of a horizon Killing vector field in that setting would prove that vacuum spacetimes with compact Cauchy horizons are non-generic, which would support the Strong Cosmic Censorship Conjecture in cosmology.During the last decades, Moncrief and Isenberg have made important progress on their conjecture, assuming that the spacetime metric is real analytic [MI83, IM85, MI20].
Remarkably, the existence of a horizon Killing vector field does often not even rely on the real analyticity of the spacetime metric.Alexakis, Ionescu and Klainerman proved in [AIK10a] (see also [IK13]) an analogue of Hawking's theorem, showing the existence of a horizon Killing vector field in a neighbourhood of any bifurcate horizon in smooth vacuum spacetimes, as opposed to real analytic.This result has been central in their approach to prove uniqueness of subextremal Kerr black holes [AIK10b,AIK14] in the smooth setting.For compact Cauchy horizons in smooth vacuum spacetimes, as opposed to real analytic, a horizon Killing vector field has been shown to exist by Petersen in [Pet21b], assuming that the surface gravity is a non-zero constant (extending [FRW99,Pet21a,PR23]).The assumption on constant surface gravity has recently been shown to be equivalent to a weak non-degeneracy assumption for compact Cauchy horizons in vacuum spacetimes, see [BR21] and [GM22].
Though the above mentioned results mainly concern vacuum spacetimes without cosmological constant, one expects them to extend to the case of positive cosmological constant and electro-vacuum spacetimes as well (c.f.[Rác00]).In conclusion, studying wave equations close to non-degenerate horizons (bifurcate or constant non-zero surface gravity), one might in quite wide generality be able to pass to modes with respect to the horizon Killing vector field and analyze the (almost) Keldysh type equation they are known to satisfy by the arguments in this paper.

Suitable coordinates near non-degenerate Killing horizons
The first step towards proving Theorem 1.7 is to define appropriate coordinates near the lightlike hypersurface H: Proposition 2.1.Assume the same as in Theorem 1.7.Then, for any p ∈ H, there is a real analytic coordinate system (x 0 , . . ., x n ), defined on an open neighborhood U ∋ p, such that • the metric g expressed in these coordinates satisfies where    is positive definite and where κ is the (nowhere vanishing) surface gravity.Example 2.3.The simplest example of a spacetime satisfying all our assumptions is M = R n+1 , equipped with the real analytic Misner metric where H = {x 1 = 0}, W = ∂ 0 and surface gravity κ = − 1 2 .
Example 2.4.In fact, even in the subextremal Kerr(-de Sitter) spacetime, one can easily choose coordinates which almost satisfy the conditions in Proposition 2.1, with one (insignificant) difference.To define these, it will be convenient to introduce an intermediate coordinate system, which will only be defined near one of the horizons.Let us start with the event horizon.In terms of Boyer-Lindquist coordinates, define near r = r e .This commonly used analytic coordinate system ( t * , r, φ * , θ) is defined near the future event horizon.Choose now the coordinates , one easily computes that the metric g * at the future event horizon is given by in these coordinates, where g * 22 | x1=0 , g * 33 | x1=0 > 0.Moreover, we have where the surface gravity κ e is given by

the computation in
Step 1 in the proof of Theorem 1.1.These coordinates coincide with the coordinates in Proposition 2.1, up to the multiplication by the positive conformal factor ψ. Since conformal changes of the geometry only reparametrize the lightlike geodesics, ψ is irrelevant for the analysis.However, it is of course natural to construct the coordinates in Proposition 2.1 without a conformal factor.This would here correspond to changing x 1 to x1 by solving the geodesic equation and changing the remaining coordinates x j to xj by demanding that In this new coordinate system, we get precisely the conditions in Proposition 2.1.One analogously constructs similar coordinates near the future cosmological horizon.
Proof of Proposition 2.1.Let first (x 0 , x 2 , . . ., x n ) be real analytic coordinates in an open neighborhood V ⊆ H of p, such that Now let L be the unique lightlike real analytic vector field (transversal to H) along for j = 2, . . ., n. Define now the real analytic coordinate x 1 in an open neighborhood U ⊂ M of p, such that V = U ∩ H, by solving the geodesic equation in direction of L, i.e. we solve and set x 1 = 0 at H. It follows that x 1 is a defining function for U ∩ H.We also extend the other coordinates to U by demanding that in U, for j = 2, . . ., n.The inverse function theorem for real analytic functions implies that this forms a coordinate system.We now show that This is equivalent to W leaving the integral curves of ∂ 1 invariant.Since W is a Killing vector field and the integral curves of ∂ 1 are geodesics, it thus suffices to prove that the initial velocity ∂ 1 | x1=0 of the geodesics are invariant under W , i.e. that Using that W is a Killing vector field, we observe that We also have Moreover, for all j = 0, 2, . . ., n, we have where we have used that g(∂ j , ∂ 1 )| x1=0 is constant by (15).This shows that Taken together, this shows our claim that It is now clear that the metric has the form (12) at x 1 = 0 and that the part (13) is positive definite.Using (10), we compute This finishes the proof.

Real analyticity near general horizons
The goal of this section is to prove Theorem 1.7.In order to explain the idea, let us start by discussing the following example: The condition (ii) in Theorem 1.7 is that Such a mode solution to u = 0 must satisfy the reduced equation This is a Keldysh type equation on the quotient space and [GZ21, Thm.1] implies that v and hence u is real analytic.
The proof of Theorem 1.7 is a generalization of the argument in Example 3.1: Proof of Theorem 1.7.Shrinking U if necessary, we can write the coordinates from Proposition 2.1 as where K ⊂ R n−1 is an open relatively compact subset and ǫ, δ > 0 are sufficiently small.Since ∂ 0 = W | U is a Killing vector field, we would like to eventually reduce P in the x 0 -direction.For this, we first set V := U/∼, where p ∼ q if and only if (x 1 (p), . . ., x n (p)) = (x 1 (q), . . ., x n (q)), i.e.only x 0 (p) and x 0 (q) may differ.The induced coordinates on the quotient space are (x 1 , . . ., x n ) : V → (−δ, δ) x1 × K x2,...,xn , i.e. we have "dropped" the x 0 -coordinate.
The complex (r, s)-tensors on U are complex linear combinations of basis elements of the form where I := (i 1 , . . ., i r , j 1 , . . ., j s ), and we of course have Let us define f := P u and write and by the assumption in Theorem 1.7, we know that It thus follows that the wave equation P u = f , restricted to the subset U, can be written as a linear system of scalar wave equations for each I := (i 1 , . . ., i r , j 1 , . . ., j s ), where the coefficients g αβ , A J I,γ , B J I are independent of x 0 .By the mode condition (ii), we note that which implies that Inserting this into (16) gives a new system of equations where the new coefficients C J I,k and D J I are independent of x 0 .Note also that the sums now exclude derivatives in x 0 .We have thus shown that the equation P u = f is equivalent to a system of equations where the principal symbol of P is where Id is the identity matrix and p(x 1 , . . ., x n , ξ 1 , . . ., ξ n ) : This is where the information about the metric g in Proposition 2.1 becomes useful.We claim that first that where N * {x 1 = 0} denotes the conormal bundle of the horizon {x 1 = 0}.In order to compute the components g ij , for i, j = 1, . . ., n, we first need to invert the full matrix of metric components.By Proposition 2.1, we conclude that for α, β = 0, . . ., n.The components appearing in (17) are given, at x 1 = 0, by is positive definite by Proposition 2.1, we have proven (18).By standard microlocal analytic hypoellipticity at elliptic points in T * V, see for example [Mar02, Thm.4.2.2 & Exe.4.6.4],we hence conclude that u I is microlocally real analytic everywhere at x 1 = 0 except potentially at the conormal bundle N * {x 1 = 0}, i.e. the analytic wave front set at x 1 = 0 is contained in the conormal bundle.We will show the real analyticity at the conormal bundle by applying [GZ21, Thm.2], which requires a computation of the Hamiltonian vector field H p at N * {x 1 = 0}.For this, we first compute ∂ 1 p| N * {x1=0} .At an arbitrary point q := (0, x 2 , . . ., x n , ξ 1 , 0, . . ., 0) ∈ N * {x 1 = 0}, using (19) and Proposition 2.1, we compute We may now compute the Hamiltonian vector field as where we recall that κ is nowhere vanishing.In particular dp| N * {x1=0}\{0} = 0 and H p | N * {x1=0}\{0} ξ • ∂ ξ , which means that the assumptions in [GZ21, Thm.2] are satisfied.Note here that [GZ21, Thm.2] is only proven for scalar valued wave equations, but the argument goes through line by line for systems of equations with a scalar principal symbol, as in our case.Indeed, [GZ21, Thm.2] relies on Haber's normal form in [Hab14] only for the principal symbol.Thus, having scalar principal symbol suffices.Hence [GZ21, Thm.2] implies that u I | t * =0 is microlocally real analytic also at the conormal bundle.It follows that u I | t * =0 is real analytic in an open subset containing {x 1 = 0}.Consequently, u I and therefore u is real analytic in an open neighborhood containing p, which completes the proof.

Joint quasinormal modes
We continue by proving the next main result of this paper: Proof of Theorem 1.1.Let us for simplicity restrict in this proof to the case of complex functions, as opposed to complex tensor fields of higher rank.This will make the proof more transparent and avoid the technical details involved with working with system of equations.All such technicalities are already present in the proof of Theorem 1.7 above.We thus consider functions of the form u(t * , r, φ * , θ) = e −i(σt * +kφ * ) v(r, θ), which are smooth on φ * ,θ ⊂ M * .We aim to prove that u is real analytic on U.
Step 1: Real analyticity near the horizons.We would like to apply Theorem 1.7 with H = H + e/c and therefore need to check that all assumptions of Theorem 1.7 are satisfied.Firstly, the Kerr(-de Sitter) spacetime is a vacuum spacetime with a non-negative cosmological constant, so the dominant energy condition is clearly satisfied.Moreover, the horizons H + e/c are real analytic lightlike hypersurfaces.Secondly, the Killing vector fields W e/c := ∂ t * + a r 2 e/c + a 2 ∂ φ * are clearly lightlike at H + e/c , respectively.Since A and B are invariant under ∂ t * and ∂ φ * , they are also invariant under W e/c .Further, the surface gravity κ e/c of the horizons is computed using the extended metric in (3) as follows: On the other hand, we have r 2 e/c + a 2 cos 2 (θ) The surface gravity of the horizons H + e/c is thus given by κ e/c = ± µ ′ (r e/c ) 2b(r 2 e/c + a 2 ) .
Since this is non-zero, we conclude that H + e/c are non-degenerate Killing horizons with respect to W e/c .Finally, we compute that e/c + a 2 k u, which shows that u is a mode with respect to W e/c .We may therefore apply Theorem 1.7 and conclude that u is real analytic in open neighborhoods of the horizons, which are invariant under W e and W c , respectively.
Step 2: Real analyticity in the domain of outer communication.We now prove real analyticity of u in the open subset φ * ,θ .(Recall that r c = ∞ in the Kerr spacetime).The Boyer-Lindquist coordinates (t, r, φ, θ) are defined on this set and are convenient to work with.Since the conditions (ii) and (iii) in Theorem 1.1 imply that u(t, r, φ, θ) = e −i(σt+kφ) w(r, θ), so we can equally well consider the modes with respect to the Boyer-Lindquist coordinates.The dual metric G of g in Boyer-Lindquist coordinates is We begin by proving real analyticity of w in the open subset (r e , r c ) r × (0, π) θ i.e. we leave out the north and the south pole of S 2 φ,θ for the moment.Since we have assumed that the coefficients of P are independent of t and φ, the function w satisfies an induced equation on (r e , r c ) r × (0, π) θ , with principal part given by 1 Since µ(r), c(θ) > 0 in this set, the induced equation for w is elliptic with real analytic coefficients.Standard analytic hypoellipticity, see for example [Mar02, Thm.4.2.2 & Exe.4.6.4],therefore implies that w is real analytic in (r e , r c ) r ×(0, π) θ and hence u is real analytic in R t × (r e , r c ) r × S 1 φ × (0, π) θ .We now turn to show that u is also real analytic at the north and south poles of S 2 φ,θ , i.e. at the limits θ = 0 and θ = π, still with r ∈ (r e , r c ).Note that the expression (21) does not extend smoothly to those points.We now write u(t, r, φ, θ) = e −iσt z(r, φ, θ), i.e. the idea is to show real analyticity of z(r, φ, θ) := e −ikφ w(r, θ), which is smooth in (r e , r c ) r × S 2 φ,θ .Since the coefficients of P are independent of t, we get an induced equation for z, with real analytic coefficients and principal part We claim that this operator is elliptic at θ = 0 and θ = π.For this, we note that µ(r)∂ 2 r + bG S 2 at the north and the south pole of S 2 .Since µ(r) > 0 for r ∈ (r e , r c ) and b > 0, we conclude that (22) is indeed elliptic there as well.Again, standard real analytic hypoellipticity, as in for example [Mar02, Thm.4.2.2 & Exe.4.6.4],implies that z, and therefore u, is real analytic also at the north and the south pole if r ∈ (r e , r c ).To sum up, we now know that u is real analytic in the domain of outer communication and slightly beyond the horizons, i.e. in a region of the form R t * × (r e − ǫ, r c + ǫ) r × S 2 φ * ,θ .
Step 3: The region beyond the horizons.It remains to prove real analyticity in the regions beyond the horizons (only the event horizon if Λ = 0).Consider first the region R t * × (r e − δ, r e ) r × S 2 φ * ,θ , beyond the event horizon.We may use Boyer-Lindquist coordinates also here, since this region does not intersect any horizon, where the coordinates would not be defined.Let us again consider z(r, φ, θ) := e −ikφ w(r, θ), which by assumption is smooth in (r e − δ, r e ) r × S 2 φ,θ .In this set, we have µ(r) < 0. Again, the coefficients of P are independent of t and the principal part of the induced equation for z can be read off from (22) to be Since the operator is elliptic on {r} × S 2 φ,θ , for all r ∈ (r e − δ, r e ), we conclude that all inextendible bicharacteristics of the induced operator pass through all hypersurfaces {r} × S 2 φ,θ precisely once.Moreover, the induced equation for z is a linear wave operator with real analytic coefficients on the (r e − δ, r e ) r × S 2 φ * ,θ .and we know by Step 1 that z is real analytic in an open subset (r e − ǫ, r e ) × S 2 φ * ,θ for some ǫ > 0. Propagation of analytic singularities, see for example [Mar02, Thm.4.3.7 & Exe.4.6.4],therefore implies that z is real analytic in (r e − δ, r e ) r × S 2 φ,θ , and hence u is real analytic in R t × (r e − δ, r e ) r × S 2 φ,θ .One similarly treats the subset φ,θ , in case Λ > 0. This finishes the proof.

Standard quasinormal modes
We finish by proving our last main result: Proof of Theorem 1.2.As in the proof of Theorem 1.1, let us for simplicity restrict to the case of complex functions, as opposed to complex tensors of higher rank.This will again make the proof more transparent and avoid technical details that are completely analogous to the corresponding part of the proof of Theorem 1.7.It is convenient to change coordinate system to one that is better suited for the quasinormal mode condition (ii) in Theorem 1.2.We introduce the new coordinate system (τ * , r, ψ * , θ), where with again r 0 ∈ (r e , r c ) is uniquely defined by µ ′ (r 0 ) = 0, and r 0 = ∞ if Λ = 0. Note that where z is smooth in (r e − δ, r c + δ) r × S 2 ψ * ,θ , where r c = ∞ if Λ = 0. Since the coefficients of P are independent of t * and φ * , and therefore of τ * , it follows that z satisfies a τ * -reduced equation in (r e − δ, r c + δ) r × S 2 ψ * ,θ , where P σ u := e iστ * P e −iστ * u .
We now further decompose into angular modes where we claim that each summand where p = (ψ * , θ), it is smooth as claimed.Since the coefficients of P are independent of t * and φ * , and therefore of ψ * , and P σ acts diagonally on the ψ * -Fourier modes, it follows that Now, if Λ > 0, then [PV, Thm.2.1] implies that the operator P σ is a Fredholm operator between appropriate function spaces containing e −ikψ * v k .Since the kernel is finite dimensional, it follows that only finitely many such terms can be non-zero.We conclude that Each term satisfies the assumption of Theorem 1.1 and are therefore analytic.Hence the finite sum is also real analytic, concluding the proof when Λ > 0.
In order to similarly proceed in the case Λ = 0, we need to instead use the Fredholm theory developed in [Vas] (remember that the above coordinate change is trivial when Λ = 0) to deduce that in fact z(r, φ * , θ) = N j=1 e −ikj φ * v kj (r, θ) is a finite sum.In this case we have r c = ∞ and the cosmological horizon is replaced by an asymptotically Euclidean end.For the analysis near the event horizon, the methods based on [Vas13] described above can be applied without changes.However, the analysis near the asymptotically Euclidean end cannot be based on [Vas13], we instead need to use a slight generalization of [Vas,Prop. 5.28].
Let us therefore briefly recall how the Fredholm problem was set up in [Vas,Prop. 5.28].We begin by bordifying the space (r e − δ, ∞) r × S 2 φ * ,θ , at r = ∞ by introducing x := 1 r , i.e. we radially compactify spacelike infinity.We thus write where R 3 is the radially compactified R 3 .On these spaces, we define where s, l are variable order differential and decay orders (as x → 0), which we will choose below.We refer to [Vas,Sec. 5.3.9]for the definition of variable order weighted Sobolev spaces H s,l (R 3 ).Note that near the spacelike hypersurface The characteristic set of P has two components, one close to the event horizon and a scattering characteristic set at x = 0, in particular, the characteristic set at fiber infinity near x = 0 is empty.By the decay assumptions on A and B, the scattering principal symbol of P at x = 0 is given by for the fixed σ with Im σ ≥ 0, and any ξ ∈ sc T * {x=0} V.
If Im σ > 0 it follows that σ 2 / ∈ [0, ∞), which implies that P is elliptic as set up in (28) and consequently a Fredholm operator for any order s, l.However, in case σ ∈ R\{0}, there is a scattering characteristic set at x = 0, given by all ξ ∈ sc T * {x=0} V with |ξ| = |σ|.As shown in [Vas,, the sets is nowhere vanishing, lightlike and tangent to H.If (11) is satisfied, then there is a constant κ ∈ R, such that Proof.Using that W | H is lightlike and tangent to H, we compute that for all vector fields X, Y , tangent to H, we have since also [X, Y ] is tangent to H. Hence ∇ X W is tangent to H and normal to H, meaning that there is a unique one-form ω on H, such that The assertion in the lemma is thus that ω(W | H ) is constant.Since W is a Killing vector field, with W | H tangent to H, it is immediate that For any X ∈ T H, we have It thus remains to show that dω(X, W | H ) = 0, for all X ∈ T H.For this, we first note that for all X, Y ∈ T H, we have We therefore conclude that for all X ∈ T H, which proves that κ is constant.

Remark 2. 2 .
These coordinates are essentially the Gaussian null coordinates introduced by Moncrief-Isenberg in [MI83], with the extra condition that ∂ 0 is the horizon Killing vector field restricted to an open neighborhood.(This is precisely what is obtained a posteriori after the construction of the horizon Killing vector field in [MI83].) * are both Killing vector fields, since a and r 0 are constant.It follows that u is a quasinormal mode if and only if u(τ * , r, ψ * , θ) = e −iστ * z(r, ψ * , θ).
sc is similar to Y s introduced above.Analogous to above, defineX s,l sc := u ∈ Y s,l sc | P u ∈ Y s−1,l+1sc and consider P : X s,l sc → Y s−1,l+1 sc .