Kerr scattering coefficients via isomonodromy

We study the scattering of a massless scalar field in a generic Kerr background. Using a particular gauge choice based on the current conservation of the radial equation, we give a generic formula for the scattering coefficient in terms of the composite monodromy parameter σ between the inner and the outer horizons. Using the isomonodromy flow, we calculate σ exactly in terms of the Painlevé V τ -function. We also show that the eigenvalue problem for the angular equation (spheroidal harmonics) can be calculated using the same techniques. We use recent developments relating the Painlevé V τ -function to Liouville irregular conformal blocks to claim that this scattering problem is solved in the combinatorial sense, with known expressions for the τ -function near the critical points.


Introduction
Black hole scattering theory is an important problem from the general relativity and string theory perspective. It serves as a way of testing linear stability [1], superradiance [2], black hole entropy [3] and relaxation times in AdS-CFT [4]. For the rich story behind it see [5,6]. Scattering is also important to study astrophysical problems related to the detection of gravitational waves and a possible mechanism for short gamma ray bursts, for example. See the introduction of [7], references therein and [8,9] for discussions. It is also a daunting problem from the mathematical side, involving unknown properties of new special functions [10]. Most of the previous study in four dimensions on the subject relied on numerical analysis, either by direct numerical integration or by expansion of the solution in terms of known special functions [11][12][13]. Also, decoupling limits like the near-extremal cases were taken directly on the metric (see [14] for applications in Kerr-CFT).
In this letter we build on previous analytical work, a reverse Hilbert-Riemann problem [15], based on isomonodromy equations [16] to relate the scattering coefficients for a massless scalar field in four-dimensional Kerr background in terms of the Painlevé V transcendent τ -function [17] and show that the problem can be solved in the combinatorial sense using the expasions of the τ -function near the critical points. We also show that the angular equation can be solved implicitly using the same techniques, in a new way of extracting eigenvalues for the spheroidal wave functions [18,19]. The results shown here can be understood as the natural conclusion to the application of the monodromy problem to solve scattering problems as applied to black holes [20][21][22][23]. The result follows from a careful consideration of the "linearization" of the Painlevé V equation in the sense given by Miwa, Jimbo et al. in [24][25][26], and of the time-reversal symmetry. These two points JHEP11(2015)144 allows us to place boundary conditions such that the two linearly independent solutions of the radial equation correspond to the two solutions involved in the connection problem of the Painlevé V equation [27]. Recent developments on the relation between Painlevé transcendents and the AGT conjecture [28,29] show the striking connection between black holes and black hole scattering to Conformal Field Theory [30].
This letter is organized as follows. In section 2 we treat the radial equation for the massless Klein-Gordon field in a generic four-dimensional Kerr background. We write the equation in a canonical confluent Heun type and discuss the time-reversal symmetry. In section 3 we discuss how the scattering problem for the radial equation is connected to the connection problem of the confluent Heun equation. We describe the (extended) monodromy data of the confluent Heun equation, including the Stokes parameters, and use the time-reversal symmetry to parametrize the scattering coefficient in terms of the monodromy parameters. In section 4 we discuss how to extract these parameters from the Painlevé V τ -function, using the isomonodromy flow. In section 5, we outline a procedure to extract the eigenvalues for the angular equation using the same structure. Finally, in section 6 we list the results of [17] which give a combinatorial expansion for the Painlevé V τ -function near t 0 = 0, which would correspond to our extremal limit. We close with a short discussion and future directions in section 7.

Klein-Gordon in Kerr
It is known that the massless Klein-Gordon equation in a Kerr background is separable [31], i. e., the solution can be written as Φ(x µ ) = e −iωt e imφ R(r)S(θ). The angular part S(θ) is solved in terms of oblate scalar spheroidal harmonics, whose eigenvalues (the separation constant C ℓ,m ) can be approximated for different asymptotic regimes, including low and high frequency limits [18,32], and will be the object of study in section 5. We will focus for now on the radial part where Q(r) and W 2 r are given by The form which will be useful for the analysis has with the parameter

JHEP11(2015)144
The radial equation is written in this parametrization as which is known as the confluent Heun equation, with parameters given by where we defined the temperatures T ± and angular velocities Ω ± for the outer and inner horizons. For future reference, we also define Also, eq. (2.1) has real coefficients, so if R(r) is a solution, so is R(r) * . We then define the radiation flux as One should bear in mind that real r means z * = −z, because of (2.3).

The Fuchsian system
Writing the solutions for (2.5) in terms of Frobenius expansions at z = 0 or z = t 0 is straightforward. On the other hand, extracting the scattering coefficients from the solutions is not straightforward in the sense that not enough information is known about the function defined by the ODE. To wit, we are interested in the so-called connection problem: how a solution of (2.5) with known behavior at one singular point, say y t 0 (z) at z = t 0 , is related to the pair of solutions with known behavior at another singular point, say, y ± ∞ (z) at z = ∞. So, for example, in C and D are the connection coefficients. If y t 0 (z) represents a purely ingoing wave at the outer horizon z = t 0 and y ± ∞ (z) ingoing and outgoing waves at infinity, then the connection coefficients are indeed the scattering coefficients. As we will see below, solving the connection problem allows us to solve the scattering problem.
In [16], the authors proposed a method based on the Schlesinger equations to extract the scattering coefficients. It relies on viewing the ODE (2.5) as the equation satisfied by the first row of the 2 × 2 matricial Fuchsian system

JHEP11(2015)144
and to extract the connection coefficients from the monodromy data of the Fuchsian system. The monodromy data for the system at hand is defined as follows: suppose A(z) is meromorphic and let {z i } be the isolated singular points of the system (3.2). The point at z = ∞ is studied by the substitution u = 1/z. The solution Φ(z) will also cease to be analytic at these points and behave singularly there. SinceΦ(z) = Φ(e 2πi (z − z i ) + z i ) is also a solution of (3.2), then we must haveΦ(z) = Φ(z)M i , with M i a constant matrix associated with a circuit around z i -the monodromy. The form of M i depends on the behavior of A(z) near the singular point z i . If A(z) diverges like (z − z i ) −1 , then the singular point is called regular and M i is conjugated to diag(e 2πiα + i , e 2πiα − i ), where α ± i are the solutions of the indicial equation near z i . This form of M i stems from the leading behavior (z − z i ) α ± i of the solution near z i , a branch point for Φ(z). We will for now on assume that the differences α + i − α − i are not integers for the radial equation. If, however, A(z) diverges like (z − z i ) −1−r i , where the integer r i > 0 is the Poincaré rank, the solution will have an essential singularity and will experience the Stokes phenomenon: the leading term in the Frobenius expansion displays an exponential behavior which diverges or converges depending on which direction the limit z → z i is taken in the complex plane. In order to describe it, we will specialize to our case. In equation (2.5), the matricial system has A(z) with the structure where there are two regular singular points, at z = 0 and z = t, and an irregular point at z = ∞ with Poincaré rank 1. We also have the usual Pauli matrix σ 3 . The elements of A 0 and A t in (3.3) are independent of z and are obtained from the ODE (2.5). This identification has extra freedom because the Fuchsian system has more free parameters than the ODE. We thus pick the gauge choice which, along with the ODE parameters (2.6), fixes the elements of A i up to conjugation by e uσ 3 . The epithet gauge choice will become clear in the next section. In order to describe the system at the irregular singular point z = ∞, one defines the sectors [27,33] S for j ∈ Z. On each S j , one defines the asymptotic behavior for These canonical solutions are connected by the Stokes matrices Ψ j+1 (z) = Ψ j (z)S j and, given that Ψ j+2 (e 2πi z) = Ψ j (z)e −πiθ∞σ 3 , we have S j+2 = e πiθ∞σ 3 S j e −πiθ∞σ 3 , and the Stokes matrices are all determined by any consecutive 2, which have the following structure and the s 2j and s 2j+1 are called Stokes parameters. Local solutions of the Fuchsian system at z = ∞ are defined in the Riemann surface of the logarithm, and these parameters describe how the solution (3.6) changes as one changes the leaf of the surface. The (extended) monodromy data for (2.5) is completely determined by the monodromy matrices M 0 and M t , at the regular singular points 0 and t, respectively, and the Stokes matrices S 0 and S 1 , which determine the (formal) monodromy matrix at ∞ in sector S j by Knowledge of the θ i and of two consecutive Stokes parameters, say, s 1 and s 2 is sufficient data for solving the connection problem. We define the connection matrices E i by and those are given in terms of θ i , σ, s 1 and s 2 in [27,33]. Note that, because M ∞ M t 0 M 0 = 1, we have that the combined monodromy parameter 2 cos πσ = Tr M 0 M t satisifies 2 cos πσ = 2 cos πθ ∞ + e iπθ∞ s 1 s 2 , (3.10) The connection matrices E i can be seen -when properly normalized by det E i = 1 -as the matrices that implement the change to the "natural basis" of solutions of the ODE at z i where the monodromy M i is diagonal. Explicit formulas for E i in terms of monodromy data are given in [33]. We use the remaining gauge symmetry A(z) → e uσ 3 A(z)e −uσ 3 so that both entries in the first row of Φ(z) have equal but opposite current density normalization (2.8), defined as j ± . So we have the entries [Φ(z)] 11 being the normalized incoming and outgoing waves near infinity. Because the j ± are conserved, this choice also means that the entries will be normalized near the outer horizon t 0 , up to a phase. Therefore, the entries of the connection matrix will give the explicit scattering coefficients, allowing us to write T and R such that and by construction the time-reversed wave is equal, up to a phase, to the entry [Φ(z)] 12 near t 0 . We notice that Det E ∞t 0 = 1 implies |R| 2 = 1 − |T | 2 , which is the condition of flux conservation. Now, by substitution of the diagonalized forms of the monodromy matrices in Tr M ∞ M t = Tr M −1 0 , and using (3.12), we obtain the formula , (3.13) which gives the coefficient in terms of the composite monodromy parameter σ. In order to compute it, we will need the structure of isomonodromy deformations, outlined in the next section. After embedding this system in the isomonodromy formulation, we will be able to give explicit coordinate dependence of the wave functions.

JHEP11(2015)144 4 The isomonodromy method
The idea of considering the Fuchsian system (3.2) instead of the scalar ODE came in the context of the Riemann-Hilbert problem, since a simple counting argument [35] shows that the generic scalar ODE does not have enough parameters to realize all possible monodromy matrices. The Fuchsian system, however, has too many parameters, so there is a family of A(z) in (3.2) with the same monodromy data. The description of this family can be understood physically: suppose we are given an invertible solution Φ(z) of (3.2) with some singular points. Then we can understand A(z) = [ d dz Φ(z)]Φ(z) −1 as a "gauge potential" describing non-abelian charges, dipoles or multipoles at each singular point. The physical configuration described by the gauge potential is invariant under generic gauge transfor- meaning that the monodromy data (non-abelian holonomies) are also invariant under these transformations.
In the case at hand (2.5), with A(z) of the form (3.3), we supply a t-component to this gauge potential, in such a way that the "field strength" The monodromy data of the Fuchsian system (3.2) will be independent of t if A 0 (t) and A t (t) in (3.3) satisfy the equations above (4.2). For the relationship between the Schlesinger equations and the theory of flat holomorphic connections we refer to [24][25][26]35]. The system (3.2) has some direct conserved charges. The gauge choice (3.4) is justified now because these are the constant of motions under Schlesinger evolution. With these constraints, the equation satisfied by any element of the first row of Φ(z) in (3.2), with A(z) given by (3.3), is of the form where C 0 , C t , λ and µ are functions of the entries of A(z). Compared to (2.5), the system above has an extra singularity at z = λ. However, this is an apparent one -its monodromy is trivial, as the indicial equation has exponents 0 and 2 and there is no logarithm behavior due to the relation between µ, λ and t, The system (4.2) is better known in another clothing, as

JHEP11(2015)144
can be checked to satisfy the Painlevé V equation. These are part of the Painlevé transcendents family of differential equations: non-linear second order differential equations whose solutions do not possess movable branch points [35]. These define new special functions, with applications in integrable systems, random matrix theory and conformal field theory -see [17,36] for references in those applications. In [27], asymptotic expressions for the Painlevé V system were given in the guise of the τ -function in terms of the monodromy data. The τ -function is the most natural isomonodromy invariant that can be defined, with a clear interpretation as a generating functional in quantum field theory applications [37][38][39].
As the problem now stands, knowledge of the τ function solves the system completely. We pick suitable initial conditions for µ and λ in (4.3), and from those recover the A 0 and A t , which then will set the initial conditions for the τ function. By choosing the initial conditions based on the parameters of the initial Heun equation (2.6) we have which can be formally inverted to yield the Stokes parameters s 1 and s 2 , and hence σ. We note that our definition (4.6) differs slightly from the one in [27] by One also notes that (4.8) allows us to interpret the τ -function as the generating function that implements the canonical transformation between the non-trivial monodromy data s 1 , s 2 and the canonically conjugated parameters of the ODE t 0 , c t 0 . The structure outlined here gives another interpretation for the τ -function. One of the most fundamental results from Miwa, Jimbo et al. on the Schlesinger system [24][25][26] -see also [35] -was the Hamiltonian structure of the isomonodromy flow. This structure stems from the symplectic structure of flat holomorphic connections (see [16] for a description), but can be intuitively understood from the algebraic constraint (4.4). One first notes that C 0 and C t are not independent, they are related to the asymptotic behavior at z = ∞, via θ ∞ , Varying t will now tie the variation of λ and µ in a symplectic system, in which C 0 can be thought of as the Hamiltonian function. Now, by "solving" the Hamiltonian system one JHEP11(2015)144 means giving a canonical transformation from the coordinates C 0 , t 0 to coordinates where the isomonodronic flow is trivial. But the latter is parametrized by the monodromy data, which by construction is invariant under the flow. The τ -function is then the generating function that implements this canonical transformation, as evidenced by (4.8). In [40] a set of Darboux coordinates were constructed for the Painlevé VI case. It would be interesting to give a parallel of that construction to the case at hand.

Asymptotic wavefunctions and normalization
Now, with the embedding of the radial equation (2.1) in the matricial system (3.2) via (4.7), we can work out the specific form of the wave functions in terms of the coordinate r. The form is not necessary for our results, but will help comparing with previous work [12,13]. The normalization of the elementary matrix Φ(z) is given by with boundary conditions at z = ∞ given by (3.6). It is straightforward to verify that the Wronskian between the two entries of the first row is given by: : we arrive at: (4.14) We remind the reader that the value k is "gauge dependent", and proceed to normalize the asymptotic from the entries of Φ(z). From (3.6), we have: where the limit is taken through the direction arg z = 1 2 π (z → +i∞). Following [27] we define the Frobenius basis at z = t 0 to be: where G t 0 is a matrix which diagonalizes A t . Its particular form will not be important for us save that for the conditions (4.7) it can be checked to be of the lower triangular form. From this we can write the Wronskian-normalized basis for wavefunctions at z = t 0 :

JHEP11(2015)144
Now, we can write the current-normalized wavefunctions in terms of the coordinate r: As anticipated in the previous section, a suitable choice of k can make the relative normalizations equal. The result (3.13) follows. In order to implement the time-reversal procedure, one has to check the relative phases and include them explicitly in the formula for the connection coefficients. At any rate, the basis at r = r + is different from the one considered by [12,13], which were mainly concerned about generic spin Teukolsky equation and did not work the normalization out explicitly near r = r + . In the non-zero spin case, one has to solve not only for the monodromy but also for the Teukolsky-Starobinski identities [41]. It is an interesting future problem to try and extend the results presented here to that case.

Angular equation
The general solution for the composite monodromy parameter σ using (4.8) (and (3.10)) can of course be applied to more generic situations. The associated equation for the angular part of the wavefunction S(θ) is also of the confluent Heun type. Following [19], the scalar spheroidal harmonics are the solutions of the angular equation: for which the solution has the following behavior near the singular points x = cos θ = ±1: for integer m. It is clear that this behavior will only happen for particular discrete values of the constant C ℓ,m , which can be fed into (2.1).
To bring this equation into the canonical confluent Heun form (2.5), we make: and the new parameters are: Now, from the discussion above, the existence of a solution of the confluent Heun equation with behavior given by (5.2) puts constraints on the connection matrices. In this particular

JHEP11(2015)144
case, one of the natural solutions at z = 0 will not "mix" with the natural solution at z = 1. This would correspond to the vanishing of one of the constants at (3.1). One should not that the second natural solution will of course mix. Also, since the parameter θ 0 and θ t are integers, there will be logarithm behavior for the other solution -in the case ω → 0 it would correspond to the associated Legendre functions of the second kind. It is not hard to see that from this condition the composite monodromy 2 cos πσ = Tr M 0 M t will have to be special: as required by the vanishing of the analogue of (3.13) for this case. Now, using (4.8), one has a -somewhat formal -solution for the eigenvalues: Since τ (t) is also a function of the Stokes parameters s 1 , s 2 , the second equation in (4.8) is also necessary to find the C ℓ .

Asymptotics of Painlevé V
In the spirit of making the analysis self-contained, we copy the relevant formulae about the Painlevé V τ -function from [17]. In the following we use their definition for the τ -function: and assume θ i , σ and θ i − θ j , σ − θ i are not integers. The expansion for the tau-function is of the formτ where the irregular conformal block B is given as a power series over the set of Young tableaux Y: with coefficients , (6.4) where λ denotes a Young tableau, λ i is the number of boxes in row i, λ ′ j is the number of boxes in column j and h λ (i, j) = λ i + λ ′ j − i − j + 1 is the hook length related to the box JHEP11(2015)144 The structure constants C are rational products of Barnes functions 5) where G(z) is defined by the functional equation G(1 + z) = Γ(z)G(z). The parameters σ and s in (6.2) are related to the "constants of integration" of the Painlevé V equation. In our treatment, they are functions of the Stokes parameters s 1 and s 2 . The parameter σ is given by (3.10), whereas the expression forŝ in (6.2) is rather long and involved. We will outline the procedure in section 10 of [33] to compute it. LetM 0 ,M t be the monodromy matrices for the hypergeometric equations. They are of the formM i =Ê −1 i e iπθ i σ 3Ê i witĥ E i given by equation (10.15) in [33]. They satisfyM tM0 = e −iπ(θ 0 +θ 1 ) diag(e −iπσ , e iπσ ). We introduce the matrix whose purpose is to match the diagonal monodromy of the hypergeometric system at infinity to the non-diagonal monodromy of the confluent hypergeometric system at the origin. It is defined up to a diagonal matrix s −σ 3 . We now have the monodromy parametrization for E i in (3.9): where R 2 = s 1 /s 2 serves to bring M ∞ to a symmetric form. With this intricate construction, the parameter s can be read implicitly from M t = M −1 ∞ M −1 0 -or more easily from its trace. The parameterŝ in (6.2) is given byŝ = − 1 2 (σ + θ ∞ )s 2 . These expressions are purely algebraic in the sense that they can be computed from the constraints M ∞ M t M 0 = 1 obeyed by the monodromy matrices. The expansion (6.2) is obtained from the AGT conjecture [28], relating the Liouville conformal blocks to the instanton partition function of N = 2 SUSY quiver theories in four dimensions via a confluence limit.
The asymptotic limit of the Painlevé V τ -function was also considered in great detail in [27,33]. For our purposes, in the limit t 0 → 0, the sum (6.2) can be approximated by three terms:τ where we assume that 2n < σ < 2n + 2. The value of n ∈ Z will be obtained below, and we will also assume that ̺ = σ − 2n < 1, which can be accomplished by shifting n, and the calculation that ensues will be analogous as below. As per (6.2) the leading parameters of the expansion of f (t) are: , (6.10)

JHEP11(2015)144
Now, (4.8) is written as: As alluded above, the second equation does not entail information about the initial value of the isomonodromy flow beyond the value of t 0 . We will see it as fixing the value ofŝ in terms of ̺. Plugging the expansion for f (t) as above, we find: (6.14) Now, back to (4.8) and isolating ̺, we find: where ǫ = min(̺, 1 − ̺). To first order in t 0 , we then have: where we see that the integer n is chosen to offset the integer part of C ℓ,m . Also, for small ω, ̺ 0 is purely imaginary, but the t 0 corrections will introduce a real part, and the expansion will have poor convergence. This is verified for the Painlevé VI expansions using formulas analogue to (6.2) in [17], as convergence of the series improves for finite t 0 . Notwithstanding these remarks, we find qualitative agreement with [12,15] even though the small t 0 expansion is not strictly speaking a low-frequency expansion, as other parameters also depend on ω. At any rate, the expression (4.8) should be valid for generic values of the parameters, and with the results of section 5 where the separation constant C ℓ,m are also computed using isomonodromy, we can claim the formal analytical solution of the scattering problem is presented. We leave for future work a survey of the values for the scattering coefficients, as well as expansions for the near extremal and near-superradiance θ 0 ≈ 0 cases.

Discussion
In this letter, we sucessfully related the scattering coefficients for a massless scalar field in a Kerr background in terms of the Painlevé V τ -function. By choosing the current normalization of the ingoing and outgoing waves, we first establish that the expression (3.13) does give the transmission coefficient directly in terms of the composite monodromy parameter σ. Also, by relating the problem of finding σ to the inverse Riemann-Hilbert problem, we are able to give analytic, if implicit solutions for the composite monodromy (4.8) and the eigenvalue of the angular equation (5.6) in terms of the Painlevé V τ -function. Also,

JHEP11(2015)144
approximate expressions for small isomonodromy parameter -related to low-frequency and/or near extremal case -can be obtained from the asymptotics of Painlevé V equation in terms of more elementary functions. It is our hope that the application of the techniques outlined here and in previous work by the authors [16] will help solve many outstanding problems in scattering for black-hole backgrounds, like Kerr-(A)dS, and higher spin perturbations. It should be stressed that the usefulness of (4.8) goes far beyond the low-energy or the near horizon expansions, -both tied to small t 0 -allowing one to study normal modes, and the onset of superradiance at θ 0 ≈ 0. One can also use essentially the same arguments given here to study scalar perturbations between the inner and outer horizon, providing yet another tool for testing for instabilities. We will leave these important directions to future work.
The result showed here relied on the conservation of current argument, which allowed us to equate the scattering problem to the (mathematically defined) connection problem of the radial equation. Although this needs proof, the same argument should work for arbitrary spin perturbations, the so-called Teukolsky Master Equation (TME), whose radial and angular equations for Kerr background are also separable and yield the same confluent Heun type considered here. This should provide a powerful analytical tool to deal with normal modes, generation of gravitational waves by black holes, and the linear stability of the inner horizon.
In a more abstract viewpoint, the results extracted here point to the striking relationship between black hole physics and Conformal Field Theory. Relations for the τfunction (6.2) were obtained because of recents results concerning Liouville Field Theories in the semi-classical limit, where the Painlevé V τ -function arises as the expansion of the irregular conformal blocks, which are obtained from a confluence limit of (normalized) 4-point functions in Liouville Field Theory (see [17] and [16] for details). The 4-point function itself arises as the Painlevé VI τ -function, and in the Black Hole case in the Kerr-de Sitter conformally coupled scalar field scattering [42]. The confluence limit of the Painlevé VI is the same as the zero cosmological constant limit. This can be seen somehow as an extension of the ideas of [30], where conformal symmetry was also used to determine Black Hole scattering in the extreme Kerr limit. As it happens, the correspondence seems exact, at least for a wide class of (vacuum and charged) Black Holes in four dimensions.
While this seems surprising at first glance, one cannot help but wonder whether this fact is yet another facet on the intricate relationship between integrable systems and Conformal Field Theory [43]. General Relativity solutions like the (static, axisymmetric) black hole systems consist an integrable sector of Einstein equations in four dimensions [44]. Requisition that the solutions of this integrable sector are singularity free basically imposes that there is a coordinate patch where the components of the metric are given in terms of rational functions. It follows then, apart from the non-trivial fact of separability of the equations, that the differential equations governing the propagation of fields in these backgrounds have its singularities easily classified, being Fuchsian or having a higher Poincaré index. Moreover, the isomonodromy flow provides a non-linear symmetry of these equations stemming from the theory of flat holomorphic connections and the Painlevé property. All in all, everything conspires to have the conformal blocks "solving" Fuchsian-type of differ-JHEP11(2015)144 ential equations, since conformal blocks are in principle extracted from the representation theory of the Virasoro algebra. We can then describe the chain of deep mathematical facts pointing to the relationship between four-dimensional black holes and Conformal Field Theory. It is a very interesting open problem to see if this holds in the generic case.