arxiv: v2 [hep-th] 1 Mar 2016

Transcription

1 KEK-TH-1867 Black rings at large D Kentaro Tanabe arxiv: v2 [hep-th 1 Mar 2016 Theory Center, Institute of Particles and Nuclear Studies, KEK, Tsukuba, Ibaraki, , Japan ktanabe@post.kek.jp Abstract We study the effective theory of slowly rotating black holes at the infinite limit of the spacetime dimension D. This large D effective theory is obtained by integrating the Einstein equation with respect to the radial direction. The effective theory gives equations for non-linear dynamical deformations of a slowly rotating black hole by effective equations. The effective equations contain the slowly rotating Myers-Perry black hole, slowly boosted black string, non-uniform black string and black ring as stationary solutions. We obtain the analytic solution of the black ring by solving effective equations. Furthermore, by perturbation analysis of effective equations, we find a quasinormal mode condition of the black ring in analytic way. As a result we confirm that thin black ring is unstable against non-axisymmetric perturbations. We also include 1/D corrections to the effective equations and discuss the effects by 1/D corrections.

2 Contents 1 Introduction 1 2 Effective equations Setup Effective equations Effective energy-momentum tensor /D corrections Black ring and its physical properties Ring coordinate embedding Black ring solution Quasinormal modes Black string limit /D corrections Phase diagram Summary 27 A Other solutions 28 A.1 Slowly rotating Myers-Perry black hole A.2 Slowly boosted black string B Trivial perturbations in ring coordinate 34 C 1/D corrections of effective equations 35

3 1 Introduction Stationary asymptotically flat black holes have only spherical horizon topology in four dimensions by the Hawking s topology theorem [1. This situation changes drastically in higher dimensions, and higher dimensional black hole has various horizon topology. As one concrete example of this fact, Emparan and Reall found the five dimensional black hole solution, black ring, whose horizon topology is S 1 S 2 [2. This is the first asymptotically flat black hole solution which has a non-spherical horizon topology. After this discovery many new five dimensional black hole solutions with multiple horizons have been constructed by the solution generating technique such as the inverse scattering method e.g., see [3 for a review. In these solutions the topology of each horizon is S 3 or S 1 S 2. We expect that black hole solutions can have much more various and non-trivial horizon topology in higher dimensions than five. To study this variety it would be the first step to extend the Emparan-Reall s discovery to higher dimensions, that is, the construction of the black ring in higher dimensions. However the solution generating technique cannot be applied to the asymptotically flat solution in higher dimensions than five. Then there are some efforts to construct the black ring in higher dimensions by numerical methods [4, 5 and by analytical approximation methods such as the blackfold method [6, 7 1. In this paper we consider the construction of the D dimensional black ring solution by using another analytical approximation method, the large D expansion method [10, 11. Recently this large D expansion method was extended to construct the stationary [12, 13 and time-dependent black hole solutions [14, 15. So it is interesting to see if we can obtain the black ring solution by the large D expansion method. The paper [13 discussed the effective theory and gave the effective equation for the large D stationary rotating black holes. This stationary solution was assumed to have the O1 horizon angular velocity at the large D limit. On the other hand the analysis by the blackfold method [6 found that the horizon angular velocity of the D dimensional thin black ring becomes Ω H = 1 1 D 3 R, 1.1 where R is a ring radius. This means that the horizon angular velocity of the thin black ring is O1/ D, not O1 at the large D limit. Thus the effective equation in [13 cannot be straightforwardly applied to the thin black ring. Hence, to construct the black ring solution by the large D expansion method, we should study the effective theory of the black hole with O1/ D horizon angular velocity anew. The purposes of this paper are to give the effective theory of black holes with O1/ D horizon angular velocity by including time-dependence and to construct the large D black ring solution by solving the effective equations. We call a black hole with O1/ D horizon angular velocity slowly rotating large D black hole in this paper. We expect that there is no black ring with O1 1 The blackfold method has achieved a great success along this direction. In fact black holes with various horizon topology, including expected and unexpected ones, has been found [8, 9. 1

4 horizon angular velocity. The horizon angular velocity of the black ring is determined by the balance condition between the string tension and centrifugal force by the boost, which is the horizon angular velocity. At large D the tension of a black string becomes small compared to its mass by OD 1 [11. As a result the horizon angular velocity should be also small at large D. Hence, generically, there is no black ring with O1 horizon angular velocity at large D. This is why we do not use the word slowly rotating black ring in our analysis. The large D effective theory for dynamical black holes has been considered in [14, 15, where it was shown that the large D expansion method gives simple effective equations for the dynamical black holes. The equations for dynamical black branes obtained in [15 can be solved very easily, and the numerical solution of the equations capture the nonlinear evolution and endpoint of the Gregory-Laflamme instability of the black brane [16 in higher dimensions than the critical dimension [17. We will study the large D effective equation for slowly rotating large D black holes in the similar manner with [15. Then we willseethatthelarged blackringisfoundasananalyticstationarysolutionoftheeffective equations. Furthermore, by perturbing effective equations, we can obtain a quasinormal mode condition of the black ring. As found numerically in [18, the quasinormal mode says that the thin black ring is unstable against non-axisymmetric perturbations, and relatively fat black ring is stable. Our large D effective equations describe the non-linear evolution of such non-axisymmetric instabilities of the black ring. We can also include the 1/D corrections to the effective equations. We find that the 1/D corrections give striking features to the instability modes of the black ring. This paper is organized as following. In section 2 we give the effective theory of slowly rotating large D black holes by the effective equations. We also discuss some general properties of the stationary solution of the equations and give formula for thermodynamic quantities such as the mass and angular momentum of the solution with 1/D corrections. In section 3 we construct the black ring solution analytically by solving the effective equations. Some physical properties such as quasinormal modes and phase diagram of the black ring will be also discussed there. We close this paper by giving discussion and outlook of this work in section 4. The appendices contain technical details and some useful byproducts of the main results in this paper. Especially the slowly rotating Myers-Perry black hole and slowly boosted black string will be rediscovered with their quasinormal mode frequencies analytically in Appendix A. Note that we perform the large D expansion of the Einstein equation to find D dimensional black hole solutions. In the expansion, defining n by D = n+4, 1.2 we use 1/n as the expansion parameter instead of 1/D in this paper. This definition of n is same with one of [6, 7 for black rings. 2

5 2 Effective equations We consider the large D effective theory for slowly rotating large D black holes. The effective equations will be given as equations for the energy and momentum density of a dynamical black hole. In this section we study general properties of solutions of the effective equations without specifying an embedding of solutions. 2.1 Setup The effective equations describe non-linear dynamical deformations of black holes. To obtain them, it is natural to use the ingoing Eddington-Finkelstein coordinate for the metric ansatz as ds 2 = Adv 2 2u v dv +u a dx a dr +r 2 G ab dx a dx b 2C a dvdx a +r 2 H 2 dω 2 n, 2.1 where X a = z,φ. A, u v, G ab, C a are functions of v,r,x a. u a and H are functions of v,x a. In this paper we use v as a time coordinate in this system and the effective theory obtained below, while t was used for it in [15. This metric ansatz partially fixes the gauge of the radial coordinate r. The residual gauge can be fixed in solving the Einstein equations. X a = z,φ are inhomogeneous spatial directions of black holes, and we take Φ as the rotational direction. We give large D scalings and boundary conditions of the metric functions. Large D scalings To solve the Einstein equation by 1/n expansion, we should specify behaviors of each metric functions at large D. In the analysis of the blackfold [6 they found that the horizon angular velocity of the D = n +4 dimensional thin black ring is given by Ω H = 1 n+1 1 R +OR 2, 2.2 where R is a ring radius. This result implies that the metric function C Φ in eq. 2.1 becomes O1/ n at large D. From this observation we define the slowly rotating large D black hole by C Φ = O1/ n 2. One may think that this assumption is restrictive for a black ring and describes the dynamics only of a thin black ring solution. However this is not correct. As we will see, solutions obtained in this paper cover wider range than one by the blackfold because we do not require the much larger ring radius than the ring thickness. In our analysis ring radius and thickness can be comparable, and our large D solution can describe the dynamics also of not-thin black ring. Inoursetting Φ is not akillingvector. Thusthemetricfunctions have Φ-dependences in general. Then, to have consistent 1/n expansions of the Einstein equations, we should 2 The stationary solution considered in [13 was assumed to have C Φ = O1. 3

6 assume Φ = O n 3. If we take Φ = O1, the Einstein equations have not only 1/n series but also 1/ n series, and the analysis becomes a bit involved. To treat the scaling Φ = O n it is useful to introduce new coordinate φ defined by Φ = φ n. 2.3 Then Φ = O n is equivalent to φ = O1. Summarizing, our metric ansatz for the slowly rotating large D black holes is ds 2 = Adv 2 2u v dv +u a dx a dr +r 2 G ab dx a dx b 2C a dvdx a +r 2 H 2 dω 2 n, 2.4 where x a = z,φ. Each metric functions have following large D expansions A = k 0 A k v,r,x a n k, C z = k 0 C k z v,r,x a n k+1, C φ = k 0 C k φ v,r,xa n k+1, 2.5 u v = k 0 u k v v,r,x a n k, u z = k 0 u k z v,x a n k+1, u φ = k 0 u k φ v,xa n k+1, 2.6 G zz = 1+ k 0 G k zz v,r,x a n k+1, G zφ = k 0 G k zφ v,r,xa n k+2, 2.7 and G φφ = Gz2 n 1+ k 0 G k φφ v,r,xa n k+1, H = Hz 2.8 Under these scaling assumptions, the Einstein equation can be consistently expanded in 1/n. φ is the rotational direction of the black hole, and we assume v and φ become the Killing vectors in a background geometry. This assumption implies that the asymptotic metric at far from the horizon does not have v and φ dependences. So the metric function Hz and the leading order of G φφ has only z dependences as seen in eq Finally we introduce new radial coordinate R to realize r = On defined by r n R =, 2.9 where a constant r 0 is a fiducial horizon size, and we set to r 0 = 1. 3 We consider the decoupled mode excitation [19. So we assume v = O1 and r = OD. r 0 4

7 Boundary conditions We give boundary conditions for metric functions. In this paper we consider non-linear dynamics of the decoupled mode excitations [19. The decoupled mode condition demands following boundary conditions in the asymptotic region R 1 A = 1+OR 1, C a = OR 1,G zφ = OR The dynamical solutions have the horizon-like surface where A vanishes at the leading order in 1/n expansion 4. We regard the surface as a horizon of the dynamical black hole. The boundary condition on the horizon is the regularity condition of the each metric functions. 2.2 Effective equations We solve the Einstein equation by the large D expansion. The leading order equations of the Einstein equations contain only R-derivatives, so we can integrate them easily. Then, the leading order solutions after imposing the boundary conditions are obtained as and A 0 = 1 p vv,x a, C a 0 = p av,x a, u 0 Hz v = 2.11 R R 1 H z 2, u 0 a = u 0 a z, G 0 zz = 0, G 0 zφ = p zp φ p v 1 R, 2.12 G 0 φφ = 2+ 2HzG zh z p 2 φ 1 Gz1 H z 2 logr+ Gz 2 p v R p v v,x a and p a v,x a are mass and momentum density of the solution. They are introduced as integration functions of R-integrations of the Einstein equations. We can see that the horizon position of this dynamical solution is at R = p v v,x a. Furthermore we found that the function Hz should satisfy the following condition 1 H z 2 +HzH z = 0, 2.14 from the boundary condition of G 0 zφ at asymptotic region R 1. In eq we already imposed the condition This condition is equivalent to the constant condition of the mean curvature of r = const. surface in the asymptotic region for static solutions [12 5. In this paper we consider small rotation Ω H = O1/ n, so the leading order solution has no effects of rotations and gives same equations with static solutions. The condition 2.14 is integrated as 1 H z 2 Hz = 2ˆκ, In this paper we consider the slowly rotating black holes. Thus the horizon position is determined mainly by g uu = A. 5 If we perform D 1+1 decomposition of the Einstein equation on R = const. surface, this condition is obtained from the momentum constraint as shown in [12. 5

8 where ˆκ is a constant related with the surface gravity of the horizon as we will see later. So we call ˆκ a reduced surface gravity. At the leading order the integration functions p v and p a are arbitrary functions of v,x a, and there is no constraint among them. u 0 a is also arbitrary function of z at the leading order. At higher order in 1/n, constraints for u 0 a appear, and we can determine their explicit forms. Solving the next-to leading order Einstein equations we obtain non-trivial conditions which p v v,x a and p a v,x a should satisfy. The non-trivial conditions can be obtained also from the momentum constraints on R = const. surface. These conditions give effective equations for the slowly rotating large D black holes. Their forms are v p v H z 2ˆκHz zp v 2 φ p v 2ˆκGz 2 + φp φ Gz 2 + H z Hz p z = 0, 2.16 v p φ H z 2ˆκHz zp φ 2 φ p [ φ 2ˆκGz p 2 φ Gz 2 φ 4ˆκ2 GzHz 2 +2G zhzh z 4ˆκ 2 GzHz 2 φ p v + H z Hz p z p φ p v + G zh z ˆκGzHz p φ = 0, 2.17 p v and v p z H z 2ˆκHz zp z 2 φ p [ z 2ˆκGz 2 + zp v + 1 p φ p z Gz 2 φ p v + H z p 2 z G z p 2 φ Hz p v Gz 3 + G z p v ˆκGz 3 φp φ + HzG z 2 H z+gzg z HzH zg z 4ˆκ 2 Gz 2 Hz 2 p v 1 2H z 2 2ˆκHz 2 p z = These equations describe dynamical non-linear deformations of mass and momentum density of the dynamical black hole. To solve these equations we should specify the embedding functions Gz and Hz. Before doing that, we study some general properties of solutions of these effective equations. Note that the effective equations for dynamical black strings obtained in [15 can be reproduced from eqs. 2.16, 2.17 and 2.18 as one example. As one simple solution of the effective equations, we consider a stationary solution p v = p v z, p a = p a z This solution has two Killing vectors v and φ. The effective equations can be solved under this ansatz. Eqs and 2.17 give p z = p vz 2ˆκ, p φz = ˆΩ H Gz 2 p v z

9 ˆΩ H is an integration constant of z-integration of eq To see the physical meaning of ˆΩ H, we rewrite the v,φ part of the leading order metric of this stationary solution as ds 2 v,φ = 1 p vz dv 2 + Gz2 dφ R n ˆΩ p v z 2 H R dv +O1/n From this expression we can see that Ω H = ˆΩ H / n gives the horizon angular velocity. Actually horizon generating Killing vector ξ is ξ = v + ˆΩ H φ = v +Ω H Φ The Killing vector ξ becomes null at the horizon R = p v z. Furthermore we can calculate the surface gravity κ of the black hole by using ξ as R=pvz κ = rξ µ ξ µ 2ξ r 1 H = n z 2 2Hz = nˆκ, 2.23 where x µ = v,x a. Here we omit O1/n terms for simplification. Thus the integration constant ˆκ represents the surface gravity of the black hole at the leading order in 1/n expansions. Finally, substituting the solutions2.20 into eq. 2.18, we obtain an equation for p v z = e Pz as P z H z Hz Pz [ G z 2 Gz 2 G z Gz + G z GzHzH z ˆΩ 2 H GzG z1 H z 2 HzH = z To solve this equation for Pz we should specify the functions Gz and Hz. The functions can be determined by embedding the solution into a background geometry. In next section we will embed the leading order solution into a flat background in the ring coordinate. Such solution describes dynamical black rings. In Appendix A we show other solutions such as the Myers-Perry black hole and boosted black strings by considering the embeddings into other backgrounds. 2.3 Effective energy-momentum tensor Next we study physical quantities of the solution. The mass and angular momentum of dynamical black holes can be evaluated by the effective energy-momentum tensor. The effective energy momentum tensor, T µν, is defined by 6 T µν = 1 [ [ K µν 8πG D Σ ḡ µν K, 2.25 Σ 6 The effective energy momentum tensor is originally defined on r = const. surface. Thus its component runs over v,x a,x I where x I is a coordinate on S n in eq However the metric components of g IJ have been taken as non-dynamical ones by the gauge choice. So we consider only x µ = v,x a components of the effective energy momentum tensor. 7

10 where x µ = v,x a in eq Σ is a r = const. surface in the asymptotic region R 1. G D is the gravitational constant in D dimensions. K µν is the extrinsic curvature of r = const. surface. The square bracket represents the background subtraction at Σ. ḡ µν is the background induced metric on Σ. The background metric ḡ µν in the asymptotic region R 1 is obtained from the leading order solution as ḡ µν dx µ dx ν = dv 2 +r 2 dz 2 + Gz2 n 1 2HzH zg zlogr ngz1 H z 2 +On 2 The background extrinsic curvature, say K µν, is calculated as K vv = O1/n, Kzz = 2ˆκ+O1/n, Kφφ = GzG zh z 2ˆκHz dφ O1/n, 2.27 in our gauge choice for the radial coordinate r. The trace part of the background extrinsic curvature, K, contains contributions from K IJ where x I is a coordinate on S n in eq It is obtained as K = n [ 2ˆκ+ 1 n 2ˆκ G zh z 2ˆκGzHz 2ˆκlogR +On By definition, the background extrinsic curvature K µν satisfies K µν K µν = O1/R, K K = O1/R Then we can compute the effective energy-momentum tensor T µν. The results are and T vv = nˆκp v 1+On 1,R 1, T va = 2ˆκp a a p v 1+On 1,R 1, πG D R 16πG D R T zz = Gz2 H zp z +Hz z p v +Hz φ p φ 1+On 1 16πG D Gz 2,R 1, HzR T zφ = 1 Gz2ˆκp z p φ p v z p φ + φ p z +2G zp φ p v 1+On 1,R 1, n 16πG D Gzp v R T φφ = 1 1 [ 2ˆκHz2ˆκp 2 φ n32πg DˆκHzp v R +p v φ p φ +Gz 2 v p v 2ˆκp v GzH zp v G zp v 2ˆκGzp z 1+On 1,R Then we can define the mass M and angular momentum J Φ of a dynamical black hole from this T µν as M = nˆκω n 8πG D J Φ = n ˆκΩ n 8πG D dz dφ GzHz n p v v,x a, 2.32 n dz dφ GzHz n p φ v,x a, 2.33 n 8

11 where Hz n term comes from the volume factor on S n. Ω n is the volume of unit S n. The n factor in JΦ comes from the relation between φ and Φ in eq For stationary solutions we can define the mass and angular momentum of the black hole by the Komar integral of each Killing vector. It can be shown easily that their definitions are equivalent. The horizon area of the solution, A H, becomes A H = Ω n dz dφ GzHz n p v v,x a, 2.34 n where we used the fact that the horizon is at R = p v. Then we find that the stationary solution satisfies the Smarr formula given by n+1 n+2 M = κ A+Ω H J Φ, πG D at the leading order in 1/n expansions. Note that since Ω H = O1/ n and κ = On, Ω H J Φ term does not contribute to the Smarr formula at the leading order. The mass and area of the solution can be calculated also for non-stationary solutions p v = p v v,x a. For this non-stationary solution, we can see that the Smarr formula 2.35 holds at the leading order in 1/n expansion. But this formula does not hold at the higher order for non-stationary solutions /D corrections Solving the next-to-next-to leading order of the Einstein equations, we can obtain 1/n corrections to the effective equations 2.16, 2.17 and The explicit forms of the 1/n corrections to the effective equations are not simple, so we show them in Appendix C. Here we show only physical effects by 1/n corrections to the solution. At first we define the mass and momentum density, p v and p a, up to 1/n corrections. The mass and momentum density are introduced as integration functions of R integrations of the Einstein equations. Thus we should specify the normalization of the integration functions. Expanding A and C a up to 1/n corrections we define the mass and momentum density, p v and p a, by coefficients of 1/R as A = 1 p vv,x a R +On 2,R 2, C a = p av,x a R +On 2,R This definition normalizes the integration functions up to O1/n by the asymptotic behavior of the metric functions. Using this definition and the next-to-next-to leading order stationary solutions, the surface gravity up to 1/n correction is obtained as [ κ n = ˆκ 1 G zh z n 4ˆκGzHz + ˆκlogp v + H z p z ˆκ p 2 φ + 2Hz p v 2Gz 2 p 2 +On v One may think that this expression is strange since the surface gravity κ seems not to be constant when the solution is stationary where p v = p v z and p a = p a z. However we 9

12 can show that [ d G zh z dz 4ˆκGzHz + ˆκlogp v + H z p z ˆκ p 2 φ + 2Hz p v 2Gz 2 p 2 = 0, 2.38 v by using the effective equations 2.16, 2.17 and 2.18 for the stationary solution. Hence the surface gravity becomes constant for the stationary solution as expected. The horizon angular velocity Ω H = ˆΩ H / n is also obtained as p φ ˆΩ H = Gz G zh z p v n2ˆκ 2 Gz 3 Hz p φ logp v p v +On The first term in r.h.s. is equivalent to ˆΩ H at the leading order as seen in eq But, at higher order in 1/n expansion, they are different. We can show that, when the solution is stationary, the horizon angular velocity is constant up to 1/n correction by using the effective equations with 1/n corrections given in Appendix C. We can calculate the mass, angular momentum and area of the solution up to O1/n. The results are complicate, so we show the results only for the stationary solutions p v = p v z and p a = p a z. The mass and angular momentum formula up to O1/n become M = nˆκω n dz dφ GzHz n Mz, πG D n J Φ = n ˆκΩ n dz dφ GzHz n J φ z, πG D n where Mz = p v + 1 n J φ z = p φ 1 n [ p v p z H z 2ˆκ Hz p v G zh z 4ˆκ 2, 2.42 GzHz [ p z p φ H z 2ˆκp v Hz + p φ G zh z 4ˆκ GzHz The mass and angular momentum can be obtained from the Komar integral for the stationary solution or from the effective energy momentum tensor for general time-dependent solution. These definitions coincide for the stationary solution. The formula of the horizon area also has the 1/n correction as A H = Ω n dz dφ n GzHz n A H v,x a 2.44 where A H = p v + 1 n [ p v logp v 1 p 2 φ 2Gz p v Then we can see that the Smarr formula, can be satisfied up to 1/n corrections. n+1 n+2 M = κ A H +Ω H J Φ, πG D 10

13 3 Black ring and its physical properties Next we solve the effective equations in an explicit embedding. Especially we consider an embedding into a flat background in the ring coordinate and find a stationary solution, that is, the black ring solution analytically. 3.1 Ring coordinate embedding The leading order metric has the following asymptotic form at R 1 ds 2 R 1 = dv 2 Hz +2 1 H u0 z 2dv a zdxa dr n +r 2 dz 2 +r 2 Gz 2 dφ 2 +r 2 Hz 2 dω 2 n. 3.1 We embed this leading order metric into a flat background in the ring coordinate. The D = n+4 dimensional flat metric in the ring coordinate is [20 7 [ ds 2 = dt 2 R 2 R 2 dr 2 + R+rcosθ 2 R 2 r 2 +R2 r 2 dφ 2 +r 2 dθ 2 +sin 2 θdω 2 n, 3.2 where 0 r R, 0 θ π and 0 Φ 2π. R is a ring radius. r = 0 is the origin of the ring coordinate. The asymptotic infinity and the axis of Φ-rotation are at r = R. In this ring coordinate r = const. surface has a topology of S 1 S n+1. θ = 0 is an inner equatorial plane, and θ = π is an outer equatorial plane. Remembering the definition of R = r/r 0 n, we can embed the leading order induced metric on R = const. surface into the flat background in this ring coordinate by r = r 0 where r 0 is a constant. As done in the previous section we set to r 0 = 1. Then, comparing eqs. 3.1 and 3.2 on r = 1 surface, we find that the embedding gives following identifications Hz = Rsinθ R+cosθ, Gz = R R 2 1 R+cosθ, dθ dz = R+cosθ. 3.3 R It is easy to confirm that this identification actually satisfies eq We can calculate the surface gravity, κ = nˆκ, by these identifications as ˆκ = 1 H z 2 2Hz = R 2 1 2R. 3.4 The embedded solution is the black ring solution because the horizon topology is now S 1 S n+1. From R r, the ring radius should be larger than unity R 1. R 1 corresponds to the fat black ring, and thin black ring is described by R 1. 7 In this note we use the ring coordinate by r,θ in [20. Of course the following analysis can be done also by using x,y coordinate in [20. 11

14 Using identifications 3.3 and the surface gravity 3.4, we obtain the leading order effective equations for the black ring from eqs. 2.16, 2.17 and 2.18 as v p v + R+y1+Ry R y p v R+y2 R 2 1 RR 2 1 3/2 2 φ p v + R+y2 R 2 R 2 1 φp φ + 1+Ry R 1 y 2p z = 0, 3.5 and v p φ + R+y1+Ry R y p R 2 φ R+y2 1 RR 2 1 3/2 2 φ p φ [ + R+y2 p 2 R 2 R 2 1 φ 1+2Ry +R2 φ R 2 φ p v 1 p v + 1+Ry R 1 y 2 p z p φ p v + 21+Ry R R 2 1 p φ = v p z + R+y1+Ry R y p z R+y2 R 2 1 RR 2 1 3/2 2 φ p z R+y 1 y 2 y p v + R+y2 R R 2 R 2 1 φ R+y2 1 y 2 R 3 R 2 1 [ p z p φ p v + 1+Ry R p 2 z 1 y 2 p v 1 y 2 p 2 φ + 2R+y2 1 y 2 p v R 2 R 2 1 3/2 φ p φ + R+y R Ry y2 +R 2 2y 2 1 R p z = R 2 11 y 2 Here we introduced a coordinate y defined by y = cosθ. 3.8 Our effective equations 3.5, 3.6 and 3.7 describe non-linear dynamical deformations of the black ring from thin R 1 to not-thin R > 1 region. 3.2 Black ring solution The black ring solution is obtained as the stationary solution of the effective equations 3.5, 3.6 and 3.7. The stationary solution is given by As done in eq. 2.20, we find from eqs. 3.5 and 3.6 p v = e Py, p a = p a y. 3.9 p z y = R+y 1 y 2 p vy, p R 2 φ = ˆΩ R 2 R 2 1 H 1 R+y 2 p vy p v 12

15 Furthermore eq. 3.7 gives an equation for p v = e Py as eq P y+ 2 R+y P y R R+y1+Ry + ˆΩ 2 R 2 R H R+y 4 1+Ry = This equation contains a pole at y = 1/R in the source term. The solution can have a singular behavior by the source term at y = 1/R, too. To obtain a regular solution at y = 1/R, ˆΩ H should be ˆΩ H = R 2 1 R If ˆΩ H doesnottakethisvalue, thefunctionpyhasalogarithmicdivergenceaty = 1/R. Thus, as seen in five dimensional black ring [2 and higher dimensional black ring by the blackfold [6, 7, the regularity condition determines the horizon angular velocity of the black ring. The solution of eq under the condition 3.12 is obtained analytically as Py = p 0 + d 0 1+Ry1+Ry +2RR+ylogR+y + R+y 2R 2 R+y 2, 3.13 where p 0 and d 0 are integration constants. They are degree of freedom associated with trivial deformations and do not affect physical properties of black rings. p 0 is the 1/n redefinition of r 0, and d 0 comes from the redefinition of φ coordinate of the ring coordinate as discussed in Appendix B. Summarizing above results, we found the black ring solution at the leading order of large D expansion by the metric ds 2 = 1 p vz dv dv 2 +2 R 2ˆκ u0 a zdxa n 2 R+y 2ˆκR +r 2 dz 2 + ˆΩ 8ˆκ 2 R H R+cosθ +r 2 R2 R 2 1 R+cosθ 2 dr p vz dz R n + ˆΩ R 2 R 2 1 H p v z dzdφ R n 2 1 2RR+cosθ R 2 1 p v z R+cosθ 2 R dφ dv n logr n + ˆΩ 2 R 2 R 2 1 p v z dφ 2 H R+cosθ 2 nr n + r2 R 2 sin 2 θ R+cosθ 2dΩ2 n, 3.14 where R = r/r 0 n with r 0 = 1. The coordinates z and θ are related by eq ˆκ and ˆΩ H are given in eqs. 3.4 and The function p v = e Py is the solution of the equation 3.11, and it is obtained in eq with y = cosθ. Our black ring solution breaks down at R = 1, and the very fat black ring cannot be captured by our large D solution 8. However, the black ring with not so large ring radius can be described by eq So we can study the properties of not only the thin black ring but also not-thin black ring R > 1 by using our solution A very fat black ring means the solution with R = 1+On 1. 13

16 3.3 Quasinormal modes We investigate quasinormal modes of the black ring solution. The quasinormal modes are obtained by perturbation analysis of the effective equations around the black ring solution. The perturbation ansatz is p v v,y,φ = e Py 1+ǫe iωv e imφ F v y, 3.15 p z v,y,φ = R+y 1 y 2 p 2ˆκR vy 1+ǫe iωv e imφ F z y, 3.16 p φ v,y,φ = RR2 1 3/2 R+y 2 p v y 1+ǫe iωv e imφ F φ y, 3.17 where we used eq There is one remark on the quantum number m associated with φ. In the ring coordinate 3.2 the coordinate Φ has the range of 0 Φ 2π. Thus the quantum number m Φ associated with Φ is quantized as m Φ = 0,±1,±2,... From the relation 2.3, these quantum numbers are related by m = m Φ n So m in the perturbations can take non-integer values in general. Perturbing the effective equations 3.5, 3.6 and 3.7 with respect to ǫ, we obtain perturbation equations for F v y, F z y and F φ y. The perturbation equations have a pole at y = 1/R again. To solve the perturbation equations we should impose regularity conditions. If we specify the behavior of the perturbation fields at the pole y = 1/R as the regularity condition by F v y 1+Ry l 1+O1+Ry, 3.19 where l is a non-negative integer, we get one non-trivial condition for the frequency ω as 1 [ R 9 ω 3 +ir 6 R 2 1 3m 2 +3imR R 2 1m 2 +imr+lr 2 ir 3 ω +3l 2R 2 ω 2 R 3 R 2 1 3m 4 +6im 3 R+23l 4m 2 R 2 +2i3l 2mR 3 +3l 1lR 4 ω ir 2 1 3/2 m 6 +3im 5 R+3l 2m 4 R 2 +6il 1m 3 R 3 4 7l+3l 2 m 2 R 4 +3ill 1mR 5 +l 2 l 1R 6 = This is the quasinormal mode condition for the black ring. One may feel strange about this derivation since we derive the quasinormal mode condition from the local condition 3.19 at y = 1/R. If the condition 3.20 is satisfied, the perturbation actually becomes regular at y = 1/R, butits global structure, e.g., of F v y, is still unknown. Furthermore the quantum number l associated with the harmonics in the ring coordinate should be given by the global solution of F v y. But the quantum number is introduced by the local 14

17 Figure 1: Plots of the leading order quasinormal modes with l = 0 and m = 2, ω l=0 ±, of the black ring. The left panel shows the real part of the frequency normalized by the reduced surface gravity ˆκ. The right panel shows the imaginary part. The black and gray line represent ω l=0 + and ω l=0 respectively. The real part of ω l=0 + and ω l=0 are same. ω l=0 + is unstable for R > 2. condition 3.19 in our derivation at the large D limit. This localization of the quantum number associated with the harmonics at the large D limit has been observed in the spherical harmonics and spheroidal harmonics in [13. So we expect that the same feature would appear in the harmonics in the ring coordinate. It is interesting to investigate this property in detail, although we do not pursue this structure in this paper. In the following we give results derived from eq for non-axisymmetric and axisymmetric modes separately. Non-axisymmetric modes m 0 The quasinormal mode condition 3.20 can be solved in a simple form for l = 0 by ω l=0 R ± = 2 1 [ ˆm±iˆm1 ˆm, ω l=0 0 = R R 2 1 R [ ˆm iˆm 2 2, 3.21 where ˆm = m/r. ω l=0 + is an instability mode when R > m. ω l=0 ± can be understood as the quasinormal mode of the boosted black string as we will see below. We regard ω l=0 0 in eq as a gauge mode in this paper 9. At large radius limit R 1 the quasinormal modes should reproduce the quasinormal modes of the boosted black string. Although ω l=0 ± has corresponding modes, ω l=0 0 does not have. So it might be natural to consider that ω l=0 0 is a gauge mode. In Figure 1 we show plots of the quasinormal mode ω l=0 ± with m = 2. ω l=0 + becomes stable in relatively fat region R < m. This plot reproduces the behavior of the numerical results in [ We could not find any 9 We have not still studied residual gauge of v,z,φ coordinate in detail. The gauge of r-coordinate was fixed in eq So we should check the gauge transformation of perturbations to check which modes are physical. Actually there are residual gauge. In Appendix B we study stationary one of them. It is a bit involved task to eliminate whole gauge modes. So instead we use another simple method to identify physical modes here. 10 The perturbation considered in [18 is the mode with m Φ = 2, not m = 2 as seen in eq Thus our result for m = 2 does not exactly correspond to the mode in [18. However, in D = 5 n = 1, we can 15

18 instabilities for l 0 and m 0 mode perturbations. Note that the instability mode ω l=0 + always saturates the superradiant condition 11 [ Re ω l=0 + mˆω H = The nature of this saturation at the large D limit is unclear. The dynamically unstable mode of the Myers-Perry black hole also shows the peculiar relation with the superradiant condition at the large D limit [13, 21. For equally spin Myers-Perry black hole, the coincidence of the onset of the superradiant and dynamically unstable regime was confirmed up to 1/D corrections [21, and it is consistent with the numerical result [22. For singly rotating Myers-Perry black hole the coincidence was also found at the large D limit [13. However it is known numerically that their onsets are at different rotation parameters for the singly rotating Myers-Perry black hole in finite dimensions [23, 24. Thus the coincidence of the onsets of superradiant and dynamically unstable regime is only the feature at the large D limit for the singly rotating Myers-Perry black hole. As we will see below, the saturation of the superradiant condition of the black ring does not always hold at higher order in 1/n expansion, and we find that onsets of the superradiant and dynamically unstable regime are same also for the black ring at the large D limit. Endpoints of instability? We found non-axisymmetric instabilities of the black ring solution. One may ask what its endpoint is. In [25 they discussed that the instability of the five dimensional thin black ring leads to the fragmentation into black holes by following two reasons. One is that the dynamical timescale to release the inhomogeneity by gravitational wave radiations is much longer than the timescale of the black ring instability. Thus the inhomogeneity by the instability of the thin black ring will grow in time. Another reason is that the fragmenting solution can be more entropic than the black ring. This can be understood by the fact that the non-uniform black string is less entropic than localized black holes in five dimensions. Then we expect that going to the fragmenting solution is preferable as the endpoint of the instability than going back to the fatter stable black ring. So the growth of the inhomogeneity of the black ring would not stop, and the fragmentation occurs. In our case this second reason of the discussion cannot hold. The non-uniform black string can be the endpoint of the Gregory-Laflamme instability as seen in [15 since the non-uniform black string becomes more entropic solution in higher dimensions than the critical dimension [17. Then, also for the black ring instability, we can say same thing. In enough higher dimensions the non-uniform solution is more entropic than the fragmenting solutions. So we expect that the black ring would not go expect that m Φ = 2 and m = 2 are not different so much. 11 The quantum numbers m and m Φ has the relation Then the superradiant factor is Re [ω m ΦΩ H = Re [ω mˆω H. 16

19 to the fragmentation by the instability. On the other hand the first reason can be applied to our case. Actually the inhomogeneity by the black ring instability and its rotation give time-dependent quadrupole moments to gravitational fields. Then the timescale of gravitational wave radiations by time-dependent quadrupole moments, t GW, is estimated as done in [25 by t GW 1 R n+1 Ω H G D M Rn+2 R n G D M en/2 R, 3.23 where the exponential factor e n/2 comes from Ω n in G D M as seen in eq We restored theblack ringthickness r 0 for usefulness. Thetimescale of the black ringinstability, t BR, is estimated by the quasinormal mode frequency as t BR = [ Im ω l=0 + r 0 1 r So t BR is O1 quantity in 1/n expansion. Then, remembering that our large D black ring solution satisfies R > r 0, the ratio of timescales is exponentially small in n at large D t BR t GW e n/2 R n 1 O r 0 This implies that inhomogeneity by the black ring instability cannot be dissipated by the gravitational wave radiation to the infinity in 1/n expansion. But the inhomogeneous black ring solution would not fragment into small black holes unlike in five dimensions as discussed above. Thus the inhomogeneity would stop at some point like the nonuniform black string observed in [15. Then we reach the conclusion for the endpoint of the black ring instability that the black ring evolves to a non-uniform black ring, NUBR, as a metastable solution at large D 12. The metastable solution at large D means that the solution is stationary in 1/n expansion, but it is not stationary in Oe n/2 by the gravitational wave emissions. The stationary solution with inhomogeneities along rotating direction is prohibited by the rigidity theorem [1, 26. Our statement for the endpoint of the black ring instability is consistent with the rigidity theorem since NUBR is not stationary in Oe n/2. Such exponentially suppressed evolutions cannot be described by our effective equations 3.5, 3.6 and 3.7 by 1/n expansions. It is interesting to investigate this possibility by solving the effective equations for the black ring 3.5, 3.6 and 3.7 directly as the investigation of the endpoint of the Gregory-Laflamme instability in [15. If there is a metastable solution such as NUBR, we can find a stationary and non-axisymmetric solution of the effective equations. This statement for the endpoint of the black ring instability is only for the leading order result in 1/n expansion. If we include 1/n corrections to the effective equations, the results show the dimensional dependences, and we can observe the critical dimension in which the stability of NUBR would change. Actually we have investigated the critical 12 The author thanks Roberto Emparan for the suggestion of this possibility. 17

20 dimension of the non-uniform black string by considering 1/n corrections [29. So our conjecture that the endpoint of the black ring instability is a stable NUBR is valid only at the leading order in 1/n expansion. Axisymmetric modes m = 0 For m = 0 modes, we can also solve eq in a simple form as ω m=0 ± = ± l 1 il Note that thenumerator of eq is acubicalgebraic equation of ω, butonesolution of the numerator equation is canceled with the denominator of eq for m = 0 modes. So there are only two modes as eq The quasinormal mode frequencies 3.26 are same with one of the Schwarzschild black hole at large D [19. This is because we consider the decoupled mode excitations. In the decoupled mode excitation, the dynamics of perturbations is determined almost locally on the horizon. For m = 0 modes, the perturbation does not have interactions along φ direction. Thus the perturbation feels the horizon as one of the D 1 dimensional Schwarzschild black hole, and its quasinormal mode becomes one of the Schwarzschild black hole. The axisymmetric perturbation does not have any instabilities at the leading order of 1/n expansion. Thus the instability found in [25, 27 for the fat black ring against axisymmetric perturbations is not contained in our setup. Such instability occurs due to interactions between horizons of the black ring, so it may be non-decoupled mode instability. However, by including 1/n corrections, we find a suggestion for the axisymmetric instabilities of the black ring as we will see later. 3.4 Black string limit As one useful observation let us see the large ring radius limit R 1 of the black ring solution found above. This limit gives the expression of the black ring as a boosted black string [6. At the large radius limit we can set the behavior of Py in eq to Py = O1/R, 3.27 by p 0 = 0 and d 0 = 0. Then we get the large radius limit of the black ring solution as ds 2 = 1 1 dv 2 +2dvdr 1 dvdx+ 1 1 dx 2 +r 2 dω 2 R nr nr n, 3.28 where we defined the black string direction dx by dx = Rdφ/ n = RdΦ. This large radius limit solution is actually regarded as the large D limit metric of the boosted string with the boost velocity sinhα = 1 n,

21 as found in [6. Using this boost relation we can reproduce the large radius limit of the quasinormal modes, ω l=0 ±, of the black ring from the quasinormal modes of the black string. The quasinormal modes of the black string, ω± BS was obtained in [11 by large D expansion. For the perturbation e iωv e ikx, the leading order result is ω BS ± = ±iˆk1 ˆk, 3.30 where ˆk = k/ n. Now the boost transformation on the black string is acting by dv coshαdv sinhαdx, dx coshαdx sinhαdv This transformation on ω BS ± gives quasinormal modes of the boosted black string, ω bbs ±, as ω± BS ωbbs ± = ksinhα+ωbs ± coshα Using eq and ˆk = k/ n, we obtain the quasinormal modes of the boosted black string with the boost velocity 3.29 as ω bbs ± = ˆk±iˆk1 ˆk This quasinormal modes precisely reproduce the large radius limit of quasinormal modes of the black ring by identifying ˆk = ˆm 13. In Appendix A we give the direct derivation of the quasinormal modes of the boosted black string from the effective equations /D corrections We can obtain 1/n corrections to results obtained above by solving the effective equations up to O1/n. In the following we set to p 0 = 0 and d 0 = 0 for Py in eq p v z and p φ z for the stationary solution are p v y = e Py [ 1+ P1 y, 3.34 n p φ y = RR2 1 3/2 R+y 2 e Py [ 1+ 1 P 1 y 21+Ry n R 2 1 Py+ logr R 1 R 2, The boost transformation is one by the exact symmetry of the black string. Thus we can obtain the quasinormal mode of the boosted black string from the black string. As for the singly rotting Myers- Perry black hole, the boost transformation comes from an approximation symmetry appearing only at the leading order of large D limit. Thus quasinormal modes of the Myers-Perry black hole do not follow boost transformation rules. 19

22 where P 1 y is given by P 1 y = p 1 + d 1 R+y + R2 1 2 R 4 R+y 2 logr R 1 1+Ry R R2 1y R1 2y 2 2R 3 R+y 2 logr+y 2 21+Ry R 3 R+y Li 2 1 [ 4R 4 R 2 1R+y Ry +2y 2 +2R 8 4 3y 2 +2R 7 y5 4y 2 R 2 5 6y 2 12R 5 y1+y 2 R y 2 +5y 4 2R 3 y10+y 4 +R y 2 5y 4 logr+y [ + R 3 R 2 1R+y 3 2+R 7 y y 2 +3R 2 1 2y 2 +R 6 2+y 2 +3R 3 y2+y 2 2Ry3+y 2 +3R 5 y1+y 2 R y 2 +y RyR 2 1R+y 2 logr R p 1 andd 1 areintegration constantsof1/ncorrections. Li 2 xisthepolylogarithmfunction. This solution has the surface gravity [ R κ = n R 2n +On 2, 3.37 and the horizon angular velocity Ω H = 1 [ R 2 1 n R n R 2 2logR R 1 2R 2 +On The horizon angular velocity up to O1/n can be obtained by imposing the regularity of p v z at the pole y = 1/R up to O1/n. We can see that this horizon angular velocity reproduces eq. 2.2 by the blackfold method at the large radius limit up to O1/n. The surface gravity has the following large radius limit [ 1 1 lim R κ = n 2 = n 2coshα 2n +On 2 1+On 2, 3.39 Then this expression implies that the surface gravity is reproduced by the boosted black string with the boost velocity 3.29 up to 1/n corrections. We can also compare the 1/n corrections of our results with the 1/R 2 corrections by the blackfold method [7 to the horizon angular velocity. However, in this comparison, there are some subtle things such as a definition of R, so the comparison is not clear. The fact that 1/n corrections do not contain 1/R contribution to Ω H is consistent with the result in [6. Finally we give 1/n corrections to quasinormal modes of the black ring. ω l=0 ± up to 1/n corrections is ω l=0 ± = [ R 2 1 R ˆm±iˆm1 ˆm+ δˆωl=0 ±, 3.40 n 20

23 Figure 2: Plots of the quasinormal mode ω l=0 + with l = 0 and m = 2 in n = 10 of the black ring up to O1/n corrections. The left panel shows the real part of the frequency normalized by the reduced surface gravity. The right panel shows the imaginary part. The dashed line represents the leading order result. The thick line is the result with O1/n corrections. where δˆω l=0 ± = 1 [ 2ˆmR 2 ˆm 2 R 2 4+4R 2ˆm 2 +iˆm2+r 2 +16ˆm 2 +8ˆm 4 +2ˆm 2 1 2iˆmlogR R 1 2ˆm2+3R 2 4ˆm 2 +i2+3r 2 +1ˆm 2 2R 2 10ˆm 4 2ˆm 2 logr R At the large radius limit R with fixed ˆm, this quasinormal mode reduces to one of the boosted black string with the boost parameter given in eq as seen in Appendix A. Thuswe could confirmthat the blackfold analysis by [6 is correct upto1/n corrections. Fromthis1/n correction ofω l=0 +, wefindthat theblack ringbecomesunstableagainst non-axisymmetric perturbations for R > R D where [ R D = m logm m 1 n 2m 2 +On At the large radius limit this threshold ring radius R D gives the marginally stable wave number ˆm D as ˆm D m R D = 1+On 2,R This corresponds to the Gregory-Laflamme mode of the boosted black string with the boost velocity 3.29 as seen in Appendix A. In Figure 2 we show the plots of the quasinormal modefrequencyω l=0 + withm = 2inn = 10 uptothe1/n corrections. Theplots infigure 2 show that the results become much closer to numerical results [18 by 1/n corrections. At R = 1 our quasinormal mode formula breaks down due to the term logr R 1. Thus the quasinormal mode frequency behavior around R = 1 is not reliable. One interesting observation on 1/n corrections to quasinormal modes is the relation of the superradiant and dynamically unstable regime. At the leading order the real part of 21

24 Figure 3: The plot of the superradiant factor is shown for the dynamically unstable mode ω l=0 + with m = 2 and n = 10. The superradiant factor can be negative and positive in dynamically unstable regime R > m, while it is always positive in stable regime. the quasinormal mode frequency ω l=0 + is marginal in the superradiant condition [ Re ω l=0 + mˆω H = On By including the 1/n corrections, we see that the real part of the quasinormal mode frequency ω l=0 + deviates from the superradiant condition. Actually the real part of the frequency does not always satisfy the superradiant condition [ Re ω l=0 + mˆω H = 2ˆκ n ˆm 12 R 2 4ˆm+2R 2ˆm 2 R 2 +On In Figure 3 we show the plot of eq for n = 10 and m = 2. The superradiant regime can be seen in the dynamically unstable regime R > m. But there is also the dynamically unstable mode which does not satisfy the superradiant condition at R > R S > m, R S = m m 2 +m m Thus the black ring has dynamically unstable modes which satisfy the superradiant condition for R S > R > m, and the superradiant condition becomes not to be satisfied by unstable modes at R > R S. This behavior of the dynamically unstable modes is the striking property of the black ring also seen in numerical results [18. In fact the all known dynamically unstable modes of black holes such as one of the Myers-Perry black holes always satisfy the superradiant condition. The physics behind this feature is not still clear. We can also obtain 1/n corrections to the quasinormal modes with m = 0 in a simple form. The result is ω m=0 ± = [ R 2 1 R ± l 1 il 1+ δˆωm=0 ± n,

25 where [ δˆω m=0 ± = ± R2 l 1 2R 2 3l 5 7l2 l R 2 l l 5 R 4 l 1 [ i R2 l 1 2R 2 2l 5 4l2 7l R l 5 l 1 R l 1 At the large radius limit R the quasinormal mode frequency ω m=0 ± deviates from one of Schwarzschild black hole [19 in 1/n corrections. This deviation can be understood as a boost effect of the black string as discussed in Appendix A. Around R = 1 the imaginary part of the 1/n correction δˆω m=0 ± becomes positive. This might mean that the axisymmetric perturbation becomes also unstable around R = 1. Actually we find that there is a marginally stable mode at R = R f as = l l 1+On 1, 3.49 R=Rf ω m=0 ± where R f is the marginal radius for the apparent axisymmetric instability of the black ring given by R f = 1 1 l n8l 1 +On The imaginary part becomes positive at 1 > R > R f, and the axisymmetric perturbation may be unstable there. As mentioned above our formula of the quasinormal mode breaks down around R = 1. Thus we cannot say immediately that the axisymmetric instability mode of the black ring is found from our quasinormal mode formula. However R f < 1 might suggest that the unstable region for the axisymmetric perturbations exists around R = 1. This may be related with the instability of the fat black ring [25, Phase diagram Let usdrawthephasediagramof theblack ringsolution obtainedbythelarge D expansion method above. To do it we collect the formula for thermodynamic quantities of the black ring. The mass M and the angular momentum J Φ formula are given in eqs and For the black ring, using the ring coordinate embedding by eq. 3.3 and the leading order stationary solution 3.14, these formula become where and M = nω n 8G D ˆM, JΦ = n Ω n 8G D ˆ JΦ, 3.51 ˆM = Jˆ Φ = dy R2 1e Py R 1 y 2 R n 1 y 2, 3.52 R+y dy R2 1 5/2 e Py R n 1 y 2 RR+y y 2 R+y 23